传送门

推一波式子：

1x+1y=1n!\frac 1 x+\frac 1 y=\frac 1 {n!}x1​+y1​=n!1​

=>xy−x∗n!−y∗n!xy-x*n!-y*n!xy−x∗n!−y∗n! = 000

=>(x−n!)(y−n!)=(n!)2(x-n!)(y-n!)=(n!)^2(x−n!)(y−n!)=(n!)2

于是把(n!)2(n!)^2(n!)2质因数分解就行了。

代码：

#include<bits/stdc++.h>
using namespace std;
const int mod=1e9+7,N=1e6+5;
int n,pri[N],tot=0,cnt[N],ans=1;
bool vis[N];
inline void solve(int x){
	for(int i=1;i<=tot&&pri[i]*pri[i]<=x;++i){
		if(x==1)return;
		if(x!=x/pri[i]*pri[i])continue;
		while(x==x/pri[i]*pri[i])++cnt[pri[i]],x=x/pri[i];
	}
	if(x^1)++cnt[x];
}
int main(){
	scanf("%d",&n);
	for(int i=2;i<=n;++i){
		if(!vis[i])pri[++tot]=i;
		for(int j=1;i*pri[j]<=n;++j){
			int k=i*pri[j];
			vis[k]=1;
			if(i==i/pri[j]*pri[j])break;
		}
		solve(i);
	}
	for(int i=1;i<=n;++i)ans=1ll*ans*(cnt[i]?(cnt[i]<<1|1):1)%mod;
	cout<<ans;
	return 0;
}