AcWing 887. 求组合数 III
原题链接
简单
作者:
我要出去乱说
,
2021-01-21 22:25:40
,
所有人可见
,
阅读 292
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long LL;
int p;
int qmi(int a, int k) //只需要两个参数是因为p已经定义为全局变量了
{
int res = 1;
while (k)
{
if (k & 1) res = (LL)res * a % p;
a = (LL)a * a % p;
k >>= 1;
}
return res;
}
int C(int a, int b)
{
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;
}
return res;
}
int lucas(LL a, LL b) //a % p与 b % p结果必定是小于p的,而a / p 和 b / p则不一定,需要继续递归
{
if (a < p && b < p) return C(a, b);
return (LL)C(a % p, b % p) * lucas(a / p, b / p) % p;
}
int main()
{
int n;
cin >> n;
while (n -- )
{
LL a, b;
cin >> a >> b >> p;
cout << lucas(a, b) << endl;
}
}