887. 求组合数 III

// 887. 求组合数 III.cpp : 此文件包含 "main" 函数。程序执行将在此处开始并结束。
//

/*

https://www.acwing.com/problem/content/889/
给定 n 组询问,每组询问给定三个整数 a,b,p,其中 p 是质数,请你输出 Cbamodp 的值。

输入格式
第一行包含整数 n。

接下来 n 行,每行包含一组 a,b,p。

输出格式
共 n 行,每行输出一个询问的解。

数据范围
1≤n≤20,
1≤b≤a≤10^18,
1≤p≤105,

输入样例:
3
5 3 7
3 1 5
6 4 13
输出样例:
3
3
2
*/

#include <iostream>

using namespace std;


long long qmi(long long a, long long b, long long M) {
	long long res = 1;
	while (b)
	{
		if (b & 1) res =  res * a % M;
		a = a * a % M;
		b >>= 1;
	}
	return res;
}


long long C(long long a, long long b, long long p) {
	if (b > a) return 0;
	long long ret = 1;
	for (long long i = 1, j = a; i <= b; i++, j--) {
		ret = ret * j % p;
		ret = ret * qmi(i, p - 2, p) % p;
	}

	return ret;
}


long long lucas(long long a, long long b, long long p) {
	if (a < p && b < p) return C(a, b, p);
	return C(a % p, b % p, p) * lucas(a / p, b / p, p)%p;
}


void solve() {
	long long a, b, p; cin >> a >> b >> p;
	 
	long long ret = lucas(a, b, p);

	cout << ret << endl;
}


int main()
{
	int n; cin >> n;
	while (n--) {
		solve();
	}

	return 0;
}

 

posted on 2024-12-09 17:00  itdef  阅读(13)  评论(0)    收藏  举报

导航