最大公约数 （数学+离散化）

题目大意：
给定 $n(n\le 3e5)$ 个正整数，$a_1,a_2,\dots ,a_n(a_i\le 1.5e7)$,求最少删去几个数，使得删去后这些数的最大公约数比原先的所有数的最大公约数大。

题目分析：
我们可以先通过把每个数除去这 $n$ 个数的 $gcd$ 来把他们的公因子剔除，我们再对处理后的数离散化一下，然后我们筛出 $\sqrt{ma{x}_{a_i}}$ 以内的质数，用一个桶记录一下这些质数是多少个 $a_i$ 因子，最后找出桶中的最大值 $maxx$
如果 $maxx=0$ 表示全是 $1$ 他们的 $gcd\equiv 1$ 无法通过删数减少其最大公约数，此时无解就输出 $-1$ ，其余情况就是有解情况，当我们把 $n-maxx$ 个数字剔除后我们剩余数字的 $gcd$ 就会增大，增大的值就是桶中数字为 $maxx$ 对应的那个质数
还有一个小细节，如果 $a_i$ 是质数且 $a_i>\sqrt{ma{x}_{a_i}}$ 我们应该将这个没筛到过的质数加入质数表中，避免遗漏

具体细节见代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<set>
#include<map>
#define ll long long
#define inf 0x3f3f3f3f
using namespace std;
int read()
{
	int res = 0,flag = 1;
	char ch = getchar();
	while(ch<'0' || ch>'9')
	{
		if(ch == '-') flag = -1;
		ch = getchar();
	}
	while(ch>='0' && ch<='9')
	{
		res = (res<<3)+(res<<1)+(ch^48);//res*10+ch-'0';
		ch = getchar();
	}
	return res*flag;
}
const int maxn = 1e6+5;
const int mod = 1e9+7;
const double pi = acos(-1);
const double eps = 1e-8;
struct node{
	int val,cnt;
}nod[maxn];
int n,a[maxn],cnt,pri[maxn],gcdd,num,ans[maxn];
bool vis[maxn];
void get_prime(int n)
{
	for(int i = 2;i <= n;i++)
	{
		if(!vis[i]) pri[++cnt] = i;
		for(int j = 1;j <= cnt && i*pri[j] <= n;j++)
		{
			vis[i*pri[j]] = true;
			if(i%pri[j] == 0) break;
		}
	}
}
int main()
{
	n = read();
	int maxx = 0;
	for(int i = 1;i <= n;i++)
	{
		a[i] = read();
		maxx = max(maxx,a[i]);
		gcdd = i==1 ? a[i] : __gcd(gcdd,a[i]);
	}
	get_prime(sqrt(maxx/gcdd));
	sort(a+1,a+n+1);
	for(int i = 1;i <= n;i++)
		if(a[i] != a[i-1])
		{
			nod[++num].val = a[i]/gcdd;
			nod[num].cnt = 1;
		}
		else nod[num].cnt++;
	for(int i = 1;i <= num;i++)
	{
		bool flag = true;
		for(int j = 1;j <= cnt;j++)
			if(nod[i].val%pri[j] == 0)
				ans[j] += nod[i].cnt,flag = false;
		if(flag && nod[i].val != 1)
		{
			pri[++cnt] = nod[i].val;
			ans[cnt] += nod[i].cnt;
		}
	}
	maxx = *max_element(ans+1,ans+cnt+1);
	if(!maxx) puts("-1");
	else printf("%d\n",n-maxx);
	return 0;
}