P1045 麦森数【快速幂】

题目训练

洛谷

P1045 麦森数【快速幂】

这是一道与幂有关、涉及高精度的题目，原理还是快速幂，分两步走，第一问计算pow(2,p)-1的位数可以参照下面的方法思考，最后得到位数为p*long10(2)+1

第二问快速幂运算的原理我是参照的下面这篇博客，讲的很详细
原文链接: https://blog.csdn.net/Sirius_han/article/details/88757132.

完整代码附上

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <math.h>

#define INF 0x3f3f3f3f
using namespace std;
typedef long long ll;
const int N = 1100;
int p,a[N];
 
void poww()
{
    int c[N];
    memset(c,0,sizeof(c));
    for(int i = 1;i <= 500;i ++)
    {
        for(int j = 1;j <= 500;j ++)
            c[i+j-1] += a[i] * a[j];
    }
    for(int i =1;i <= 500;i ++)
        if(c[i]/10 !=0)
        {
            c[i+1] += c[i] / 10;
            c[i] %= 10;
        }
    for(int i =0 ;i <= 500;i ++)
        a[i] = c[i];
}
void add_one()
{
    for(int i = 1;i <= 500;i ++)
        a[i] *= 2;
    for(int i =1;i <= 500;i ++)
        if(a[i] / 10 !=0)
        {
            a[i+1] += a[i] / 10;
            a[i] %= 10;
        }
}

//快速幂，只不过乘法的时候需要十分的精确
void get_num(int x)
{
    if(x/2 != 1)       //每次都去除以2，知道结果得到一的时候，就不再进行递归的操作
        get_num(x/2);
    poww();
    if(x %2 == 1)      //如果最后还余下一个1的话，这样就需要再乘以一个2，和快速幂的方法类似
        add_one();
}
 
int main()
{
    scanf("%d",&p);
    printf("%d\n",ceil(p*log10(2)));  //第一问求位数
    a[1] = 2;
    get_num(p);
    a[1] -- ;
    for(int i = 500;i >=1;i --)
    {
        printf("%d",a[i]);
        if((i-1)%50==0)
            printf("\n");
    }
    return 0;
}