卢卡斯定理求组合数 C ( N, M ) % P,当 P 是素数时适用,这里可以当一个模板用 HUD 3037
#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int SIZE = 100231;
LL F[SIZE];
void init(LL Mod) {
F[0] = 1;
for (int pos = 1; pos <= Mod; pos++)
F[pos] = F[pos - 1] * pos % Mod;
}
LL inv(LL A, LL M) {
return A == 1 ? 1 : inv(M % A, M) * (M - M / A) % M;
}
LL Lucas(LL N, LL M, LL Mod) {
init(Mod);
LL ans = 1;
while (N && M) {
LL A = N % Mod, B = M % Mod;
if (A < B)
return 0;
ans = ans * F[A] % Mod * inv(F[B] * F[A - B] % Mod, Mod) % Mod;
N /= Mod; M /= Mod;
}
return ans;
}
int main() {
ios::sync_with_stdio(false);
int quan; cin >> quan;
LL N, M, Mod;
while (quan--) {
cin >> N >> M >> Mod;
cout << Lucas(N + M, N, Mod) << endl;
}
}0.0分
0 人评分
C语言网提供由在职研发工程师或ACM蓝桥杯竞赛优秀选手录制的视频教程,并配有习题和答疑,点击了解:
一点编程也不会写的:零基础C语言学练课程
解决困扰你多年的C语言疑难杂症特性的C语言进阶课程
从零到写出一个爬虫的Python编程课程
只会语法写不出代码?手把手带你写100个编程真题的编程百练课程
信息学奥赛或C++选手的 必学C++课程
蓝桥杯ACM、信息学奥赛的必学课程:算法竞赛课入门课程
手把手讲解近五年真题的蓝桥杯辅导课程
发表评论 取消回复