丢番图

题意：
对于一个给定的正整数n ，有多少种本质不同的解满足方程
$\frac{1}{x}$ + $\frac{1}{y}$ = $\frac{1}{n}$ (x,y,n∈N*)

$n\le 1e14$

分析：
n = $\prod\limits_{i=1}^kp^{c_i}_i$
包括本质相同的解得个数
ans = $\prod\limits_{i=1}^k(2 * c_i + 1)$
证明：

$$\begin{aligned} \frac{x + y}{xy} & =\frac{1}{n} \\ n(x +y) & = xy \\ nx &= (x - n)y \\ \frac{nx}{x-n} &=y \end{aligned}$$
令k=x-n，则x = k + n
$\therefore$有 $\frac{n^2+kn}{k}$= y
即 $\frac{n^2}{k}$+n=y
$\because$y,n,k∈N*
$\therefore$ 原式解的个数为 $n^2$的约数个数
原题要求本质不相同的解的个数为(ans + 1) / 2 (x,y为2n的解只计算了一次)

#include <cstdio>
#include <cstring>

typedef long long LL;
const int N = 1e7 + 10,M = 1e6;
int prime[M],tot;
bool p[N];

void work() {
    for (int i = 2;i <= N - 10;i ++) {
        if (!p[i]) prime[++ tot] = i;
        for (int j = 1;j <= tot;j ++) if (prime[j] * i > N - 10) break;
        else p[prime[j] * i] = true;
    }
}

int main() {
    work();
    LL n,ans = 1;
    scanf("%lld",&n);
    for (int i = 1;i <= tot;i ++) if (prime[i] > n) break;else if (n % prime[i] == 0){
        LL t = 0;
        while (n % prime[i] == 0) {
            t ++;
            n /= prime[i];
        }
        ans *= 2 * t + 1;
    }
    if (n != 1) ans *= 3;
    printf("%lld",(ans + 1) >> 1);
}

考试时脑抽打了个不必要的线筛….