C. Reducing Fractions

该题的思路来自：http://blog.csdn.net/acm_ted/article/details/7966780

C. Reducing Fractions

比较好的素因子分解的方法，就是每个数记录该数的最大素因子，这样只要预处理后每个数只要最多十几次运算就可以将10^7以内的数素因子分解了

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>

using namespace std;
const int maxn=10000010;
int prime[maxn],a[100001],b[100001],n,m;
int up[maxn],down[maxn];

void init()
{
    for(int i=2;i<maxn;i++)
      if(!prime[i]){
         prime[i]=i;
         for(int j=2*i;j<maxn;j+=i) prime[j]=i;
      }
}

int main()
{
    init();
    while(scanf("%d%d",&n,&m)==2)
    {
        for(int i=0;i<n;i++) scanf("%d",&a[i]);
        for(int i=0;i<m;i++) scanf("%d",&b[i]);
        memset(up,0,sizeof(up));
        memset(down,0,sizeof(down));
        for(int i=0;i<n;i++)
        for(int j=a[i];j>1;j/=prime[j])  up[prime[j]]++;

        for(int i=0;i<m;i++)
        for(int j=b[i];j>1;j/=prime[j])  down[prime[j]]++;

        //cout<<up[2]<<" "<<up[5]<<endl;
        //cout<<down[2]<<" "<<down[5]<<endl;
        printf("%d %d\n",n,m);
        for(int i=0;i<n;i++)
        {
            int tmp=1;
            for(int j=a[i];j>1;j/=prime[j])
               if(down[prime[j]]>0) down[prime[j]]--;
               else tmp*=prime[j];
            if(i) putchar(' ');
            printf("%d",tmp);
        }
        puts("");
        for(int i=0;i<m;i++)
        {
            int tmp=1;
            for(int j=b[i];j>1;j/=prime[j])
              if(up[prime[j]]>0) up[prime[j]]--;
              else tmp*=prime[j];
            if(i) putchar(' ');
            printf("%d",tmp);
        }
        puts("");
    }
    return 0;
}