【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\n",j,res);
        }
        //printf("%d %d\n",i,res);
        ans=min(ans,res);
    }
    printf("%I64d\n",ans);
    return 0;
}

summary

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

转载于:https://www.cnblogs.com/lxyyyy/p/11439085.html