[bzoj2186]: 莎拉公主的困惑(欧拉函数+乘法逆元)
来源:互联网 发布:vip解析接口源码 编辑:程序博客网 时间:2024/04/25 22:36
2186: [Sdoi2008]沙拉公主的困惑
Time Limit: 10 Sec Memory Limit: 259 MBSubmit: 4233 Solved: 1476
[Submit][Status][Discuss]
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
4 2
Sample Output
数据范围:
对于100%的数据,1 < = N , M < = 10000000
HINT
Source
题解
看到此题觉得好难啊 其实也好难 于是乎去翻题解 找了好多都没有解释清楚 终于找到下面这个写的不错
参考http://blog.csdn.net/u013486414/article/details/48878353
题意:求1~n!中与m!互质的数的个数,且m<=n(1 < = n , m < = 10000000)
思路:因为m<=n,所以m!能被n!整除,那么我们很容易推出下面结论:
对于两个正整数n和m,,如果n是m的倍数,那么1~n中与m互质的数的个数为n/m*Eular(m)。
因为欧拉函数的通式为:φ(x)=x(1-1/p1)(1-1/p2)(1-1/p3)(1-1/p4)…..(1-1/pn)。所以对于本题目来说,答案如下:
其中pi为<=m的所有素数,因为对于1-1/pi=(pi-1)/pi=(pi-1)*inv(pi),所以这里要用到逆元,对于1~p中模p的逆元,如果p很大,那么有一个递推式可以求逆元:inv[i]=(m-m/i)*inv[m%i]%m,初始化inv[1]=1。
以上引用
但是本人不会这个求逆元的递推式 用扩展欧几里得做也同样可以求出逆元
此题的关键是预处理出 阶乘 逆元 和质数 和ans=(1-1/p1)(1-1/p2).....(1-1/pk)
//参考:http://blog.csdn.net/u013486414/article/details/48878353#include<iostream>#include<cstdio>#include<cstdlib>#include<cstring>#include<algorithm>#include<cmath>#include<queue>#define N 10000000#define maxk 1000000#define ll long longusing namespace std;ll prime[maxk],ine[N+100],fac[N+100],tot,R,ans[N+100];bool flag[N+100];int T,n,m;int x,y;void exgcd(int a,int b){if(b==0) {x=1;y=0;}else {exgcd(b,a%b);int t=x;x=y;y=t-(a/b)*y;}}void pre(){fac[1]=1; for(int i=2;i<=N;i++)fac[i]=fac[i-1]*i%R;for(int i=2;i<=N;i++){if(!flag[i]){prime[++tot]=i;exgcd(i,R);ine[i]=(x+R)%R;}for(int j=1;j<=tot;j++){int to=i*prime[j];if(to>N) break;flag[to]=1;if(i%prime[j]==0) break;}}ans[1]=1;for(int i=2;i<=N;i++){ans[i]=ans[i-1];//不是素数就传递 if(!flag[i])ans[i]=(ll)ans[i]*(i-1)%R*ine[i]%R;}}int main(){scanf("%d%d",&T,&R);pre(); while(T--){scanf("%d%d",&n,&m); printf("%lld\n",fac[n]*ans[m]%R);}return 0;}
- [bzoj2186]: 莎拉公主的困惑(欧拉函数+乘法逆元)
- bzoj2186 莎拉公主的困惑【线性筛+欧拉函数+逆元】
- Bzoj2186:[Sdoi2008]沙拉公主的困惑:欧拉函数+乘法逆元
- bzoj2186[沙拉公主的困惑] 欧拉函数 线性筛 乘法逆元
- BZOJ2186 [Sdoi2008]沙拉公主的困惑 【数论,欧拉函数,线性筛,乘法逆元】
- [bzoj2186]沙拉公主的困惑 欧拉函数+逆元
- [BZOJ2186][SDOI2008]沙拉公主的困惑 欧拉函数
- 2186: [Sdoi2008]沙拉公主的困惑 线性筛素数+欧拉函数+乘法逆元
- 欧拉函数+逆元 [Sdoi2008]沙拉公主的困惑
- bzoj2186: [Sdoi2008]沙拉公主的困惑 逆元
- bzoj2186(欧拉函数+逆元)
- BZOJ 2186 沙拉公主的困惑 逆元+素数筛+欧拉函数
- BZOJ 2186: [Sdoi2008]沙拉公主的困惑 [欧拉函数][逆元]
- BZOJ2186沙拉公主的困惑
- bzoj2186沙拉公主的困惑
- 【Bzoj2186】沙拉公主的困惑
- BZOJ2186 [Sdoi2008]沙拉公主的困惑 数论:递推求逆元
- 【codevs2301】【BZOJ2186】沙拉公主的困惑,数论练习之逆元与φ
- 求两个链表的第一个公共子结点
- 日志组件log4j详解
- 拥塞控制
- 修改UISearchBar的TextField的背景颜色
- PHP类型
- [bzoj2186]: 莎拉公主的困惑(欧拉函数+乘法逆元)
- MaterialDesign(一)
- Hibernate批量处理操作(批量查询 批量删除)
- Activity的启动模式(二)
- java基础总结(2)----语言基础
- 面向对象的选项卡
- 九度1078:二叉树遍历(已经中序,先序遍历求后序遍历)
- shell脚本
- python笔记1_标识符_保留字_注释_字符串