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

1

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

HINT

Source

代码

#include<cstdio>
#include<cmath>
using namespace std;
typedef long long ll;
const int N=1e7+10;
//const int M=N/10;
bool check[N]={1,1};
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=2;i<=m;i++){
        if(!check[i]){
            for(int j=i*i;j<N;j+=i){
                check[j]=1;
            }
        }
    }
}
int main(){
    //freopen("sh.txt","r",stdin);
    get_prime();
    scanf("%d%d",&T,&mod);
    ans1[0]=1;
    for(int i=1;i<N;i++) ans1[i]=ans1[i-1]*i%mod;
    inv[1]=1;
    for(int i=2;i<N;i++){
        if(i>=mod) break;
        inv[i]=(mod-mod/i)*inv[mod%i]%mod;
    }
    ans2[1]=1;
    for(int i=2;i<N;i++){
        if(!check[i]){
            ans2[i]=ans2[i-1]*(i-1)%mod;
            ans2[i]=ans2[i]*inv[i%mod]%mod;
        }
        else ans2[i]=ans2[i-1];
    }
    while(T--){
        scanf("%d%d",&n,&m);
        ll ans=ans1[n]*ans2[m]%mod;
        printf("%lld\n",ans);
    }
    return 0;
}

__EOF__

本文作者：shenben
本文链接：https://www.cnblogs.com/shenben/p/6264347.html
关于博主：评论和私信会在第一时间回复。或者直接私信我。
版权声明：本博客所有文章除特别声明外，均采用 BY-NC-SA 许可协议。转载请注明出处！
声援博主：如果您觉得文章对您有帮助，可以点击文章右下角【推荐】一下。您的鼓励是博主的最大动力！