【C++】 排列与组合算法详解(进阶篇)

这篇具有很好参考价值的文章主要介绍了【C++】 排列与组合算法详解(进阶篇)。希望对大家有所帮助。如果存在错误或未考虑完全的地方,请大家不吝赐教,您也可以点击"举报违法"按钮提交疑问。

c++排列组合,笔记,c++,lucas定理,费马小定理,逆元,递推

写在前面

我上次发了一篇题解:C++排列与组合算法详解
最开始,我是抱着水题解的想法写的,但却成为了阅读量 最高 的文章,没有之一。

所以今天咱们来重制一篇文章,介绍几个进阶优化版的算法。


算法1:朴素算法

思路

具体见 C++排列与组合算法详解

缺点

不能将结果取模,因为朴素的组合公式在取模意义下没用。


算法2:递推预处理

思路

我们发现:
C a 0 = 1 C a b = C a − 1 b + C a − 1 b − 1 ( a , b > 0 ) C_a^0 = 1\\ C_a^b=C_{a-1}^b+C_{a-1}^{b-1}(a,b>0) Ca0=1Cab=Ca1b+Ca1b1(a,b>0)

所以我们可以写一个递推函数(部分非主要内容已省略):

void init_C()
{
    for (int i = 0; i < N; i ++ ) // N 表示预处理最大的下标
        for (int j = 0; j <= i; j ++ )
            if (!j) c[i][j] = 1;
            else c[i][j] = (c[i - 1][j] + c[i - 1][j - 1]) % P;
}

再预处理阶乘:

void f(int n)
{
    f[0] = 1;
    for (int i = 1; i <= n; i ++ )
        f[i] = (LL)f[i - 1] * i % P;
}

需要排列的话还可以预处理排列:

void init_A(int n)
{
    for (int i = 0; i <= n; i ++ )
        for (int j = 0; j <= i; j ++ )
            a[i][j] = (LL)f[i - j] * c[i][j] % P;
}
时间复杂度: O ( n 2 ) O(n^2) O(n2)

可以处理 5000 5000 5000 以内规模的数据


算法3:阶乘逆元

思路

根据费马小定理可得,当 p p p 为质数时, a p − 1 ≡ 1 ( m o d p ) a^{p-1} \equiv 1\pmod p ap11(modp)
∴ a p − 2 ≡ 1 a ( m o d p ) \therefore a^{p-2} \equiv \frac{1}{a}\pmod p ap2a1(modp)
这就是乘法逆元,通常使用在需要除法取模的情况。

这里再次提一下排列、组合公式: C a b = a ! b ! ( a − b ) ! ,    A a b = a ! b ! C_a^b=\frac{a!}{b!(a-b)!},\ \ A_a^b=\frac{a!}{b!} Cab=b!(ab)!a!,  Aab=b!a!

求逆元需要用到快速幂:

LL qpow(LL a, LL b, LL p)
{
	LL res = 1;
	while (b)
	{
		if (b & 1) res = res * a % p;
		b >>= 1;
		a = a * a % p;
	}
    return res;
}

然后预处理阶乘和阶乘逆元:

f[0] = uf[0] = 1;
for (int i = 1; i < N; i ++ )
{
	f[i] = (LL)f[i - 1] * i % mod;
	uf[i] = (LL)uf[i - 1] * qpow(i, mod - 2, mod) % mod;
}

同样的,如果输出排列、组合结果的话需要利用公式。

时间复杂度: O ( n log ⁡ n ) O(n \log n) O(nlogn)

可以处理 1 0 5 10^5 105 以内规模的数据

思考:读者也可以尝试写 O ( n ) O(n) O(n) 预处理阶乘逆元。

算法4:Lucas 定理

思路

由 Lucas 定理可得:当 p p p 为质数时,

C a b = C a p b p × C a   m o d   p b   m o d   p \large{C_a^b=C_{\frac{a}{p}}^{\frac{b}{p}} \times C_{a \bmod p}^{b \bmod p}} Cab=Cpapb×Camodpbmodp

因此,我们可以写一个递归函数 LL lucas(int a, int b),递归出口是 a k < p , b k < p a_k<p, b_k<p ak<p,bk<p

递归的过程相当于自上向下将 C a 1 b 1 , C a 2 b 2 , … , C a k b k C_{a_1}^{b_1},C_{a_2}^{b_2},…,C_{a_k}^{b_k} Ca1b1,Ca2b2,,Cakbk 添加到乘式里。

LL 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;
}

这里面的 C(a, b) 是指 算法3 ,即用阶乘和阶乘逆元求组合数。

LL qpow(LL a, LL b, LL p)
{
	int res = 1;
	while (b)
	{
		if (b & 1) res = res * a % p;
		b >>= 1;
		a = a * a % p;
	}
	return res;
}

LL C(LL a, LL b)
{
	LL res = 1;
	for (int i = 1, j = a; i <= b; i ++ , j -- )
	{
		res = (LL)res * j % p;
		res = (LL) res * qpow(i, p - 2, p) % p;
	}
	return res;
}

同样的,如果输出排列、组合结果的话需要利用公式。

时间复杂度: O ( p × log ⁡ p n ) O(p \times \log_p n) O(p×logpn)

可以处理 a , b ≤ 1 0 18 , p ≤ 1 0 5 a,b \le 10^{18},p \le 10^5 a,b1018,p105 以内规模的数据


c++排列组合,笔记,c++,lucas定理,费马小定理,逆元,递推

最后,如果觉得对您有帮助的话,点个赞再走吧!文章来源地址https://www.toymoban.com/news/detail-643645.html

到了这里,关于【C++】 排列与组合算法详解(进阶篇)的文章就介绍完了。如果您还想了解更多内容,请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章,希望大家以后多多支持TOY模板网!

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处: 如若内容造成侵权/违法违规/事实不符,请点击违法举报进行投诉反馈,一经查实,立即删除!

领支付宝红包 赞助服务器费用

相关文章

  • java实现排列组合算法

    我这里只写了组合的算法。         假设现有 M=4 个数据 a,b,c,d。从中随机抽取n个数,n为1—4个数据进行组合。那么数学中的计算组合方式为C(4,1) + C(4,2) + C(4,3) + C(4,4)  = 4 + 6 + 4 + 1 = 15。那么共有15种组合方式。 方案一:此方法容易理解但是效率慢         我的做

    2024年02月13日
    浏览(47)
  • 算法--排列组合

    2024年02月16日
    浏览(36)
  • 数论——卢卡斯定理、求组合数 学习笔记

    温馨提示:组合数一般较大,下面的示范代码均无视数据范围,如果爆 int 请自行开 long long 或高精度处理。 从 (n) 个不同元素中,任取 (m) 个元素组成一个集合,叫做从 (n) 个不同元素中取出 (m) 个元素的一个组合;从 (n) 个不同元素中取出 (m leq n) 个元素的所有组合

    2024年02月08日
    浏览(39)
  • 算法第二期——排列组合(Python)

    目录 排列函数permutations( ) 易错点 组合函数combinations( ) 手写排列和组合代码  1.暴力法/

    2023年04月18日
    浏览(34)
  • 【算法教程】排列与组合的实现

    在讲排列与组合之前,我们先定义数据元素类型Fruit 对N个不同元素进行排序,总共有多少不同的排列方式? 例子:某水果店有以下水果,请对所有水果进行全排列,请输出所有排列 排列算法的javascript实现模板(DSF,最优解in-place) 测试结果 对N个不同元素进行排序,总共有多少不

    2024年02月07日
    浏览(43)
  • C++前缀和算法的应用:DI序列的有效排列的原理、源码及测试用例

    C++算法:前缀和、前缀乘积、前缀异或的原理、源码及测试用例 包括课程视频 给定一个长度为 n 的字符串 s ,其中 s[i] 是: “D” 意味着减少,或者 “I” 意味着增加 有效排列 是对有 n + 1 个在 [0, n] 范围内的整数的一个排列 perm ,使得对所有的 i: 如果 s[i] == ‘D’,那么

    2024年02月08日
    浏览(55)
  • 详解Python中的排列组合生成器

    在实际的开发场景中,经常需要遍历多个数组中的元素,将它们组合在一起使用。要取完所有可能的组合,最基本的方法是使用嵌套的循环,有多少个数组就嵌套多少层循环。嵌套循环是基本的原理,但不够简洁,Python中有更优雅的方式来实现这种功能。 在Python的内置模块

    2024年02月10日
    浏览(44)
  • 【C/C++练习】经典的排列组合问题(回溯算法)——电话号码的字母组合

    📖题目描述 题目出处 :电话号码的字母组合 示例: 📖题解  这是一道典型的排列组合问题,根据输入,我们需要找到所有的组合。下面以输入字符串 digits = \\\"23\\\" 为例来讲解这道题目。 图解: 分析:  首先要知道输入的字符串 \\\"23\\\" 中的数字字符分别对应哪些字符串,其中

    2024年02月16日
    浏览(39)
  • (3)回溯算法团灭子集、组合、排列问题 -- 元素无重可复选

    回溯算法团灭子集、组合、排列问题 根据元素是否重复和是否可复选,分为以下三种变体: 1、元素无重不可复选 2、元素有重不可复选 3、元素无重可复选 该篇给出第三种情况的代码,另外两种的实现见上方变体链接。

    2024年02月09日
    浏览(30)
  • 【华为OD机试真题 】1011 - 第K个排列 (JAVA C++ Python JS) | 机试题+算法思路+考点+代码解析

    🍂个人博客首页: KJ.JK   🍂专栏介绍: 华为OD机试真题汇总,定期更新华为OD各个时间阶段的机试真题,每日定时更新,本专栏将使用Java语言进行更新解答,包含真题,思路分析,代码参考,欢迎大家订阅学习 🎃题目描述

    2023年04月25日
    浏览(71)

觉得文章有用就打赏一下文章作者

支付宝扫一扫打赏

博客赞助

微信扫一扫打赏

请作者喝杯咖啡吧~博客赞助

支付宝扫一扫领取红包,优惠每天领

二维码1

领取红包

二维码2

领红包