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\n",j,res);
}
//printf("%d %d\n",i,res);
ans=min(ans,res);
}
printf("%I64d\n",ans);
return 0;
}
summary
- ==还是不要骚想
- 要动手! 不能瓜想 浪费时间
- 多积累! 数学好重要 QAQ!