[高精度] [快速幂] [洛谷] P1045 麦森数

最新推荐文章于 2024-08-08 12:45:00 发布

Zeolim

最新推荐文章于 2024-08-08 12:45:00 发布

阅读量318

点赞数

CC 4.0 BY-SA版权

分类专栏：题目高精度文章标签：洛谷 P1045 麦森数

本文链接：https://blog.youkuaiyun.com/Zeolim/article/details/81318654

题目同时被 2 个专栏收录

18 篇文章

订阅专栏

高精度

3 篇文章

订阅专栏

本文深入探讨了高精度计算中快速幂算法的应用，通过C++实现了一个能够处理大规模数值运算的程序。文章详细介绍了算法的基本原理，以及如何通过数组进行大数乘法和幂运算，适用于需要高精度计算的场景。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

基础的快速幂和高精度应用

数位不会求

只能百度

#include <iostream>
#include <cstring>
#include <cmath>
using namespace std;

typedef long long ll;

const int MAXN = 1e3 + 10;

int ans[MAXN] = {1};

int brr[MAXN] = {2};

void muti(int arr[], int brr[], int crr[])
{
    int ret[MAXN] = {0};

    for(int i = 0; i < 550; i++)
    {
        for(int j = 0; j < 550 && i + j < 600; j++)
        {
            ret[i + j] += (arr[i] * brr[j]);
        }
    }

    for(int i = 0; i < 550; i++)
    {
        ret[i + 1] += (ret[i] / 10);
        ret[i] %= 10;
    }

    memcpy(crr, ret, sizeof(ret));
}

void qpow(ll n)
{
    while(n)
    {
        if(n&1)
            muti(ans, brr, ans);
        muti(brr, brr, brr);
        n >>= 1;
    }
}

ll sisze(ll x)
{
    return int(log10(2) * x + 1);
}

void puts(ll n)
{
    cout<<sisze(n)<<endl;

    ans[0]--;

    for(int i = 499; i >= 0; i--)
    {
        cout<<ans[i];
        if(i != 499 && i % 50 == 0)
            cout<<endl;
    }
}

int main()
{
    int n;

    cin>>n;

    qpow(n);

    puts(n);

    return 0;
}