麦森数

形如2^P−1的素数称为麦森数，这时P一定也是个素数。

但反过来不一定，即如果P是个素数，2^P−1不一定也是素数。

到1998年底，人们已找到了37个麦森数。

最大的一个是P=3021377，它有909526位。

麦森数有许多重要应用，它与完全数密切相关。 

任务：从文件中输入P，计算2^P−1的位数和最后500位数字（用十进制高精度数表示）。

输入格式

文件中只包含一个整数P。

输出格式

第一行：十进制高精度数2^P−1的位数。 

第2-11行：十进制高精度数2^P−1的最后500位数字。（每行输出50位，共输出10行，不足500位时高位补0） 

不必验证2^P−1与P是否为素数。

数据范围

1000<P<3100000

代码：

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <cmath>

using namespace std;

const int N = 510;

int a[N], c[N], res[N];

void mul(int a[N], int b[N], int res[N])
{
    memset(c, 0, N * 4);
    for (int i = 0; i < 500; i ++ )
        for (int j = 0; j < 500; j ++ )
            if (i + j < 500)
                c[i + j] += a[i] * b[j];
    for (int i = 0, t = 0; i < 500; i ++ )
    {
        t += c[i];
        res[i] = t % 10;
        t /= 10;
    }
}

void qmi(int p)
{
    res[0] = 1;
    a[0] = 2;
    while (p)
    {
        if (p & 1) mul(a, res, res);
        mul(a, a, a);
        p >>= 1;
    }
    res[0] -- ;
}


int main()
{
    int p;
    scanf("%d", &p);
    printf("%d\n", int(log10(2) * p) + 1);

    qmi(p);
    for (int i = 499, j = 0; j < 10; j ++ )
    {
        for (int k = 0; k < 50; k ++, i -- )
            printf("%d", res[i]);
        puts("");
    }
    return 0;
}