题意:
给定n个数,问最少要去掉几个数,使得剩下的数gcd 大于原来n个数的gcd值。
思路:
这道题可以枚举素数x,对于每个x,找到所有(a【i】/gcd(all)) 是x倍数的个数,就是一个次数。找这个次数的过程正好与素数筛的过程一致。
分解质因数会T
#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <string>
#include <cmath>
#include <queue>
#include <map>
using namespace std;
int n;
int a[300005];
int isp[15000005];
long long primes[10000005];
int cnt;
int ct[15000005];
map<int ,int>mp;
int gcd(int a,int b)
{
if(b==0)
return a;
return gcd(b,a%b);
}
void init()
{
cnt=0;
for(int i=2;i<=15000000;i++)
{
if(isp[i]==0)
primes[cnt++]=i,mp[i]=cnt-1;
for(int j=0;j<cnt;j++)
{
if(1ll*i*primes[j]>15000000)
break;
isp[i*primes[j]]=1;
if(i%primes[j]==0)
break;
}
}
}
int main() {
init();
while(~scanf("%d",&n))
{
for(int i=1;i<=n;i++)
scanf("%d",&a[i]);
int gcdd=a[1];
for(int i=2;i<=n;i++)
gcdd=gcd(a[i],gcdd);
for(int i=1;i<=n;i++)
a[i]/=gcdd;
memset(ct,0,sizeof(ct));
int ans=0;
for(int i=1;i<=n;i++)
{
ct[a[i]]++;
}
for(int i=0;i<cnt;i++)
{
int tct=0;
for(int j=primes[i];j<=15000000;j+=primes[i])
{
tct+=ct[j];
}
ans=max(ans,tct);
}
ans=n-ans;
if(ans==n)
printf("-1\n");
else
printf("%d\n",ans);
}
return 0;
}