扩展卢卡斯 (板子)

求c(m,n)%a

不知道原理，不过套板子就对了：

#include<bits/stdc++.h>using namespace std;typedef long long ll;ll exgcd(ll a,ll b,ll &x,ll &y)
{if(!b){x=1;y=0;return a;}ll res=exgcd(b,a%b,x,y),t;t=x;x=y;y=t-a/b*y;return res;
}ll p;inline ll power(ll a,ll b,ll mod)
{ll sm;for(sm=1;b;b>>=1,a=a*a%mod)if(b&1)sm=sm*a%mod;return sm;
}ll fac(ll n,ll pi,ll pk)
{if(!n)return 1;ll res=1;for(register ll i=2;i<=pk;++i)if(i%pi)(res*=i)%=pk;res=power(res,n/pk,pk);for(register ll i=2;i<=n%pk;++i)if(i%pi)(res*=i)%=pk;return res*fac(n/pi,pi,pk)%pk;
}inline ll inv(ll n,ll mod)
{ll x,y;exgcd(n,mod,x,y);return (x+=mod)>mod?x-mod:x;
}inline ll CRT(ll b,ll mod){return b*inv(p/mod,mod)%p*(p/mod)%p;}const int MAXN=11;static ll n,m;static ll w[MAXN];inline ll C(ll n,ll m,ll pi,ll pk)
{ll up=fac(n,pi,pk),d1=fac(m,pi,pk),d2=fac(n-m,pi,pk);ll k=0;for(register ll i=n;i;i/=pi)k+=i/pi;for(register ll i=m;i;i/=pi)k-=i/pi;for(register ll i=n-m;i;i/=pi)k-=i/pi;return up*inv(d1,pk)%pk*inv(d2,pk)%pk*power(pi,k,pk)%pk;
}inline ll exlucus(ll n,ll m)
{ll res=0,tmp=p,pk;static int lim=sqrt(p)+5;for(register int i=2;i<=lim;++i)if(tmp%i==0){pk=1;while(tmp%i==0)pk*=i,tmp/=i;(res+=CRT(C(n,m,i,pk),pk))%=p;}if(tmp>1)(res+=CRT(C(n,m,tmp,tmp),tmp))%=p;return res;
}int main()
{scanf("%lld%lld%d",&n,&m,&p);printf("%d\n",exlucus(n,m));return 0;
}

扩展卢卡斯 (板子)相关推荐

P4720 【模板】扩展卢卡斯定理/exLucas(无讲解，纯记录模板)
P4720 [模板]扩展卢卡斯定理/exLucas 题意: CnmmodpC_{n}^{m}\bmod pCnmmodp 对于 100% 的数据,1≤m≤n≤1018,2≤p≤106,不保证 p 是 ...
『数学』--数论--组合数+卢卡斯定理+扩展卢卡斯定理
组合数: 在N个数中选取M个数,问选的方式有几种? 直接递归暴力简单 #include<cstdio> const int N = 2000 + 5; const int MOD = (i ...
P2183 [国家集训队]礼物（扩展卢卡斯）
P2183 [国家集训队]礼物题意: 有n个礼物,分给m个人,分给第i个人的礼物数量是wi,问送礼物的方案数. 题解: 扩展卢卡斯模板题很容易看出和组合数有关的题目,对于总方案,完美可以将其分解为 ...
组合数学 —— 组合数取模 —— 卢卡斯定理与扩展卢卡斯定理
[卢卡斯定理] 1.要求:p 是质数,m.n 很大但 p 很小或者 n.m 不大但大于 p 2.定理内容其中, 3.推论当将 n 写成 p 进制:,将 m 写成 p 进制: 时,有: 4.实现 ...
洛谷P4720 【模板】扩展卢卡斯
P4720 [模板]扩展卢卡斯题目背景这是一道模板题. 题目描述求 C(n,m)%P 其中 C 为组合数. 输入输出格式输入格式: 一行三个整数 n,m,p ,含义由题所述. 输出格式: 一行 ...
洛谷 P4720 【模板】扩展卢卡斯定理/exLucas
[模板]扩展卢卡斯定理/exLucas 题目背景这是一道模板题. 题目描述求 Cnmmodp{\mathrm{C}}_n^m \bmod{p}Cnmmodp 其中 C\mathrm{C}C 为组 ...
卢卡斯定理扩展卢卡斯
卢卡斯定理&扩展卢卡斯 Lucas EXlucas 例题 Lucas 卢卡斯定理: ( m n ) = ( m p n p ) ∗ ( m ( m o d p ) n ( m o d p ) ...
【知识总结】扩展卢卡斯定理（exLucas）
扩展卢卡斯定理用于求如下式子(其中 p p p不一定是质数): C n m m o d p C_n^m\ mod\ p Cnm mod p 我们将这个问题由总体到局部地分为三个层次解决. 层次一:原 ...
【数学】扩展卢卡斯定理
Description 求 ( n m ) m o d p \dbinom{n}{m}\bmod p (mn)modp 其中 p p p 较小且不保证 p p p 是质数. Method 前置芝士 ...
扩展卢卡斯定理(Exlucas)
题目链接戳我前置知识中国剩余定理(crt)或扩展中国剩余定理(excrt) 乘法逆元组合数的基本运用扩展欧几里得(exgcd) 说实话Lucas真的和这个没有什么太大的关系,但是Lucas还 ...

扩展卢卡斯 (板子)

扩展卢卡斯 (板子)相关推荐

最新文章

热门文章