【NOIP】麦森数

【NOIP】麦森数

题目

题目描述

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

任务：从文件中输入P（1000<P<3100000），计算2 P −1的位数和最后500位数字（用十进制高精度数表示）

输入格式

文件中只包含一个整数PP（1000<P<31000001000<P<3100000）

输出格式

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

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

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

输入输出样例

输入 #1 复制
1279
输出 #1 复制
386
00000000000000000000000000000000000000000000000000
00000000000000000000000000000000000000000000000000
00000000000000104079321946643990819252403273640855
38615262247266704805319112350403608059673360298012
23944173232418484242161395428100779138356624832346
49081399066056773207629241295093892203457731833496
61583550472959420547689811211693677147548478866962
50138443826029173234888531116082853841658502825560
46662248318909188018470682222031405210266984354887
32958028878050869736186900714720710555703168729087

分析

第一个问题：需要用的 log10（）。。（ 头文件 #include< cmath > ）
2^p-1 和 2 ^p 的位数是相同的 比如我们求 x= 2 ^ p .
（根据 10^log10 =1 ，）可以转换成 x = ( 10^log10*2 ) ^p 。

第二个问题：求最后的500个数。
这个用高精度乘法就ok了。 （这个可以直接做模版，基本都能通用）

代码

#include<iostream>
#include<cstring> 
#include<cmath>
using namespace std;


int p;
int f[1001],res[1001],sav[1001];


void result_1(){
	memset(sav,0,sizeof(sav));
	
	//求每位上的值 
	for(int i=1;i<=500;i++){
		for(int j=1;j<=500;j++){
			sav[i+j-1] += res[i]*f[j];
		}
	}
	//进位 
	for(int i=1;i<=500;i++){
		sav[i+1] += sav[i]/10;
		sav[i] %= 10;
	}
	memcpy(res,sav,sizeof(res));
}

void result_2(){
	memset(sav,0,sizeof(sav));
	
	//求每位上的值 
	for(int i=1;i<=500;i++){
		for(int j=1;j<=500;j++){
			sav[i+j-1] += f[i]*f[j];
		}
	}
	//进位 
	for(int i=1;i<=500;i++){
		sav[i+1] += sav[i]/10;
		sav[i] %= 10;
	}
	memcpy(f,sav,sizeof(res));
}
int main(){
	cin>>p;
	cout<<(int)(log10(2)*p +1) <<endl;
	
	res[1] = 1;
	f[1] = 2;
	while(p!=0){
		if(p%2==1)result_1();
		
		p /= 2;
		result_2();
	}
	res[1] -= 1;       //我们求的是2^p,,题目要求 2^p -1； 
	for(int i=500;i>=1;i--){
		if(i!=500 && i%50 ==0) {
			cout<<endl<<res[i];
		}else{
			cout<<res[i];
		}
	}
	return 0;
}