高精度快速幂

原题：

洛谷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;
	}
}