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

线性求逆元后好神的数论啊。。。。

#include<cstdio>
#include<cstring>
#include<cctype>
#include<algorithm>
using namespace std;
#define rep(i,s,t) for(register int i=s;i<=t;i++)
#define dwn(i,s,t) for(register int i=s;i>=t;i--)
#define clr(x,c) memset(x,c,sizeof(x))
#define ll long long
int read(){
    int x=0;char c=getchar();
    while(!isdigit(c)) c=getchar();
    while(isdigit(c)) x=x*10+c-'0',c=getchar();
    return x;
}
const int nmax=1e7+5;
int prime[nmax>>3],ans[nmax],fac[nmax],inv[nmax];bool vis[nmax];
int n,m,mod,T;
void init(){
    int cnt=0,tmp;
    rep(i,2,nmax-1){
        if(!vis[i]) prime[++cnt]=i;
        rep(j,1,cnt){
            tmp=prime[j];
            if((ll)i*tmp>=nmax) break;
            vis[i*tmp]=1;
            if(i%tmp==0) break;
        }
    }
     
    fac[1]=1;
    rep(i,2,nmax-1) fac[i]=(ll)fac[i-1]*i%mod;
     
    inv[1]=1;
    rep(i,2,nmax-1) inv[i]=(ll)(mod-mod/i)*inv[mod%i]%mod;
    ans[1]=1;
    rep(i,2,nmax-1){
        if(!vis[i]) ans[i]=(ll)ans[i-1]*(i-1)%mod*inv[i%mod]%mod;
        else ans[i]=ans[i-1];
    }
}
int main(){
    T=read(),mod=read();
    init();
    while(T--){
        int n=read(),m=read();
        printf("%d\n",(ll)fac[n]*ans[m]%mod);
    }
    return 0;
}

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

Time Limit: 10 Sec Memory Limit: 259 MB
Submit: 3190 Solved: 1089
[Submit][Status][Discuss]

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

[ Submit][ Status][ Discuss]

秒客网