BZOJ2186 [Sdoi2008]沙拉公主的困惑 【数论,欧拉函数,线性筛,乘法逆元】
来源:互联网 发布:cf手游龙爪技能数据 编辑:程序博客网 时间:2024/04/25 19:23
2186: [Sdoi2008]沙拉公主的困惑
Time Limit: 10 Sec Memory Limit: 259 MBSubmit: 5003 Solved: 1725
[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
1 11
4 2
4 2
Sample Output
1
数据范围:
对于100%的数据,1 < = N , M < = 10000000
数据范围:
对于100%的数据,1 < = N , M < = 10000000
首先我们由gcd的性质知道互质的数的个数以M!为周期循环
即gcd(a,b) = gcd(a,b + a)
所以我们只要求出M!以内与N!互质的个数,乘上N!/M!就可以了
ans = M! phi(M!) * N! / M! = N! * phi(M!)
M!包含M以内的所有质数,所以我们用线性筛求出所有质数就可以预处理出所有phi
求phi时用到除法,可以用逆元递推预处理:inv[i] = P - P / i * inv[P % i] % P
8s压线过QAQ
#include<cstdio>#include<bitset>#include<cstring>#include<algorithm>#define LL long long int#define REP(i,n) for (register long long int i = 1; i <= (n); i++)using namespace std;const int maxn = 10000005,Max = 10000000,maxm = 2000000,INF = 1000000000;int fac[maxn],prime[maxm],inv[maxn],sum[maxn],primei = 0,P,N,M,ans;bitset<maxn> isn;void Pre(){fac[0] = 1;REP(i,Max) fac[i] = (LL)fac[i - 1] * i % P;inv[1] = 1;for (int i = 2 ; i <= Max && i <= P; i++) inv[i] = P - (LL)(P / i) * inv[P % i] % P;sum[0] = 1;for (LL i = 2; i <= Max; i++){if (!isn[i]) prime[++primei] = i;for (int j = 1; j <= primei && (LL)i * prime[j] <= Max; j++){isn[i * prime[j]] = true;if (i % prime[j] == 0) break;}}sum[1] = 1;for (LL i = 2; i <= Max; i++){sum[i] = sum[i - 1];if (!isn[i]) sum[i] = (LL)sum[i] * (i - 1) % P * inv[i] % P;}}int main(){int T;scanf("%d%d",&T,&P);Pre();while (T--){scanf("%d%d",&N,&M);printf("%lld\n",(LL)fac[N] * sum[M] % P);}return 0;}
阅读全文
0 0
- BZOJ2186 [Sdoi2008]沙拉公主的困惑 【数论,欧拉函数,线性筛,乘法逆元】
- bzoj2186[沙拉公主的困惑] 欧拉函数 线性筛 乘法逆元
- Bzoj2186:[Sdoi2008]沙拉公主的困惑:欧拉函数+乘法逆元
- 2186: [Sdoi2008]沙拉公主的困惑 线性筛素数+欧拉函数+乘法逆元
- [bzoj2186]沙拉公主的困惑 欧拉函数+逆元
- bzoj2186【Sdoi2008】沙拉公主的困惑 数论:欧拉筛法,线性逆元
- [BZOJ2186][SDOI2008]沙拉公主的困惑 欧拉函数
- 【bzoj2186】[Sdoi2008]沙拉公主的困惑 数论 线性筛逆元
- 欧拉函数+逆元 [Sdoi2008]沙拉公主的困惑
- 【bzoj2186】【sdoi2008】【沙拉公主的困惑】【数论】
- bzoj2186: [Sdoi2008]沙拉公主的困惑 逆元
- BZOJ2186 [Sdoi2008]沙拉公主的困惑 数论:递推求逆元
- [SDOI2008]沙拉公主的困惑 线性筛 素数+欧拉
- bzoj2186 莎拉公主的困惑【线性筛+欧拉函数+逆元】
- 【SDOI2008】【BZOJ2186】【沙拉公主的困惑】【题解】【数论】
- [BZOJ2186][Sdoi2008]沙拉公主的困惑(数论)
- [bzoj2186][Sdoi2008]沙拉公主的困惑(数论)
- [BZOJ2186][SDOI2008]沙拉公主的困惑(数论)
- 《Linux内核设计与实现》读书笔记(四)- 进程的调度
- 面向对象--简单多态例子
- iOS-92-获取网络图片尺寸
- 项目部署:同一个ip,不同端口号,session失效
- 清澄A1026. 字符统计
- BZOJ2186 [Sdoi2008]沙拉公主的困惑 【数论,欧拉函数,线性筛,乘法逆元】
- js里==和===有什么区别
- [Office] 稿纸设置
- python接口自动化1-发送get请求
- VMware虚拟机的安装
- MVP设计模式
- mongdb数据库学习
- 解决VM中安装CentOS 7.4后无法联网的问题
- Linux DISPLAY 变量设置