codeforces 1043 F. Make It One 数学？随机？

https://codeforces.com/contest/1043/problem/F

没搜到中文题解，看了一下官方题解，翻译解释一发。

题意，从n个数中选最少个数的数使得取出的数gcd为1。数据范围是n最大3e5，每个数最大3e5.

思路：我们很容易观察出，如果存在可行答案，那么取出数量不超过7，因为最小的7个质数相乘都大于3e5了。

然后定义dp[ i ][ j ] 的意义为，取出i个不同的数，他们的gcd为 j 。

通过容斥原理可得  

其中cntj表示有多少个能整除 j 的数，整个式子的意思就是 从能 整除 j 的数中选出 i 个，其中有一些的gcd不为 j ，后边一项就是减去不符合的。整个式子从后往前求一下就得出了。

最后答案要判断 dp[ i ][ 1 ] 是否不为0，由于数太大，我们选取一个大质数在运算过程中取模。

虽然有一定几率求错，但概率应该比较小。

#include <algorithm>
#include  <iostream>
#include   <cstring>
#include    <string>
#include    <cstdio>
#include    <vector>
#include    <bitset>
#include     <stack>
#include     <queue>
#include     <cmath>
#include       <set>
#include       <map>
using namespace std;
#define rep(i,a,b) for(int i=a;i<=b;i++)
#define per(i,a,b) for(int i=a;i>=b;i--)
#define clr(a,x) memset(a,x,sizeof(a))
#define pb push_back
#define all(x) (x).begin(),(x).end()
#define fi first
#define se second
#define mp make_pair
#define SZ(x) ((int)(x).size())
typedef unsigned long long ull;
typedef long long ll;
typedef vector<int> vi;
typedef pair<int,int> pii;
/*************head******************/
const int MAXN=3e5+10;
const ll MOD=1e9+7;
ll dp[10][MAXN];
ll a[MAXN];
ll cnt[MAXN];
ll fac[MAXN];
ll inv[MAXN];
ll pow_mod(ll a,ll b){
    ll res=1;
    while(b){
        if(b&1)res=res*a%MOD;
        a=a*a%MOD;
        b>>=1;
    }
    return res;
}
ll C(ll n,ll m){
    return fac[m]*inv[n]%MOD*inv[m-n]%MOD;
}
int main(){
    fac[0]=1;
    rep(i,1,300000)fac[i]=fac[i-1]*i%MOD;
    inv[300000]=pow_mod(fac[300000],MOD-2);
    per(i,300000,1)inv[i-1]=inv[i]*i%MOD;
    int n;
    scanf("%d",&n);
    rep(i,1,n){
        ll x;
        scanf("%I64d",&x);
        a[x]++;
    }
    rep(i,1,300000){
        for(int j=i;j<=300000;j+=i){
            cnt[i]+=a[j];
        }
    }
    int ans=-1;
    for(int i=7;i>=1;i--){
        for(int j=300000;j>=1;j--){
            if(cnt[j]<i)continue;
            dp[i][j]+=C(i,cnt[j]);
            if(dp[i][j]>MOD)dp[i][j]-=MOD;
            for(ll k=2*j;k<=300000;k+=j){
                dp[i][j]-=dp[i][k];
                if(dp[i][j]<0)dp[i][j]+=MOD;
            }
        }
        if(dp[i][1]>0)ans=i;
    }
    printf("%d",ans);
    return 0;
}