洛谷P1445 [Violet]樱花

题意

给定一个整数 $\mathrm{n}$，求方程 $\frac{1}{x}+\frac{1}{y}=\frac{1}{n!}$ 的正整数解的组数。
其中$1\le n\le 10^6$。

---

思路

原方程左边通分得 $\frac{x+y}{xy}=\frac{1}{n!}$
交叉相乘，得 $(x+y)n!=xy$
即 $-(x+y)n!+xy=0$
等式两边同时加上$(n!{)}^2$，得$(n!{)}^2-(x+y)n!+xy=(n!{)}^2$
因式分解得 $(x-n!)(y-n!)=(n!{)}^2$
令 $a=x-n!$，$b=y-n!$，则 $ab=(n!{)}^2$
由唯一分解定理知 $n!=\prod _{i=1}^k{p}_i^{c_i}$
则 $(n!{)}^2=\prod _{i=1}^k{p}_i^{2c_i}$
即 $ab=\prod _{i=1}^k{p}_i^{2c_i}$
$\because n!$是确定的
$\therefore$ 确定了$a,b$，就能确定$x,y$。
\quad 而且只要确定了$a$，就能确定$b$。
$\because a$ 是 $(n!{)}^2$ 的因式
$\therefore a$有 $\prod _{i=1}^k(2c_i+1)$ 种取值
那么答案就是 $\prod _{i=1}^k(2c_i+1)\mathrm{m}\mathrm{o}\mathrm{d}(10^9+7)$

（以上推导过程转载自M_sea神仙的题解，如果 $\mathrm{M}\_\mathrm{s}\mathrm{e}\mathrm{a}$ 神仙不同意我的行为，请告诉我，我将很快删除这一段qwq）

---

代码

const
	nn=1000000;
	mo=1000000007;
var
	a,s,v:array[0..nn+10] of longint;
	i,j,cnt,n,t:longint;
	ans:int64; //安全起见，防止爆炸
begin
	read(n);
	fillchar(s,sizeof(s),0);
	cnt:=0;
	fillchar(v,sizeof(v),0);
	for i:=2 to n do
	begin
		if v[i]=0 then
		begin
			inc(cnt);
			a[cnt]:=i;
			v[i]:=i;
		end;
 		for j:=1 to cnt do
		begin
			t:=a[j];
			if (t>v[i]) or (t*i>n) then break;
			v[t*i]:=a[j];
		end; 
	end; //线性筛法求素数
	for i:=2 to n do 对n!进行质因数分解
	begin
		t:=i;
		while t<>1 do
		begin
			inc(s[v[t]]);
			t:=t div v[t]; //因为v[t]必然是t的因子
		end;
	end;
	ans:=1;
	for i:=1 to cnt do 
		ans:=ans*(s[a[i]] shl 1+1) mod mo; //计算答案
	writeln(ans);
end.

时间复杂度为 $\mathcal{O}(n\log n)$