AcWing算法基础课笔记——求组合数3

求组合数Ⅲ

20万组数据, 1 ≤ b ≤ a ≤ 1 0 18 , 1 ≤ p ≤ 1 0 5 1 \le b \le a \le 10^{18}, 1\le p \le 10 ^5 1≤b≤a≤1018,1≤p≤105,使用卢卡斯定理。

卢卡斯定理:
C a b ≡ C a m o d p b m o d p C a / p b / p ( m o d p ) C_a^b \equiv C_{a mod p}^{b mod p}C_{a /p}^{b / p} (mod p) Cab≡CamodpbmodpCa/pb/p(modp)

时间复杂度为 O ( l o g p N ⋅ p l o g p ) O(log _p N \cdot plogp) O(logpN⋅plogp)

题目

代码

cpp 复制代码
#include<iostream>
#include<algorithm>

using namespace std;

typedef long long LL;

int p;

int qmi(int a, int k){
	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) {
	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;
	}
	return 0;
}
相关推荐
生锈的键盘5 分钟前
推荐算法实践:movielens数据集
算法
董董灿是个攻城狮5 分钟前
Transformer 通关秘籍9:词向量的数值实际上是特征
算法
nuise_13 分钟前
李宏毅机器学习笔记06 | 鱼和熊掌可以兼得的机器学习 - 内容接宝可梦
人工智能·笔记·机器学习
林泽毅14 分钟前
SwanLab x EasyR1:多模态LLM强化学习后训练组合拳,让模型进化更高效
算法·llm·强化学习
小林熬夜学编程16 分钟前
【高并发内存池】第八弹---脱离new的定长内存池与多线程malloc测试
c语言·开发语言·数据结构·c++·算法·哈希算法
刚入门的大一新生22 分钟前
归并排序延伸-非递归版本
算法·排序算法
独好紫罗兰27 分钟前
洛谷题单3-P1980 [NOIP 2013 普及组] 计数问题-python-流程图重构
开发语言·python·算法
独好紫罗兰32 分钟前
洛谷题单3-P1009 [NOIP 1998 普及组] 阶乘之和-python-流程图重构
开发语言·python·算法
曦月逸霜43 分钟前
蓝桥杯高频考点——高精度(含C++源码)
c++·算法·蓝桥杯
ゞ 正在缓冲99%…1 小时前
leetcode152.乘积最大子数组
数据结构·算法·leetcode