大视野2186: 沙拉公主的困惑(求逆元)
2016-04-09 13:00
309 查看
2186: [Sdoi2008]沙拉公主的困惑
Time Limit: 10 Sec Memory Limit: 259 MBSubmit: 2616 Solved: 880
[Submit][Status][Discuss]
Description
大富翁国因为通货膨胀,以及假钞泛滥,政府决定推出一项新的政策:现有钞票编号范围为1到N的阶乘,但是,政府只发行编号与M!互质的钞票。房地产第一大户沙拉公主决定预测一下大富翁国现在所有真钞票的数量。现在,请你帮助沙拉公主解决这个问题,由于可能张数非常大,你只需计算出对R取模后的答案即可。R是一个质数。Input
第一行为两个整数T,R。R<=10^9+10,T<=10000,表示该组中测试数据数目,R为模后面T行,每行一对整数N,M,见题目描述 m<=nOutput
共T行,对于每一对N,M,输出1至N!中与M!素质的数的数量对R取模后的值Sample Input
1 114 2
Sample Output
1数据范围:
对于100%的数据,1 < = N , M < = 10000000
HINT
Source
[Submit][Status][Discuss]
HOME Back
题意:求
中
互质的数的个数,其中
。
分析:因为
,所以
,我们很容易知道如下结论
对于两个正整数
和
,如果
是
的倍数,那么
中与
互素的数的个数为
本结论是很好证明的,因为
中与
互素的个数为
,又知道
,所以
结论成立。那么对于本题,答案就是
其中
为小于等于
的所有素数,先筛选出来即可。由于最终答案对一个质数取模,所以要用逆元,这里
求逆元就有技巧了,用刚刚介绍的递推法预处理,否则会TLE的。
代码:
#include <iostream> #include <string.h> #include <stdio.h> #include <bitset> using namespace std; typedef long long LL; const int N = 10000005; bitset<N> prime; void isprime() { prime.set(); for(int i=2; i<N; i++) { if(prime[i]) { for(int j=i+i; j<N; j+=i) prime[j] = false; } } } LL ans1 ,ans2 ; LL inv ; int main() { isprime(); int MOD,m,n,T; scanf("%d%d",&T,&MOD); ans1[0] = 1; for(int i=1; i<N; i++) ans1[i] = ans1[i-1] * i % MOD; inv[1] = 1; for(int i=2;i<N;i++) { if(i >= MOD) break; inv[i] = (MOD - MOD / i) * inv[MOD % i] % MOD; } ans2[1] = 1; for(int i=2; i<N; i++) { if(prime[i]) { ans2[i] = ans2[i-1] * (i - 1) % MOD; ans2[i] = ans2[i] * inv[i % MOD] % MOD; } else { ans2[i] = ans2[i-1]; } } while(T--) { scanf("%d%d",&n,&m); LL ans = ans1 * ans2[m] % MOD; printf("%lld\n",ans); } return 0; }
接下来还有一个关于逆元的有意思的题目,描述如下
证明:由
其中
所以只需要证明
,而我们知道
模
的逆元对应全部
中的所有数,既是单射也是满射。
所以进一步得到
证明完毕!
相关文章推荐
- IDEA插件配置之Eclipse Code Formatter
- [算法学习]20150409.1.插入排序
- Scrapy设置之Performance概览
- LeetCode 17 Letter Combinations of a Phone Number
- CSP最大的矩形
- 结对编程项目——电梯调度
- java动态数组 2.0
- Java之XMemcached使用及源码详解
- 将单向链表按某值划分成左边小、中间相等、右边大的形式
- B树
- 最快速建立Android项目的JavaBean实体类
- HDFS详解
- 关于hibernate的 Disabling contextual LOB creation as JDBC driver reported JDBC version [3] less than 4
- java动态数组
- getBackground().setAlpha(0)设置透明度连下层控件背景色都被改变
- 第四次程序设计作业 C++计算器计算及命令行的使用 前缀表达式方法实现
- [fluentd学习]安装使用
- Couchbase:更好的Cache系统
- 257. Binary Tree Paths
- java中volatile关键字的含义