CodeForces-1205B Weakened Common Divisor

先找出最小的一对数字，把他们所有的因数找出来，保存在数组里

接下来过一遍数字对，每次都把能够整除至少一个数字的因数保留，最后剩余的都可以作为答案

如果因数特别多会超时，所以在开始时尽量对因数分解，分解为质数，能够大大减少因数数量

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long ll;
const int N=150000+10;
ll prime[100010];
void get_prime()
{
	memset(prime,0,sizeof(prime));
	for(ll i=2;i<=100000;i++)
	{
		if(!prime[i]) prime[++prime[0]]=i;
		for(ll j=1;j<=prime[0]&&prime[j]<=100000/i;j++)
		{
			prime[prime[j]*i]=1;
			if(i%prime[j]==0) break;
		}
	}
}
ll a[N],b[N];
ll g[200000];
bool vis[100000];
int cnt;
int main()
{
	int n;
	ll m=0;
	int p=0;
	get_prime();
	scanf("%d",&n);
	for(int i=0;i<n;i++)
	{
		scanf("%I64d%I64d",&a[i],&b[i]);
		if(a[i]+b[i]>m)
		{
			m=a[i]+b[i];
			p=i;
		}
	}
	cnt=0;
	ll e=sqrt(a[p]);
	for(ll i=2;i<=e;i++)
		if(a[p]%i==0)
		{
			g[cnt++]=a[p]/i;
			g[cnt++]=i;
		}
	e=sqrt(b[p]);
	for(ll i=2;i<=e;i++)
		if(b[p]%i==0)
		{
			g[cnt++]=b[p]/i;
			g[cnt++]=i;
		}
	g[cnt++]=a[p];
	if(b[p]!=a[p]) g[cnt++]=b[p];
	sort(g,g+cnt);
	int t=1;
	for(int i=1;i<cnt;i++) if(g[i]!=g[i-1]) g[t++]=g[i];
	cnt=t;
	memset(vis,false,sizeof(vis));
	t=0;
	for(int i=0;i<cnt;i++)
	{
		for(int j=1;j<prime[0];j++)
		{
			if(g[i]%prime[j]==0) vis[prime[j]]=true;
			while(g[i]%prime[j]==0)
				g[i]/=prime[j];
		}
		if(g[i]>1) g[t++]=g[i];
	}
	for(int i=2;i<100000;i++)
		if(vis[i]) g[t++]=i;
	cnt=t;
	for(int i=0;i<n;i++)
	{
		int t=0;
		for(int j=0;j<cnt;j++)
			if(a[i]%g[j]==0||b[i]%g[j]==0)
			g[t++]=g[j];
		cnt=t;
	}
	if(cnt) printf("%I64d\n",g[cnt-1]);
	else printf("-1\n");
	return 0;
}