CodeForces1043 F. Make It One(容斥)

题意：

给定长度为n的数组a，问你最少选择其中多少个数
满足这些数的gcd为1
如果无解输出-1

数据范围：n<=3e5，a(i)<=3e5

解法：

因为a(i)<=3e5
2×3×5×7×11×13×17>3e5
因为最小的7个质因子超过3e5
所以每个a(i)的质因子数量不超过7个,
那么答案最大为7
因此枚举选择的数的个数即可
关键在于如何check

设d[i][j]表示选择i个数,gcd为j的方案数
d[i][j]=C(cnt[j],i)-sigma[d[i][k]],其中k是j的倍数
原理是选择从j的倍数中选择i个出来的方案数
然后容斥掉gcd为j的倍数的方案,最后判断cnt[i][1]是否不为0即可

因为计算过程中数字很大,所以还要要定一个模数

code：

#include<bits/stdc++.h>
using namespace std;
#define int long long
int ppow(int a,int b,int mod){
    int ans=1%mod;a%=mod;
    while(b){
        if(b&1)ans=ans*a%mod;
        a=a*a%mod;
        b>>=1;
    }
    return ans;
}
//
const int maxm=3e5+5;
const int mod=1e9+7;
int fac[maxm],inv[maxm];
void init(){
    fac[0]=1;
    for(int i=1;i<maxm;i++){
        fac[i]=fac[i-1]*i%mod;
    }
    inv[maxm-1]=ppow(fac[maxm-1],mod-2,mod);
    for(int i=maxm-2;i>=0;i--){
        inv[i]=inv[i+1]*(i+1)%mod;
    }
}
int C(int n,int m){
    if(m<0||m>n)return 0;
    return fac[n]*inv[m]%mod*inv[n-m]%mod;
}
//
int a[maxm];
int cnt[maxm];
int d[10][maxm];
signed main(){
    init();
    int n;cin>>n;
    for(int i=1;i<=n;i++){
        cin>>a[i];cnt[a[i]]++;
    }
    for(int i=1;i<maxm;i++){
        for(int j=i+i;j<maxm;j+=i){
            cnt[i]+=cnt[j];
        }
    }
    for(int i=1;i<10;i++){
        for(int j=maxm-1;j>=1;j--){
            d[i][j]=C(cnt[j],i);
            for(int k=j+j;k<maxm;k+=j){
                d[i][j]=(d[i][j]-d[i][k]+mod)%mod;
            }
        }
        if(d[i][1]){
            cout<<i<<endl;
            return 0;
        }
    }
    cout<<-1<<endl;
    return 0;
}

---