(快速幂算法+高精度）洛谷P1045 麦森数

本文链接：https://blog.csdn.net/a1845613403/article/details/121863403

本文解析了遇到高精度运算时如何使用快速幂算法优化，介绍了`qpower()`和`grow()`函数，以及如何利用对数法则计算2的p次方的位数。重点讲解了如何避免低精度乘法模板限制，通过模拟人工运算解决新运算类型。适合对算法深入理解及高精度运算实践的同学。

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

前言

故事的最后，让我们以一道十分经典的题目——《麦森数》来结尾。接受现实吧，总会有我们没准备过的高精度运算出现。我们固然可以提前把高精度的快速幂模板也准备好，但是总会有百密一疏的时候，即使出现过一种没见过的高精度运算类型，也要做出合理的应对。

题目概述

在这里插入图片描述

AC代码

#include<iostream>
#include<cstring>
#include<cmath>
using namespace std;
#define maxn 1002
int mid[maxn],base[maxn],ans[maxn];
void qpower()//替代低精度快速幂算法中的相乘部分 
{
	memset(mid,0,sizeof(mid));
	for(int i=1;i<=500;++i)
	  for(int j=1;j<=500;++j)
	  mid[i+j-1]+=ans[i]*base[j];
 for(int i=1;i<=500;++i)
   mid[i+1]+=mid[i]/10,mid[i]%=10;
 memcpy(ans,mid,sizeof(mid));
}
void grow()//替代低精度快速幂算法中的base自平方部分 
{
	memset(mid,0,sizeof(mid));
	for(int i=1;i<=500;++i)
	  for(int j=1;j<=500;++j)
	  	mid[i+j-1]+=base[i]*base[j];
	for(int i=1;i<=500;++i)
	mid[i+1]+=mid[i]/10,mid[i]%=10;
	memcpy(base,mid,sizeof(base));
}
int main()
{
	int p;
	scanf("%d",&p);
	printf("%d\n",(int)(log10(2)*p+1));
	ans[1]=1,base[1]=2;//ans存的是答案，base存的是快速幂算法里已经积累的次数 
	while(p>0)//快速幂 
	{
		if(p&1)
		qpower();
		grow();
		p>>=1;
	}
	int sum=0;
	ans[1]-=1;//2的次方数末尾不会是0，可以直接减去1不用考虑借位的问题
	//(因为如果要出来0前一个末尾得是5，这显然不可能） 
	for(int i=500;i>=1;--i)
	{
		printf("%d",ans[i]);
		++sum;
		if(sum%50==0)
		printf("\n"); 
	}
  return 0;
}

分析思路

1.首先可能会想到用朴素的做法，因为大整数类里是有存储长度的，然后把2乘p次，但是这样有两个问题，一个是时间复杂度达到O(n),必然要TLE，另一个是题目暗示最后的位数可能会很大，数组也没有那么多位数，所以这题不能再使用Bigint模板了，而要自己实现一个更高效的乘方运算。并且只需要后五百位数字。所以啊，刚开始学要背好模板，但是题目越往后做，反而不能再被模板限制，有点“无招胜有招”那味。

2.快速幂算法是如何实现的？比如 $2^{10}$ ,朴素的做法是乘10次，但是我们需要吗？想想我们自己计算的时候，我们是直接用 $2^5*2^5$ ,如果 $2^5$ 也不会呢？那就考虑用更小的次方。快速幂其实类似于这样的思想，从小的次方开始指数扩大，这样就大大降低了复杂度。事实上，快速幂算法的时间复杂度是O( $log_n$ )的。这个在很多题目里都可以使用来做优化。一般地，对于一个低精度快速幂 $a^b$ (因为幂运算会很大，所以模板还加了取余k的操作)，有以下的做法：

#define ll long long
ll qpower(ll a,ll b)
{
	ll base=a,ans=1;
	while(b>0)
	{
		if(b&1)
		{
		ans*=base;
	  	ans%=k;
	  }
	  base*=base;
	  base%=k;
	  b>>=1;
	}
	return ans;
}

其中"b>>1"实际是通过二进制来检查有没有当前位可用。也就是说实际上是这样做的，先检查有没有 $a^1$ ,再检查有没有 $a^2$ ,再检查有没有 $a^4...$ 有的话就把他们相乘。结束的条件是b在二进制下全都是0了，也就是全部分解完毕。那么对于高精度的呢？这就需要用到高精度乘高精度做法了，对应的替换写在注释行里了，这里不再赘述。

3.那么位数如何计算呢(不用高精度Bigint存储长度的情况下)?首先， $2^p-1$ 的位数是等于 $2^p$ 次方的位数的。为什么呢？因为 $2^p$ 的末尾不会是0，所以不可能出现说减去1把位数给减少了。那为什么最后一位不可能是0？因为最后一位是0，那么 $2^{p-1}$ 末尾就得是5,这与它是偶数相矛盾。其次我们知道， $10^n$ 有 $n + 1$ 位。那么，有如下的变形：(用到了高中的指对变换方法)
$2^p={10^{log_{10}2}}^p$
所以，它的位数就是：
$log_{10}2^p+1=p*log_{10}2+1$