BZOJ 2186 沙拉公主的困惑

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

Time Limit: 10 Sec Memory Limit: 259 MB

Submit: 3397 Solved: 1164
[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

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

HINT

Source

代码

#include<cstdio>

#include<cmath>

using namespace std;

typedef long long ll;

const int N=1e7+;

//const int M=N/10;

bool check[N]={,};

int n,m,mod,tot,T;

//int prime[M];

ll ans1[N],ans2[N],inv[N];

/*void get_prime(){

    for(int i=2;i<N;i++){

        if(!check[i]) prime[++tot]=i;

        for(int j=1;j<=tot&&i*prime[j]<=n;j++){

            check[i*prime[j]]=1;

            if(i%prime[j]==0) break;

        }

    }

}*/

//面向std编程，不要问我为什么

void get_prime(){

    int m=sqrt(N+0.5);

    for(int i=;i<=m;i++){

        if(!check[i]){

            for(int j=i*i;j<N;j+=i){

                check[j]=;

            }

        }

    }

}

int main(){

    //freopen("sh.txt","r",stdin);

    get_prime();

    scanf("%d%d",&T,&mod);

    ans1[]=;

    for(int i=;i<N;i++) ans1[i]=ans1[i-]*i%mod;

    inv[]=;

    for(int i=;i<N;i++){

        if(i>=mod) break;

        inv[i]=(mod-mod/i)*inv[mod%i]%mod;

    }

    ans2[]=;

    for(int i=;i<N;i++){

        if(!check[i]){

            ans2[i]=ans2[i-]*(i-)%mod;

            ans2[i]=ans2[i]*inv[i%mod]%mod;

        }

        else ans2[i]=ans2[i-];

    }

    while(T--){

        scanf("%d%d",&n,&m);

        ll ans=ans1[n]*ans2[m]%mod;

        printf("%lld\n",ans);

    }

    return ;

}

巴特西