BZOJ2186[Sdoi2008]沙拉公主的困惑 线性筛素数
Posted CQzhangyu
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了BZOJ2186[Sdoi2008]沙拉公主的困惑 线性筛素数相关的知识,希望对你有一定的参考价值。
【BZOJ2186】[Sdoi2008]沙拉公主的困惑
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
4 2
Sample Output
1
数据范围:
对于100%的数据,1 < = N , M < = 10000000
数据范围:
对于100%的数据,1 < = N , M < = 10000000
题解:跟欧拉函数的求法类似,就是找出1-m中所有的素数pi,再将n!分别乘上(pi-1)/pi就行了
#include <cstdio> #include <cstring> #include <iostream> using namespace std; typedef long long ll; ll n,m,mod,T,num; int p[10010]; int pri[1000000]; int jc[10000010],jpc[10000010]; bool np[10000010]; ll pm(ll x,ll y) { ll z=1; while(y) { if(y&1) z=z*x%mod; x=x*x%mod,y>>=1; } return z; } int main() { scanf("%d%d",&T,&mod); int i,j; jc[1]=jpc[1]=1; for(i=2;i<=10000000;i++) { jc[i]=1ll*jc[i-1]*i%mod,jpc[i]=jpc[i-1]; if(!np[i]) pri[++num]=i,jpc[i]=1ll*jpc[i]*(i-1)%mod*pm(i,mod-2)%mod; for(j=1;j<=num&&1ll*i*pri[j]<=10000000;j++) { np[i*pri[j]]=1; if(i%pri[j]==0) break; } } while(T--) { scanf("%d%d",&n,&m); printf("%lld\n",1ll*jc[n]*jpc[m]%mod); } return 0; }
以上是关于BZOJ2186[Sdoi2008]沙拉公主的困惑 线性筛素数的主要内容,如果未能解决你的问题,请参考以下文章
BZOJ2186[Sdoi2008]沙拉公主的困惑 线性筛素数