题目大意:给出n 求1/x + 1/y = 1/n满足x <= y这样的对数.
从题目可以看出x>=n, y>=n
设y = n + k 带入原式化简移项得到x=n*n/k + n
题目就变成了求n*n/k结果为整数的个数,也就是求n*n的因子个数.
如果一个数字 n = p1^r1 * p2^r2 * ... pk^rk ,那么 n*n = p1^r1 * p2^r2 * ... pk^rk * p1^r1 * p2^r2 * ... pk^rk ,
它的因子的个数就是 (2*r1+1)*(2*r2+1)*...*(2*rk+1).(+1是因为1是每个数的因子)
结果要除以2因为存在x和y相同的情况,取一种就可以了.
#include <memory.h>
#include <cstdio>
#include <cmath>
const int MAX = 32000;
int prime[3500], prime_idx;
bool is_prime[MAX];
void init_prime(){
memset(is_prime + 2, true, sizeof(is_prime));
for(int i = 2; i < MAX; ++i){
if(is_prime[i]){
prime[prime_idx++] = i;
for(int j = i + i; j < MAX; j += i){
is_prime[j] = false;
}
}
}
}
int main(){
init_prime();
int T, caseno = 0;
scanf("%d", &T);
while(T--){
int n;
long long ans = 1, boundary;
scanf("%d", &n);
boundary = sqrt(n);
for(int i = 0; i < prime_idx && prime[i] <= boundary + 1; ++i){
int expo = 0;
while(n % prime[i] == 0){
++expo;
n /= prime[i];
}
ans *= expo * 2 + 1;
}
if(n > 1)ans *= 3;
printf("Scenario #%d:\n%I64d\n\n", ++caseno, (ans + 1) / 2);
}
return 0;
}