[BZOJ2186][SDOI2008]沙拉公主的困惑 欧拉函数
来源:互联网 发布:linux c一站式编程pdf 编辑:程序博客网 时间:2024/04/30 01:09
为什么一开始我再用容斥+莫比乌斯函数乱推。。。
答案即为φ(m!)*n!/m!,考虑gcd(x,y)==gcd(x+y,y),所以关于某个数x互质的数是以x为长度的循环节,因为m!|n!,所以恰好n!/m!个完整的循环节。
预处理阶乘,和φ(m!)*m!(用φ的公式计算),其实后者也可O(n)预处理。
代码:
#include<iostream>#include<cstdio>#define ll long longusing namespace std;int mod,ca,num,pri[1000010],n,m;ll a[10000010],jie[10000010];ll ksm(ll a,int b){ll re=1; while(b){if(b&1)re=re*a%mod;a=a*a%mod;b>>=1;} return re;}void getpri(int n){ a[1]=1; for(int i=2;i<=n;i++) { if(!a[i]) pri[++num]=i; for(int j=1;j<=num&&i*pri[j]<=n;j++) { a[i*pri[j]]=1; if(i%pri[j]==0) break; } }}int main(){ scanf("%d%d",&ca,&mod); getpri(1e7); a[0]=jie[0]=1; for(int i=1;i<=1e7;i++) { if(a[i]==1) a[i]=a[i-1]; else a[i]=a[i-1]*(i-1)%mod*ksm(i,mod-2)%mod; jie[i]=jie[i-1]*i%mod; } while(ca--) { scanf("%d%d",&n,&m); printf("%lld\n",jie[n]*a[m]%mod); } return 0; }
阅读全文
0 0
- [BZOJ2186][SDOI2008]沙拉公主的困惑 欧拉函数
- Bzoj2186:[Sdoi2008]沙拉公主的困惑:欧拉函数+乘法逆元
- BZOJ2186 [Sdoi2008]沙拉公主的困惑 【数论,欧拉函数,线性筛,乘法逆元】
- [bzoj2186]沙拉公主的困惑 欧拉函数+逆元
- 【bzoj2186】【sdoi2008】【沙拉公主的困惑】【数论】
- BZOJ2186: [Sdoi2008]沙拉公主的困惑
- bzoj2186: [Sdoi2008]沙拉公主的困惑
- bzoj2186【SDOI2008】沙拉公主的困惑
- 【SDOI2008】【BZOJ2186】沙拉公主的困惑
- [BZOJ2186][Sdoi2008]沙拉公主的困惑
- [Sdoi2008](bzoj2186)沙拉公主的困惑
- 【SDOI2008】【BZOJ2186】沙拉公主的困惑
- [BZOJ2186][Sdoi2008]沙拉公主的困惑
- bzoj2186: [Sdoi2008]沙拉公主的困惑
- BZOJ2186:[Sdoi2008]沙拉公主的困惑
- bzoj2186 [Sdoi2008]沙拉公主的困惑
- bzoj2186: [Sdoi2008]沙拉公主的困惑
- BZOJ2186 [Sdoi2008]沙拉公主的困惑
- git看不到别人创建的远程分支
- spring事务管理几种方式
- 第3课:在IDE下开发一个Scala程序,透彻解析及Scala控制结构
- 【Mono】C# + Mono + Linux 编程初体验之 HelloWorld
- java设计模式之代理模式
- [BZOJ2186][SDOI2008]沙拉公主的困惑 欧拉函数
- 基于日志的同步数据一致性和实时抽取
- Django自定义标签实现多级评论
- 图解25匹马的选马问题
- pandas把所有大于0的数设置为1
- CF 454B Little Pony and Sort by Shift KMP算法简单运用
- vector扩容
- python爬虫OS X环境搭建
- eclipse clean 之后无法运行的问题