LG4720 【模板】扩展卢卡斯定理
扩展卢卡斯定理
求 \(C_n^m \bmod{p}\),其中 \(C\) 为组合数。
\(1≤m≤n≤10^{18},2≤p≤1000000\) ,不保证 \(p\) 是质数。
Fading的题解
设
\]
求出
C_n^m & \mod & {p_1^{\alpha_1}} \\
C_n^m & \mod & {p_2^{\alpha_2}}\\
\vdots\\
C_n^m & \mod & {p_k^{\alpha_k}}
\end{align*} \right.
\]
最后用中国剩余定理合并即可。
假设我们现在要求
=\frac {n!}{m!(n-m)!}\bmod{P^k}
\]
我们无法求得\(m!,(n-m)!\)的逆元,理由很简单,不一定有解。
怎么办呢?发现\(a\)对\(p\)有逆元的充要条件为\(gcd(a,p)=1\)。所以我们换一个式子:
\]
其中\(x\)为\(n!\)中\(P\)因子的个数(包含多少\(P\)因子)。其他都是同理的。
所以如果我们对于一个\(n\),可以求出\(\frac {n!}{P^x}\),那不就完事了吗?(这样我们就可以求逆元了)
问题等价于求\(\frac {n!}{P^x}\bmod{P^k}\)。我们对\(n!\)进行变形:
\]
左边是\(1\sim n\)中所有\(\leq n\)的\(P\)的倍数,右边是其它剩余的。
因为\(1\sim n\)中有\(\lfloor \frac nP\rfloor\)个\(P\)的倍数,
=P^{\lfloor \frac nP\rfloor}(\lfloor \frac nP\rfloor)!\prod_{i=1,i\not\equiv 0\pmod {P}}^ni
\]
显然后面这个\(\bmod\ P\)是有循环节的。
\]
其中\((\prod_{i=1,i\not\equiv 0\pmod {P}}^{P^k}i)^{\lfloor \frac n{P^k}\rfloor}\)是循环节\(1\sim P\)中所有无\(P\)因子数的乘积,\((\prod_{i=P^k\lfloor \frac n{P^k}\rfloor,i\not\equiv 0\pmod {P}}^ni)\)是余项的乘积。
比如
\equiv(1\cdot 2\cdot 4\cdot 5\cdot 7\cdot 8)^2(19\cdot 20\cdot 22)3^7(1\cdot 2\cdot 3\cdot 4\cdot 5\cdot 6\cdot 7)\pmod {3^2}\\
\equiv3^77!(1\cdot 2\cdot 4\cdot 5\cdot 7\cdot 8)^2(19\cdot 20\cdot 22)\pmod {3^2}
\]
正好是一样的,$\lfloor \frac {22}{3^2}\rfloor=2 $。理解了吧...
发现前面这个\(P^{\lfloor \frac nP\rfloor}\)是要除掉的,然而\((\lfloor \frac nP\rfloor)!\)里可能还包含\(P\)。
所以,我们定义函数:
=f(\lfloor \frac nP\rfloor)(\prod_{i=1,i\not\equiv 0\pmod {P}}^{P^k}i)^{\lfloor \frac n{P^k}\rfloor}(\prod_{i=P^k\lfloor \frac n{P^k}\rfloor,i\not\equiv 0\pmod {P}}^ni))
\]
这样就好了。时间复杂度是\(O(\log_{P}n)\)。边界\(f(0)=1\)。
看看原式
=\frac {f(n)}{f(m)f(n-m)}P^{x-y-z}\bmod{P^k}
\]
由于\(f(m)\)一定与\(P^k\)互质,所以我们可以直接用\(exgcd\)求逆元啦。
最后我们还有一个问题,\(P^{x-y-z}\)咋求?其实很好求。
比如说,我们要求\(f(n)=\frac {n!}{P^x}\)中的\(x\)。
设\(g(n)=x\)(上述的\(x\)),再看看阶乘的式子
\]
很显然我们要的就是\(P^{\lfloor \frac nP\rfloor}\),而且\((\lfloor \frac nP\rfloor)!\)可能还有\(P\)因子,所以我们可以得到以下递推式
\]
时间复杂度是\(O(\log_{P}n)\)。边界\(g(n)=0(n<P)\)
所以答案就是
\]
最后用中国剩余定理合并答案即可得到\(C_n^m\)。总时间复杂度\(O(P \log P)\)。
#include<bits/stdc++.h>
#define co const
#define il inline
template<class T> T read(){
T x=0,w=1;char c=getchar();
for(;!isdigit(c);c=getchar())if(c=='-') w=-w;
for(;isdigit(c);c=getchar()) x=x*10+c-'0';
return x*w;
}
template<class T> il T read(T&x){
return x=read<T>();
}
using namespace std;
#define int long long
#define gcd __gcd
int pow(int a,int b,int mod){
a=(a%mod+mod)%mod;
int ans=1;
for(;b;b>>=1,a=a*a%mod)
if(b&1) ans=ans*a%mod;
return ans;
}
void exgcd(int a,int b,int&x,int&y){
if(!b) x=1,y=0;
else exgcd(b,a%b,y,x),y-=a/b*x;
}
int inv(int a,int mod){
if(!a) return 0;
if(gcd(a,mod)!=1) return -1;
int x,y;
exgcd(a,mod,x,y);
return (x%mod+mod)%mod;
}
int fac(int n,int pi,int pk){ //n! % pi^ki pk=pi^ki
if(n<=1) return 1;
int ans=1;
for(int i=2;i<pk;++i)
if(i%pi) ans=ans*i%pk;
ans=pow(ans,n/pk,pk);
for(int i=2;i<=n%pk;++i)
if(i%pi) ans=ans*i%pk;
return ans*fac(n/pi,pi,pk)%pk;
}
int binom(int n,int m,int pi,int pk){
int up=fac(n,pi,pk),dl=fac(m,pi,pk),dr=fac(n-m,pi,pk);
int ex=0;
for(int i=n;i;i/=pi) ex+=i/pi;
for(int i=m;i;i/=pi) ex-=i/pi;
for(int i=n-m;i;i/=pi) ex-=i/pi;
return up*inv(dl,pk)%pk*inv(dr,pk)%pk*pow(pi,ex,pk)%pk;
}
signed main(){
int n=read<int>(),m=read<int>(),P=read<int>();
int ans=0,p=P;
for(int i=2;i<=p;++i)if(p%i==0){
int pi=i,pk=1;
while(p%i==0) p/=i,pk*=i;
ans=(ans+binom(n,m,pi,pk)*(P/pk)%P*inv(P/pk,pk)%P)%P;
}
printf("%lld\n",ans);
return 0;
}
LG4720 【模板】扩展卢卡斯定理的更多相关文章
- bzoj2142 礼物——扩展卢卡斯定理
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=2142 前几天学了扩展卢卡斯定理,今天来磕模板! 这道题式子挺好推的(连我都自己推出来了) , ...
- 【知识总结】扩展卢卡斯定理(exLucas)
扩展卢卡斯定理用于求如下式子(其中\(p\)不一定是质数): \[C_n^m\ mod\ p\] 我们将这个问题由总体到局部地分为三个层次解决. 层次一:原问题 首先对\(p\)进行质因数分解: \[ ...
- 【Luogu3807】【模板】卢卡斯定理(数论)
题目描述 给定\(n,m,p(1≤n,m,p≤10^5)\) 求 \(C_{n+m}^m mod p\) 保证\(P\)为\(prime\) \(C\)表示组合数. 一个测试点内包含多组数据. 输入输 ...
- 卢卡斯定理&扩展卢卡斯定理
卢卡斯定理 求\(C_m^n~mod~p\) 设\(m={a_0}^{p_0}+{a_1}^{p_1}+\cdots+{a_k}^{p_k},n={b_0}^{p_0}+{b_1}^{p_1}+\cd ...
- P3807 【模板】卢卡斯定理
P3807 [模板]卢卡斯定理 求 \(C_{m + n}^{m} \% p\) ( \(1\le n,m,p\le 10^5\) ) 错误日志: 数组开小(哇啊啊啊洼地hi阿偶我姑父阿贺佛奥UFO爱 ...
- 洛谷 P3807 【模板】卢卡斯定理
P3807 [模板]卢卡斯定理 题目背景 这是一道模板题. 题目描述 给定n,m,p(1\le n,m,p\le 10^51≤n,m,p≤105) 求 C_{n+m}^{m}\ mod\ pCn+mm ...
- 洛谷——P3807 【模板】卢卡斯定理
P3807 [模板]卢卡斯定理 洛谷智推模板题,qwq,还是太弱啦,组合数基础模板题还没做过... 给定n,m,p($1\le n,m,p\le 10^5$) 求 $C_{n+m}^{m}\ mod\ ...
- 【luogu P3807】【模板】卢卡斯定理/Lucas 定理(含 Lucas 定理证明)
[模板]卢卡斯定理/Lucas 定理 题目链接:luogu P3807 题目大意 求 C(n,n+m)%p 的值. p 保证是质数. 思路 Lucas 定理内容 对于非负整数 \(n\),\(m\), ...
- 【刷题】洛谷 P3807 【模板】卢卡斯定理
题目背景 这是一道模板题. 题目描述 给定\(n,m,p( 1\le n,m,p\le 10^5)\) 求 \(C_{n+m}^{m}\ mod\ p\) 保证 \(p\) 为prime \(C\) ...
随机推荐
- 第2/7Beta冲刺
1.团队成员 成员姓名 成员学号 秦裕航 201731062432(组长) 刘东 201731062227 张旭 201731062129 王伟 201731062214 2.SCRU部分 2.1各成 ...
- GridView树状结构显示
下面的树形结构代码需要GridVIew中的数据要求是按照上下级关系已经排列好的顺序,比如: GridView ID ParentID Name 1 0 父1 2 1 父1子1 3 1 父1子2 4 3 ...
- java8新特性1--Lambda表达式
一.Lambda表达式是什么? Lambda表达式有两个特点 一是匿名函数,二是可传递. 匿名函数的应用场景是 通常是在需要一个函数,但是又不想费神去命名一个函数的场合下使用Lambda表达式.lam ...
- LeetCode 653. 两数之和 IV - 输入 BST(Two Sum IV - Input is a BST)
653. 两数之和 IV - 输入 BST 653. Two Sum IV - Input is a BST 题目描述 给定一个二叉搜索树和一个目标结果,如果 BST 中存在两个元素且它们的和等于给定 ...
- 09 Servlet相关知识点---学习笔记
1.概念:运行在服务器端的小程序 Servlet就是一个接口,定义了Java类被浏览器访问到(tomcat识别)的规则.将来我们自定义一个类,实现Servlet接口,复写方法. 2.快速入门:(1)创 ...
- python中将已有链接的视频进行下载
使用python爬取视频网站时,会得到一系列的视频链接,比如MP4文件.得到视频文件之后需要对视频进行下载,本文写出下载视频文件的函数. 首先导入requests库,安装库使用pip install ...
- quartz2.3.0(十)xml配置方式定义quartz定时任务
1.新增pom依赖 除了按照<quartz2.3.0系列目录——带您由浅入深全面掌握quartz2.3.0>添加依赖之外,pom.xml里新增加依赖: <dependency> ...
- CF627E Orchestra [矩阵计数]
也许更好的阅读体验 \(\mathcal{Description}\) 题目大意 有一个\(r * c\)的矩阵上有\(n\)个点,问有多少个子矩阵里包含至少\(k\)个点 输入格式 第一行四个数\( ...
- [LOJ#3119][Luogu5400][CTS2019]随机立方体(容斥+DP)
https://www.cnblogs.com/cjyyb/p/10900993.html #include<cstdio> #include<algorithm> #defi ...
- python多线程一些知识点梳理
学习python的进程和线程以来,对这两个概念一直都处于模糊状态,所以决定花点时间好好学习一下这块知识.以下是我自己在学习过程中形成的一些疑问以及搜集的一些相应的比较好的答案,整理如下,方便复习自查. ...