链接:https://www.nowcoder.com/acm/contest/90/F
来源:牛客网
题目描述
给定n,求1/x + 1/y = 1/n (x<=y)的解数。(x、y、n均为正整数)
输入描述:
在第一行输入一个正整数T。 接下来有T行,每行输入一个正整数n,请求出符合该方程要求的解数。 (1<=n<=1e9)
输出描述:
输出符合该方程要求的解数。
输入
3 1 20180101 1000000000
输出
1 5 181
解 :易得 x=n*y/(y-n) 设 y=y1+n; 则 x=n*n/y1 +n 因为 n,x,y1都为整数则y1是n*n的因子
那么问题变成了求n^2的因子个数,比如 n=4 ,n*n=16;
y1 y x
1 5 20
2 6 12
4 8 8
8 12 6
16 20 5
显然 答案为ans=((n^2)的因子数+1)/2;
设 p1^a1 *p2^a2 * ..*pn ^ an=n
那么 (a1+1)*(a2+1)*..*(an+1)为n的因子数
显然 n*n的因子数为 (2*a1+1)*(2*a2+1)*..*(2*an+1)
#include <bits/stdc++.h>
using namespace std;
int main()
{
int T,n;
scanf("%d",&T);
while(T--)
{
scanf("%d",&n);
int res=1;
for(int i=2;i*i<=n;i++)
if(n%i==0)
{
int cnt=0;
while(n%i==0)
{
cnt++;
n/=i;
}
res*=(cnt*2+1);
}
if(n>1) res*=3;
printf("%d\n",(res+1)/2);
}
return 0;
}