【SDOI2008】【BZOJ2186】沙拉公主的困惑

Description

　　大富翁国因为通货膨胀，以及假钞泛滥，政府决定推出一项新的政策：现有钞票编号范围为1到N的阶乘，但是，政府只发行编号与M!互质的钞票。房地产第一大户沙拉公主决定预测一下大富翁国现在所有真钞票的数量。现在，请你帮助沙拉公主解决这个问题，由于可能张数非常大，你只需计算出对R取模后的答案即可。R是一个质数。

Input

第一行为两个整数T，R。R<=10^9+10，T<=10000，表示该组中测试数据数目，R为模后面T行，每行一对整数N，M，见题目描述 m<=n

Output

共T行，对于每一对N，M，输出1至N！中与M！素质的数的数量对R取模后的值

Sample Input

1 11

4 2

Sample Output

1

数据范围：

对于100%的数据，1 < = N , M < = 10000000

HINT

Source

答案是 ϕ(m!)∗n!/m!
考察线性求逆元(以前的数论板子现在也有用了..

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define GET (ch>='0'&&ch<='9')
#define MAXN 10000010
using namespace std;
int T,R,n,m;
bool not_prime[MAXN];
int prime[MAXN],top;
int fac[MAXN],inv[MAXN],ret[MAXN];
void in(int &x)
{
    char ch=getchar();x=0;
    while (!GET)    ch=getchar();
    while (GET) x=x*10+ch-'0',ch=getchar();
}
void init()
{
    fac[1]=1;inv[1]=1;ret[1]=1;
    for (int i=2;i<=MAXN-10;i++)    fac[i]=1ll*fac[i-1]*i%R;
    for (int i=2;i<=MAXN-10;i++)
    {
        if (!not_prime[i])  prime[++top]=i;
        for (int j=1;j<=top&&i*prime[j]<=MAXN-10;j++)
        {
            not_prime[i*prime[j]]=1;
            if (i%prime[j]==0)  break;
        }
    }
    for (int i=2;i<=MAXN-10&&i<R;i++)   inv[i]=1ll*(R-R/i)*inv[R%i]%R;
    for (int i=2;i<=MAXN-10;i++)
    {
        ret[i]=ret[i-1];
        if (!not_prime[i])  ret[i]=1ll*ret[i]*(i-1)%R*inv[i]%R%R;
    }
}
int main()
{
    for (in(T),in(R),init();T;T--)  in(n),in(m),printf("%d\n",(1ll*fac[n]*ret[m]+R)%R);
}