洛谷P1045题解

一、问题引出

[NOIP2003 普及组] 麦森数

题目描述

形如 $2^P-1$ 的素数称为麦森数，这时 $P$ 一定也是个素数。但反过来不一定，即如果 $P$ 是个素数，$2^P-1$ 不一定也是素数。到 1998 年底，人们已找到了 37 个麦森数。最大的一个是 $P=3021377$，它有 909526 位。麦森数有许多重要应用，它与完全数密切相关。

任务：输入 $P(1000<P<3100000)$，计算 $2^P-1$ 的位数和最后 $500$ 位数字（用十进制高精度数表示）

输入格式

文件中只包含一个整数 $P(1000<P<3100000)$

输出格式

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

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

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

样例 #1

样例输入 #1

1279

样例输出 #1

386
00000000000000000000000000000000000000000000000000
00000000000000000000000000000000000000000000000000
00000000000000104079321946643990819252403273640855
38615262247266704805319112350403608059673360298012
23944173232418484242161395428100779138356624832346
49081399066056773207629241295093892203457731833496
61583550472959420547689811211693677147548478866962
50138443826029173234888531116082853841658502825560
46662248318909188018470682222031405210266984354887
32958028878050869736186900714720710555703168729087

提示

【题目来源】

NOIP 2003 普及组第四题

二、原始思路

我一开始以为这就是一道简单的高精度问题，于是我直接拿一个超大数组去存计算过程当中的每个数，完全没有考虑到超时的问题。事实证明，计算过程完全可以再简化。

#include <iostream>
#include <cmath>
using namespace std;
int a[10000000]={0};
int main()
{
    int p,c=1;
    cin>>p;
    a[c]=1;
    for (int i = 0; i < p; i++)
    {
        int jw=0,j;
        for (j = 1; j <= min(c,500); j++)
        {
            int sum=a[j]*2+jw;
            a[j]=sum%10;
            jw=sum/10;
        }        
        if (c<=500)
        {
            while (jw>0)
            {
                a[j]=jw%10;
                jw/=10;
                j++;
            }
        }
        c=j-1;
    }  
    a[1]-=1;
    cout<<c<<endl;
    for (int i = 500; i >= 1; i--)
    {
        printf("%d",a[i]);
        if ((501-i)%50==0)
        {
            printf("\n");
        }
    }
}

三、正确思路

总结来说，此题有以下几个点需要重点注意：

• 如何求位数？

首先我们知道$2^p$与$2^p-1$有着相同的位数，因为2的次方满足了最后一位不为零的要求，所以减一后位数并不会改变，那么我们可以直接求$2^p$的位数。那么怎么求位数呢？我们不妨设$k=2^p$，根据$10^n$的位数为$n+1$，我们只要想办法把$k=2^p$中的底数2改为10，指数加一就是位数了。由此想到用10的几次方来代替2，那么就不难想到$10^{lo{g}_{10}2}=2$，这样便可以把$k=2^p$中的2代换掉，变为$k=(10^{lo{g}_{10}2}{)}^p$。根据乘方的原理，将p乘进去，原式便可化为我们最终想要的形式$k=10^{p\ast lo{g}_{10}2}$了，所以位数就是$p\ast lo{g}_{10}2+1$。（提醒一下，C++中cmath库自带log10()函数…）

• 求最后500位数

其实上面的题解是正确的，也就是说我们只用去管最后500位，其它的位我们是没有必要去管的

#include <iostream>
#include <cmath>
using namespace std;
long long a[501]={0};
long long cp(int n)
{
    long long result=1;
    for (int i = 1; i <= n; i++)
    {
        result*=2;
    }
    return result;
}
int main()
{
    int p;
    cin>>p;
    int num1=p/32;
    int num2=p-num1*32;
    long long x=cp(32);
    a[500]=1;
    /**
     * @brief 
     * 为了加快速度，我们在前num1轮对每个数乘以2^32，由于p/32并不一定刚好整除，因此我们在剩下的num1轮当中对每个数每次乘2
     */
    for (int i = 1; i <= num1; i++)
    {
        for (int j = 500; j >= 1; j--)
        {
            a[j]*=x;
        }
        for (int j = 500; j >= 1; j--)
        {
            a[j-1]+=a[j]/10;
            a[j]%=10;
        }
    }
    for (int i = 1; i <= num2; i++)
    {
        for (int j = 500; j >= 1; j--)
        {
            a[j]*=2;
        }
        for (int j = 500; j >= 1; j--)
        {
            a[j-1]+=a[j]/10;
            a[j]%=10;
        }
    }
    a[500]--;
    cout<<(int)(p*log10(2)+1)<<endl;
    for (int i = 1; i <= 500; i++)
    {
        cout<<a[i];
        if (i%50==0)
        {
            cout<<endl;
        }
    }
}