[模板] 数学基础:快速幂/乘/逆元/exGCD/(ex)CRT/(ex)Lucas定理
方便复制
快速乘/幂
时间复杂度 \(O(\log n)\).
ll nmod;
//快速乘
ll qmul(ll a,ll b){
ll l=a*(b>>hb)%nmod*(1ll<<hb)%nmod;
ll r=a*(b&((1<<hb)-1))%nmod;
return (l+r)%nmod;
}
//快速幂
ll qpow(ll a,ll b){
ll res=1;
while(b){
if(b&1)res=res*a%nmod;
a=a*a%nmod;
b>>=1;
}
return res;
}
exgcd
内容
解不定方程 $ ax+by = c $
时间复杂度 \(O(\log n)\).
void exgcd(ll a,ll b,ll& x,ll& y,ll& d){ //a&b should > 0
b==0?(x=1,y=0,d=a):(exgcd(b,a%b,y,x,d),y-=x*(a/b));
}
//use
ll a,b,c;
ll m,x,y;
exgcd(a,b,x,y,m);
if(c%m!=0)cout<<"No\n"; //无解
else{
c/=m,x*=c,y*=c,a/=m,b/=m;
//令x取最小非负整数
x1=x%b;if(x1<0)x1+=abs(b);
y1=(c-a*x1)/b;
//令y取最小非负整数
y1=y%a;if(y1<0)y1+=abs(a);//(y>=0?y%a:y%a+abs(a));
x1=(c-b*y1)/a;
}
逆元
内容
求\(n * x \equiv 1 (mod m)\) 最小正整数解.
单个数
时间复杂度 \(O(\log n)\).
//1:qp(n,nmod-2)
//2
ll getv(ll n){return n<0?n+nmod:n;}
ll inv(ll n){
ll x,y,d;
exgcd(getv(n%nmod),p,x,y,d);
return x%p+(x<0?p:0);
}
线性求逆元
时间复杂度 \(O(n)\).
- 公式法(并不能记住板子)
- 阶乘法
利用下面的公式:
\]
\]
代码
ll fac[nsz],ifac[nsz];
void init(int bnd){
fac[0]=1;
rep(i,1,bnd)fac[i]=i*fac[i-1]%nmod;
ifac[bnd]=inv(fac[bnd]);
repdo(i,bnd-1,0){
ifac[i]=ifac[i+1]*(i+1)%nmod;
}
}
CRT
中国剩余定理 && 扩展中国剩余定理 - niiick - CSDN博客
内容
解线性同余方程组 \(x \equiv a_i \pmod{m_i}, \forall i \in \{1, 2, \cdots , n\}\). 其中\(m_i\)两两互质.
设\(M=\prod_{i=1}^nm_i\), \(M_i=\frac M{m_i}\);
\(M_i^{-1}\) 为 \(M_i\) 关于 \(\bmod m_i\)的逆元,
则可以构造出通解
时间复杂度 $O(n \cdot 逆元)$, 通常为 $O(n \log n)$.
### 代码
```
ll crt(ll a,ll m,ll m0){//m0 | m; (m0,m/m0)=1
return m/m0*inv(m/m0,m0)%m*a%m;
}
```
## excrt
模数不互质.
利用合并的思想求解.
时间复杂度$O(n \log n)$.
### 代码
```
ll excrt(ll *a,ll *m,ll n){
ll a0=a[1],m0=m[1],x,y,g;
rep(i,2,n){
g=exgcd(m0,m[i],x,y);
if((a[i]-a0)%g!=0)return -1;
x=(a[i]-a0)/g*x%(m[i]/g);
a0+=x*m0;
m0=m0/g*m[i];
a0%=m0;
}
return a0<0?(a0%m0+m0):(a0%m0);
}
```
## Lucas定理
[Lucas定理 - permui - 博客园](https://www.cnblogs.com/owenyu/p/6724560.html)
### 内容
求 $\binom n m \bmod p$, 保证$p \in \{ prime \}$.
设$n=(a_0a_1\dots a_k)_p$, $m=(b_0b_1\dots b_k)_p$, 有
$$ \binom n m\equiv \prod _{i=0}^k\binom {a_i} {b_i} \pmod p \]
也即
\]
递归求解.
时间复杂度 \(O(p \log_p n)\) , 或者 \(O(p)\) 预处理, \(O(\log_p n)\) 单次询问.
代码
ll c(ll n,ll m){
if(m<0||m>n)return 0;
return fact(n)*inv(fact(m)*fact(n-m)%p)%p;
}
ll lucas(int n,int m){//c(n,m)%p
return m?lucas(n/p,m/p)*c(n%p,m%p)%p:1;
}
exLucas
【知识总结】扩展卢卡斯定理(exLucas) - Inspector_Javert - CSDN博客
p为合数.
分解质因数+阶乘取模+组合数+excrt
注意如果计算 \(n! \bmod p^k\) 时如果计算 \(p^x\) 对结果的贡献, 将无法求逆元. 因此需要求 \(\frac{n!}{p^x} \bmod p^k\), 即忽略p的幂, 然后在求组合数时再乘回来.
码量++
复杂度太长... 当它是 \(O(p \log p)\) 好惹
代码
ll qpow(ll a,ll b,ll nmod){
ll res=1;
while(b){
if(b&1)res=res*a%nmod;
a=a*a%nmod;
b>>=1;
}
return res;
}
void exgcd(ll a,ll b,ll &x,ll &y,ll &d){
b?(exgcd(b,a%b,y,x,d),y-=x*(a/b)):(x=1,y=0,d=a);
}
ll inv(ll a,ll m){
ll x,y,d;
exgcd(a,m,x,y,d);
return x>=0?x%m:x%m+m;
}
ll crt(ll a,ll m,ll m0){//m0 | m; (m0,m/m0)=1
return m/m0*inv(m/m0,m0)%m*a%m;
}
ll fact(ll n,ll p,ll pk){//(n!/p^x)%(p^k)
if(n<=1)return 1;
ll ans=1,tmp=n%pk;
rep(i,1,pk){
if(i%p)ans=ans*i%pk;
}
ans=qpow(ans,n/pk,pk);
rep(i,1,tmp){
if(i%p)ans=ans*i%pk;
}
return ans*fact(n/p,p,pk)%pk;
}
ll c(ll n,ll m,ll p,ll pk){//c(n,m)%(p^k)
ll sum=0;
for(ll i=n;i;i/=p)sum+=i/p;
for(ll i=m;i;i/=p)sum-=i/p;
for(ll i=n-m;i;i/=p)sum-=i/p;
return qpow(p,sum,pk)*fact(n,p,pk)%pk*inv(fact(m,p,pk),pk)%pk*inv(fact(n-m,p,pk),pk)%pk;
}
ll fac[40][2],pf; //0 p; 1 pk
void getfac(ll n){
ll tmp=sqrt(n);
rep(i,2,tmp){
if(n%i==0){
fac[++pf][0]=i,fac[pf][1]=1;
while(n%i==0)n/=i,fac[pf][1]*=i;
}
}
if(n>1)fac[++pf][0]=n,fac[pf][1]=n;
}
ll exlucas(ll n,ll m,ll p){
ll ans=0;
getfac(p);
rep(i,1,pf){
ans=(ans+crt(c(n,m,fac[i][0],fac[i][1]),p,fac[i][1]))%p;
}
return ans;
}
[模板] 数学基础:快速幂/乘/逆元/exGCD/(ex)CRT/(ex)Lucas定理的更多相关文章
- HDU4869:Turn the pokers(快速幂求逆元+组合数)
题意: 给出n次翻转和m张牌,牌相同且一开始背面向上,输入n个数xi,表示xi张牌翻转,问最后得到的牌的情况的总数. 思路: 首先我们可以假设一开始牌背面状态为0,正面则为1,最后即是求ΣC(m,k) ...
- 洛谷 P1226 【模板】快速幂||取余运算
题目链接 https://www.luogu.org/problemnew/show/P1226 题目描述 输入b,p,k的值,求b^p mod k的值.其中b,p,k*k为长整型数. 输入输出格式 ...
- NOIP2011多项式系数[快速幂|组合数|逆元]
题目描述 给定一个多项式(by+ax)^k,请求出多项式展开后x^n*y^m 项的系数. 输入输出格式 输入格式: 输入文件名为factor.in. 共一行,包含5 个整数,分别为 a ,b ,k , ...
- ACM-ICPC 2018 焦作赛区网络预赛- L:Poor God Water(BM模板/矩阵快速幂)
God Water likes to eat meat, fish and chocolate very much, but unfortunately, the doctor tells him t ...
- XTU 1260 - Determinant - [2017湘潭邀请赛A题(江苏省赛)][高斯消元法][快速幂和逆元]
是2017江苏省赛的第一题,当时在场上没做出来(废话,那个时候又不懂高斯消元怎么写……而且数论也学得一塌糊涂,现在回来补了) 省赛结束之后,题解pdf就出来了,一看题解,嗯……加一行再求逆矩阵从而得到 ...
- Codeforces1106F 【BSGS】【矩阵快速幂】【exgcd】
首先矩阵快速幂可以算出来第k项的指数,然后可以利用原根的性质,用bsgs和exgcd把答案解出来 #include<bits/stdc++.h> using namespace std; ...
- hdu-4990 Reading comprehension(快速幂+乘法逆元)
题目链接: Reading comprehension Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K ( ...
- 洛谷 P1226 【模板】快速幂||取余运算 题解
Analysis 快速幂模板,注意在最后输出时也要取模. 快速幂模板 inline ll ksm(ll x,ll y) { ll ans=; ) { ) { ans*=x; ans%=k; } x*= ...
- Codeforces 963 A. Alternating Sum(快速幂,逆元)
Codeforces 963 A. Alternating Sum 题目大意:给出一组长度为n+1且元素为1或者-1的数组S(0~n),数组每k个元素为一周期,保证n+1可以被k整除.给a和b,求对1 ...
随机推荐
- Asp.Net MVC @Html.TextBox 只允许输入数字问题
转到Html.TextBox()看可以看出 Html.TextBox()方法,创建文本框[<input type="text"/>],并且可以带上name,value和 ...
- Dynamics 365 POA表记录的查询
微软动态CRM专家罗勇 ,回复313或者20190311可方便获取本文,同时可以在第一间得到我发布的最新博文信息,follow me!我的网站是 www.luoyong.me . PrincipalO ...
- 驰骋工作流引擎JFlow与activiti的对比之2种结构化模式
1. 任意循环(Arbitrary Cycles) ACTIVITI : 某一个或多个活动可以反复执行. 例子:用户买了瓶汽水,拿到汽水后,中了一瓶,又去兑换了一瓶汽水,如果又中了,再去兑换一瓶汽水- ...
- WPF:间接支持虚拟化的ListBox
/// <summary> /// 间接实现了虚拟化的ListBox /// 子项必须实现IVisible接口 /// 你可以在IsVisible发生改变时实现一系列自定义动作 /// 比 ...
- 南京邮电大学java程序设计作业在线编程第八次作业
程序设计类课程作业平台 王利国 主页 教学资源 我的作业列表 程序设计课 账户 王利国的"Java语言程序设计第8次作业(2018)"详细 主页 我的作业列表 作业结果详细 总分: ...
- 程序员利器Tmux使用手册
转:https://blog.csdn.net/chenqiuge1984/article/details/80132042
- 【原】Java学习笔记003 - 数据类型
package cn.temptation; public class Sample01 { public static void main(String[] args) { System.out.p ...
- c/c++ 网络编程 文件传输
网络编程 文件传输 1,文件发送端 2,文件接收端 文件发送端: #include <iostream> #include <string.h> #include <sy ...
- 持续代码质量管理-SonarQube-7.3简单使用
安装了SonarQube以及Sonar Scanner之后,就需要那代码检测了.当然为了方便我们使用已有现成的demo,知道到对应的git地址下载即可. 1. sonar-examples下载 htt ...
- LeetCode算法题-Island Perimeter(Java实现)
这是悦乐书的第238次更新,第251篇原创 01 看题和准备 今天介绍的是LeetCode算法题中Easy级别的第105题(顺位题号是463).您将获得一个二维整数网格形式的地图,其中1代表土地,0代 ...