1. 程式人生 > >51nod 1616 最小集合(枚舉倍數)

51nod 1616 最小集合(枚舉倍數)

ima 個數 return 復雜度 amp ems ios 比較 esp

技術分享

分析:也就是取任意多個數,它們的最大公約數都在這個集合裏。考慮到ai比較小,可以枚舉小於a中最大值的所有數,判斷是否為其中若幹個數的gcd。記c[k]為a中k的倍數的個數,然後枚舉k的倍數i*k,c[i]<2直接跳過,如果c[i*k]==c[k],說明k的那些倍數也同時是i*k的倍數,k就可以不在集合中,反之,如果任意i,c[i*k]<c[k],說明存在一個倍數和其它k的倍數的gcd是k,所以k一定在集合中。兩次枚舉倍數,復雜度為O(nlogn)。

 1 #include<iostream>
 2 #include<cstring>
 3 #include<cstdio>
 4
using namespace std; 5 const int maxn=1e6+5; 6 bool In_set[maxn]; 7 int count_multiple[maxn]; 8 int main(){ 9 int n,ans=0,a; 10 scanf("%d",&n); 11 ans=n; 12 memset(In_set,0,sizeof(In_set)); 13 memset(count_multiple,0,sizeof(count_multiple)); 14 for(int i=0;i<n;i++){ 15
scanf("%d",&a); 16 if(In_set[a])ans--; 17 In_set[a]=true; 18 } 19 for(int i=1;i<maxn;i++){ 20 //if(In_set[i])continue; 21 for(int j=i;j<maxn;j+=i){ 22 if(In_set[j]) 23 count_multiple[i]++; 24 } 25 } 26 for
(int i=1;i<maxn;i++){ 27 if(In_set[i]||count_multiple[i]<2)continue; 28 bool ok=true; 29 for(int j=2*i;j<maxn;j+=i){ 30 if(count_multiple[j]==count_multiple[i]){ 31 ok=false;break; 32 } 33 } 34 if(ok) 35 ans++; 36 } 37 cout<<ans<<endl; 38 return 0; 39 }

51nod 1616 最小集合(枚舉倍數)