bzoj 2186（欧拉函数）

传送门
题解：
引理1：若gcd（x,m）==1，则gcd（x+k*m,m）==1(k>0)（反证法）
引理2：由引理1，又因为1~b中与b互质的数有φ(b)个，有以下事实成立：若b|a，则1~a中与b互质的数有a/b*φ(b)个。
原问题：求1~n！中与m！互质的数的个数，设nn=n!，mm=m!。
则ans=nn/mm*φ(mm)
又因为φ(mm)=mm*∏((pi-1)/pi)，其中pi为所有不大于mm的质数。
所以ans=nn*∏((pi-1)/pi)。
然后线性筛预处理出阶乘fac[ ]，逆元rev[ ]，还有上面那个ans[ ]。
最后对于每组询问O(1)查询。

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int MAXM=1e7+5;
typedef long long ll;
bool vis[MAXM];
ll prime[MAXM/10],tot=0,ans[MAXM],fac[MAXM],rev[MAXM];
int n,m,kase,p;
void init() {
    memset(vis,false,sizeof(vis));
    fac[1]=1,rev[1]=1,ans[1]=1;
    for (register int i=2;i<MAXM;++i) {
        fac[i]=fac[i-1]*i%p;
        if (i<p) rev[i]=(p-p/i)*rev[p%i]%p;
        ans[i]=!vis[i]?ans[i-1]*(i-1)%p*rev[i]%p:ans[i-1];
        if (!vis[i]) prime[++tot]=i;
        for (int j=1;j<=tot&&i*prime[j]<MAXM;++j) {
            vis[i*prime[j]]=true;
            if (i%prime[j]==0) break;
        }
    }       
}
int main() {
    scanf("%d%d",&kase,&p);
    init();
    while (kase--) {
        scanf("%d%d",&n,&m);
        printf("%lld\n",fac[n]*ans[m]%p);
    }
    return 0;
}