[Sdoi2008]沙拉公主的困惑
来源:互联网 发布:高性能mysql第四版pdf 编辑:程序博客网 时间:2024/05/16 02:05
【问题描述】
大富翁国因为通货膨胀,以及假钞泛滥,政府决定推出一项新的政策:现有钞票编号范围为1到N的阶乘,但是,政府只发行编号与M!互质的钞票。房地产第一大户沙拉公主决定预测一下大富翁国现在所有真钞票的数量。现在,请你帮助沙拉公主解决这个问题,由于可能张数非常大,你只需计算出对R取模后的答案即可。R是一个质数。
【问题分析】
若gcd(a,b)=1,则易证gcd(a+b,b)=1
证明:gcd(a,b)=1 =>a,b没有公共约数 假设gcd(a+b,b)=k(k!=1)
则k是b的约数,又k是(a+b)的约数,因此k是a的约数,与假设矛盾。
所以我们可以讲题目所求即n!以内与m!互质的数的个数转换为
φ(m!)*(n!)/(m!) %p(φ是欧拉函数,φ(x)表示比x小的与x互质的数的个数)化简整理后,(φ(x)=∏(pi-1)/pi pi为x的质因数)即为
n!*∏(pi-1)/pi
所以我们要预处理出质数,阶乘还有逆元。
乘法逆元
定义:对于任意x∈[1,p),若(x,p)=1,则存在唯一的正整数y∈[1,p)满足x*y≡1(mod p),我们称y为x在模p意义下的逆元。
a/x=a*1/x=a*x*inv[x]/x=a*inv[x] (mod p)
这说明,当我们需要除掉x的时候,我们只需要乘上x的逆元就行了。
以上 出自Po神
//强行膜Po神#include <iostream>#include <cstdio>using namespace std;const int N=10000001;int p,t,n,m,tot; int prime[N/2];bool mark[N]; long long rev[N],jie[N],ans[N];void chuli(){ int i,j; for(i=2;i<=N;i++) { if(!mark[i]) prime[++tot]=i; for(j=1;j<=tot;j++) { if(i*prime[j]>N) break; mark[i*prime[j]]=1; if(i%prime[j]==0) break; } } jie[1]=1; for (i=2;i<N;i++) jie[i]=jie[i-1]*i%p; rev[1]=1; for(i=2;i<=N&&i<p;i++) rev[i]=(p-p/i)*rev[p%i]%p; ans[1]=1; for(i=2;i<=N;i++) { if(!mark[i]) ans[i]=ans[i-1]*(i-1)%p*rev[i%p]%p; else ans[i]=ans[i-1]; } return;}int main(){ int i; long long anss; scanf("%d%d",&t,&p); chuli(); for (i=1;i<=t;i++) { scanf("%d%d",&n,&m); anss=jie[n]*ans[m]%p; printf("%lld\n",anss); } return 0;}
0 0
- [Sdoi2008]沙拉公主的困惑
- [SDOI2008]沙拉公主的困惑
- bzoj 2186: [Sdoi2008]沙拉公主的困惑
- 【bzoj 2186】: [Sdoi2008]沙拉公主的困惑
- 【bzoj2186】【sdoi2008】【沙拉公主的困惑】【数论】
- BZOJ2186: [Sdoi2008]沙拉公主的困惑
- bzoj2186: [Sdoi2008]沙拉公主的困惑
- bzoj2186【SDOI2008】沙拉公主的困惑
- 【SDOI2008】【BZOJ2186】沙拉公主的困惑
- bozj2186 [Sdoi2008]沙拉公主的困惑
- [BZOJ2186][Sdoi2008]沙拉公主的困惑
- [Sdoi2008](bzoj2186)沙拉公主的困惑
- 【SDOI2008】【BZOJ2186】沙拉公主的困惑
- 2186: [Sdoi2008]沙拉公主的困惑
- BZOJ 2186 [Sdoi2008]沙拉公主的困惑
- [BZOJ2186][Sdoi2008]沙拉公主的困惑
- bzoj2186: [Sdoi2008]沙拉公主的困惑
- 2186: [Sdoi2008]沙拉公主的困惑
- 使用Oracle绿色客户端(InstantClient)连接远程Oracle的配置方法
- div验证样式
- HTML学习2-符号标签
- PyCUDA教程(二):教程
- 安装visualsvn的注意事项
- [Sdoi2008]沙拉公主的困惑
- 自定义Notification
- imx6q搭建NFS
- Linux安装VMware Tools教程
- JAVA 数组的排序方法冒泡、选择、快速、数组递增
- RabbitMQ 基础教程(1) - Hello World
- RSA加密
- HTML学习0-head
- C++11 in Qt