麦森数（数组模拟大数+快速幂）

麦森数

1. 题目描述
   求2^p-1的位数和后五百位数
2. 思路
   求位数的话可以直接推断
   12345的位数为5
   1000010位数为7
   多写几组就可以知道位数是lg(n)+1；对于2^p-1的位数直接代入即可。计算时因为只有100，1000，10000减一与本身位数差一，而2的幂次不会出现这样的值，所以直接求2的p次方的幂次即可。
   求后五百位数值
   对于我来说是个难题，尤其还是高精度问题，直接求肯定会超时
   求2的p次方可以想到快速幂，用快速幂解决了时间复杂度，但高精度还是个问题
   解决高精度，可以用数组模拟大数的方法，用数组储存数

 for(int i=0;i<500;i++)
        for(int j=0;j<500-i;j++)
        c[i+j]+=x[i]*y[j];
    for(int i=0;i<500;i++)
    {
        c[i+1]+=c[i]/10;
        c[i]%=10;
    }

3. 代码

#include<bits/stdc++.h>
using namespace std;
int a[510],b[510];//数组a储存后五百位数，数组b储存2的幂次
int p;
//数组模拟大数乘法
void mul(int x[],int y[])
{
    int c[510];
    memset(c,0,sizeof(c));
    for(int i=0;i<500;i++)
        for(int j=0;j<500-i;j++)
        c[i+j]+=x[i]*y[j];
    for(int i=0;i<500;i++)
    {
        c[i+1]+=c[i]/10;
        c[i]%=10;
    }
    for(int i=0;i<500;i++)
        x[i]=c[i];
}
void solve()
{
    memset(a,0,sizeof(a));
    memset(b,0,sizeof(b));
    a[0]=1;//乘法初始化为1
    b[0]=2;
    //快速幂实现求值
    while(p)
    {
        if(p&1)
            mul(a,b);
        p/=2;
        mul(b,b);
    }
    a[0]--;//2^p-1
}
int main()
{
    cin>>p;
    //公式直接求位数
    int n=p*(double)(log(2)/log(10))+1;
    cout<<n<<endl;
    solve();
    int ans=0;
    for(int i=499;i>=0;i--)
    {
        ans++;
        if(ans==50)//50位数一行
        {
            ans=0;
            cout<<a[i];
            cout<<endl;
        }
        else
            cout<<a[i];
    }
    //cout<<endl；//这个换行没有，还pe了两次
    return 0;
}