题意:
给长度为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;
}