hdu 3758 Factorial Simplification

这题主要是质因数分解！！

求出每个因子的幂，如果有负数，则输出-1；

如果2的幂数为0，这输出0；

最后就是开始凑阶乘了……

#include<iostream>
#include<stdio.h>
#include<algorithm>
#include<iomanip>
#include<cmath>
#include<string>
#include<cstdlib>
#include<vector>
using namespace std;
int prime[1300],cnt,p[1300],q[1300],num[1300],len;
int s[1300],r[1300],a[1300];
bool f[10010];
int min(int a,int b)
{
    return a>b?b:a;
}
int max(int a,int b)
{
    return a>b?a:b;
}
void init()
{
    int i,j;
    cnt=0;
    for(i=2;i<=10010;i++)
    {
        if(f[i]==0) prime[cnt++]=i;
        for(j=0;j<cnt&&prime[j]*i<=10010;j++)
        {
            f[i*prime[j]]=1;
            if(i%prime[j]==0) break;
        }
    }
}
int count(int n,int p)//计算n!中p因子的个数
{
    int sum=0;
    while(n>=p)
    {
        sum+=(n/=p);
    }
    return sum;
}
int divs(int n,int m)//分解因子
{
    int i,j;
    len=0;
    memset(num,0,sizeof(num));
    for(i=0;i<n;i++)
    {
        for(j=0;prime[j]<=p[i];j++)
            num[j]+=count(p[i],prime[j]);
        len=max(len,j);
    }
    for(i=0;i<m;i++)
    {
        for(j=0;prime[j]<=q[i];j++)
        {
            if(j>=len) return -1;
            num[j]-=count(q[i],prime[j]);
            if(num[j]<0) return -1;
        }
    }
    if(num[0]==0) return 0;
    return 1;
}
int main()
{
    init();
    int i,j,k,ans,n,m,temp,t,cur,mul;
    bool flag;
    cin>>t;
    while(t--)
    {
        cin>>n>>m;
        for(i=0;i<n;i++) cin>>p[i];
        for(i=0;i<m;i++) cin>>q[i];
        ans=divs(n,m);
        if(ans<=0)
        {
            cout<<ans<<endl;
            continue;
        }
        for(i=0;i<len&&num[i]>0;i++) ;
        len=min(len,i);
        cur=prime[i]-1;
        cnt=0;
        while(num[0]>0&&cur>=2)
        {
            flag=true;
            a[0]=count(cur,prime[0]);
            if(a[0]>num[0]) flag=false;
            else mul=num[0]/a[0];
            for(i=1;flag&&i<len&&prime[i]<=cur;i++)
            {
                a[i]=count(cur,prime[i]);
                if(a[i]>num[i]) flag=false;
                else mul=min(mul,num[i]/a[i]);
            }
            if(flag && mul>0)
            {
                if(i<len) len=i;
                r[cnt]=cur;
                s[cnt++]=mul;
                for(i=0;i<len;i++)
                {
                    num[i]-=(mul*a[i]);
                    if(num[i]==0) break;
                }
                if(len>i) len=i,cur=prime[len]-1;
                else cur--;
            }
            else cur--;
        }
        cout<<cnt<<endl;
        for(i=0;i<cnt;i++)
            cout<<r[i]<<' '<<s[i]<<endl;
    }
    return 0;
}