洛谷 P3807 【模板】卢卡斯定理

题目背景

这是一道模板题。

题目描述

给定n,m,p(1\le n,m,p\le 10^51≤n,m,p≤105)

求 C_{n+m}^{m}\ mod\ pCn+mm mod p

保证P为prime

C表示组合数。

一个测试点内包含多组数据。

输入输出格式

输入格式：

第一行一个整数T(T\le 10T≤10)，表示数据组数

第二行开始共T行，每行三个数n m p，意义如上

输出格式：

共T行，每行一个整数表示答案。

输入输出样例

输入样例#1：复制

2
1 2 5
2 1 5

输出样例#1：复制

3
3

卢卡斯定理

$C(n,m) \% p = C(n \% p, m \% p) * C(n / p, m / p)$

对于这道题来说，p是素数，解逆元的时候用快速幂

 1 #include<cstdio>
 2 #include<cstring>
 3 #include<cmath>
 4 #include<algorithm>
 5 #define LL long long
 6 using namespace std;
 7 const LL MAXN=1e6+10;
 8 const LL INF=0x7fffff;
 9 inline LL read()
10 {
11     char c=getchar();LL flag=1,x=0;
12     while(c<'0'||c>'9')    {if(c=='-')    flag=-1;c=getchar();}
13     while(c>='0'&&c<='9')    x=x*10+c-48,c=getchar();return x*flag;
14 }
15 LL js[MAXN];
16 LL fastpow(LL a,LL p,LL mod)
17 {
18     LL base=1;
19     while(p)
20     {
21         if(p&1)    base=(base*a)%mod;
22         a=(a*a)%mod;
23         p>>=1;
24     }
25     return base;
26 }
27 LL C(LL n,LL m,LL mod)
28 {
29     if(m>n)    return 0;
30     return js[n]*fastpow(js[m],mod-2,mod)*fastpow(js[n-m],mod-2,mod)%mod;
31 }
32 LL Lucas(LL n,LL m,LL mod)
33 {
34     if(m==0)    return 1;
35     else return C(n%mod,m%mod,mod)*(Lucas(n/mod,m/mod,mod))%mod;
36 }
37 int main()
38 {
39     LL T=read();
40     js[0]=1;
41     while(T--)
42     {
43         LL n=read(),m=read(),mod=read();
44         for(LL i=1;i<=mod;i++)    js[i]=(js[i-1]*i)%mod;
45         printf("%lld\n",Lucas(n+m,m,mod)%mod);
46     }
47     return 0;
48 }

转载于:https://www.cnblogs.com/zwfymqz/p/7788177.html

洛谷 P3807 【模板】卢卡斯定理相关推荐

洛谷——P3807 【模板】卢卡斯定理
P3807 [模板]卢卡斯定理洛谷智推模板题,qwq,还是太弱啦,组合数基础模板题还没做过... 给定n,m,p($1\le n,m,p\le 10^5$) 求 $C_{n+m}^{m}\ mod\ ...
专题·树链剖分【including 洛谷·【模板】树链剖分
初见安~~~终于学会了树剖~~~ [兴奋]当初机房的大佬在学树剖的时候我反复强调过:"学树剖没有前途的!!!" 恩.真香. 一.重链与重儿子所谓树剖--树链剖分,就是赋予一个链的 ...
洛谷·【模板】点分树 | 震波【including 点分树
初见安-这里是传送门:洛谷P6329 [模板]点分树 | 震波一.点分树其实你会点分治的话,点分树就是把点分治时的重心提出来重新连城一棵树. 比如当前点是u,求出子树v的重心root后将root与 ...
OI模板卢卡斯定理
卢卡斯定理 / Lucas\text{Lucas}Lucas定理现有问题: 给定整数 n,m,pn,m,pn,m,p,求 (nm)modp\dbinom{n}{m}\bmod p(mn)modp, ...
洛谷P3807 【模板】卢卡斯定理exgcd
题目背景这是一道模板题. 题目描述给定n,m,p(1\le n,m,p\le 10^51≤n,m,p≤105 ) 求 C_{n+m}^{m}\ mod\ pCn+mm mod p 保证P为pri ...
[洛谷P3807]【模板】卢卡斯定理
题目大意:给你$n,m,p(p \in \rm prime)$,求出$C_{n + m}^m\bmod p(可能p\leqslant n,m)$ 题解:卢卡斯$Lucas$定理,$C_B^A\bmod ...
洛谷—— P3807 【模板】卢卡斯定理
https://www.luogu.org/problemnew/show/3807 题目背景这是一道模板题. 题目描述给定n,m,p(1\le n,m,p\le 10^51≤n,m,p≤105) ...
洛谷.4245.[模板]任意模数NTT(MTT/三模数NTT)
题目链接三模数$NTT$: 就是多模数$NTT$最后$CRT$一下...下面两篇讲的都挺明白的. https://blog.csdn.net/kscla/article/details/ ...
洛谷.4897.[模板]最小割树(Dinic)
题目链接最小割树模板.具体见:https://www.cnblogs.com/SovietPower/p/9734013.html. ISAP不知为啥T成0分了.. Dinic: //1566ms ...