题目
n(n<=2e3)个数,第i个数ai(1<=ai<=1e9),每次你可以选两个不同的数x,y
并把它们合并为gcd(x,y)或者min(x,y),问最后剩下的值有哪些可能
思路来源
官方题解
题解
注意到gcd(x,y)<=min(x,y),所以最后的值只可能在[1,min(a1,..,an)]里,
考虑这里面的每个值v,n个值里v的倍数的值就直接取gcd,这样显然是最优的,
若干个数k1*v,k2*v,...,ki*v,ki越多,gcd(k1,...ki)越可能为1,
首先,只能选择v的倍数,否则gcd一定不为v;那如果某个子序列的gcd为v,选上其他的v的倍数,gcd仍为v
写法随意,为了卡常写的比较阴间,其实可以考虑1e7以下用数组,1e7以上用map,能卡的更短
代码
#include<bits/stdc++.h>
using namespace std;
const int N=2e3+10,M=2e6+10;
int n,a[N],mn,ans[M],c,res;
bool vis[M];
int main(){
scanf("%d",&n);
for(int i=1;i<=n;++i){
scanf("%d",&a[i]);
}
sort(a+1,a+n+1);
for(int i=1;i<=n;++i){
for(int j=1;j<=sqrt(a[i]);++j){
if(a[i]%j==0){
if(!vis[j])ans[++c]=j,vis[j]=1;
if((a[i]/j)!=j){
if((a[i]/j)>=M)ans[++c]=a[i]/j;
else if(!vis[a[i]/j])ans[++c]=a[i]/j,vis[a[i]/j]=1;
}
}
}
}
sort(ans+1,ans+c+1);
c=unique(ans+1,ans+c+1)-(ans+1);
for(int i=1;i<=c;++i){
if(ans[i]>=a[1])break;
int x=0;
for(int j=1;j<=n;++j){
if(a[j]%ans[i]==0){
x=(!x)?a[j]:__gcd(x,a[j]);
}
}
res+=(x==ans[i]);
}
printf("%d\n",res+1);
return 0;
}