light oj 1236 - Pairs Forming LCM & uva 12546 - LCM Pair Sum

第一题给定一个大数，分解质因数，每个质因子的个数为e1,e2,e3,……em，

则结果为((1+2*e1)*(1+2*e2)……(1+2*em)+1)/2.

代码如下：

 #include <stdio.h>
 #include <string.h>
 #include <algorithm>
 #include <iostream>
 #define M 10000005
 #define mod 1000000007
 #define ll unsigned long long
 using namespace std;
 ll prime[M/],cnt;
 bool f[M];
 void init()
 {
     int i,j;
     cnt=;
     for(i=;i<M;i++){
         if(!f[i]) prime[cnt++]=i;
         for(j=;j<cnt&&i*prime[j]<M;j++){
             f[i*prime[j]]=;
             if(i%prime[j]==) break;
         }
     }
 }
 ll solve(ll n)
 {
     ll ans=;
     for(ll i=;i<cnt&&(ll)prime[i]*prime[i]<=n;i++){
         if(n%prime[i]==){
             ll t=;
             n/=prime[i];
             while(n%prime[i]==){
                 t++;
                 n/=prime[i];
             }
             ans*=(ll)(+t*);
         }
     }
     if(n>) ans=*ans;
     return (ans+)/;
 }
 int main()
 {
     int t,ca=;
     ll n;
     init();
     scanf("%d",&t);
     while(t--){
         scanf("%llu",&n);
         printf("Case %d: %llu\n",++ca,solve(n));
     }
     return ;
 }

第二题给定一个数N，其质因子pi及个数ei给定，

则结果为(1+p1+p1^2+……+p1^ei+ei*p1^ei)*……*(1+pm+pm^2+……+pm^em+em*pm^em)+N.

代码如下：

 #include <stdio.h>
 #include <string.h>
 #include <algorithm>
 #include <iostream>
 #define mod 1000000007
 #define ll long long
 using namespace std;
 ll Pow(ll a,int b)
 {
     ll ans=;
     while(b){
         if(b&) ans=(ans*a)%mod;
         b>>=;
         a=(a*a)%mod;
     }
     return ans;
 }
 ll sum(ll a,int b)
 {
     if(b==) return ;
     if(b&) return (+Pow(a,(b+)/))*sum(a,b/)%mod;
     else return ((+Pow(a,b/+))*sum(a,b/-)%mod+Pow(a,b/))%mod;
 }
 int main()
 {
     int t,ca=,n,p,e;
     scanf("%d",&t);
     while(t--){
         scanf("%d",&n);
         ll ans=,xx,num=;
         for(int i=;i<n;i++){
             scanf("%d%d",&p,&e);
             xx=Pow(p,e);
             num=(num*xx)%mod;
             xx=(xx*e)%mod;
             xx=(xx+sum(p,e))%mod;
             ans=(ans*xx)%mod;
         }
         printf("Case %d: %lld\n",++ca,(ans+num)%mod);
     }
     return ;
 }