恭介的法则(rule)

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

Input
一行，一个正整数n

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

Sample Input
6

Sample Output
135

Data Constraint
对于30%的数据 n<=6
对于60%的数据 n<=50
对于100%的数据 n<=700000

.
.
.
.
.
.
分析
对原式瞎搞得：


我们将y表示为一般的形式，即y=n!+k(k代表任意符合条件的数)，带入前面的式子中，得

合并同类项化简得

因为x是整数，所以n!n!+kn!必定为k的倍数，因为kn!是k的倍数，故不必考虑第二个多项式，我们只要考虑第一个多项式，即满足如下条件

这样，答案就转化成求n!n!的约数个数了

设n!=p1^e1 * p2^e2 * p3^e3 * … * pk^ek，则 n! * n!= p1^(2 * e1) * p2^(2 * e2) * … * pk^(2 * ek) 。 则因子个数sum=(2 * e1+1) * (2 * e2+1) * … * (2 * ek+1)。
对于每个质数p，它的出现次数等于[N/p]+[N/p^2]+…，这样对于每个p是O(log_p(N))的，而1～N之间的质数有O(N/lnN)个，均摊一下就是O(n)的。
至于怎么算每个质数p出现的次数，p的次数等于[n/p]+[n/p ^ 2]+[n/p ^ 3]+….这个可以看成是把1 ~ n里p的倍数的数字都除掉一个p，这样的数有[n/p]个，除完之后这[n/p]个数会变成1 ~ [n/p]，然后再除掉一个p，就是[[n/p]/p]个，然后变成1~ [n/p^2]，一直做下去。

.
.
.
.
.
.
.
程序：

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;

int len=1,mo=100000000;
long long ans[30000],k=1;
bool bz[700010];

void mul(int x)
{
    for (int i=1;i<=len;i++) 
		ans[i]*=x;
    for (int i=1;i<=len;i++)
    {   
        ans[i+1]+=ans[i]/mo;
        ans[i]%=mo;
    }
    while (ans[len+1]>0)
    {
        len++;
        ans[len+1]+=ans[len]/mo;
        ans[len]%=mo;
    }
}

int main()
{
    int n;
    scanf("%d",&n);
    ans[1]=1;
    memset(bz,false,sizeof(bz));
    for (int i=2;i<=n;i++)
    	if (bz[i]==false)
    	{
    		int x=n,x1=0;
    		while (x>=i)
    		{
    			x/=i;
    			x1+=x;
    		}
    		x1=x1*2+1;
    		if (x1*k>mo)
    		{
    			mul(k);
    			k=1;
    		}
    		k*=x1;
    		for (int j=i;j<=n;j+=i)
    			bz[j]=true;
    	}
    if (k>1) mul(k);
    printf("%lld",ans[len]);
    for (int i=len-1;i>=1;i--)
		printf("%08lld",ans[i]);
	return 0;
}