P1045 [NOIP2003 普及组] 麦森数

题目部分思路来源自洛谷用户：憧憬未来，sqrt_7

[NOIP2003 普及组] 麦森数

题目描述

形如 2^P-1的素数称为麦森数，这时P一定也是个素数。但反过来不一定，即如果P是个素数，2^P-1不一定也是素数。到 1998 年底，人们已找到了 37 个麦森数。最大的一个是 P=3021377，它有909526位。麦森数有许多重要应用，它与完全数密切相关。

任务：输入P(1000<P<3100000)计算2^P-1的位数和最后500位数字（用十进制高精度数表示）

输入格式

文件中只包含一个整数P(1000<P<3100000)

输出格式

第一行：十进制高精度数2^P-1的位数。

第2行：十进制高精度数2^P-1的最后500位数字。（每行输出50位，共输出10行，不足500位时高位补0）

不必验证2^P-1与P是否为素数。

样例 #1

样例输入 #1

1279

样例输出 #1

386
00000000000000000000000000000000000000000000000000
00000000000000000000000000000000000000000000000000
00000000000000104079321946643990819252403273640855
38615262247266704805319112350403608059673360298012
23944173232418484242161395428100779138356624832346
49081399066056773207629241295093892203457731833496
61583550472959420547689811211693677147548478866962
50138443826029173234888531116082853841658502825560
46662248318909188018470682222031405210266984354887
32958028878050869736186900714720710555703168729087

思路：

整体上有两个工作，一个是输出位数，另一个是输出后500位；

1.一个数的位数通常与10的幂指数有关，具体为n+1（n为幂指数），因此如果把2^P-1变为10^n，即可求出位数。2的幂指数保证最后一位不是0，因此-1不会造成退位，则求2^P的位数即可，而2=10^( log10(2) )，所以2^P =10^( P*log10(2) )，位数即为P*log10(2)+1；

2.500位会牵扯到高精度运算，如果依次乘2，数量过于庞大，一定会超时，对于这种情况，采取以下两个方法：快速幂，压位；

#include<bits/stdc++.h>
using namespace std;
int base[1000],n,ans[1000],sav[1000];
void doublebase(){//将2^n平方
	memset(sav,0,sizeof(sav));
	for(register int i=0;i<500;i++){//乘操作
		for(register int j=0;j<500;j++){
			sav[i+j]+=base[i]*base[j];
		}
	}
	for(register int i=0;i<1000;i++){//进位操作
		sav[i+1]+=sav[i]/10;
		sav[i]%=10;
	}
	memcpy(base,sav,sizeof(base));
}
void multi(){//进行乘法运算
	memset(sav,0,sizeof(sav));
	for(register int i=0;i<500;i++){
		for(register int j=0;j<500;j++){
			sav[i+j]+=ans[i]*base[j];
		}
	}
	for(register int i=0;i<1000;i++){
		sav[i+1]+=sav[i]/10;
		sav[i]%=10;
	}
	memcpy(ans,sav,sizeof(ans));
}
int main(){
	cin>>n;
	ans[0]=1,base[0]=2;
	cout<<(int)(log10(2)*n+1)<<endl;
	while(n!=0){//快速幂运算
		if(n&1)multi();
		doublebase();
		n>>=1;
	}
	ans[0]-=1;
	for(int i=499;i>=0;i--){
		cout<<ans[i];
		if(i%50==0)cout<<endl;
	}
}

快速幂：

1.快速幂基础详见：快速幂算法（全网最详细地带你从零开始一步一步优化）；

2.输出时注意要倒序输出；

#include<bits/stdc++.h>
using namespace std;
int main(){
	int n;
	unsigned long long num[500]={1};
	cin>>n;
	cout<<(int)(n*log10(2))+1<<endl;
	for(;n>0;n-=60){
		unsigned long long in=0;
		for(int i=0;i<500;i++){
			if(n>=60)num[i]<<=60;
			else num[i]<<=n;
			num[i]+=in;
			in=num[i]/10;
			num[i]%=10;
		}
	}
	num[0]--;
	for(int i=499;i>=0;i--){
		cout<<num[i];
		if(i%50==0)cout<<endl;
	}
}

压位：

1.压位的思想就是一次性乘2^60以减少运算规模；

小技巧：

1.在第一段代码中，采用了关键字“register”，目的是将变量储存在寄存器里，方便CPU运算，常用在定义大量使用的数据时，减少运算时间；

2.第二段代码中，采用了“<<=”的运算，意为“左移”，和其意义等同于*2^n，使代码更加简便。