Codeforces Round #511 (Div. 1) A. Enlarge GCD（数论）

Codeforces Round #511 (Div. 1) A. Enlarge GCD（数论）

题目大意：要你在一组数组中找最少的元素的个数，使得删去这些个数的元素后，这组数组元素的最大公因数会增大。
题解：设这组数组的原本最大公因数为$D$。考虑让改变后的数组增大最大公因数，可以找出一些元素，这组元素的最大公因数比原本的$D$要大，说明这组元素中均存在着一些除D之外的公因数$d$，而这些公因数$d$就是这组元素全部被$D$整除后的得到的最大公因数，即$gcd(a_i/D,a_{i+1}/D,\cdots )=d$。那么我们只要找出这组数组的元素全部被$D$整除后，有哪些元素相互之间的公因数大于$1$，而数组元素的总量减去此类元素的数量的最大值，就是我们要求的答案。因为删去此类元素之外的元素，这组数组的最大公因数就变成了$Dd$。
下面是AC代码。

#include<bits/stdc++.h>
using namespace std;
struct fastio{
    fastio(){
        ios::sync_with_stdio(0);
        cin.tie(0);
        cout.tie(0);
    }
}fio;
#define endl '\n'
const int N=1e7+5e6+5;
vector<int> prime,part(N,0);
int gcd(int a,int b){
    if(!b) return a;
    return gcd(b,a%b);
}
void init(){
    part[0]=part[1]=-1;
    for(int i=2;i<N;++i){
        if(!part[i]){
            prime.push_back(i);
            part[i]=i;
        }
        for(auto j:prime){
            if(i*j>=N) break;
            part[i*j]=j;
            if(i%j==0) break;
        }
    }
}
unordered_map<int,int> prif;
void factor(int n){
    while(n>1){
        int p=part[n];
        while(part[n]==p) n/=p;
        ++prif[p];
    }
}
int main(){
    init();
    vector<int> data;
    int n,m,t,ans=0x3f3f3f3f,d=0;
    cin>>n;
    t=n;
    while(t--){
        cin>>m;
        data.push_back(m);
        d=gcd(d,m);
    }
    for(auto i:data) factor(i/d);
    for(auto i:prif) ans=min(ans,n-i.second);
    if(ans==0x3f3f3f3f) cout<<-1<<endl;
    else cout<<ans<<endl;
}