P1445 [Violet]樱花

纪念一下，独立做出来蓝色数学题不容易。

看到阶乘不自觉地想到递推，然后发现好像不太行，然后突然想起数学书上一题~~（其实是作业）~~，问对于任意正整数 $n$，$a,b$ 满足
$$\frac{1}{n+a}+\frac{1}{n+b}=\frac{1}{n}$$
时 $a,b$ 与 $n$ 的关系，显然
$$\begin{gathered} \frac{1}{n+a}+\frac{1}{n+b}=\frac{1}{n} \\ \frac{(n+a)+(n+b)}{(n+a)(n+b)}=\frac{1}{n} \\ \frac{2n+a+b}{n^2+an+bn+ab}=\frac{1}{n} \\ n(2n+a+b)=n^2+an+bn+ab \\ ab=n^2 \end{gathered}$$
考虑把原来的问题转化成上述形式，我们发现 $a,b$ 同号，当 $a,b>0$，$a,b$ 的解数是 $n^2$ 的因数个数，考虑到要求正整数解，需要满足 $n+a,n+b>0$，即 $0>a,b>-n$，但是这样显然 $ab<n^2$ 所以 $a,b>0$ 一定成立。

而我们原来的问题也就转化为了求 $(n!{)}^2$ 的因数个数，考虑到 $n\le 10^6$，我们预处理出来 $\le 10^3$ 的质数然后质因数分解，最后根据因数个数定理就完了，需要注意的一点是对于有大于自己的的平方根的因子的数，这个因子要单独处理。

代码：

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define mem(a, x) memset(a, x, sizeof(a))
#define rps(i, b, e) for(int i=(b);i<=(e);i++)
#define prs(i, e, b) for(int i=(e);i>=(b);i--)
#define rp(i, e) rps(i, 1, e)
#define pr(i, e) prs(i, e, 1)
#define rp0(i, e) rps(i, 0, (e)-1)
#define pr0(i, e) prs(i, (e)-1, 0)
#define rpg(i, x) for(int i=head[x];i;i=e[i].nxt)
#define opf(s) freopen(s".in", "r", stdin), freopen(s".out", "w", stdout)
using namespace std;
const int XR=1100, NR=1e6+10;
int p[XR], c;
bool f[XR+10];
void getprime() {
	f[0]=f[1]=true;
	rps(i, 2, XR) {
		if(!f[i])p[++c]=i;
		for(int j=1;j<=c && i*p[j]<=XR;j++) {
			f[i*p[j]]=true;
			if(i%p[j]==0)break;
		}
	}
}
typedef long long LL;
const LL MOD=1e9+7;
LL ans=1, n, ct[XR], sxd[NR], cnt, t[NR];
int main()
{
	cin>>n;
	getprime();
	rps(i, 2, n) {
		int x=i;
		for(int i=1;i<=c && x>1;i++)
			while(x%p[i]==0)x/=p[i], ct[i]++;
		if(x>1)sxd[++cnt]=x;
	}
	sort(sxd+1, sxd+cnt+1);
	LL sum=0;
	rp(i, cnt+1) {
		if(sxd[i]!=sxd[i-1])ans=ans*(sum*2+1)%MOD, sum=1;
		else sum++;
	}
	rp(i, c)ans=ans*(ct[i]*2+1)%MOD;
	cout<<ans;
	return 0;
}

一遍过，开心。

（建议删掉题目背景，一点也不酸