[NOIP2003 普及组] 麦森数

#include <algorithm>
#include <cmath>
#include <cstdio>
#include <cstring>
#include <iomanip>
#include <iostream>
#include <map>
#include <sstream>
#include <string>
#include <vector>
using namespace std;
vector<int> mul(vector<int> &A, vector<int> &B) { //高精度乘法函数
    int i, j;
    vector<int> C(A.size() + B.size());
    for (i = 0; i < A.size(); i++) {
        for (j = 0; j < B.size(); j++) {
            C[i + j] += A[i] * B[j];
        }
    }
    int temp = 0;
    for (i = 0; i < C.size(); i++) {
        temp += C[i];
        C[i] = temp % 10;
        temp /= 10;
    }
    while (C.size() < 500) C.push_back(0); //不足500位的补零
    if (C.size() > 500) {    //高于500位直接用resize函数删除后面的位。
        C.resize(500);        //如果不处理会导致超时，因为结果只需考虑后面的500位，超出的部位可以直接删掉，不会影响后面500位存储的数字。
    }
    return C;
}
vector<int> poww(int b) {
    vector<int> A, B;
    string a = "1";
    string base = "2";
    int i;
    for (i = a.size() - 1; i >= 0; i--) A.push_back(a[i] - '0');
    for (i = base.size() - 1; i >= 0; i--) B.push_back(base[i] - '0');
    while (b != 0) {        //快速幂处理
        if (b & 1 != 0) A = mul(A, B);
        B = mul(B, B);
        b >>= 1;
    }
    return A;
}

int main() {
    vector<int> C;
    int P, N, i;
    cin >> P;
    N = P * log10(2) + 1;
    cout << N << endl;
    C = poww(P);
    C[0]--;
    int cnt = 0;
    for (i = C.size() - 1; i >= 0; i--) {
        cout << C[i];
        cnt++;
        if (cnt / 50) {
            cout << endl;
            cnt = 0;
        }
    }
    return 0;
}