POJ-2917 / HDU-1299 Diophantus of Alexandria

题目链接：

$\quad$ POJ-2917 $\quad$ HDU-1299

题目大意：

$\quad$ 给一个 $n$ ，求满足等式 $\frac{1}{x} + \frac{1}{y} = \frac{1}{n}$ 且 $x \leq y$ 的解的对数？

数据范围：

$\quad$ $1 \leq n \leq 10^9$

解题思路：

$\quad$ 显然 $x \gt n, y \gt n$ ；所以，令 $x = n + a, y = n + b$；
$\quad$ 化简可得 : $n ^ 2 = a b$ ；
$\quad$ 那么，问题就跟求 $n^2$ 的约数有关系了。
$\quad$ 由唯一分解定理，$n = {p_1}^{r_1} \cdot {p_2}^{r_2} \ldots {p_k}^{r_k}$ ；$n$ 的约数个数 ${Factor_n} = (1 + {r_1}) (1 + {r_2}) \ldots (1 + {r_k})$ 个。
$\quad$ 同理：${n^2} = {p_1}^{2 \cdot r_1} \cdot {p_2}^{2 \cdot r_2} \ldots {p_k}^{2 \cdot r_k}$ ; $n^2$ 的约数 ${Factor_{n^2}} = (1 + 2{r_1}) (1 + 2{r_2}) \ldots (1 + 2{r_k})$;
$\quad$ 所以，对 $n$ 分解质因数就可以了。

$\quad$ PS: 大于 $sqrt(n)$ 的质因子如果存在，它的指数必然是 $1$ 。
$\quad$ PPS: 对于 $x = y$ (即 $a = b$) 的情况，这种情况只会被计算 $1$ 次，其他情况都是成对出现的；也就是说，$Factor_{n^2}$ 一定是一个奇数。当然，从上面的分解情况来看，这个也是成立的。

$\quad$ 最后既然要求 $x \leq y$ 的对数，那么把 $Factor_{n^2}$ 砍半就可以了。

AC代码：

/**********************************************
 *Author*        :XzzF
 *Created Time*  : 2018/3/28 12:28:03
 *Ended  Time*  : 2018/3/28 12:32:04
*********************************************/

#include <cstdio>
#include <cmath>
#include <cstdlib>
#include <cstring>
#include <algorithm>
using namespace std;
typedef long long LL;
const int inf = 1 << 30;
const LL INF = 1LL << 60;

int n, T;

int main()
{
    scanf("%d", &T);
    for (int cas = 1; cas <= T; cas++) {
        scanf("%d", &n);
        LL ans = 1;
        for (int i = 2; i * i <= n; i++) {
            if (n % i == 0) {
                int r = 0;
                while (n % i == 0) {
                    r++;
                    n /= i;
                }
                ans = (ans * (2 * (LL)r + 1));
            }
        }
        if (n > 1) ans = (ans * (2 * 1 + 1));   //存在大于sqrt(n)的质数
        printf("Scenario #%d:\n", cas);
        //printf("%lld\n", (ans + 1) / 2);   //砍半
        printf("%lld\n\n", ans / 2 + 1);
    }
    return 0;
}