当前位置: 代码迷 >> 综合 >> HDU 1299 Diophantus of Alexandria (数论+思维)
  详细解决方案

HDU 1299 Diophantus of Alexandria (数论+思维)

热度:82   发布时间:2023-11-15 12:16:31.0

题目链接: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;
}