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

数据范围

对于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);}
}

转载于:https://www.cnblogs.com/Anoxiacxy/p/6934937.html

[bzoj2186]沙拉公主的困惑相关推荐

bzoj2186【SDOI2008】沙拉公主的困惑
2186: [Sdoi2008]沙拉公主的困惑 Time Limit: 10 Sec Memory Limit: 259 MB Submit: 2363 Solved: 779 [ Submi ...
[SDOI2008]沙拉公主的困惑线性筛素数+欧拉
本文为博主原创文章,欢迎转载,请注明出处 www.cnblogs.com/yangyaojia [SDOI2008]沙拉公主的困惑线性筛素数+欧拉题目大意给定n,m,求在1到n!内与m!互质的 ...
φ(1)=1 [Sdoi2008]沙拉公主的困惑
问题 F: [Sdoi2008]沙拉公主的困惑时间限制: 1 Sec 内存限制: 259 MB 题目描述大富翁国因为通货膨胀,以及假钞泛滥,政府决定推出一项新的政策:现有钞票编号范围为1到N的阶乘 ...
沙拉公主的困惑 cash
一.沙拉公主的困惑 (cash.pas/c/cpp, 限时3 秒,内存256M) [题目描述] 大富翁国因为通货膨胀,以及假钞泛滥,政府决定推出一项新的政策:现有钞票编号范围为1 到N 的阶乘,但是 ...
【bzoj2186】[Sdoi2008]沙拉公主的困惑欧拉函数
题目描述大富翁国因为通货膨胀,以及假钞泛滥,政府决定推出一项新的政策:现有钞票编号范围为1到N的阶乘,但是,政府只发行编号与M!互质的钞票.房地产第一大户沙拉公主决定预测一下大富翁国现在所有真钞票的 ...
[Sdoi2008]沙拉公主的困惑
Description 大富翁国因为通货膨胀,以及假钞泛滥,政府决定推出一项新的政策:现有钞票编号范围为1到N的阶乘,但是,政府只发行编号与M!互质的钞票.房地产第一大户沙拉公主决定预测一下大富翁国现 ...
bzoj 2186: [Sdoi2008]沙拉公主的困惑
Time Limit: 10 Sec Memory Limit: 259 MB Submit: 2463 Solved: 820 [Submit][Status][Discuss] Descrip ...
[欧拉函数] Bzoj P2186 沙拉公主的困惑
Description 大富翁国因为通货膨胀,以及假钞泛滥,政府决定推出一项新的政策:现有钞票编号范围为1到N的阶乘,但是,政府只发行编号与M!互质的钞票.房地产第一大户沙拉公主决定预测一下大富翁国现 ...
【codevs2301】【BZOJ2186】沙拉公主的困惑，数论练习之逆元与φ
传送门1 传送门2 写在前面:数论!数论!数论! 思路: 1.分析出题目本意就是求phi(m!)*n!/m!%r(这一步大概挺难想的,想出来这个后面就简单多了) 2.n!%r可以O(10^7)预处理 ...

[bzoj2186]沙拉公主的困惑