P1045麦森数

P1045麦森数

#include<iostream>
#include <cmath>
#include <cstring>

const int maxn = 1e5 + ;
#define int long long
using namespace std;

//第一个表示位数，之后的数字表示每个位数上的数
int ans[maxn] = {, }, k[maxn] = {, }, c[maxn];

void muti(int a[], int b[]) {
    memset(c, , sizeof(c));
    c[] = a[] + b[];
    if (c[] > ) c[] = ;
    for (int i = ; i <= b[]; i++) {//模拟乘法 32 * 23
//          3 2
//        * 2 3
//        ———————
//          9  6
//       6  4
//       ————————
//       7  3  6
        for (int j = ; j <= a[]; j++) {
            c[j + i - ] += b[i] * a[j];
            if (c[j + i - ] > ) {
                c[i + j] += c[i + j - ] / ;
                c[i + j - ] %= ;
            }
        }
    }
    for (int i = ; i <= c[]; i++) a[i] = c[i];
}

void quick(int x)//快速幂
{
    while (x) {
        if (x & ) muti(ans, k);
        muti(k, k);
        x >>= ;
    }
}

signed main() {
    int x;
    cin >> x;
    cout << (int) (x * log10() + );
    quick(x);
    ans[]--;//-1操作
    for (int i = ; i >= ; i--) {
        if (i %  == ) cout << endl;
        cout << ans[i];
    }
    return ;
}