原题:
洛谷P1045 [NOIP2003 普及组] 麦森数
输入格式
文件中只包含一个整数 P(1000<P<3100000)
输出格式
第一行:十进制高精度数 2^P−1的位数。
第 2∼11行:十进制高精度数 2^P−1的最后 500位数字。(每行输出 50位,共输出 10行,不足 500位时高位补0)
不必验证 2^P−1与 P是否为素数。
这道题就是高精度乘法与快速幂的结合,没有想的那么复杂
高进度乘法:
void mul(int* a, int* b,int* c)
{
int t[1000];
memset(t, 0, sizeof(t));
for (int i = 0; i < 500; i++)
{
for (int j = 0; j < 500; j++)
{
t[i + j] += a[i] * b[j];
t[i + j + 1] += t[i + j] / 10;//进位
t[i + j] %= 10;
}
}
memcpy(c, t, sizeof(t));
}
快速幂:
void quickpow(int p)
{
a[0] = 2;
b[0] = 1;
while (p)
{
if (p & 1)mul(b,a,b);
mul(a,a,a);
p >>= 1;
}
b[0]--;//这里是2^p-1中的-1
}
只是在这里要注意,因为在高精度乘法的那个函数里面,将 t memcpy给了 c ,而 c 在下面快速幂的位置,是a,b为实参,所以a,b数组的大小要设为2*500
当然,我们也可以用vector,然后高精度函数返回vector类型的 t ,这样a,b的大小设为500,t的大小设为1000就可以了,注意一下就行。
pass代码:
#include<iostream>
using namespace std;
int a[1000], b[1000];
void mul(int* a, int* b,int* c)
{
int t[1000];
memset(t, 0, sizeof(t));
for (int i = 0; i < 500; i++)
{
for (int j = 0; j < 500; j++)
{
t[i + j] += a[i] * b[j];
t[i + j + 1] += t[i + j] / 10;
t[i + j] %= 10;
}
}
memcpy(c, t, sizeof(t));
}
void quickpow(int p)
{
a[0] = 2;
b[0] = 1;
while (p)
{
if (p & 1)mul(b,a,b);
mul(a,a,a);
p >>= 1;
}
b[0]--;
}
int main()
{
int p;
cin >> p;
cout << (int)(log10(2) * p) + 1 << endl;
quickpow(p);
for (int i = 0, k = 499; i < 10; i++)
{
for (int j = 0; j < 50; j++, k--)
{
cout << b[k];
}
cout << endl;
}
}