快速幂（递归与非递归）

最新推荐文章于 2024-09-08 15:12:17 发布

x_miracle

最新推荐文章于 2024-09-08 15:12:17 发布

阅读量207

点赞数

分类专栏：分治文章标签：算法 c++

本文链接：https://blog.csdn.net/qq_42783346/article/details/126238732

版权

分治专栏收录该内容

2 篇文章 0 订阅

订阅专栏

本文介绍了快速幂(Binary Exponentiation)算法，这是一种在Ο(logn)时间内计算乘方的高效方法。通过递归和位操作两种方式详细阐述了算法的实现过程，包括将指数转化为二进制并逐位进行乘法运算的优化策略。

摘要由CSDN通过智能技术生成

快速幂(Binary Exponentiation)

它可以以 $\Theta (log⁡n)$ 的时间复杂度计算乘方。

法一：递归

原理：

若 $n$ 为偶数：
$\begin{equation*} \begin{aligned} a^n &= a^{n/2} * a^{n/2} \\ &= a^{n/4} *a^{n/4}*a^{n/4}*a^{n/4} \\ \end{aligned} \end{equation*}$
运用了二分的思想，使求幂的过程尽可能少。

递归边界为： $a^0=1$

long long binaryPow(long long a, long long b)
{
	if(b == 0)	return 1;
	else if(b%2==1)	return a*binaryPow(a,b-1);
	else
	{
		long long num = binaryPow(a, b/2);
		return num * num;
	}
	// 不写成return binaryPow(a, b/2) * binaryPow(a, b/2)
}

long long binaryPow(long long a, long long b)
{
  if (b == 0) return 1;
  long long res = binaryPow(a, b / 2);
  if (b % 2)
    return res * res * a;
  else
    return res * res;
}

法二：

我们把幂数 $b$ 写成二进制 $b=(b_nb_{n-1}...b_2b_1)_2$ 。

有： $b=(b_n*2^n+b_{n-1}*2^{n-1}+...+b_2*2+b_1*1)$ 。其中，$b_i = 0,1 $。

所以， $a^b=a^{(b_nb_{n-1}...b_2b_1)_2}$

而 $n$ 有 $\lfloor log_2n \rfloor +1$ 个二进制位，由此，幂运算得到优化：

long long binaryPow(long long a, long long b) 
{
  long long res = 1;
  while (b > 0) 
  {
    if (b & 1) res = res * a;//如果b的二进制末位为1。
    a = a * a;//对应的二进制位倍数。
    b >>= 1;//b右移1位。
  }
  return res;
}