http://www.cnblogs.com/BLADEVIL/p/3490321.html

http://www.cnblogs.com/zyfzyf/p/3997986.html

翻了翻题解,这两个合起来比较明白……

题意:求1~n!中与m!互质的数的数量(mod R)。

∵由欧几里得算法得gcd(a,b)=gcd(b,a%b)

∴gcd(a+b,b)=gcd(b,(a+b)%b)=gcd(b,a) 即 gcd(a,b)=gcd(a+b,b)

推广:gcd(a,b)=gcd(a+k*b,b)

∴若a与b互质,则a+k*b也与b互质

令n=5,m=4,则n!=120,m!=24

符合题意的数有以下:

5 29 53 77 101

7 31 55 79 103

11 35 59 83 107

13 37 61 85 109

17 41 65 89 113

19 43 67 91 115

23 47 71 95 119

∴若x(x<n!)与m!互质,则x+m!*k也和m!互质,根据我们刚才推出的性质,仔细观察,发现 x 在 mod m!的剩余系中恰好有n!/m!个数符合题意。

∴观察上表的第一列,我们只需考虑小于m!的数中与n!互质的数。

∴这样的数恰好有φ(m!)个

∴答案为φ(m!)*(n!/m!)

但是这样仍然超出了我们能承受的时间复杂度。

∴ans=φ(m!)*(n!/m!)

=m!*((p1-1)/p1)+((p2-1)/p2)+......+((pk-1)/pk)*n!/m!  (p1......pk为m!的质因子,根据阶乘的定义,易证,m!的质因子恰好为<=m的所有质数)

=((p1-1)/p1)+((p2-1)/p2)+......+((pk-1)/pk)*n!

∴我们递推预处理出((p1-1)/p1)+((p2-1)/p2)+......+((pk-1)/pk)和n!即可O(1)地回答每个询问。

在预处理连续乘积的时候,我们要用到快速幂求乘法逆元。

综上,要预处理的是:

①maxm以内的素数

②maxm以内的素数的逆元

③1!~maxn!

④maxm以内的素数的((p1-1)/p1)+((p2-1)/p2)+......+((pk-1)/pk)

 #include<cstdio>
#include<algorithm>
#include<iostream>
using namespace std;
typedef long long ll;
ll MOD;
int T,n[],m[],maxn,maxm,Ni[],Fac[],Muls[];
bool unPrime[];
void Shai_Prime()
{
unPrime[]=true;
for(ll i=;i<=maxm;i++)
for(ll j=i*i;j<=maxm;j+=i) unPrime[j]=true;
}
ll pow_mod(ll a,ll p,ll MOD)
{
if(!p) return ;
ll ans=pow_mod(a,p>>,MOD);
ans=ans*ans%MOD;
if((p&)==) ans=ans*a%MOD;
return ans;
}
void Init_Ni()
{
for(int i=;i<=maxm;i++)
if(!unPrime[i])
Ni[i]=(int)pow_mod(i,MOD-,MOD);
}
void Init_Fac()
{
Fac[]=;
for(int i=;i<=maxn;i++) Fac[i]=(int)((ll)Fac[i-]*(ll)i%MOD);
}
void Init_Muls()
{
ll res=;
for(int i=;i<=maxm;i++)
{
if(!unPrime[i]) res=((res*(ll)(i-))%MOD)*(ll)Ni[i]%MOD;
Muls[i]=res;
}
}
int main()
{
scanf("%d",&T); cin>>MOD;
for(int i=;i<=T;i++)
{
scanf("%d%d",&n[i],&m[i]);
maxn=max(n[i],maxn);
maxm=max(m[i],maxm);
}
Shai_Prime(); Init_Ni(); Init_Fac(); Init_Muls();
for(int i=;i<=T;i++)
printf("%d\n",(int)((ll)Muls[m[i]]*(ll)Fac[n[i]]%MOD));
return ;
}

【数论】【欧拉函数】【筛法求素数】【乘法逆元】【快速幂取模】bzoj2186 [Sdoi2008]沙拉公主的困惑的更多相关文章

  1. 数学--数论--HDU 4675 GCD of Sequence(莫比乌斯反演+卢卡斯定理求组合数+乘法逆元+快速幂取模)

    先放知识点: 莫比乌斯反演 卢卡斯定理求组合数 乘法逆元 快速幂取模 GCD of Sequence Alice is playing a game with Bob. Alice shows N i ...

  2. POJ 1845 Sumdiv [素数分解 快速幂取模 二分求和等比数列]

    传送门:http://poj.org/problem?id=1845 大致题意: 求A^B的所有约数(即因子)之和,并对其取模 9901再输出. 解题基础: 1) 整数的唯一分解定理: 任意正整数都有 ...

  3. 九度OJ 1085 求root(N, k) -- 二分求幂及快速幂取模

    题目地址:http://ac.jobdu.com/problem.php?pid=1085 题目描述: N<k时,root(N,k) = N,否则,root(N,k) = root(N',k). ...

  4. The 2018 ACM-ICPC China JiangSu Provincial Programming Contest快速幂取模及求逆元

    题目来源 The 2018 ACM-ICPC China JiangSu Provincial Programming Contest 35.4% 1000ms 65536K Persona5 Per ...

  5. 二分求幂/快速幂取模运算——root(N,k)

    二分求幂 int getMi(int a,int b) { ; ) { //当二进制位k位为1时,需要累乘a的2^k次方,然后用ans保存 == ) { ans *= a; } a *= a; b / ...

  6. BZOJ2186 [Sdoi2008]沙拉公主的困惑 【数论,欧拉函数,线性筛,乘法逆元】

    2186: [Sdoi2008]沙拉公主的困惑 Time Limit: 10 Sec  Memory Limit: 259 MB Submit: 5003  Solved: 1725 [Submit] ...

  7. UVA 10006 - Carmichael Numbers 数论(快速幂取模 + 筛法求素数)

      Carmichael Numbers  An important topic nowadays in computer science is cryptography. Some people e ...

  8. BZOJ2186: [Sdoi2008]沙拉公主的困惑(求[1,N!]与M!互素的个数)(线性筛)

    2186: [Sdoi2008]沙拉公主的困惑 Time Limit: 10 Sec  Memory Limit: 259 MBSubmit: 6103  Solved: 2060[Submit][S ...

  9. 欧拉函数phic以及超大数的快速幂

    题目:求a^b*c%mod; 其中b<=10^100000; 是不是很大..... /*当你要计算 A^B%C的时候 因为此题中的B很大,达到10^100000,所以我们应该联想到降幂公式. 降 ...

随机推荐

  1. CSS3中transform属性的用法

    有时候网站也要愚弄一下访客,比如愚人节.下面我给大家推荐个效果,就是整个页面左右颠倒了.css3 很强大,简单的几行代码就可以帮我们实现这个效果. view source   print? 01 &l ...

  2. JavaScript的大括号的语义

    Javascript中大括号"{}"有四种语义作用: 语义1. 组织复合语句,这是最常见的: view source   print? 1 if( condition ) { 2 ...

  3. Dilworth定理证明

    命题:偏序集能划分成的最少的全序集的个数与最大反链的元素个数相等. (离散数学结构第六版课本P245:把一个偏序集划分成具有全序的子集所需要的最少子集个数与元素在偏序下都是不可比的最大集合的基数之间有 ...

  4. [bzoj 2844]线性基+高斯消元

    题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=2844 又用到线性基+高斯消元的套路题了,因为经过高斯消元以后的线性基有非常好的序关系,所以 ...

  5. UVA10480:Sabotage(最小割+输出)

    Sabotage 题目链接:https://vjudge.net/problem/UVA-10480 Description: The regime of a small but wealthy di ...

  6. 设置eclipse控制台上的信息输入到某个文件

    转摘自:http://binary.duapp.com/2013/09/1511.html Run->Run Configurations->Common->File

  7. Windows下安装Mycat-web

    Mycat-web是基于Mycat的一个性能监控工具,如:sql性能监控等. 在安装Mycat-web之前需要先安装Zookeeper: 可参考: http://blog.csdn.net/tlk20 ...

  8. Struts2 内建的验证规则

    Struts2 内建的验证规则 conversion validator:转换验证器 date validator:日期验证器 double validator:浮点验证器 email validat ...

  9. RPC-Thrift(四)

    Client Thrift客户端有两种:同步客户端和异步客户端. 同步客户端 同步客户端比较简单,以RPC-Thrift(一)中的的例子为基础进行研究源码,先看一下类图. TServiceClient ...

  10. Cpython解释器支持的线程

    因为Python解释器帮你自动定期进行内存回收,你可以理解为python解释器里有一个独立的线程,每过一段时间它起wake up做一次全局轮询看看哪些内存数据是可以被清空的,此时你自己的程序 里的线程 ...