1. 基本知识
GCD、欧拉定理、因式分解
最方便因式分解
scanf("%d",&a);
for(j=2;j*j<=a;j++)
{
while(a%j==0)
a/=j,A[j]++;
}
if(a>1)
A[a]++;
根据欧拉公式,我们可以将任何一个数表示成如下形式:
n=p1^x1*p2^x2*p3^x3.......*pm^xm;
如果将A和B分别表示成
A = p1^a1 * p2^a2 * … * pn^an
B = p1^b1 * p2^b2 * … * pn^bn
其中p1,p2....都是素数,a1,a2...都是整数
那么
GCD( A, B ) = p1^min(a1,b1) * p2^min(a2,b2) * … * pn^min(an,bn)对于这样由N个数乘积组成的数,我们可以将N个数的拆分结果全部保存在同一个A数组中,最终的结果就是超大数自身的拆分。
还有一个需要注意的是,给出的a<=1000000000,直接开数组内存上和时间上都很吃力,所以推荐用map<int,int>容器。
#include<iostream>
#include<string>
#include<map>
using namespace std;
#define mod 1000000000
map<int,int>A,B;
int n,m,a;
__int64 ans;
int min(int a,int b)
{
return a<b?a:b;
}
int main()
{
int i,j;
freopen("D:\\in.txt","r",stdin);
while(scanf("%d",&n)==1)
{
A.clear();B.clear();
for(i=0;i<n;i++)
{
scanf("%d",&a);
for(j=2;j*j<=a;j++)
{
while(a%j==0)
a/=j,A[j]++;
}
if(a>1)
A[a]++;
}
scanf("%d",&m);
for(i=0;i<m;i++)
{
scanf("%d",&a);
for(j=2;j*j<=a;j++)
{
while(a%j==0)
a/=j,B[j]++;
}
if(a>1)
B[a]++;
}
ans=1;
int flag=0;
map<int,int>::iterator it;
for(it=A.begin();it!=A.end();it++)
{
if(B.count(it->first))
{
int p=it->first;
int r=min(A[p],B[p]);
for(i=0;i<r;i++)
{
ans*=p;
if(ans>=mod)
{
ans%=mod;
flag=1;
}
}
}
}
if(flag)
printf("%09d\n",(int)ans);
else
printf("%I64d\n",ans);
}
return 0;
}