卢卡斯定理求组合数  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分

0 人评分

C语言网提供由在职研发工程师或ACM蓝桥杯竞赛优秀选手录制的视频教程,并配有习题和答疑,点击了解:

一点编程也不会写的:零基础C语言学练课程

解决困扰你多年的C语言疑难杂症特性的C语言进阶课程

从零到写出一个爬虫的Python编程课程

只会语法写不出代码?手把手带你写100个编程真题的编程百练课程

信息学奥赛或C++选手的 必学C++课程

蓝桥杯ACM、信息学奥赛的必学课程:算法竞赛课入门课程

手把手讲解近五年真题的蓝桥杯辅导课程

评论列表 共有 0 条评论

暂无评论