洛谷P1045 [NOIP2003 普及组] 麦森数 C语言 详解 简单易懂

题目描述

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

任务：从文件中输入PP（1000<P<31000001000<P<3100000），计算2^{P}-12P−1的位数和最后500位数字（用十进制高精度数表示）

输入格式

文件中只包含一个整数PP（1000<P<31000001000<P<3100000）

输出格式

第一行：十进制高精度数2^{P}-12P−1的位数。

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

不必验证2^{P}-12P−1与PP是否为素数。

输入输出样例

输入 #1复制

1279

输出 #1复制

386
00000000000000000000000000000000000000000000000000
00000000000000000000000000000000000000000000000000
00000000000000104079321946643990819252403273640855
38615262247266704805319112350403608059673360298012
23944173232418484242161395428100779138356624832346
49081399066056773207629241295093892203457731833496
61583550472959420547689811211693677147548478866962
50138443826029173234888531116082853841658502825560
46662248318909188018470682222031405210266984354887
32958028878050869736186900714720710555703168729087

说明/提示

【题目来源】

NOIP 2003 普及组第四题

AC代码如下

#define _CRT_SECURE_NO_WARNINGS 1
#include <stdio.h>
#include <math.h>
#include <string.h>
/***********************************************
假设2^p有n位，即2^p = 10^n + a  （a<10)
n = log10(2^p-a)=log10(2^p)=p*log(2)
输出的位数由math.h中的log10函数计算
p = a(0)*2^0 + a(1)*2^1 +....+ 2^n (其中a(n)是p的各二进制位)
2^p = 2^(a(0)*2^0) * ... * 2^(2^n)
为了避免超时，大整数以10000进制存储
************************************************/
/***********************************************
大整数a，b相乘，末500位存放在a中
************************************************/
void mul(int* a, int* b)//高精度  A*B Problem
{
    int i, j, res[125];
    memset(res, 0, sizeof(res));//初始化数组  res暂存a*b结果
    for (i = 0; i < 125; i++)
        for (j = 0; j < 125 - i; j++)
        {
            res[i + j] += a[i] * b[j];
            res[i + j + 1] += res[i + j] / 10000;
            res[i + j] %= 10000;
        }
    memcpy(a, res, sizeof(res));//把结果从res赋值到a数组中
}
int p, result[125], base[125]/**用来存放2^(2^n)**/;  //result是存结果的，p是指数，base是底数
int main()
{
    int i;
    scanf("%d", &p);
    printf("%d\n", (int)(p * log10(2.0)) + 1);//+1是因为n = log10(2^p-a)=log10(2^p)=p*log(2) 时省去了a，如3*10^4 4+1位数
    result[0] = 1;//存结果
    base[0] = 2;//存底数，起初是2
    while (p)//高精度快速幂 最优版
    {
        if (p & 1)//相当于p%2==1  如果是奇数次幂，如2^5 
            mul(result, base); //高精度result = result * base ;  单独运算  result*底数^1 即resule*2^1                
        p >>= 1;//指数右移实现指数折半，p/=2;
        mul(base, base); //base = (base * base) ; 底数变大成原来的平方 mul是高精度乘法函数
    }
    result[0]--;
    //result一共125位，每一位都是四位数，125*4=500
    for (i = 124; i >= 0; i--)
        if (i % 25 == 12)//
            printf("%02d\n%02d", result[i] / 100, result[i] % 100);
        else
        {
            printf("%04d", result[i]);
            if (i % 25 == 0)printf("\n");//要求输出500位，每行五十位
        }
    return 0;
}

关于高精度算法和快速幂，详可参见 洛谷P1303 A*B Problem （高精）C语言（代码详解和注释）_Prudento的博客-CSDN博客https://blog.csdn.net/Prudento/article/details/122227017?spm=1001.2014.3001.5501

快速幂算法（全网最详细地带你从零开始一步一步优化）_扬俊的小屋-CSDN博客_快速幂算法https://blog.csdn.net/qq_19782019/article/details/85621386