bzoj2186 [Sdoi2008]沙拉公主的困惑（数论）

求1~n!中与m！互质的数的个数。
最后结论就是$n!\prod \frac{p_i-1}{p_i}$
有点像欧拉函数那样，我们线性筛预处理素数，阶乘，逆元即可。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 10000010
inline char gc(){
    static char buf[1<<16],*S,*T;
    if(S==T){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=gc();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x*f;
}
int mod,n,m,inv[N],prime[N/10],tot=0,fac[N],mul[N];
bool notprime[N];
inline void init(){
    notprime[1]=1;inv[1]=1;mul[1]=1;fac[0]=1;
    for(int i=2;i<=1e7;++i){
        if(!notprime[i]) prime[++tot]=i;
        for(int j=1;prime[j]*i<=1e7;++j){
            notprime[prime[j]*i]=1;
            if(i%prime[j]==0) break;
        }
    }for(int i=2;i<=prime[tot];++i){
        if(i>=mod) break;
        inv[i]=(ll)inv[mod%i]*(mod-mod/i)%mod;
    }for(int i=1;i<=1e7;++i) fac[i]=(ll)fac[i-1]*i%mod;
    for(int i=2;i<=1e7;++i){
        if(notprime[i]) mul[i]=mul[i-1];
        else mul[i]=(ll)mul[i-1]*(i-1)%mod*inv[i%mod]%mod;
    }
}
int main(){
//  freopen("a.in","r",stdin);
    int tst=read();mod=read();init();
    while(tst--){
        n=read();m=read();
        printf("%d\n",(ll)fac[n]*mul[m]%mod);
    }return 0;
}