【Ybtoj 第28章例3】不定方程【质数和约数】

---

解题思路

由原题的式子可推出：$y=\frac{xn!}{x-n!}$，我们设$t=x-n!$，则$x=t+n!$，得到$y=n!+\frac{(n!{)}^2}{t}$，因为n和y是整数，所以t一段要是$(n!{)}^2$的约数，题目就变成了求$(n!{)}^2$的约数个数，即求${\prod }_{i=1}^m(c_i+1)$
对于一个合数n^2，它的约数个数为${\prod }_{i=1}^m(c_i\ast 2+1)$，所以首先求n中的质数个数，然后求出所有的质数在n！中出现的次数，

最后通过${\prod }_{i=1}^m(c_i\ast 2+1)$得到$(n!{)}^2$的约数个数。

---

代码

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<iomanip>
#include<cstring>
#include<cmath>
#include<map>
#include<queue>
#define ll long long
#define ldb long double
using namespace std;

const ll mod=1e9+7;
ll n,p,cnt,pri[10000010],v[10000010],c[10000010],ans=1;

int main(){
	scanf("%lld",&n);
	for(ll i=2;i<=n;i++)
	{
		if(!v[i])
		{
			v[i]=i;
			pri[++cnt]=i;
		}
		for(ll j=1;j<=cnt;j++)
		{
			if(pri[j]>v[i]||pri[j]>n/i)break;
			v[pri[j]*i]=pri[j];
		}
	}
	for(ll i=1;i<=cnt;i++)
	{
		for(ll j=pri[i];j<=n;j*=pri[i])
		{
			c[i]+=n/j;
			c[i]%=mod;
		}
	}
	for(int i=1;i<=cnt;i++)
		ans=ans*((2*c[i]+1)%mod)%mod;
	printf("%lld",ans); 
}