YBTOJ C. 【例题3】不定方程

最新推荐文章于 2024-07-24 10:46:51 发布

龑钧

最新推荐文章于 2024-07-24 10:46:51 发布

阅读量105

点赞数

分类专栏：质数与约数文章标签： c语言算法 c++

本文链接：https://blog.csdn.net/DL24yanjun/article/details/126444304

版权

质数与约数专栏收录该内容

3 篇文章 0 订阅

订阅专栏

在这里插入图片描述
方程移项可得 $y=\frac{xn!}{x-n!}$ 设 $x - n! = t$ 得到 $y=n!+\frac{n!^2}{t}$ 为题就转换为 $n!^2$ 的约数，所以先算n以内的质数，再计算每个质数出现的次数，即可
公式 $\prod_{i}^{cnt}(c[i]*2)+1$
注意：记得取余！
Code

#include<bits/stdc++.h>
#define re register

#define inl inline

#define int long long

using namespace std;
int read(){
	int sum=0,f=1;char c=getchar();
	while(!isdigit(c)){if(c=='-') f=-1;c=getchar();}
	while(isdigit(c)){sum=(sum<<3)+(sum<<1)+(c^48);c=getchar();}
	return f*sum;
}
int c[10000010],cnt;
const int mod=1e9+7;
int pri[10000010];
int v[10000010];
inl void init(int n){
	for(re int i=2;i<=n;i++){
		if(!v[i]) v[i]=i,pri[++cnt]=i;
		for(int j=1;j<=cnt;j++){
			if(pri[j]>v[i]||pri[j]>n/i) break;
			v[pri[j]*i]=pri[j];
		}
	}
}
int ans=1;
signed main(){
	int n=read();
	init(n);
	for(re int i=1;i<=cnt;i++){
		for(re int j=pri[i];j<=n;j*=pri[i]){
			c[i]+=n/j;
			c[i]%=mod;
		}
	}
	for(re int i=1;i<=cnt;i++){
		ans=(ans*(c[i]*2+1)%mod+mod)%mod;
	}
	printf("%lld",ans);
}