[bzoj2186]沙拉公主的困惑

2186 沙拉公主的困惑

Time Limit: 10 Sec  
Memory Limit: 259 MB

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

源代码

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#define LL long long
using namespace std;
int T,R;
const int N=10000010;
int fac[N],prefix[N];
bool flag[N];
int prime[N],cnt;
int exgcd(int a,int b,int&x,int&y){
    if(b==0){x=1;y=0;return a;}
    int xx,yy;int g=exgcd(b,a%b,xx,yy);
    x=yy;y=xx-(a/b)*yy;
    return g;
}
int getinv(int a,int m){
    int x,y;
    exgcd(a,m,x,y);
    return (x%m+m)%m;
}
int now=2; 
void pre(int lmt){
    fac[1]=1;prefix[0]=1;
    for(int i=now;i<=lmt;i++,now++){
        fac[i]=fac[i-1]*(LL)i%R;
        if(!flag[i]){
            prime[++cnt]=i; 
            prefix[cnt]=(int)((LL)prefix[cnt-1]*(i-1)%R)*(LL)getinv(i,R)%R;
        }
        for(int j=1;j<=cnt&&i*prime[j]<=N-1;j++){
            flag[i*prime[j]]=true;
            if(i%prime[j]==0)break;
        }
    }
}
 
int getphifac(int n){
    int rtn=fac[n];
    int l=0,r=cnt;
    while(l!=r){
        int mid=(l+r)/2+1;
        if(prime[mid]<=n)l=mid;
        else r=mid-1;
    }
    rtn=(LL)rtn*prefix[l]%R;
    return rtn;
}
int main(){
    scanf("%d%d",&T,&R);
    pre(N-1);
    while(T--){
        int n,m;
        scanf("%d%d",&n,&m);
        int p=getphifac(m);
        int ans=(p*(LL)fac[n])%R*(LL)getinv(fac[m],R)%R;
        printf("%d\n",ans);
    }
}