【2019.8.31】

FFF团

先缩点成DAG 然后在拓扑排序过程中有两个点的话就退出

#include<bits/stdc++.h>
using namespace std;
#define Max(x,y) ((x)>(y)?(x):(y))
#define Min(x,y) ((x)>(y)?(y):(x))
#define Abs(x) ((x)<0?-(x):(x))
const int N=1000+5,M=6000+5,inf=0x3f3f3f3f,P=19650827;
int n,m,in[N];
template <class t>void rd(t &x){
    x=0;int w=0;char ch=0;
    while(!isdigit(ch)) w|=ch=='-',ch=getchar();
    while(isdigit(ch)) x=(x<<1)+(x<<3)+(ch^48),ch=getchar();
    x=w?-x:x;
}

int head[N],tot,hd[N],tt;
struct edge{int u,v,nxt;}e[M],E[M];
void add(int u,int v){
	e[++tot]=(edge){u,v,head[u]},head[u]=tot;
}
void Add(int u,int v){
	E[++tt]=(edge){u,v,hd[u]},hd[u]=tt;
}

int idx,Bcnt,dfn[N],low[N],bl[N];
stack<int>s;bool inst[N];
void tarjan(int u){
	dfn[u]=low[u]=++idx;
	s.push(u),inst[u]=1;
	for(int i=head[u],v;i;i=e[i].nxt){
		if(!dfn[v=e[i].v]) tarjan(v),low[u]=Min(low[u],low[v]);
		else if(inst[v]&&dfn[v]<low[u]) low[u]=dfn[v];
	}
	if(low[u]==dfn[u]){
		int v;++Bcnt;
		do{
			v=s.top();s.pop();
			bl[v]=Bcnt,inst[v]=0;
		}while(u!=v);
	}
}

queue<int>q;
bool topsort(){
	while(!q.empty()) q.pop();
	for(int i=1;i<=Bcnt;++i) if(!in[i]) q.push(i);
	if(q.size()>1) return 0;
	while(!q.empty()){
		int u=q.front();q.pop();
		for(int i=hd[u];i;i=E[i].nxt)
		if(!(--in[E[i].v])) q.push(E[i].v);
		if(q.size()>1) return 0;
	}
	return 1;
}

int main(){
	freopen("FFF.in","r",stdin);
	freopen("FFF.out","w",stdout);
	int T;rd(T);
	while(T--){
		memset(head,0,sizeof(head));
		memset(inst,0,sizeof(inst));
		memset(dfn,0,sizeof(dfn));
		memset(low,0,sizeof(low));
		memset(hd,0,sizeof(hd));
		memset(in,0,sizeof(in));
		idx=Bcnt=tot=tt=0;
		rd(n),rd(m);
		for(int i=1,u,v;i<=m;++i) rd(u),rd(v),add(u,v);
		for(int i=1;i<=n;++i) if(!dfn[i]) tarjan(i);
		for(int i=1;i<=m;++i)
		if(bl[e[i].u]!=bl[e[i].v]) ++in[bl[e[i].v]],Add(bl[e[i].u],bl[e[i].v]);
		if(!topsort()) puts("Light my fire!");
		else puts("I love you my love and our love save us!");
	}
    return 0;
}

maple做数学题

有一天，maple上数学课，他学会了如何求一个数的因子，于是老师给他布置了一道题，在区间[L,R]里面，找出所有满足下列条件的数字：这个数的第二小的因子为K。找到这些数字以后，maple还要去把这些数加起来，请问最终结果是什么？

==素数筛 筛出来的v数组里面存的是i的最小素因数

前面加一个特判就好了

70昏

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define Max(x,y) ((x)>(y)?(x):(y))
#define Min(x,y) ((x)>(y)?(y):(x))
#define Abs(x) ((x)<0?-(x):(x))
const int N=100000+5,M=1e7+5,inf=0x3f3f3f3f,P=1e9+7;
ll l,r,k,n,sum=0;
template <class t>void rd(t &x){
    x=0;int w=0;char ch=0;
    while(!isdigit(ch)) w|=ch=='-',ch=getchar();
    while(isdigit(ch)) x=(x<<1)+(x<<3)+(ch^48),ch=getchar();
    x=w?-x:x;
}

int cnt=0,prime[M],v[M];
void Prime(){
	for(int i=2;i<=r;++i){
		if(!v[i]) v[i]=i,prime[++cnt]=i;
		for(int j=1;j<=cnt&&prime[j]*i<=r;++j){
			v[i*prime[j]]=prime[j];
			if(!(i%prime[j])) break;
		}
	}
}

int main(){
	freopen("math.in","r",stdin);
	freopen("math.out","w",stdout);
	rd(l),rd(r),rd(k);n=r/k+1;
	if(r/k<k){
		if(l<=k) return printf("%lld",k%P),0;
		return puts("0"),0;
	}
	Prime();
 	for(int i=Max(l,k);i<=r;++i)
	if(v[i]==k) sum=(sum+i)%P;
	printf("%lld",sum);
    return 0;
}

#include <bits/stdc++.h>

#define LL long long
#define LD double
#define FOR(i,a,b) for (int i = (a);i <= (b); i++)
#define DFOR(i,a,b) for (int i = (a);i >= (b); i--)
#define debug(x) cerr << "debug: " << (#x) << " = " << (x) <<endl;
#define PI acos(-1)
#define mp make_pair
#define pb push_back
#define itr iterator
#define bit(x) (1LL<<(x))
#define lb(x) ((x)&(-x))
#define sqr(x) ((x)*(x))
#define gn 3
#define l(x) ch[x][0]
#define r(x) ch[x][1]
#define y0 Y0
#define y1 Y1
#define y2 Y2
#define fir first
#define sec second

using namespace std;

const int N = 320000, NN = 10010, MM = 110;
const LL P = 1000000007LL;

int tot,prime[N],tim;
bool v[N];
LL f[NN][MM],inv2;

LL qpow(LL x,LL n){
	LL ans=1;
	for(;n;n>>=1,x = x*x%P) if(n&1) ans=ans*x%P;
	return ans;
}

LL S(LL n,int m){
	if(n==1) return 1ll;
	if(n==0) return 0ll;
	if(m==0) return n%P*(n%P+1)%P * inv2%P;
	if(n<NN && m<MM) return f[n][m];
	if(m && prime[m]>n){
		while(m && prime[m]>n) m--;
		return S(n,m);
    }
	return (S(n,m-1) + P - prime[m]*(LL)S(n/prime[m], m-1)%P) % P;
}

bool isprime(LL n){
	if(n==2) return 1;
	for(int i=2;1ll*i*i<=n;i++) if(n%i==0) return 0;
	return 1;
}

int main(){
	//freopen("math9.in","r",stdin);
	//freopen("math9.out","w",stdout);
	FOR(i,2,N-1) if(!v[i]){
		prime[++tot] = i;
		for(int j=i+i;j<N;j+=i) v[j] = 1;
	}
	FOR(i,1,NN-1){
		f[i][0] = (1ll*i*(i+1)/2) % P;
		FOR(j,1,MM-1) f[i][j] = (f[i][j-1] + P - prime[j]*(LL)f[i/prime[j]][j-1]%P) % P;
	}
	inv2 = qpow(2,P-2);
	int T;
	T=1;
	while(T--){
		LL L,R,K;
		cin>>L>>R>>K;
		
		if(!isprime(K)) cout<<0<<endl;
		else if(K>N){
			if(K<=R && K>=L) cout<<K%P<<endl;
			else cout<<0<<endl;
		}else{
			int tmp = 0;
			while(prime[tmp+1]<K) tmp++;
			cout << (S(R/K,tmp)*K%P - S((L-1)/K,tmp)*K%P + P) % P << endl;
		}
	}
	return 0;
}

数字

有一个n个数字的列表，你可以对列表进行两种操作

1.用x的代价删除一个数字

2.用y的代价选择一个数字增加1

两种操作没有次数限制，问使得这个列表所有数字的gcd（最大公因数）大于1的最少花费是多少？

枚举可能的gcd 然后就暴力...

改了之后的30昏==

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define Max(x,y) ((x)>(y)?(x):(y))
#define Min(x,y) ((x)>(y)?(y):(x))
#define Abs(x) ((x)<0?-(x):(x))
const int N=500000+5,M=1000000+5,inf=0x3f3f3f3f,P=19650827;
int n,x,y,a[N],cnt=0,mx,pri[M];
bool v[M];
template <class t>void rd(t &x){
    x=0;int w=0;char ch=0;
    while(!isdigit(ch)) w|=ch=='-',ch=getchar();
    while(isdigit(ch)) x=(x<<1)+(x<<3)+(ch^48),ch=getchar();
    x=w?-x:x;
}

int gcd(int a,int b){return !b?a:gcd(b,a%b);}

int main(){
	freopen("T3.txt","r",stdin);
//	freopen("numbers.out","w",stdout);
	rd(n),rd(x),rd(y);
	ll sum,ans=1000000000000000000ll;
	for(int i=1;i<=n;++i) rd(a[i]),mx=Max(mx,a[i]);
	for(int i=2;i<=50;i++){
		if(!v[i]){
			pri[++cnt]=i;
			sum=0;
			bool ye=1;
			for(int j=1;j<=n;j++){
				int del=i-(a[j]%i);
				if(del==i) del=0;
				sum+=Min((ll)x,(ll)y*del);
				if(sum>=ans) {ye=0;break;}
			}
			if(ye) ans=sum;
		}
		for(int j=1;j<=cnt&&i*pri[j]<=mx;j++){
			v[i*pri[j]]=1;
			if(i%pri[j]==0) break;
		}
	}
	printf("%lld",ans);
    return 0;
}

数据太水 约数到50就能过了==

约数50

#include
using namespace std;
#define ll long long
#define Max(x,y) ((x)>(y)?(x):(y))
#define Min(x,y) ((x)>(y)?(y):(x))
const int N=500000+5,M=1000000+5,inf=0x3f3f3f3f,P=19650827;
int n,x,y,a[N],cnt=0,mx,pri[M];
bool v[M];
template void rd(t &x){
    x=0;int w=0;char ch=0;
    while(!isdigit(ch)) w|=ch=='-',ch=getchar();
    while(isdigit(ch)) x=(x<<1)+(x<<3)+(ch^48),ch=getchar();
    x=w?-x:x;
}
int main(){
	freopen("T3.txt","r",stdin);
//	freopen("numbers.out","w",stdout);
	rd(n),rd(x),rd(y);
	ll sum,ans=1000000000000000000ll;
	for(int i=1;i<=n;++i) rd(a[i]),mx=Max(mx,a[i]);
	for(int i=2;i<=50;i++){
		if(!v[i]){
			pri[++cnt]=i;
			sum=0;
			bool ye=1;
			for(int j=1;j<=n;j++){
				int del=i-(a[j]%i);
				if(del==i) del=0;
				sum+=Min((ll)x,(ll)y*del);
				if(sum>=ans) {ye=0;break;}
			}
			if(ye) ans=sum;
		}
		for(int j=1;j<=cnt&&i*pri[j]<=mx;j++){
			v[i*pri[j]]=1;
			if(i%pri[j]==0) break;
		}
	}
	printf("%lld",ans);
    return 0;
}
    

==先把标程放这 放弃思考

#include<bits/stdc++.h>
using namespace std;
const int N=2000000;

long long a[N+100],s[N+100];
int n,x,y;
int main()
{
	freopen("numbers8.in","r",stdin);
	freopen("numbers8.out","w",stdout);
	scanf("%d%d%d",&n,&x,&y);
	for (int i=0;i<n;i++)
	{
		int t;
		scanf("%d",&t);
		a[t]++;
	}
	for (int i=N;i>0;i--)		 //Don因为要i到j的值（i<j）必须计算后缀和； 
		s[i]=s[i+1]+a[i]*(N-i);  //此时a表示的是数字在数轴上的分布； 
	for (int i=1;i<=N;i++)      //个数的前缀和； 
		a[i]=a[i-1]+a[i];
	long long ans=1LL*n*x;
	for (int i=2;i<=N;i++)   //枚举因子； 
	{
		long long res=0;
		for (int j=i;j<=N;j+=i)  //O（N） 
		{
			int l=j-i+1,r=max(l,j-x/y);
			long long tot=a[j-1]-a[r-1];
			res+=(s[r]-s[j]-tot*(N-j))*y+(a[r-1]-a[l-1])*x;  //分两部分去掉和加上；
			//由x，y的意义可知：最多加x/y个1，即增加的范围从max（l，j-x/y开始）； 
			//printf("  %d %d
",j,res);
		}
		//printf("%d %d
",i,res);
		ans=min(ans,res);
	}
	printf("%I64d
",ans);
	return 0;
}

summary

• ==还是不要骚想
• 要动手！ 不能瓜想 浪费时间
• 多积累！ 数学好重要 QAQ!