算法训练 麦森数 二分+高精度

题目链接                     蓝桥杯 算法训练---------题解

锦囊1

二分，高精度计算。

锦囊2

使用数组来保存答案的最后500位，实现乘法运算。计算幂时使用二分，则计算a^b先算a^(floor(b/2))，再平方一下，根据需要看是不是再乘a。

问题描述

　　形如2P-1的素数称为麦森数，这时P一定也是个素数。但反过来不一定，即如果P是个素数，2P-1不一定也是素数。到1998年底，人们已找到了37个麦森数。最大的一个是P=3021377，它有909526位。麦森数有许多重要应用，它与完全数密切相关。
　　任务：从文件中输入P（1000<P<3100000），计算2P-1的位数和最后500位数字（用十进制高精度数表示）

 分析：二分快速幂运算。位数：log10(2^p-1)向上取整 = p*log(2)+1.（2^p的结果个位数都是大于1的-1不会影响计算结果位数）

#include <cstdio>
#include <cmath>
#include <cstring>
using namespace std;
const int N = 500;
int a[N], b[N],res[N*N];
void multAarry(int *a,int *b){
	memset(res, 0, sizeof(res));	
	for(int i = 0; i < N; i++)
		for(int j = 0; i+j < N; j++){ // 超过N都不需要 
			res[i+j] += a[i]*b[j];
		}			
	for(int i = 0; i < N; i++) //最后计算规整，减少了重复判断 
		if(res[i] >= 10){
			res[i+1] += res[i] / 10;
			res[i] = res[i] % 10;
		}
	for(int i = 0;  i < N; i++) 
		a[i] = res[i];
}

void minusOne(int *a, int i){
	if(a[i] != 0) a[i]--;
	else{
		a[i] = 9;
		minusOne(a,i+1);	
	} 
}

int main(int argc, char** argv) {
	int n;
	scanf("%d",&n);
	printf("%d", int(n*log10(2)) + 1);
	a[0] = 1; //initialize;
	b[0] = 2;
	while(n){
		if(n%2 == 1) multAarry(a, b);
		multAarry(b,b);
		n /= 2;
	}
	minusOne(a,0); //-1
	for(int i = 500-1, cnt = 0; i >= 0; i--, cnt++){
		if((i+1)%50 == 0) printf("\n");
		printf("%d",a[i]);
	}
	return 0;
}

输入格式

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

输出格式

　　第一行：十进制高精度数2P-1的位数。
　　第2-11行：十进制高精度数2P-1的最后500位数字。（每行输出50位，共输出10行，不足500位时高位补0）
　　不必验证2P-1与P是否为素数。

样例输入

1279

样例输出

386
00000000000000000000000000000000000000000000000000
00000000000000000000000000000000000000000000000000
00000000000000104079321946643990819252403273640855
38615262247266704805319112350403608059673360298012
23944173232418484242161395428100779138356624832346
49081399066056773207629241295093892203457731833496
61583550472959420547689811211693677147548478866962
50138443826029173234888531116082853841658502825560
46662248318909188018470682222031405210266984354887
32958028878050869736186900714720710555703168729087