洛谷 CF1043F Make It One 解题报告

CF1043F Make It One

题意

从一堆数中选择最少的数，使它们的(gcd=1)

输入输出格式

输入格式

第一行：一个正整数(n)。
第二行：(n)个正整数，给出了这个数列。

输出格式

一行，(-1)（如果任意选择都不能得到(1)分）或一个正整数（表示选择的数的数量的最小值）

数据范围

(1leq nleq 300,000 , 1leq a_i leq 300,000).

---

第二次遇到此类题了，反演理解的不好，不清楚怎么用反演理解这个题。

不过容斥的方法还挺神的（虽然可能感觉有点套路？

发现选择的数不会超过(7)个，因为(2 imes 3 imes 5 imes 7 imes 11 imes 13 imes 17 >3 imes 10^5)

直接暴力枚举选几个数

然后设(dp_{i,j})表示选(i)了(i)个数且这(i)个数(gcd=j)的方案数

(dp_{i,j}=inom{cnt_j}{i}-sum_{j|k}f_{i,k})

(cnt_j)表示有多少个数是(j)的倍数

---

Code:

#include <cstdio>
#include <cstring>
#define ll long long
const int N=3e5+10;
const ll mod=1e9+7;
ll qp(ll d,ll k)
{
    ll f=1;
    while(k)
    {
        if(k&1) f=f*d%mod;
        d=d*d%mod;
        k>>=1;
    }
    return f;
}
ll fac[N],inv[N],dp[N];
int n,cnt[N],mx;
int main()
{
    scanf("%d",&n);
    fac[0]=1;
    for(int i=1;i<=n;i++) fac[i]=fac[i-1]*i%mod;
    inv[n]=qp(fac[n],mod-2);
    for(int i=n-1;~i;i--) inv[i]=inv[i+1]*(i+1)%mod;
    for(int a,i=1;i<=n;i++) scanf("%d",&a),++cnt[a],mx=mx>a?mx:a;
    for(int i=1;i<=mx;i++)
        for(int j=i<<1;j<=mx;j+=i)
            cnt[i]+=cnt[j];
    for(int i=1;i<=7;i++)
    {
        memset(dp,0,sizeof(dp));
        for(int j=mx;j;j--)
        {
            dp[j]=cnt[j]>=i?fac[cnt[j]]*inv[cnt[j]-i]%mod*inv[i]%mod:0;
            if(dp[j])
                for(int k=j<<1;k<=mx;k+=j)
                    (dp[j]-=dp[k])%=mod;
            (dp[j]+=mod)%=mod;
        }
        if(dp[1]) return printf("%d
",i),0;
    }
    puts("-1");
    return 0;
}

---

2018.11.5