HDU 1299 Diophantus of Alexandria （数论+思维）

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=1299

题目大意

要求求解满足1/x+1/y=1/n的
二元组(x,y)解的个数。

题目分析 

令y=n+k,那么式子变为

我们只需要确定k的个数就可以得到解的个数，
当然要去个重。问题转化为n*n的因子个数，
这个个数有一半是重复的，当然n和n这一对不算。

#include<bits/stdc++.h>
using namespace std;

#define debug puts("YES");
#define rep(x,y,z) for(int (x)=(y);(x)<(z);(x)++)
#define ll long long

#define lrt int l,int r,int rt
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define root l,r,rt
#define mst(a,b) memset((a),(b),sizeof(a))
#define pii pair<int,int>
#define fi first
#define se second
#define mk(x,y) make_pair(x,y)
int mod=20100501;
const int maxn=5e3+10;
const int ub=1e6;
ll powmod(ll x,ll y){ll t; for(t=1;y;y>>=1,x=x*x%mod) if(y&1) t=t*x%mod; return t;}
ll gcd(ll x,ll y){return y?gcd(y,x%y):x;}
/*
题目大意：
要求求解满足1/x+1/y=1/n的
二元组(x,y)解的个数。

题目分析：
(xy)/(x+y)=n。
xy=(x+y)*n，
令y=n+k,那么式子变为
x=(n*n)/k+n,
y=n+k,我们只需要确定k的个数就可以得到解的个数，
当然要去个重。问题转化为n*n的因子个数，
这个个数有一半是重复的，当然n和n这一对不算。
*/
int n;
ll ans;
int main(){
    int t;scanf("%d",&t);
    for(int i=0;i<t;i++){
        scanf("%lld",&n);
        ans=1;
        for(int j=2;1LL*j*j<=n;j++) if(n%j==0){
            int cnt=0;
            while(n%j==0) n/=j,cnt++;
            ans=ans*(2*cnt+1);
        }
        if(n>1) ans*=3;
        printf("Scenario #%d:\n",i+1);
        printf("%lld\n\n",(ans+1)/2);
    }
    return 0;
}