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;
}