洛谷P5253 丢番图

题目大意

求$\frac{1}{x}+\frac{1}{y}=\frac{1}{n}$的解的个数

题目分析

其实就是化(luan)简(gao)
$\frac{1}{x}+\frac{1}{y}=\frac{1}{n}$
$\frac{y}{xy}+\frac{x}{xy}=\frac{1}{n}$
两边同乘xy得
$x+y=\frac{xy}{n}$
再同乘 $-n$得
$-n(x+y)=-xy$
同时加上$n^2$得
$n^2-n(x+y)=n^2-xy$
去括号得
$n^2-nx-ny=n^2-xy$
移项得
$n^2-nx-ny+xy=n^2$
根据多项式的乘法化简左边的式子得
$(x-n)(y-n)=n^2$
化(luan)简(gao)完毕
  对于每一个可能的(x-n)，都必定会有一个(y-n)与其对应；而(x-n)是$n^2$的约数，去除重合的情况，问题转化成求$n^2$的约数的个数向上取整。
  思考到这里，开心的你看了看题目数据，然后发现用(sang)心(xin)良(bing)苦(kuang)的出题人把数据出到了$10^{14}$。于是，O(n)求约数个数的办法被无情地否决了。
  那我们应该怎么办呢？这时我们应该祭出一个定理——“唯一分解定理”！！！那么它是什么呢？

$任意大于1的自然数都能被质数的指数幂的乘积表示$

然并卵……
好了其实我们用到的其实是它的推论，准确地来说是

• $若可以将n分解为p{1}^{c1}\ast p{2}^{c2}\ast \dots \ast p{k}^{ck}(n>1)$
• $则n的正约数个数为(c1+1)\ast (c2+1)\ast \dots \ast (ck+1)$

它的原理就不在这里介绍了，我认为度娘会解释得很好的。

Code

#include<iostream>
#include<cstdio>
using namespace std;
long long ans=1,n;
int main(){
	scanf("%lld",&n);//输入n
	for(long long i=2,cnt=0;i*i<=n;i++){//唯一分解定理的推论
		if(n%i==0){
			long long cnt=0;
			while(n%i==0)n/=i,++cnt;//求c的值
			ans*=cnt*2+1;//直接乘上去
		}
	}
	if(n>1) ans*=3;//假如n没有被除完，ans*=1*2+1,这里的代码是化简后的
	ans=(ans+1)/2;//ans除以2向上取整
	printf("%lld",ans);//输出
	return 0;
}