Enlarge GCD CodeForces - 1047C 线性筛加合数分解

题目大意：
给你n个数，他们有一个最大公约数，问你删掉最少几个数，可以让他们的最大公约数变大。
分析：
首先想到的是分解者n个数的因子，然后对于2到maxn的数字，选择一个出现次数最多的因子，n-这个数目就是答案，但是分解因子的复杂度是n*sqrt（a[i]）,这样会超时，所以换一种求因子的方法，对于每一个数，找他的最小质数因子，每次除它，统计因子

#include <iostream>
#include <cstring>
#include <cstdio>
#include <queue>
#include <set>
#include <cmath>

using namespace std;
const int maxn = 1e7+5e6+10;
int a[maxn];
int b[maxn],prime[maxn],div1[maxn];
void init()
{
   for(int i=2;i<maxn;i++)
   {
       if(!prime[i])
       {
           prime[++prime[0]]=i;
           div1[i]=i;
       }
       for(int j=1;j<=prime[0]&&prime[j]<=maxn/i;j++)
       {
           prime[prime[j]*i]=1;
           div1[prime[j]*i]=prime[j];
           if(i%prime[j]==0) break;
       }
   }
}
//void solve(int x)
//{
//    for(int i=1;i<=sqrt(x);i++)
//    {
//        if(x%i==0)
//        {
//            a[i]++;
//            if(x/i!=i)
//            a[x/i]++;
//        }
//    }
//}
void solve(int x)
{
    if(x==1||prime[x]==0)
    {
        a[x]++;
    }
    else
    {
       // a[x]++;
        while(x>1)
        {
            int t=div1[x];
            //cout<<t<<endl;
            a[t]++;
            while(x%t==0&&x!=1)
            {
                x/=t;
            }
        }
    }
}
int gcd(int a,int b)
{
    if(b==0) return a;
    return gcd(b,a%b);
}
int main()
{
    init();
    int n;
    scanf("%d",&n);
    int g;int x;
    for(int i=0;i<n;i++)
    {
      scanf("%d",&b[i]);
      if(i==0) g=b[i];
      else
      g=gcd(b[i],g);
      //cout<<g<<endl;
      //solve();
    }
    for(int i=0;i<n;i++)
    {
        solve(b[i]/g);
    }
    int ans=maxn;
   // cout<<g<<endl;
    for(int i=2;i<maxn;i++)
    {
        if(a[i]) ans=min(ans,n-a[i]);
        //cout<<n-a[i]<<endl;
    }
    if(ans==maxn) printf("-1\n");
    else         printf("%d\n",ans);


    return 0;
}