题目描述
形如 2 P − 1 2^{P}-1 2P−1 的素数称为麦森数,这时 P P P 一定也是个素数。但反过来不一定,即如果 P P P 是个素数, 2 P − 1 2^{P}-1 2P−1 不一定也是素数。到 1998 年底,人们已找到了 37 个麦森数。最大的一个是 P = 3021377 P=3021377 P=3021377,它有 909526 位。麦森数有许多重要应用,它与完全数密切相关。
任务:输入 P ( 1000 < P < 3100000 ) P(1000<P<3100000) P(1000<P<3100000),计算 2 P − 1 2^{P}-1 2P−1 的位数和最后 500 500 500 位数字(用十进制高精度数表示)
输入格式
文件中只包含一个整数 P ( 1000 < P < 3100000 ) P(1000<P<3100000) P(1000<P<3100000)
输出格式
第一行:十进制高精度数 2 P − 1 2^{P}-1 2P−1 的位数。
第 2 ∼ 11 2\sim 11 2∼11 行:十进制高精度数 2 P − 1 2^{P}-1 2P−1 的最后 500 500 500 位数字。(每行输出 50 50 50 位,共输出 10 10 10 行,不足 500 500 500 位时高位补 0 0 0)
不必验证 2 P − 1 2^{P}-1 2P−1 与 P P P 是否为素数。
样例 #1
样例输入 #1
1279
样例输出 #1
386
00000000000000000000000000000000000000000000000000
00000000000000000000000000000000000000000000000000
00000000000000104079321946643990819252403273640855
38615262247266704805319112350403608059673360298012
23944173232418484242161395428100779138356624832346
49081399066056773207629241295093892203457731833496
61583550472959420547689811211693677147548478866962
50138443826029173234888531116082853841658502825560
46662248318909188018470682222031405210266984354887
32958028878050869736186900714720710555703168729087
提示
【题目来源】
NOIP 2003 普及组第四题
思路
求位数的推导:参考大佬题解(同时这也是另一种做法)。
输出后500位: 参考大佬题解,真的是太妙了。
具体见代码注释。
代码
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cmath>// 这个不要漏了,log要用
#define endl '\n'
using namespace std;
typedef unsigned long long ULL;// 范围是0 ~ 2^64 - 1
const int N = 510;
int p;
ULL a[N];
int main() {
ios::sync_with_stdio(false), cin.tie(nullptr), cout.tie(nullptr);
cin >> p;
// 输出位数
cout << (int)(p * log10(2) + 1) << endl;
// 初始化乘积为1,即最低位为1
a[0] = 1;
for (; p > 0; p -= 60) {// 每次乘上2^60次方
ULL t = 0;// 进位
// 每次都先乘到最低位上,在不断往前进位,这样就可以只算出最低500位
// 2^60 * 个位数 < 2^64 - 1,不会溢出
for (int i = 0; i < 500; i++) {
if (p >= 60) {
a[i] <<= 60;
}
else {// 如果剩不到60次方可以乘
a[i] <<= p;
}
// 处理进位
a[i] += t;
t = a[i] / 10;
a[i] %= 10;
}
}
//算的是 2^p - 1
a[0] -= 1;
int cnt = 500;
for (int i = 0; i < 10; i++) {
for (int j = 0; j < 50; j++) {
cout << a[--cnt];
}
cout << endl;
}
return 0;
}