Poj.Grids 2951 浮点数求高精度幂

2951:浮点数求高精度幂

总时间限制:

1000ms

内存限制:

65536kB

描述

有一个实数 R ( 0.0 < R < 99.999 ) ,要求写程序精确计算 R 的 n 次方。n 是整数并且 0 < n <= 25。

输入

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

输出

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

样例输入

95.123 12
0.4321 20
5.1234 15
6.7592  9
98.999 10
1.0100 12

样例输出

548815620517731830194541.899025343415715973535967221869852721
.00000005148554641076956121994511276767154838481760200726351203835429763013462401
43992025569.928573701266488041146654993318703707511666295476720493953024
29448126.764121021618164430206909037173276672
90429072743629540498.107596019456651774561044010001
1.126825030131969720661201

 //@author langx
 //G++4.4
 #include <cstdio>
 #include <cstring>
 #define SIZE 126
 #define LEN 5

 int a[SIZE],b[LEN];

 int mul( int alen,int blen){
     int tmp[SIZE], i, j, k;
     memset(tmp,,sizeof(tmp));
     for( i = ; i < alen; ++i){
         for( j = ; j < blen; ++j){
             tmp[i + j] += a[i] * b[j];
         }
     }
     k = alen + blen - ;
     for( i = ;i < k; ++i){
         if( tmp[i] >= ){
             tmp[ i + ] += tmp[i] / ;
             tmp[i] = tmp[i] % ;
         }
         a[i] = tmp[i];
     }
     a[i] = tmp[i];
     if( tmp[k] != )k++;
     return k;
 }

 int main(){
     char r[LEN + ];
     int n, i, j, len, alen, mark, formark, ok;
     while( scanf("%s%d",r, &n) ==  ){
         len = strlen(r);
         mark = ;
         for( i = len - , j =  ; i >=  ; --i){
             if(r[i]=='.'){
                 mark =i;continue; //记录小数点位置
             }
             a[j] = b[j++] = r[i]-'';
         }
         --len;
         mark = len - mark;
         alen = len;
         for( i = ;i < n; i++){
             alen = mul( alen,len);
         }
         for( i = ; i < alen; ++i){
             if( a[i] != )break;
         }
         mark *= n;
         formark = i; //记录后缀，如70.10000
         ok = ;
         for( i = alen-; i >= formark; --i){
             if( a[i] ==  && i >= mark && !ok )continue;
             if( i == mark - ) printf(".");
             printf("%d", a[i]);
             ok = ;
         }
         for(formark--; formark >= mark; --formark){
             printf("");
         }
         printf("\n");
     }
     return ;
 }