洛谷——P3807 【模板】盧卡斯定理
阿新 • • 發佈:2017-08-14
|| turn thml 數據 mod text clu -h eset
P3807 【模板】盧卡斯定理
題目背景
這是一道模板題。
題目描述
給定n,m,p(1\le n,m,p\le 10^51≤n,m,p≤10?5??)
求 C_{n+m}^{m}\ mod\ pC?n+m?m?? mod p
C表示組合數。
一個測試點內包含多組數據。
輸入輸出格式
輸入格式:
第一行一個整數T(T\le 10T≤10),表示數據組數
第二行開始共T行,每行三個數n m p,意義如上
輸出格式:
共T行,每行一個整數表示答案。
輸入輸出樣例
輸入樣例#1:2 1 2 5 2 1 5輸出樣例#1:
3 3
模板:
#include<cstdio> #include<cstdlib> #include<cstring> #include<iostream> #include<algorithm> #define N 100001 #define ll long long using namespace std; ll t,n,m,p,ans,f[N]; ll read() { ll x=0,f=1; char ch=getchar(); while(ch<‘0‘||ch>‘9‘){if(ch==‘-‘)f=-1; ch=getchar();} while(ch>=‘0‘&&ch<=‘9‘){x=x*10+ch-‘0‘; ch=getchar();} return x*f; } ll mi(ll n,ll k,ll p) { ll res=1; while(k) { if(k&1) res=n*res%p; n=n*n%p; k>>=1; } return res; } ll c(ll n,ll m,ll p) {if(m>n) return 0; return f[n]*mi(f[n-m]*f[m],p-2,p)%p; } ll lus(ll n,ll m,ll p) { if(m==0) return 1; return c(n%p,m%p,p)*lus(n/p,m/p,p)%p; } int main() { t=read(); while(t--) { f[0]=1; n=read(),m=read(),p=read(); for(ll i=1;i<=p;i++) f[i]=f[i-1]*i%p; ans=lus(m+n,m,p); printf("%d\n",ans); } return 0; }
洛谷——P3807 【模板】盧卡斯定理