1034A. Enlarge GCD

题目链接：http://codeforces.com/problemset/problem/1034/A
题目大意：给出n个数，他们有一个最大公因子x，现在可以剔除n个数中的k个数，求一个最小的k，使得剔除之后的（n-k）个数的最大公因子大于x。

分析：首先求出n个数的最大公因子x，再把这个最大公因子从这n个数中剔除（n/x），余数y=（n/x），将y所在位置标记cnt++，表示这个位置有cnt个y，和素数筛原理一样，接下来筛数，将i倍数地方的cnt全部相加，表示n个数中有因子i的个数，保存最大个数。

代码：

#include<cstdio>
#include<iostream>
#include<cmath>
#include<algorithm>
#include<cstring>
using namespace std;
typedef long long ll;
const int MAXN = 2e7+7;
int a[MAXN],cnt[MAXN];
bool vis[MAXN];


int main()
{
    int n;
    memset(vis,true,sizeof(vis));
    memset(cnt,0,sizeof(cnt));
    scanf("%d", &n);
    int x,ans;
    scanf("%d", &x);
    a[0] = x;

    for(int i = 1; i < n; i++)//求所有数的最大公约数
    {
        scanf("%d", &a[i]);
        x = __gcd(x,a[i]);
    }
    for(int i = 0; i < n; i++)
    {
        cnt[a[i]/x]++;//缩小x倍
    }
    vis[0] = vis[1] = false;
    ans = 0;
    for(int i = 2; i < MAXN; i++)//筛子,同素数筛原理，找出含有公共因子个数最多的组合，记为ans
    {
        int res = cnt[i];
        if(vis[i])
        {
            for(int j = 2*i; j < MAXN; j+=i)
            {
                vis[j] = false;
                res+=cnt[j];
            }
        }
        ans = max(ans,res);
    }//含有相同单个因子的数最多有ans个，输出需要剔除的
    printf("%d\n", ans==0?-1:n-ans);


    return 0;
}