2023.2.26【模板】扩展Lucas定理
2023.2.26【模板】扩展Lucas定理
题目概述
求\(\binom {n}{m} mod\) \(p\) 的值,不保证\(p\)为质数
算法流程
(扩展和普通算法毫无关系)
由于\(p\)不是质数,我们考虑[SDOI2010]古代猪文 - 洛谷中的处理方法:将\(p\)质因数分解得:
\]
所以我们考虑计算$\binom nm mod $ \({p_i}^{c_i}\)的值,再用CRT合并即可
展开上式:
\]
我们发现由于\(m!(n - m)!\)中可能含有因数p,我们无法求出\(m!(n - m)!\)模\({p_i}^{c_i}\)意义下的逆元,所以我们考虑除去三个数中所有的p因子,假设\(p^x | n!\)且\(p^{x+1} \nmid n!\),即x是\(n!\)中p因子的个数(对于\(m!\)和\((n - m)!\)同理)
\]
由于\(\frac{n!}{p^x}、\frac{m!}{p^y}、\frac{(n - m)!}{p^z}\)三式同构,我们考虑计算其中一个式子(以下用\(p\)替换\(p_i\))
\]
展开为
\]
提出p的倍数
\]
即
\]
如果暴力计算\(\Pi_{i = 1;i \not\equiv 0}^{n}\)复杂度过高,不难发现其有一个循环节,即每过p个数就会少乘上第p个数,又因为\({p_i}^{c_i}+ r \equiv r\ mod\ {p_i}^{c_i}\),所以我们以\({p_i}^{c_i}\)作为这个循环节
\]
对于\(\Pi_{i = 1;i \not\equiv 0}^{p^{c_i}}\)和\(\Pi_{i = {p^{c_i}}\lfloor\frac{n}{p^{c_i}}\rfloor;i \not\equiv 0}^{n}\),暴力计算即可
不管\(x\)取何值,最终p因子都会消除,所以计算时去掉\(p^{\lfloor \frac np \rfloor}\)
因为\(\lfloor \frac np \rfloor!\)中可能含有p因子,所以我们将其进行递归:
设\(f(n) = \frac {n!}{p^x}\ mod \ {p}^{c_i}\),则:
\]
根据此式递推即可,时间复杂度为\(O(log_pn)\),不会证明qwq
对于外面的\(p^{x - y - z}\),只要求出\(x、y、z\)的值就可以计算了
观察以上函数可知,每次在\(f(n)\)这一层就会去掉\(\lfloor \frac np \rfloor\)个p因子
定义\(g(n)\)为\(n!\)中p因子的个数,则:
\]
此结论对于其他题目也同样有效
所以原始式子就转化成了
\]
因为去掉了p因子,所以\(f(m)\)和\(f(n - m)\)与\(p^{c_i}\)互质,可以求逆元
因为\(p^{c_i}\)不是质数,不能直接用费马小定理计算,所以我们采用\(exgcd\)求逆元
最后进行CRT合并答案
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll res[101],d[101],zs[101],tot = 0,M[101];
inline ll g(ll n,ll p)
{
if(n == 0) return 0;
return g(n / p,p) + n / p;
}
inline ll ksm(ll base,ll pts,ll mod)
{
ll ret = 1;
for(;pts > 0;pts >>= 1,base = base * base % mod)
if(pts & 1)
ret = ret * base % mod;
return ret;
}
inline ll F(ll n,ll p,ll k)
{
if(n == 0) return 1;
ll P = ksm(p,k,1e18 + 1);
ll mul = 1;
for(ll i = 1;i <= P;i++)
if(i % p)
mul = mul * i % P;
mul = ksm(mul,n / P,P);
for(ll i = P * (n / P);i <= n;i++)
if(i % p)
mul = mul * (i % P) % P;
return F(n / p,p,k) * mul % P;
}
inline void exgcd(ll a,ll b,ll &x,ll &y)
{
if(b == 0)
{
x = 1;
y = 0;
return;
}
ll tmp;
exgcd(b,a % b,x,y);
tmp = y;
y = x - (a / b) * y;
x = tmp;
}
inline ll exlucas(ll n,ll m,ll p)
{
ll tmp = p;
for(ll i = 2;i <= sqrt(p);i++)
{
if(tmp % i == 0)
{
++tot;
d[tot] = i;
while(tmp % i == 0)
{
tmp /= i;
++zs[tot];
}
}
}
if(tmp != 1)
{
++tot;
d[tot] = tmp;
zs[tot] = 1;
}
for(int i = 1;i <= tot;i++)
{
ll P = ksm(d[i],zs[i],1e18 + 1);
ll inv1,inv2,yy;
exgcd(F(m,d[i],zs[i]),P,inv1,yy);
exgcd(F(n - m,d[i],zs[i]),P,inv2,yy);
inv1 = (inv1 % P + P) % P;
inv2 = (inv2 % P + P) % P;
res[i] = F(n,d[i],zs[i]) * inv1 % P * inv2 % P * ksm(d[i],g(n,d[i]) - g(m,d[i]) - g(n - m,d[i]),P) % P;
M[i] = P;
}
ll ans = 0;
for(int i = 1;i <= tot;i++)
{
ll inv,yy;
exgcd(p / M[i],M[i],inv,yy);
inv = (inv % M[i] + M[i]) % M[i];
ans = (ans + res[i] * (p / M[i]) % p * inv % p) % p;
}
return ans;
}
int main()
{
ll n,m,p;
cin>>n>>m>>p;
cout<<exlucas(n,m,p);
return 0;
}
2023.2.26【模板】扩展Lucas定理的更多相关文章
- 2015 ICL, Finals, Div. 1 Ceizenpok’s formula(组合数取模,扩展lucas定理)
J. Ceizenpok’s formula time limit per test 2 seconds memory limit per test 256 megabytes input stand ...
- BZOJ - 2142 礼物 (扩展Lucas定理)
扩展Lucas定理模板题(貌似这玩意也只能出模板题了吧~~本菜鸡见识鄙薄,有待指正) 原理: https://blog.csdn.net/hqddm1253679098/article/details ...
- 【learning】 扩展lucas定理
首先说下啥是lucas定理: $\binom n m \equiv \binom {n\%P} {m\%P} \times \binom{n/P}{m/P} \pmod P$ 借助这个定理,求$\bi ...
- [bzoj2142]礼物(扩展lucas定理+中国剩余定理)
题意:n件礼物,送给m个人,每人的礼物数确定,求方案数. 解题关键:由于模数不是质数,所以由唯一分解定理, $\bmod = p_1^{{k_1}}p_2^{{k_2}}......p_s^{{k_ ...
- Lucas定理和扩展Lucas定理
1.Lucas定理 首先给出式子:\(C_n^m\%p = C_{\lfloor\frac{n}{p}\rfloor}^{\lfloor\frac{m}{p}\rfloor} * C_{n\%p}^{ ...
- Ceizenpok’s formula Gym - 100633J 扩展Lucas定理 + 中国剩余定理
http://codeforces.com/gym/100633/problem/J 其实这个解法不难学的,不需要太多的数学.但是证明的话,我可能给不了严格的证明.可以看看这篇文章 http://ww ...
- [笔记] 扩展Lucas定理
[笔记] 扩展\(Lucas\)定理 \(Lucas\)定理:\(\binom{n}{m} \equiv \binom{n/P}{m/P} \binom{n \% P}{m \% P}\pmod{P} ...
- [学习笔记]扩展LUCAS定理
可以先做这个题[SDOI2010]古代猪文 此算法和LUCAS定理没有半毛钱关系. [模板]扩展卢卡斯 不保证P是质数. $C_n^m=\frac{n!}{m!(n-m)!}$ 麻烦的是分母. 如果互 ...
- BZOJ3129/洛谷P3301方程(SDOI2013)容斥原理+扩展Lucas定理
题意:给定方程x1+x2+....xn=m,每个x是正整数.但是对前n1个数做了限制x1<=a1,x2<=a2...xn1<=an1,同时对第n1+1到n1+n2个数也做了限制xn1 ...
- 扩展Lucas定理
(1)Lucas定理:p为素数,则有: (2)证明: n=(ak...a2,a1,a0)p = (ak...a2,a1)p*p + a0 = [n/p]*p+a0,m=[m/p]*p+b0其次,我们 ...
随机推荐
- Nmap安装
Nmap(Network Mapper,网络映射器)是一款开放源代码的网络探测和安全审核工具.它被设计用来快速扫描大型网络,包括主机探测与发现.开放的端口情况.操作系统与应用服务指纹识别.WAF识别及 ...
- Vue element 自定义表单验证(验证手机号)
<el-form :model="ruleForm" status-icon :rules="rules" ref="ruleForm" ...
- CH432,CH438,CH9434串口扩展芯片常见问题
目前WCH有三款串口扩展芯片CH432,CH438以及CH9434. 型号 CH432 CH438 CH9434 扩展串口数量 2 8 4 通讯接口 并口/SPI(具体需要看芯片封装) 并口 SPI ...
- [Computer Networks]一个http请求的完成的全过程
摘要 本文主要讲述了一个 http request 请求从发出到收到 response 的整个生命周期,希望可以通过对整个流程的一个描述来梳理清楚五层网络协议的定义以及各层之间是如何协作的. 对于后端 ...
- Python自动化操作sqlite数据库
你好,我是悦创. 原文首发:https://bornforthis.cn/column/pyauto/ 1. 什么是数据库 数据库是"按照数据结构来组织.存储和管理数据的仓库",是 ...
- idea 函数名灰色
idea被引用的方法名突然全部灰掉了 idea被引用的方法名突然全部灰掉了[已解决]_weixin_42554373的博客-CSDN博客_idea方法名灰色
- BBS项目 未完待续
项目开发基本流程 1.需求分析 2.架构设计 3.分组开发 4.提交测试 5.交付上线 创建项目配置 环境配置 TEMPLATES = [ { 'BACKEND': 'django.template. ...
- 简述HashSet的扩容机制以及我们在重写equals()的时候为何会重写hashcode()
简述HashSet的扩容机制以及我们在重写equals()的时候为何会重写hashcode() 摘要:在背面试知识点的时候存在这样一条著名的面试题:我们重写equals()的时候为什么要重写has ...
- SPOJ GCDMAT - GCD OF MATRIX
简要题意 给出三个整数 \(T,n,m\),\(T\) 组询问,每组询问给出四个整数 \(i_1,j_1,i_2,j_2\)(数据保证 \(i_1,j_1\leq n\ \ i_2,j_2\leq m ...
- do while 出口條件循環