JZOJsenior3456.【NOIP2013模拟联考3】恭介的法则(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

思路：
比赛时根本弃疗颓废了。这题思路巧妙
首先题目说要求满足1/x+1/y=1/n！的正整数对（x，y）的对数
（说在前面，x和y是大于等于n才成立的，我不会证明，手玩一下可知）
接下来通过简单的性质解决这道题（下面的步骤最好在纸上面演算）

---

将上式通分，得到（x+y）/xy=1/n！

左右各乘n！，得到n！（x+y）=xy

拆项移项，得到x（y-n！）=yn！——>x=yn！/（y-n！）
这里是很重要的一步，我们设y=n！+k（k>0且k∈N）
那么代入得到x=n！（n！+k）/k——>x=（（n！）^2+kn！）/k

因为x＞0且x∈N，所以k |（（n！）^2+kn！）
而kn！必为k的倍数，满足x为正整数的必要条件是k |（n！）^2
所以原问题转化为求（n！）^2的因数个数

---

拿样例举个栗子
（6！）^2=518400=2^8×3^4×5^2
所以因数个数为（8+1）×（4+1）×（2+1）=135，没有错

求（n！）^2的因数个数有个简单的方法：
把1~n的每个数的质因数求出来，所有质因数的指数平方即可
这道题就解决了，记得用高精度加压位

代码（这里高精度压八位）：

#include<bits/stdc++.h>
#define MAXN 700001
#define mod 100000000

using namespace std;

long long a[MAXN];
bool bz[MAXN];
int n,len;

void work(int x)
{
     int k=0;
     for (int i=1;i<=len;i++)
     {
         a[i]=a[i]*x+k;
         k=a[i]/mod;
         a[i]%=mod;
     }
     if (k)
     {
         a[++len]=k;
     }
     while(a[len]>=mod)
     {
         a[len+1]=a[len]/mod;
         a[len]=a[len]%mod;
         len++;
     }
}

int main()
{
    scanf("%d",&n);
    len=1;
    a[1]=1;
    long long x=1;
    memset(bz,1,sizeof(bz));
    for (int i=2;i<=n;i++)
    {
        if (bz[i])
        {
            long long y=n,z=0;
            while(y>=i)
            {
                y/=i;
                z+=y;
            }
            z=2*z+1;
            if (x*z>mod)
            {
                work(x);
                x=1;
            }
            x*=z;
            for (int j=1;j<=n/i;j++)bz[i*j]=0;
        }
    }
    if (x>1)work(x);
    printf("%d",a[len]);
    for (int i=len-1;i>=1;i--)printf("%08lld",a[i]);
}