倍增思想——理解和使用

倍增思想是一种通过成倍增长的方式高效解决问题的算法思想，核心是利用 “预先计算” 和 “跳跃式查询” 来减少重复计算，从而将时间复杂度从线性级（如 O (n)）优化到对数级（如 O (log n)）。

核心原理
倍增思想的本质是利用二进制的特性，将一个复杂问题拆解为若干个规模成倍增长的子问题，从而快速得到结果。

1.P1226 【模板】快速幂

链接: link

题目描述
给你三个整数 $a,b,p$，求 $a^b\mathrm{mod}p$。
输入格式
输入只有一行三个整数，分别代表 $a,b,p$。
输出格式
输出一行一个字符串 a^b mod p=s，其中 $a,b,p$ 分别为题目给定的值， $s$ 为运算结果。
输入输出样例
输入

2 10 9

输出

2^10 mod 9=7

说明/提示
样例解释
$2^{10}=1024$，$1024\mathrm{mod}9=7$。
数据规模与约定
对于 $100\%$ 的数据，保证 $0\le a,b<2^{31}$，$a+b>0$，$2\le p<2^{31}$。

【解法】
利⽤「二进制」以及「倍增」的思想，通过⼀个具体的例子说明，比如 3¹¹ ：
幂运算的运算法则：a^(b+c) = a^b × a^c;
可以通过将指数转化成其二进制形式，写成若干数相加。
11 = (1011)₂ = 1 × 2³ + 0 × 2² + 1 × 2¹ + 1 × 2⁰;
两者结合： 3 =
11 = 3⁽¹⁰¹¹⁾₂ = 3^(1 × 2³ + 0 × 2² + 1 × 2¹ + 1 × 2⁰)^ = 3⁸ × 3⁰ × 3² × 3¹ ;
如何快速算出 3¹ , 3² , 3³ , …, 3^log₂n
。其实很简单，从前往后看，后⼀个数是前⼀个数的平方：
3¹ = 3；
3² = 3¹ × 3¹ = 9；
3⁴ = 3² × 3² = 81;
因此，计算 3¹¹ ，我们只需将 11 的二进制表示中 1 所对应的幂乘起来即可。
如何实现这个算法呢，以 ab 为例：
• 一边提取 b 的二进制表示中的每一位；
• 一边让 a = a ∗ a ，不断变成之前的平方（倍增的思想）；
• 在提取 b 的二进制表示时，如果这一位是 1 ，就乘上对应位置的权值。
代码：

#include <iostream>
using namespace std;
typedef long long LL;
LL a, b, p;

// a^b^ % p 的值
LL quickpow(LL a, LL b, LL p)
{
	LL ret = 1;
	while (b)
	{
		if (b & 1) ret = ret * a % p;
		a = a * a % p;
		b >>= 1;//提取b的二进制位
	}
	return ret;
}
int main()
{
	cin >> a >> b >> p;
	printf("%lld^%lld mod %lld=%lld\n", a, b, p, quickpow(a, b, p));
	return 0;
}

P10446 64位整数乘法

链接: link

题目描述
求 $a$ 乘 $b$ 对 $p$ 取模的值。
输入格式
第一行输入整数 $a$，第二行输入整数 $b$，第三行输入整数 $p$。
输出格式
输出一个整数，表示 a*b mod p 的值。
输入输出样例
输入

3
4
5

输出

2

说明/提示
$1\le a,b,p\le 10^{18}$

【解法】
思路同第一题
如何实现这个算法呢，以 a × b 为例：
• 一边提取 b 的二进制表示中的每一位；
• 一边让 a = a + a ，不断变成之前的两倍（倍增的思想）；
• 在提取 b 的二进制表示时，如果这一位是 1 ，就加上对应位置的权值。

#include <iostream>
#include <cstdio>
using namespace std;

typedef long long LL;
LL a, b, p;
//防止溢出
LL qmul(LL a,LL b,LL p)
{
	LL sum = 0;
	while(b)//枚举b的二进制位
	{
		if(b & 1) sum = (sum + a) % p;
		a = (a + a) % p;//计算下一个权值
		b >>= 1;
	}
	return sum;
}
int main()
{
	cin >> a >> b >> p;
	cout << qmul(a,b,p) << endl;
	return 0;
}