求两个超级大数的最大公约数

由于给出的数太大了，所以我们将两个数A，B拆成了N个数相乘和M个数相乘的形式。N,M<=1000，拆成的数<=1000000000。是不是够大？

最终的结果最多保留9位输出。

 

例如：

3
358572 83391967 82
3
50229961 1091444 8863

输出为：

000012028

 

根据欧拉公式，我们可以将任何一个数表示成如下形式：

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)。

 

对于任何一个整数，我们可以快速的将它因式分解

scanf("%d",&a);
for(j=2;j*j<=a;j++)
{
    while(a%j==0)
        a/=j,A[j]++;
}
if(a>1)
    A[a]++;

对于这样由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;
}