求高精度幂 POJ1001

最新推荐文章于 2023-08-29 17:03:16 发布

原创最新推荐文章于 2023-08-29 17:03:16 发布 · 462 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#algorithm #算法 #数据结构 #c++

本文介绍了如何解决POJ1001问题，即实现高精度幂运算。通过讲解算法思路和数据结构的应用，详细阐述了如何处理大整数的乘方计算，同时给出了输入输出格式及样例分析。

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

求高精度幂

Description

对数值很大、精度很高的数进行高精度计算是一类十分常见的问题。比如，对国债进行计算就是属于这类问题。

现在要你解决的问题是：对一个实数R( 0.0 < R < 99.999 )，要求写程序精确计算 R 的 n 次方(Rⁿ)，其中n 是整数并且 0 < n <= 25。

Input

T输入包括多组 R 和 n。 R 的值占第 1 到第 6 列，n 的值占第 8 和第 9 列。

Output

对于每组输入，要求输出一行，该行包含精确的 R 的 n 次方。输出需要去掉前导的 0 后不要的 0 。如果输出是整数，不要输出小数点。

Sample Input

Sample Output

548815620517731830194541.899025343415715973535967221869852721
.00000005148554641076956121994511276767154838481760200726351203835429763013462401
43992025569.928573701266488041146654993318703707511666295476720493953024
29448126.764121021618164430206909037173276672
90429072743629540498.107596019456651774561044010001
1.126825030131969720661201

水题，大数乘法，开两个大数组，分别存储中间值和最终结果。需注意：输出需要去掉前导的0后不要的0（整数部分为0也不保留）。如果输出是整数，不要输出小数点。

#include <stdio.h>
#include<iostream>
#include<cstring>
using namespace std;

const int LEN=5;
const int MAX_POWER=25;
int result[200];
int r[LEN+1];
int n;
int main(){
    char str[LEN+2];
    while(scanf("%s %d",str,&n)!=EOF){
        str[LEN+2]=0;
        int flag=LEN+1;

        if(n>MAX_POWER||n<=0)
            return 0;

        int j=LEN;
        for(int i=LEN;i>=0;i--){
            if(str[i]!='.')
            {
                r[j]=str[i]-'0';
                result[j]=r[j];
                j--;
            }
            else
                flag=LEN-i;
        }

        int len_result=LEN;
        int len_r=LEN;
        for(int k=1;k<=n-1;k++){
            int tmp[200]={0};
            for(int i=len_result;i>=1;i--){
                for(int j=len_r;j>=1;j--){
                    tmp[i+j]+=r[j]*result[i];
                    if(tmp[i+j]>9){
                        tmp[i+j-1]+=tmp[i+j]/10;
                        tmp[i+j]%=10;
                    }
                }
            }
            len_result+=len_r;
            for(int i=1;i<=len_result;i++)
                result[i]=tmp[i];
        }
        flag*=n;
        int flag1=0;
        for(int i=0;i<=len_result-flag;i++)
        {
            if(result[i]==0&&flag1==0)
                continue;
            else{
                flag1=1;
                cout<<result[i];
            }
        }
        int flag2=0;
        for(int j=len_result-flag+1;j<=len_result;j++){
            if(result[j]!=0)
                flag2=1;
        }
        if(flag2==1){
            int flag3=0;
            for(int j=len_result;j>=len_result-flag+1;j--){
                if(result[j]!=0){
                    flag3=j;
                    break;
                }
            }
            cout<<'.';
            for(int j=len_result-flag+1;j<=flag3;j++)
                cout<<result[j];
        }
        cout<<endl;
    }
    return 0;
}