Codeforces Round #511 (Div. 2) C

这是我在比赛的时候写的代码，好复杂啊。然后看了一下网上其他人，思路是一模一样的，但感觉对代码的精炼能力还是有差距

题目大意：给出n个数，他们有一个最大公因子x，现在可以剔除n个数中的k个数，求一个最小的k，使得剔除之后的（n-k）个数的最大公因子大于x。

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

两份代码见下。

总结一下写法上的不同之处

我写的其实就是欧拉函数的写法，因为每个数都是可以分解成质因子相乘而得，于是就先求出所有质因子，然后每个数一个个遍历，将其分解，如果有这个质因子，那么那个质因子就+1.最后取数量最大的那个质因子即可。

网上写法：更为简洁，直接与素数筛一起结合了。他是开了一个10^7的数组（当时真没意识到codeforce上给的内存竟然都这么大），直接记录这些数。然后进行素数筛，如果有这个因子，那前面进行的时候自然会+1.

我写的：

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<set>
#include<map>
#include<vector>
#include<iterator>
using namespace std;
int s[300009],prm[4000];
int cnt;
map<int,int> mp;
map<int,int>::iterator it;
int gcd(int a,int b)
{
	if(b==0) return a;
	else return gcd(b,a%b);
}
void isprime()
{
	bool prime[4000];
	memset(prime,true,sizeof(prime));
	for(int i=2;i<3880;i++)
	{
		//cout<<i<<" ";
		if(prime[i]){
			prm[cnt++]=i;
			mp[i]=0;
			for(int j=i+i;j<3880;j+=i)
			{
				prime[j]=false;
			}
		}
	}
}
int main()
{
	int n,t;
	//freopen("t.txt","r",stdin);
	mp.clear();
	cnt=0;
	isprime();
	scanf("%d",&n);
	scanf("%d",&s[0]);
	t=s[0];
	for(int i=1;i<n;i++)
	{
		scanf("%d",&s[i]);
		t=gcd(t,s[i]);
	}
	for(int i=0;i<n;i++)
	{
		s[i]/=t;
	}

	for(int i=0;i<n;i++)
	{
		for(int j=0;j<cnt;j++)
		{
			if(s[i]%prm[j]==0)
			{
				//cout<<s[i]<<" "<<prm[j]<<endl;
				//cout<<prm[j]<<" ";
				mp[prm[j]]++;
				while(s[i]%prm[j]==0)
				{
					//cout<<s[i]<<" ";
					s[i]/=prm[j];
				}
			}
		}
		//cout<<s[i]<<" ";
		if(s[i]>1)
		{
			if(mp.count(s[i])==0)
			{
				mp[s[i]]=1;
			}else
			{
				mp[s[i]]++;
			}
		}
	}
	int ans=-1;
	//cout<<n<<" ";
	for(it=mp.begin();it!=mp.end();it++)
	{
		//cout<<it->first<<" "<<it->second<<endl;
		ans=max(ans,it->second);
	}
	if(ans!=0)
	printf("%d\n",n-ans);
	else
	printf("-1\n");
	return 0;	
} 

网上公布的代码

#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;
}