恭介的法则(rule)

$$恭介的法则(rule)$$

题目链接：$jzoj3456$

题目

终于，在众亲们的奋斗下，最终 $boss$ 恭介被关进了库特设计的密室。正当她们松了一口气时，这个世界却发生了天翻覆地的变化：地面开始下沉，天空开始变成血红色，海水沸腾……一幅世界末日的图景。美鱼从她手中的古籍《若山牧水诗歌集》中发现了原因：$白鸟はかなしからずや空の青海のあをにも染まずただよふ$ 。大$（xia）$意$（shuo）$就是狡猾的恭介在创造这个世界的时候就篡改了法则。而这个法则的起源，就是一只生死之间的猫。这个猫被关在一个黑盒子里，盒子里有两个毒气罐，如果有任意一个毒气罐被打开那么猫将会被杀死，法则也能得到纠正。然而外界能控制的仅仅是这两个毒气罐被打开的概率。假设第一个毒气罐被打开的概率为$1/x$，第二个毒气罐为$1/y$（$x$,$y$ 为正整数），那么当两个概率和为$1/(n!)$时，猫将会被莫名其妙地杀死。现在美鱼想知道，有多少对(x,y)可以让猫被莫名其妙杀死。

输入

一行，一个正整数$n$

输出

一行，满足题意的$(x,y)$对数。

样例输入

6

样例输出

135

数据范围

对于$30\%$的数据 $n<=6$

对于$60\%$的数据 $n<=50$

对于$100\%$的数据 $n<=700000$

思路

这道题是一道数论题。
（下面的推论与解法是看了某位大佬的博客做出来的）

由题意可以得出这个式子：
$$\frac{1}{x}+\frac{1}{y}=\frac{1}{n!}$$
两边都乘一个$x$，就得到这个：
$$1+\frac{x}{y}=\frac{x}{n!}$$
移项得到这个：
$$1=\frac{x}{n!}-\frac{x}{y}$$
然后通分得到这个：
$$1=\frac{x\ast (y-n!)}{y\ast n!}$$
把$y\ast n!$移到左边，就是这个：
$$y\ast n!=x\ast (y-n!)$$
再把右边的$(y-n!)$移过去，就变成了这一个：
$$\frac{y\ast n!}{y-n!}=x$$
即：
$$x=\frac{y\ast n!}{y-n!}$$

又因为这三个数（$x$，$y$和$n$）都是正整数，就可以得到$y\ast n!$是$y-n!$的倍数，而且$y>n!$。
那我们再设一个正整数$k$，让$k=y-n!$，那$y=k+n!$，那么我们前面的式子就变成了这个：$$x=\frac{(n!+k)\ast n!}{(n!+k)-n!}$$
分解一下就是：
$$x=\frac{n{!}^2+k\ast n!}{k}$$
那我们要让$x$是正整数，就要让$\frac{n{!}^2}{k}$都是正整数（因为$x=\frac{k\ast n!}{k}$一定是正整数）

那就变成求$n{!}^2$有多少个因数。

那我们让$k^{zhi}=n$，那我们通过证明可以得到有$zhi\ast 2+1$种方案。
（因为平方就是$n\ast n$，就是$k^{zhi}\ast k^{zhi}$即$k^{2\ast zhi}$，那从不选到选$2\ast zhi$个，就一共有$zhi\ast 2+1$种方案了）

那就素数筛选法，求$n$素数做。
（但是因为答案有“一”点大，要压位高精，我这里压了$8$位）

代码

#include<cstdio>
#define ll long long
#define mo 100000000

using namespace std;

ll n, len, ans[700001];
bool in[700001];

void cheng(ll now) {//高精乘
	ll jin = 0;
	for (ll i = 1; i <= len; i++) {
		ans[i] = ans[i] * now + jin;
		jin = ans[i] / mo;
		ans[i] %= mo;
	}
	while (jin) {
		ans[++len] = jin;
		jin = ans[len] / mo;
		ans[len] %= mo;
	}
}

int main() {
	len = 1;
	ans[1] = 1;//初始化
	
	scanf("%lld", &n);//读入
	
	ll an = 1;
	for (ll i = 2; i <= n; i++)
		if (!in[i]) {//找到质数
			ll num = n, zhi = 0, kind;
			while (num >= i) {//求出指数
				num /= i;
				zhi += num;
			}
			kind = 2 * zhi + 1;
			if (an * kind > mo) {//需要压位
				cheng(an);
				an = kind;
			}
			else an *= kind;
			for (ll j = 1; j * i <= n; j++)
				in[j * i] = 1;
		}
	
	if (an != 1) cheng(an);//剩下的
	
	printf("%lld", ans[len]);//输出
	for (ll i = len - 1; i >= 1; i--)
		printf("%08lld", ans[i]);//输出
	
	return 0;
}