AcWing 887. 求组合数 III
原题链接
简单
作者:
术
,
2021-03-15 12:56:17
,
所有人可见
,
阅读 268
#include <iostream>
using namespace std;
const int mod=1e9+7;
typedef long long LL;
int qmi(int a,int b,int p){
int res=1%p;
while(b){
if(b&1) res=(LL)res*a%p;
a=(LL)a*a%p;
b=b>>1;
}
return res;
}
//用定义求C
int C(int a,int b,int p){
int res=1;
for(int i=1,j=a;i<=b;i++,j--){
res=(LL)res*j%p;
res=(LL)res*qmi(i,p-2,p)%p;
}
return res;
}
int lucas(LL a,LL b,int p){//注意ab的类型
if(a<p&&b<p) return C(a,b,p);
return (LL)C(a%p,b%p,p)*lucas(a/p,b/p,p)%p;
}
int main()
{
int n;
cin>>n;
while(n--){
LL a,b,p;
cin>>a>>b>>p;
cout<<lucas(a,b,p)<<endl;
}
//cout << "Hello world!" << endl;
return 0;
}