【ybt高效进阶6-2-3】不定方程（数学）

最新推荐文章于 2021-08-08 00:27:21 发布

SSL_TJH

最新推荐文章于 2021-08-08 00:27:21 发布

阅读量91

点赞数

分类专栏： # 数学或数论文章标签：质数因子数学

要转发的话记得附上链接哦~（不过也不会有人转发的啦）

本文链接：https://blog.csdn.net/weixin_43346722/article/details/119450739

版权

数学或数论专栏收录该内容

285 篇文章 3 订阅

订阅专栏

不定方程

题目链接：ybt高效进阶6-2-3

题目大意

给你 n，问你 1/x+1/y=1/n! 这个不定方程有多少个正整数解。

思路

看到这个就开始化式子：
$1 / x + 1 / y = 1 / n!$
$x n! + y n! = x y$
$y (x - n!) = x n!$
$y = x n! / (x - n!)$

然后发现分母不太好搞的样子，直接用值代替。
$t = x - n!$
$x = n! + t$
$y = (n! + t) n! / t$
$y=(n!)^2/t+n!$

然后你分析一下就会发现，你就是要 $t$ 是 $n!)^2$ 的因子。
那就把 $n!$ 质因数分解，答案就是 $\prod(2\times c_i+1)$ 。（要乘二是因为平方这个质数的个数翻倍了）

代码

#include<cstdio>
#define ll long long
#define mo 1000000007

using namespace std;

int n, prime[1000001];
ll ans;
bool np[1000001];

void get_prime(int n) {
	for (int i = 2; i <= n; i++) {
		if (!np[i]) prime[++prime[0]] = i;
		for (int j = 1; j <= prime[0] && 1ll * i * prime[j] <= n; j++) {
			np[i * prime[j]] = 1;
			if (i % prime[j] == 0) break;
		}
	}
}

int main() {
	scanf("%d", &n);
	
	get_prime(n);
	
	ans = 1;
	for (int i = 1; i <= prime[0]; i++) {
		ll num = 0;
		for (ll j = prime[i]; j <= n; j *= prime[i]) {
			num = (num + n / j) % mo;
		}
		ans = (ans * (2 * num + 1)) % mo;
	}
	
	printf("%lld", ans);
	
	return 0;
}