【C++题解】【NOIP2003PJ】麦森数

题目链接：https://gmoj.net/junior/#main/show/1238

题目描述

输入

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

输出

更新：此题输出两行
1. 第一行输出位数
2. 第二行输出结果的最后五百位

样例输入

1279

样例输出

386

00000000000000000000000000000000000000000000000000

00000000000000000000000000000000000000000000000000

00000000000000104079321946643990819252403273640855

38615262247266704805319112350403608059673360298012

23944173232418484242161395428100779138356624832346

49081399066056773207629241295093892203457731833496

61583550472959420547689811211693677147548478866962

50138443826029173234888531116082853841658502825560

46662248318909188018470682222031405210266984354887

32958028878050869736186900714720710555703168729087

数据范围限制

如题

做法

这里给出 $3$ 个求位数的公式：

• $n/lo{g}_2(10)+1$

• $nlo{g}_{10}(2)+1$

• $nlog(2)/log(10)+1$

直接暴力模拟，肯定会时间超限， 但可以用压位的方法来解决。

定义一个数组 a 来存麦森数后 $500$ 位，会因为进行 $p\times 500$ 次乘 $2$ 的操作从而时间超限。

我们可以利用压位的方法让 $500$ 位数分成 $50$ 个部分，即数组 $a_i$ 存 $10$ 位，使得程序最多只进行 $p\times 50$ 次乘 $2$ 的操作，从而大大地提高了时间的效率。

#include<cmath>
#include<cstdio>
long long a[55]={0,1},x;
int p,l=1;
int main()
{
	scanf("%d",&p);
	for(int i=1;i<=p;i++)
	{
		x=0;
		for(int j=1;j<=50;j++)
		{
			a[j]=a[j]*2+x;
			//进行压位的操作
			x=a[j]/10000000000;
			a[j]%=10000000000;
		}
	}
	a[1]--;
	printf("%d\n",(int)(p/log2(10)+1));
	for(int i=50;i>=1;i--) printf("%010lld",a[i]);//注意输出
	return 0;
}

写文章不易，望大家给 JG_DF_发个评论！