CodeForces-1205B Weakened Common Divisor

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

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

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

#include
#include
#include
#include
#include
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;im)
		{
			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;i1) g[t++]=g[i];
	}
	for(int i=2;i<100000;i++)
		if(vis[i]) g[t++]=i;
	cnt=t;
	for(int i=0;i

 

你可能感兴趣的:(CodeForces-1205B Weakened Common Divisor)