bzoj2186: [Sdoi2008]沙拉公主的困惑

显然答案等于$\phi(m!)*(n!/m!)$
n!很容易求。
$\phi(x)=x*\frac{p_1-1}{p_1}*\frac{p_2-1}{p_2}*...$
所以$\frac{\phi(m!)}{m!}=\frac{p_1-1}{p_1}*\frac{p_2-1}{p_2}*...$pi为<=m的质数。
定义其为f(i).
则有

$$\begin{cases} &f(i)=f(i-1) \text{ } x 为合数\\ & f(i)=f(i-1)*(i-1)/m\text{ } x为质数 \end{cases}$$

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int read()
{
    char ch=getchar();int f=0;
    while(ch<'0'||ch>'9')ch=getchar();
    while(ch>='0'&&ch<='9'){f=(f<<1)+(f<<3)+ch-'0';ch=getchar();}
    return f;
}
bool check[10000005];long long prime[3000005],fac[10000005],n,mod,tot,f[10000005];
inline void exgcd(int a,int b,int &x,int &y)  
{  
    if (!b){x=1;y=0;return;}  
    exgcd(b,a%b,x,y);  
    int t=x;x=y;y=t-a/b*x;  
}  
inline int getinv(int a)  
{  
    int x=0,y=0;  
    exgcd(a,mod,x,y);  
    return (x%mod+mod)%mod;  
}  
int main()
{
    n=read(),mod=read();
    for(int i=2;i<=10000000;i++)
    {
        if(!check[i])
        {
            prime[tot++]=i;
        }
        for(int j=0;j<tot;j++)
        {
            if(i*prime[j]>10000000)
            break;
            check[i*prime[j]]=1;
            if(i%prime[j]==0)
            break;
        }
    }
    f[1]=1;fac[1]=1;
    for(int i=2;i<=10000000;i++)
    {
        if(!check[i])
        f[i]=1LL*(1LL*f[i-1]*(i-1))%mod*getinv(i)%mod;
        else f[i]=f[i-1];
        fac[i]=(1LL*fac[i-1]*i)%mod;
    }
    for(int i=1;i<=n;i++)
    {
        int x=read(),y=read();
        printf("%lld\n",1LL*f[y]*fac[x]%mod);
    }
}