題干：

F先生有n個正整數(shù)，a1，a2，...，an

他認為這些整數(shù)的最大公約數(shù)太小了,所以他想通過刪除一些整數(shù)來擴大它

您的任務是計算需要刪除的最小整數(shù)數(shù),以便剩余整數(shù)的最大公約數(shù)大于所有整數(shù)的公約數(shù).

Input

3
1 2 4

Output

1

Input

4
6 9 15 30

Output

2

Input

3
1 1 1

Output

-1

題目大意：

n個數(shù)，要你刪掉最少的數(shù)，使得剩下的數(shù)的gcd大于所有n個數(shù)的gcd。
n<=3*10^5, 1<=ai<=1.5*10^7
解題報告：

從gcd入手，考慮gcd的本質(zhì)是什么，不過就是所有n個數(shù)的公因子。所以可以想到可以從枚舉公因子的角度出發(fā)。

把所有的數(shù)都除以它們的gcd，那么這n個新數(shù)的gcd一定為1（常見套路），則問題變?yōu)?#xff1a;刪掉最少的數(shù)，使得它們的gcd不等于1。假設所有數(shù)都等于1，顯然無解。假設某一數(shù)不等于1，則只保留最大的這一個數(shù) 就是一個解， 所以一定有解。

我們先用cnt[i]表示a里面有多少個數(shù)等于i。

我們枚舉保留下來的數(shù)的gcd為p，則可以保留下來的個數(shù)為

可以發(fā)現(xiàn)，因此我們枚舉gcd的時候可以只枚舉質(zhì)數(shù)。可以通過線性篩將所有的質(zhì)數(shù)O(n)時間復雜度篩出來

最后那個枚舉的復雜度分析：

AC代碼：

#include<cstdio> #include<iostream> #include<algorithm> #include<queue> #include<map> #include<vector> #include<set> #include<string> #include<cmath> #include<cstring> #define F first #define S second #define ll long long #define pb push_back #define pm make_pair using namespace std; typedef pair<int,int> PII; const int MAX = 3e5 + 5; const int MAXX = 1e7 + 5e6 + 5; int n,a[MAX],tot,p[MAXX],cnt[MAXX]; bool is[MAXX]; void prime() {memset(is,1,sizeof is);is[1]=is[0]=0;for(int i = 2; i<MAXX; i++) {if(is[i] == 0) continue;p[++tot] = i;for(int j = i+i; j<MAXX; j+=i) {is[j] = 0;}} } int main() {int g=0,ans = 3e5;cin>>n;for(int i = 1; i<=n; i++) scanf("%d",a+i),g = __gcd(g,a[i]);for(int i = 1; i<=n; i++) cnt[a[i]/g]++;if(cnt[1] == n) {cout << "-1" << endl;return 0 ;}prime();for(int i = 1; i<=tot; i++) {int sum = 0;for(int j = p[i]; j<MAXX; j+=p[i]) {sum += cnt[j];}ans = min(ans,n-sum);}cout << ans << endl;return 0 ; }

?

總結(jié)

以上是生活随笔為你收集整理的【CodeForces - 1047C】Enlarge GCD（数学，枚举，预处理打表，思维）的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。