洛谷：麦森数（P1045）C++

1. 题意

麦森数，其实就是题目中给定n，计算2的n次方的最后500位数，题意的理解上是非常简单的。但是如果用我们平常的循环去做，肯定会超时，此题需要用到两个算法：快速幂，即可解决，不过500位，还需用到高精度来解决。

2. 思路

• 首先输入数n，然后计算位数，有固定公式，记住就好。

printf("%d\n",(int)(1.0*n*log10(2))+1);

• 由于要高精度，因此定义数组，然后用快速幂的算法来解决，下面给大家写一下快速幂的代码。

//一般几次方都是非常大的，因此要用到取余，否则会爆
# include <iostream>
using namespace std;
int main ()
{
    int n;scanf("%d",&n);

    //主要模板
    while(n--)
    {
        long long a,b,p;scanf("%lld%lld%lld",&a,&b,&p);
        long long num=1;
        while(b)
        {
            if(b%2==1)  num=num*a%p;
            a=a*a%p;
            b=b/2;
        }


        cout<<num<<endl;
    }
    return 0;
}

• 在计算之中还有一些难点，主要是高精度，不过结果要求500位，因此进位等问题并没有去考虑，大家结合代码就知道我所讲的内容。然后就是利用两个库函数，memcpy：复制数组；memset：将所有变为0，两者的头文件都是<cstring>。

• 最后就是输出，还要注意题目给出的，按照题目所给输出。

3.源码

# include <iostream>
# include <cstring>
# include <math.h>
using namespace std;
int a[1010],b[1010],c[1010];


void jisuan1(int a[],int b[])
{
    memset(c,0,sizeof(c));
    for(int i=0;i<500;i++)
    {
        for(int j=0;j<500;j++)
        {
            c[i+j]+=a[i]*b[j];
            if(c[i+j]>9)
            {
                c[i+j+1]+=c[i+j]/10;
                c[i+j]%=10;
            }
        }
    }
    memcpy(a,c,sizeof(c));
}


void jisuan2(int b[])
{
    memset(c,0,sizeof(c));
    for(int i=0;i<500;i++)
    {
        for(int j=0;j<500;j++)
        {
            c[i+j]+=b[i]*b[j];
            if(c[i+j]>9)
            {
                c[i+j+1]+=c[i+j]/10;
                c[i+j]%=10;
            }
        }
    }
    memcpy(b,c,sizeof(c));
}


int main ()
{
    int n;cin>>n;
    b[0]=2;a[0]=1;//初始化
    printf("%d\n",(int)(1.0*n*log10(2))+1);
    //主要模板，计算
    while(n)
    {
        if(n%2==1)   jisuan1(a,b);
        jisuan2(b);
        n=n/2;
    }
    a[0]=a[0]-1;


    for(int i=499;i>=0;i--)
    {
        printf("%d",a[i]);
        if(i%50==0) cout<<endl;
    }
    return 0;
}