同:http://hao666.info/blog/?p=31
题目
题目详情:
我们定义 F(x)是满足 x mod(a*b) == 0这样的a,b的组数。现在给你一个n,你需要求出 F(n)
输入格式:
多组数据,每组第一行有一个整数n, 0 < n <= 10^11。
输出格式:
每组输出一行,满足条件的(a,b)对数
答题说明:
输入样例
1
2
3
4
输出样例:
1
3
3
6
解释
第一组: (1,1)
第二组: (1,1) (1,2) (2, 1)
第三组: (1,1) (1,3) (3,1)
第四组: (1,1) (1,2) (1, 4) (2, 1) (2,2) (4,1)
解法
很简单:
求得x的所有质因子,并统计每个质因子的个数factors[],factors[i]代表第i个质因子的个数。
那么 multi( (factors[i]+1+1)*(factors[i]+1) /2) 即为所求
代码
#include <iostream>
#include <cmath>
#include <cstdlib>
#include <map>
void getFactors(long long i,std::map<long long,int>&ret){
if(i<2)
return ;
long long s=2;
long long e=(long long)sqrt((double)i);
int x=0;
while(0==i%s){
++ret[s];
i/=s;
}
++s;
while(s<=e){
while(0==i%s){
++ret[s];
i/=s;
}
s+=2;
}
if(i!=1)
++ret[i];
}
int main(){
long long i=0;
while(std::cin>>i){
std::map<long long,int>fs;
getFactors(i,fs);
long long ret=1;
for(std::map<long long,int>::iterator it=fs.begin();it!=fs.end();++it)
ret*=(it->second+2)*(it->second+1)/2;
std::cout<<ret<<std::endl;
}
}