C. Reducing Fractions

来源：互联网发布：淘宝软件编辑：程序博客网时间：2024/06/01 21:56

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

比较好的素因子分解的方法，就是每个数记录该数的最大素因子，这样只要预处理后每个数只要最多十几次运算就可以将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;}