CodeForces 1047 C.Enlarge GCD (gcd+质因子分解+思维)

题意：

给长度为n的数组，假设当前数组的gcd为k
问最少删除多少个数能使得数组的gcd变大

思路:

先求出初始数组的gcd=k
把所有数除以k，然后对每个数质因子分解，记录每种质因子出现的次数
出现次数最多的质因子一定是新gcd的因子，假设出现p次，则出现新gcd最多保留p个数，因此答案就是n-p
如果p为0说明无解

需要素数线性筛

质因子分解的时候如果prime[j]*prime[j]>a[i]的时候就停止，否则会t！
因为质因子分解只需要到根号a[i]就行了！

code:

#include<bits/stdc++.h>
using namespace std;
int gcd(int a,int b){
    return b==0?a:gcd(b,a%b);
}
const int maxm=1.5e7+5;
int mark[maxm];
int prime[maxm],cnt;
void init(){
    for(int i=2;i<maxm;i++){
        if(!mark[i]){
            prime[cnt++]=i;
        }
        for(int j=0;j<cnt;j++){
            if(1ll*prime[j]*i>=maxm)break;
            mark[prime[j]*i]=1;
            if(i%prime[j]==0)break;
        }
    }
}
int a[maxm];
int num[maxm];
signed main(){
    init();
    int n;
    scanf("%d",&n);
    int ans=0;
    for(int i=1;i<=n;i++){
        scanf("%d",&a[i]);
        ans=gcd(ans,a[i]);
    }
    for(int i=1;i<=n;i++){
        a[i]/=ans;
        for(int j=0;j<cnt;j++){
            if(prime[j]*prime[j]>a[i])break;//这一步要计算prime[j]*prime[j],否则超时
            if(a[i]%prime[j]==0){
                num[prime[j]]++;
                while(a[i]%prime[j]==0)a[i]/=prime[j];
            }
        }
        if(a[i]>1)num[a[i]]++;
    }
    ans=0;
    for(int i=2;i<maxm;i++){
        ans=max(ans,num[i]);
    }
    if(ans){
        printf("%d\n",n-ans);
    }else{
        puts("-1");
    }
    return 0;
}