ccsu 1546 数字游戏

Description

自从小黑上次在和大宝的游戏中取得了信心之后，这次他又想到了一个问题，他要去计算两个数A和B的最大公约数。 
由于这两个数非常的大，于是大宝把数字A拆成了N个数，这N个数相乘的结果是A，把B拆成了M个数，同样这M个数的乘积等于B，小黑迫不及待地想去计算出这两个数的最大公约数，这次你能帮帮他吗？如果结果超过了9位数，输出最后的9位数。 

Input

第一行包含一个正整数N ( 1 <= N <= 1000 )。 
第二行包含N个以空格隔开的小于 1 000 000 000的正整数，它们的乘积是A。 
第三行包含一个正整数M ( 1 <= M <= 1000 )。 
第四行包含M个以空格隔开的小于 1 000 000 000的正整数，它们的乘积是B。 

Output

输出A和B的最大公约数，如果结果超过了9位数，输出后9位。

Sample Input

3 
2 3 5 
2 
4 5 
4 
6 2 3 4 
1 
1 

Sample Output

10 
1

思路：

根据任意一个整数都可以划分为多个素数的积，因为所给数十分大，判断素数时是看它是否有小于sqrt(x)的因数，所以这里只用求到1~~~ sqrt(999999999)的素数，再算出A和B的每一个因数分解为不同值素数的个数，因为为防止其中有因子为大于sqrt(999999999)的素数，所以必须也要储存不能被已经求出的素数表中的数所整除的数，在将A 与 B的这个数组相比较，若有相同，则乘以它。

代码：

#include<iostream>
#include<string.h>
#include<algorithm> 
#include<math.h>
#define N 10000
#define Max 1000000000
using namespace std;
int a[6010] , b[6010] ;
int c[6010], d[6010];
int prime[6010]={2} ,top;
int isprime[N+1]={0};

void prim()
{
     int i,j;
     for (i=3;i<N;i+=2)
     {
         if (!isprime[i])
         {
            prime[top++]=i;
            for (j=i*i;j<N;j+=i) isprime[j]=1;                
         }    
     }    
}
int main()
{
    int n , m ,flag;
    int x, y ,t; 
    __int64 ans;   
    top=1;
    prim();
    while(scanf("%d",&n)!=EOF)
    {
        int lc=-1, ld=-1;
        ans=1;
        flag=0;
        memset(a,0,sizeof(a));
        memset(b,0,sizeof(b));
        memset(c,0,sizeof(c));
        memset(d,0,sizeof(d)); 
        for(int i=0;i<n;i++)
        {
            scanf("%d",&t);
            for(int j=0;j<top && t>1;j++)
            {//  system("pause");
                while(t%prime[j]==0)
                {
                    a[j]++; // cout<<"a["<<j<<"]= "<<a[j]<<endl;
                    t/=prime[j];
                }
            } 
            if(t>1) c[++lc]=t;
        }
        scanf("%d",&m) ;
        for(int i=0;i<m;i++)
        {
            scanf("%d",&t); 
            for(int j=0;j<top && t>1;j++)
            {
                while(t%prime[j]==0)
                {
                    b[j]++;  //cout<<"b["<<j<<"]= "<<b[j]<<endl;
                    t/=prime[j];
                }
            }
           if(t>1)  d[++ld]=t;
        }  
        for(int i=0;i<top;i++)
        {
             if(a[i]>0 && b[i]>0)
             {
                 for(int j=0;j<min(a[i], b[i]);j++)
                 {
                     ans*=prime[i];
                     if(ans>=Max) flag=1;
                     ans=ans%Max;
                 }
             }
        }//cout<<"&&7"<<endl;
       // cout<<"***********"<<endl;
      
        sort(c,c+lc+1);
        sort(d,d+ld+1);
        x=y=0;
        while(x<=lc && y<=ld)
        {//cout<<"c["<<x<<"]= "<<c[x]<<"   d["<<y<<"]= "<<d[y]<<endl;
            if(c[x]==d[y])
            {
                ans*=c[x];   // cout<<"ans= "<<ans<<endl;
                if(ans>=Max) flag=1;
                ans%=Max;
                x++; y++;
            }
            else if(c[x]<d[y]) x++;
            else if(c[x]>d[y])  y++; 
        } 
        int c=-1; 
        if(flag==1)
        {
             printf("%09I64d\n",ans);
        }
        else printf("%I64d\n",ans);
        
    }
    return 0;
}