大数快速幂

最新推荐文章于 2024-11-18 15:59:13 发布

转载最新推荐文章于 2024-11-18 15:59:13 发布 · 399 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

原文链接：https://www.cnblogs.com/NKHTH/p/14372982.html

标签

#c++ #算法

本文介绍了如何使用C++实现快速幂算法，通过优化计算过程解决大数幂运算问题，包括operate1和operate2函数的详细实现。适用于任意a^b的计算需求。

转自11.C++：快速幂 - 黄石forever - 博客园

在快速幂的基础上实现大数求幂方法，可以满足任意a^b的求法

#include<iostream>
#include<algorithm>
#include<iomanip>
#include<math.h>
#include<cstring>
#include<cstdio>
#include<cctype>
using namespace std;
int index,base_number[1001], answer[1001],temp[1001];
void operate1()
{
    memset(temp, 0, sizeof(temp));
    for (int i = 1; i <= 500; i++)
        for (int j = 1; j <= 500; j++)
            temp[i + j - 1] += answer[i] * base_number[j];
    for (int i = 1; i <= 500; i++)
    {
        temp[i + 1] += temp[i] / 10;
        temp[i] %= 10;
    }
    memcpy(answer, temp, sizeof(temp));
}
void operate2()
{
    memset(temp, 0, sizeof(temp));
    for (int i = 1; i <= 500; i++)
        for (int j = 1; j <= 500; j++)
            temp[i + j - 1] += base_number[i] * base_number[j];
    for (int i = 1; i <= 500; i++)
    {
        temp[i + 1] += temp[i] / 10;
        temp[i] %= 10;
    }
    memcpy(base_number, temp, sizeof(temp));
}
int main()
{
    //输入底数a与指数b
    int base;
    cin >> base >> index;
    for (int i = 1; base; i++)
    {
        base_number[i] = base % 10;
        base /= 10;
    }
    answer[1] = 1;
    while (index)
    {
        if (index % 2 == 1)
            operate1();
        index /= 2;
        operate2();
    }
    int flag = 0;
    for (int i = 1000; i >= 1; i--)
    {
        if (answer[i])flag = 1;
        if (flag)cout << answer[i];
    }
}