【LG1445】樱花

【LG1445】樱花

题面

洛谷

题解

\[\frac 1x+\frac 1y=\frac 1{n!}\\
\frac{x+y}{xy}=\frac 1{n!}\\
n!(x+y)=xy\\
xy-n!(x+y)=0\\
(n!)^2+xy-n!(x+y)=(n!)^2\\
(x-n!)(y-n!)=(n!)^2
\]

设\(a=x-n!,b=y-n!\)

\[ab=(n!)^2
\]

枚举约数就好啦

代码

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;
const int Mod = 1e9 + 7;
const int MAX_N = 1e6 + 5;
void pls(int &x, int y) { x += y; if (x >= Mod) x -= Mod; }
bool is_prime[MAX_N];
int N = 1e6, prime[MAX_N], num, ans = 1;
void sieve() {
	for (int i = 1; i <= N; i++) is_prime[i] = 1;
	is_prime[1] = 0;
	for (int i = 2; i <= N; i++) {
		if (is_prime[i]) prime[++num] = i;
		for (int j = 1; prime[j] * i <= N && j <= num; j++) {
			is_prime[i * prime[j]] = 0;
			if (!(i % prime[j])) break;
		}
	}
}
int main () {
	sieve();
	cin >> N;
	for (int i = 1; i <= num; i++) {
		if (prime[i] > N) break;
		int cnt = 0;
		for (long long j = prime[i]; j <= N; j = j * prime[i]) pls(cnt, 1ll * N / j);
		pls(cnt, cnt), pls(cnt, 1);
		ans = 1ll * ans * cnt % Mod;
	}
	printf("%d\n", ans);
	return 0;
}