3158. 【JSOI2013】丢番图

Problem

题目传送门

题目描述

丢番图是亚历山大时期埃及著名的数学家。他是最早研究整数系数不定方程的数学家之一。

为了纪念他，这些方程一般被称作丢番图方程。最著名的丢番图方程之一是 $x^n$+$y^n$=$z^n$。费马提出，对于n≥2，$x^n$+$y^n$=$z^n$没有正整数解。这被称为“费马大定理”，它的证明直到最近才被安德鲁·怀尔斯（Andrew Wiles）证明。

考虑如下的丢番图方程：$\frac{1}{x}+\frac{1}{y}=\frac{1}{n}$ $(x,y,n\in N^{+})$

小G 对下面这个问题十分感兴趣：对于一个给定的正整数n ，有多少种本质不同的解满足该方程？例如n=4，有三种本质不同(x≤y)的解：

$\frac{1}{5}+\frac{1}{20}=\frac{1}{4}$

$\frac{1}{6}+\frac{1}{12}=\frac{1}{4}$

$\frac{1}{8}+\frac{1}{8}=\frac{1}{4}$
  ~  

样例输入

4

样例输出

3

数据范围

对于30%的数据， n≤$10^5$
对于50%的数据， n≤$2^{31}$
对于100%的数据，n≤$10^{14}$

Solution

思路

看到这道题目要求求解$\frac{1}{x}$+$\frac{1}{y}$=$\frac{1}{n}$，直接这么求肯定是不好求的，我们考虑将这个式子变形：
通分得：$nx+ny=xy$
移项得：$xy-nx-ny=0$
整理，得：$xy-n(x+y)=0$
$\because 等式两边同时加上一个数，等式仍然成立$
$\therefore 配方，得：xy-n(x+y)+n^2=n^2$
用十字相乘法因式分解，得：$(x-n)(y-n)=n^2$
由此可以推出，只要(x-n)或者(y-n)中有一个是$n^2$的因数，那么方程就有一组解。
所以，我们只需要求出$n^2$的因数个数，将其除以2向上取正就好了。用数学的语言表述就是：令$d(x)$表示$x$的因数个数，答案即为$\lceil \frac{d(n^2)}{2}\rceil$

代码

#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
int  main()
{
//	freopen(".in","r",stdin);
//	freopen(".out","w",stdout);
	long long n;
	cin>>n;
	long long ans=1,t=0;
	for(int i=2;i<=sqrt(n);i++)
	{
		t=0;
		if(n%i==0)
		{
			while(n%i==0)
			{
				n/=i;
				t++;
			}
			ans*=(t*2+1);
		}
		if(n==1)
		{
			break;
		}
	}
	if(n!=1)
	{
		ans*=3;
	}
	cout<<(ans+1)/2;
}