https://www.luogu.org/problemnew/show/3807

題目背景

這是一道模板題。

題目描述

給定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表示組合數(shù)。

一個測試點內(nèi)包含多組數(shù)據(jù)。

輸入輸出格式

輸入格式：

?

第一行一個整數(shù)T(T\le 10T≤10)，表示數(shù)據(jù)組數(shù)

第二行開始共T行，每行三個數(shù)n m p，意義如上

?

輸出格式：

?

共T行，每行一個整數(shù)表示答案。

?

輸入輸出樣例

輸入樣例#1：?復制 2 1 2 5 2 1 5 輸出樣例#1：?復制 3 3


1 #include <cstdio> 2 3 #define LL long long 4 inline void read(int &x) 5 { 6 x=0; register char ch=getchar(); 7 for(; ch>'9'||ch<'0'; ) ch=getchar(); 8 for(; ch>='0'&&ch<='9'; ch=getchar()) x=x*10+ch-'0'; 9 } 10 const int N(1e5+5); 11 LL fac[N]; 12 13 inline LL Pow(LL a,int b,int p) 14 { 15 LL ret=1; 16 for(; b; b>>=1,a*=1ll*a,a%=p) 17 if(b&1) ret*=1ll*a,ret%=p; 18 return ret; 19 } 20 21 inline LL C(LL n,LL m,LL p) 22 { 23 if(n<m) return 0; 24 return fac[n]%p*Pow(fac[m],p-2,p)%p*Pow(fac[n-m],p-2,p)%p; 25 } 26 27 inline LL lus(LL n,LL m,LL p) 28 { 29 if(m==0) return 1; 30 return C(n%p,m%p,p)*lus(n/p,m/p,p)%p; 31 } 32 33 int Presist() 34 { 35 int t; read(t); fac[0]=1; 36 for(int n,m,p; t--; ) 37 { 38 read(n),read(m),read(p); 39 for(int i=1; i<=n+m; ++i) 40 fac[i]=1ll*fac[i-1]%p*i%p; 41 printf("%lld\n",lus(n+m,m,p)); 42 } 43 return 0; 44 } 45 46 int Aptal=Presist(); 47 int main(int argc,char**argv){;}

?

轉(zhuǎn)載于:https://www.cnblogs.com/Shy-key/p/7793228.html

總結(jié)

以上是生活随笔為你收集整理的洛谷—— P3807 【模板】卢卡斯定理的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。