【SDOI2008】【BZOJ2186】沙拉公主的困惑
2016-02-20 14:33
399 查看
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
Sample Output
1
数据范围:
对于100%的数据,1 < = N , M < = 10000000
HINT
Source
答案是ϕ(m!)∗n!/m!
考察线性求逆元(以前的数论板子现在也有用了..
大富翁国因为通货膨胀,以及假钞泛滥,政府决定推出一项新的政策:现有钞票编号范围为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
Sample Output
1
数据范围:
对于100%的数据,1 < = N , M < = 10000000
HINT
Source
答案是ϕ(m!)∗n!/m!
考察线性求逆元(以前的数论板子现在也有用了..
#include<iostream> #include<cstdio> #include<cstring> #include<cmath> #include<algorithm> #define GET (ch>='0'&&ch<='9') #define MAXN 10000010 using namespace std; int T,R,n,m; bool not_prime[MAXN]; int prime[MAXN],top; int fac[MAXN],inv[MAXN],ret[MAXN]; void in(int &x) { char ch=getchar();x=0; while (!GET) ch=getchar(); while (GET) x=x*10+ch-'0',ch=getchar(); } void init() { fac[1]=1;inv[1]=1;ret[1]=1; for (int i=2;i<=MAXN-10;i++) fac[i]=1ll*fac[i-1]*i%R; for (int i=2;i<=MAXN-10;i++) { if (!not_prime[i]) prime[++top]=i; for (int j=1;j<=top&&i*prime[j]<=MAXN-10;j++) { not_prime[i*prime[j]]=1; if (i%prime[j]==0) break; } } for (int i=2;i<=MAXN-10&&i<R;i++) inv[i]=1ll*(R-R/i)*inv[R%i]%R; for (int i=2;i<=MAXN-10;i++) { ret[i]=ret[i-1]; if (!not_prime[i]) ret[i]=1ll*ret[i]*(i-1)%R*inv[i]%R%R; } } int main() { for (in(T),in(R),init();T;T--) in(n),in(m),printf("%d\n",(1ll*fac *ret[m]+R)%R); }
相关文章推荐
- 1.m分解阶乘之和
- 2.几种递推数
- 3.欧拉函数
- 4.快速幂模m算法
- 5.扩展欧几里得&&中国剩余定理
- 6.数论_web
- Project Euler Problem 387 - Harshad Numbers - 深度优先
- 编程之美2015初赛A
- 数论题集
- 原根
- 阶与原根学习笔记
- HDU 1299 Diophantus of Alexandria
- Leftmost Digit(HDU 1060)
- Rightmost Digit(HDU 1061)
- Python-在奇数中寻找素数
- ZOJ 2674 Strange Limit 欧拉定理
- LeetCode-Palindrome Number
- 欧拉函数
- 组合数求模总结
- 【数论】组合数求模