hdu 3758 Factorial Simplification

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

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

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

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

 

  1 #include<iostream>
  2 #include<stdio.h>
  3 #include<algorithm>
  4 #include<iomanip>
  5 #include<cmath>
  6 #include<string>
  7 #include<cstdlib>
  8 #include<vector>
  9 using namespace std;
 10 int prime[1300],cnt,p[1300],q[1300],num[1300],len;
 11 int s[1300],r[1300],a[1300];
 12 bool f[10010];
 13 int min(int a,int b)
 14 {
 15     return a>b?b:a;
 16 }
 17 int max(int a,int b)
 18 {
 19     return a>b?a:b;
 20 }
 21 void init()
 22 {
 23     int i,j;
 24     cnt=0;
 25     for(i=2;i<=10010;i++)
 26     {
 27         if(f[i]==0) prime[cnt++]=i;
 28         for(j=0;j<cnt&&prime[j]*i<=10010;j++)
 29         {
 30             f[i*prime[j]]=1;
 31             if(i%prime[j]==0) break;
 32         }
 33     }
 34 }
 35 int count(int n,int p)//计算n!中p因子的个数
 36 {
 37     int sum=0;
 38     while(n>=p)
 39     {
 40         sum+=(n/=p);
 41     }
 42     return sum;
 43 }
 44 int divs(int n,int m)//分解因子
 45 {
 46     int i,j;
 47     len=0;
 48     memset(num,0,sizeof(num));
 49     for(i=0;i<n;i++)
 50     {
 51         for(j=0;prime[j]<=p[i];j++)
 52             num[j]+=count(p[i],prime[j]);
 53         len=max(len,j);
 54     }
 55     for(i=0;i<m;i++)
 56     {
 57         for(j=0;prime[j]<=q[i];j++)
 58         {
 59             if(j>=len) return -1;
 60             num[j]-=count(q[i],prime[j]);
 61             if(num[j]<0) return -1;
 62         }
 63     }
 64     if(num[0]==0) return 0;
 65     return 1;
 66 }
 67 int main()
 68 {
 69     init();
 70     int i,j,k,ans,n,m,temp,t,cur,mul;
 71     bool flag;
 72     cin>>t;
 73     while(t--)
 74     {
 75         cin>>n>>m;
 76         for(i=0;i<n;i++) cin>>p[i];
 77         for(i=0;i<m;i++) cin>>q[i];
 78         ans=divs(n,m);
 79         if(ans<=0)
 80         {
 81             cout<<ans<<endl;
 82             continue;
 83         }
 84         for(i=0;i<len&&num[i]>0;i++) ;
 85         len=min(len,i);
 86         cur=prime[i]-1;
 87         cnt=0;
 88         while(num[0]>0&&cur>=2)
 89         {
 90             flag=true;
 91             a[0]=count(cur,prime[0]);
 92             if(a[0]>num[0]) flag=false;
 93             else mul=num[0]/a[0];
 94             for(i=1;flag&&i<len&&prime[i]<=cur;i++)
 95             {
 96                 a[i]=count(cur,prime[i]);
 97                 if(a[i]>num[i]) flag=false;
 98                 else mul=min(mul,num[i]/a[i]);
 99             }
100             if(flag && mul>0)
101             {
102                 if(i<len) len=i;
103                 r[cnt]=cur;
104                 s[cnt++]=mul;
105                 for(i=0;i<len;i++)
106                 {
107                     num[i]-=(mul*a[i]);
108                     if(num[i]==0) break;
109                 }
110                 if(len>i) len=i,cur=prime[len]-1;
111                 else cur--;
112             }
113             else cur--;
114         }
115         cout<<cnt<<endl;
116         for(i=0;i<cnt;i++)
117             cout<<r[i]<<' '<<s[i]<<endl;
118     }
119     return 0;
120 }
View Code

 

 

 

posted @ 2013-07-22 15:14  _随心所欲_  阅读(239)  评论(0)  编辑  收藏  举报