xdu2017校赛F

Problem F Dogs of Qwordance Senior Backend R&D Engineers
问题描述
那年夏天，锘爷和杰师傅漫步在知春公园的小道上。他们的妻子、孩子牵 着狗在前面嬉戏，二人笑语盈盈，他们不深究一个小的编程问题，而是对整个 Qwordance （四字舞蹈）公司的发展前景加以描绘。这样的场景，想想就觉得好 美，想想就好向往，想想就好激动。然而，他们的狗觉得这非常的无聊，决定自 己去玩。 杰师傅的狗非常挑剔。它希望找到一块面积为 x 的长方形广场，还要求广 场的长和宽都是整数。锘爷的狗不屑地说：“这还不简单，总共有 σ0(x) 种方案 呢。”杰师傅的狗却摇了摇头说：“找到这个广场后，我们要在上面画上水平和垂 直的网格，使得水平、垂直网格之间的间距分别是相同的整数。算了这太难了， 我们去听他们讨论 Qwordance 未来发展大方向吧。”锘爷的狗说：“这也不难啊， 总共有，咦，多少方案啊？”杰师傅的狗鄙视道：“你看你连个方案数都算不出， 我不要和你玩了。像你这样的狗在朝鲜是会被做成狗肉火锅的。你走吧。”锘爷 的狗想去问锘爷的儿子（中关村小学生信息学竞赛冠军），然而怕被小主人嘲笑， 请你帮帮它吧。
输入格式 输入包含多组数据（最多 450 组），请处理到文件结束。
每组数据包含一个整数，即 x 。 对于所有数据有 1 ≤ x ≤ 1012。
输出格式 对于每组数据输出 1 行，包含面积为 x 的广场的方案数。
只要广场的长或 宽不同，或网格划分方案不同，就是不同的方案（见样例解释）。
输入输出样例 输入样例 输出样例
1 4
1 10

关键：令$\zeta (x)$=一个数的质因子的个数，$f(x)$为解的个数，

　　　因为任何一个数可以按照唯一分解定理，分解为$x = \prod\limits_{i = 1}^k {{p_i}^{{a_i}}} $,${p_i}$为质数，

　　　　因此，$\zeta (x) = \prod\limits_{i = 1}^k {({a_i} + 1)}$，

　　　　$f(x) = \sum\limits_{d|n} {\prod\limits_{i = 1}^k {({a_i} + 1)*({c_i} - {a_i} + 1)} } $

　　　　$f(x) = \prod\limits_{i = 1}^k {\sum\limits_{j = 1}^{{a_i}} {(j + 1)({a_i} - j + 1)} } $(乘法原理)

注意：（1）分解质因数时，只需到$\sqrt n $即可，因为大于等于$\sqrt n $的质数最多只有一个，这时候，只需将结果乘4即可。

　　　（2）分解素数的两种方法，埃氏筛法、线性筛法，其中线性筛法的复杂度更低，埃氏筛法可以舍弃了。

　　　（3）时常注意优化复杂度（dalao的教导）

　　　（4）ios::sync_with_stdio(false);可用于加快cin与cout的速度。

线性筛：

 #include<bits/stdc++.h>
 using namespace std;
 typedef long long ll;
 #define N 10000001
 ll prime[N+],k;
 bool is_not_prime[N+]={true,true};
 void sieve(){
     for(ll i=;i<N;i++){
         if(!is_not_prime[i]){
             prime[k++]=i;
         }
         for(ll j=;j<k&&i*prime[j]<N;j++){
             is_not_prime[i*prime[j]]=true;
             if(i%prime[j]==) break;
         }
     }
 }

 #include<bits/stdc++.h>
 using namespace std;
 typedef long long ll;
 #define N 1000002
 int prime[N+],p;
 ll n;
 bool is_prime[N+];
 int index[N+];
 void sieve(){
     p=;
     for(int i=;i<N;i++){
         is_prime[i]=true;
     }
     is_prime[]=is_prime[]=false;
     for(int i=;i<N;i++){
         if(is_prime[i]){
             prime[p++]=i;
             for(int j=*i;j<N;j+=i){
                 is_prime[i]=false;
             }
         }
     }
 }//埃筛
 /*void sieve(){
     is_prime[1]=1;
     for(ll i=2;i<N;++i){
         if(!is_prime[i])prime[p++]=i;
         for(ll j=0;j<p&&prime[j]*i<N;++j){
             is_prime[prime[j]*i]=1;
             if(i%prime[j]==0)break;
         }
     }
 }*/
 //线性筛 

 int main(){
     ios::sync_with_stdio(false);
     sieve();
     while(cin>>n){
         ll ans=;
         for(int i=;1ll*prime[i]*prime[i]<=n;i++){
             if(n%prime[i]==){
                 ll index=,part=;
                 while(n%prime[i]==){
                     n/=prime[i];
                     index++;
                 }
                 for(ll j=;j<=index;j++) part+=(j+)*(index-j+);
                 ans*=part;
             }
         }if(n!=) ans*=;
         cout<<ans<<endl;
     }
     return ;
 }