noip2003普及组 麦森数

题目

麦森数

算法

(高精度,快速幂,数学) $O(500^{2}logP)$

首先考虑如何求位数。

我们发现在 $10^k\sim 10^{k+1}-1$,$(k\ge 0)$之间的数均有 $k+1$ 位。因此对于任意正整数 $x$，它的位数是 $\lfloor log10x\rfloor +1$。

而由于2的整次幂的末位数字不为0，因此 $2^P-1$的位数和 $2^P$ 的位数相同，所以 $2^P-1$ 的位数是 $\lfloor lo{g}_{10}{2}^P\rfloor +1=\lfloor Plo{g}_{10}2\rfloor +1$。

cmath库中有函数 log10()，直接使用即可。

然后考虑如何求最后500位数。

• 直接乘 $P$ 次，时间复杂度是 $500P=1.5\times 10^9$
• 用快速幂，时间复杂度是 $500^{2}logP=5\times 10^6$

因此我们选择第二种方法。

时间复杂度

算法瓶颈在于求用快速幂 + 高精度，计算量是 $O(500^{2}logP)$。

C++代码

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

using namespace std;

const int N = 500;

int p;
int a[N], tmp[N], res[N];

void mul(int c[], int a[], int b[])
{
    memset(tmp, 0, sizeof tmp);

    for(int i = 0; i < N; i++)
        for(int j = 0; i + j < N; j++)
            tmp[i + j] += a[i] * b[j];

    for(int i = 0, t = 0; i < N; i++)
    {
        t += tmp[i];
        c[i] = t % 10;
        t /= 10;
    }
}

void qui(int k)
{
    res[0] = 1;
    a[0] = 2;

    while(k)
    {
        if(k & 1)
            mul(res, res, a);
        mul(a, a, a);
        k >>= 1;
    }

    res[0]--;
}

int main()
{

    cin >> p;

    cout << int(p * log10(2)) + 1 << endl;

    qui(p);

    for(int i = 0, k = N - 1; i < 10; i++)
    {
        for(int j = 0; j < 50; j++, k--)
            cout << res[k];
        cout << endl;
    }

    return 0;
}