【poj 1001】高精密度幂运算

思路就是将每位数字分别存入一个整数数组中，并记录小数点的位置，循环相乘，按列相加因子，并进位去权，和我们草稿纸上算乘法差不多，循环n次

#include <fstream>
#include <cstring>
using namespace std;

int main()
{

    int g = 0 , h , i , j , k , p , q , n , B[6] , BW , P , C[500] , CW , D[500] , SSS[500][500] , WWW[500] , PPP[500] , NUM , ZERO ;
    int table[10][10] ={0 ,  0 ,  0 ,  0 ,  0 ,  0 ,  0 ,  0 ,  0 ,  0 ,
                        0 ,  1 ,  2 ,  3 ,  4 ,  5 ,  6 ,  7 ,  8 ,  9 ,
                        0 ,  2 ,  4 ,  6 ,  8 , 10 , 12 , 14 , 16 , 18 ,
                        0 ,  3 ,  6 ,  9 , 12 , 15 , 18 , 21 , 24 , 27 ,
                        0 ,  4 ,  8 , 12 , 16 , 20 , 24 , 28 , 32 , 36 ,
                        0 ,  5 , 10 , 15 , 20 , 25 , 30 , 35 , 40 , 45 ,
                        0 ,  6 , 12 , 18 , 24 , 30 , 36 , 42 , 48 , 54 ,
                        0 ,  7 , 14 , 21 , 28 , 35 , 42 , 49 , 56 , 63 ,
                        0 ,  8 , 16 , 24 , 32 , 40 , 48 , 56 , 64 , 72 ,
                        0 ,  9 , 18 , 27 , 36 , 45 , 54 , 63 , 72 , 81 }; //九九乘法表
    char buf[7] ;

    while ( scanf( "%s%d" , buf , &n ) == 2 )
    {
        memset( B , 0 , sizeof(B) );                    //初始化数组
        memset( C , 0 , sizeof(C) );
        C[1] = 1 ; CW = 1 ;
        memset( D , 0 , sizeof(D) );


        i = BW = P = 0 ; j= 5 ;
        while ( buf[j] == '0' ) j-- ;       //去掉小数点后无效的0
        while ( j >= 0 )
        {
            if ( buf[j] == '.' ) { P = i ; j-- ;}       //记录小数点的位置p，跳过小数点
            B[++i] = buf[j--] - '0' ;                   //倒序的方式将读入的字符串中的数字逐个存入数组B，B[1]个位、B[2]十位，以此类推
        }
        BW = i ;
        while ( B[BW] == 0 ) BW-- ;     //得到的整数B有BW位（去掉前导的0）

        for ( h = 1 ; h <= n ; h++ )
        {
            for ( i = 1 ; i <= BW ; i++ )
            {
                k = i ;
                for ( j = 1 ; j <= CW ; j++ )
                {
                    D[k++] += table[B[i]][C[j]] ;
                }
            }
            CW += BW ;                      //位的增加
            C[1] = 0 ;
            for ( i = 1 ; i < CW ; i++ )
            {
                C[i] += D[i] ;              //累加
                D[i] = 0 ;
                C[i+1] = C[i] / 10 ;        //进位
                C[i] %= 10 ;                //去权
            }
            if ( C[CW] == 0 ) CW-- ;        //去首位0
        }


        P *= n ;    //小数点的移动
        g++ ;

        for ( i = 1 ; i <= CW ; i++ ) SSS[g][i] = C[i] ;
        WWW[g] = CW ; PPP[g] = P ;

    }


    for ( h = 1 ; h <= g ; h++ )
    {
        j = 1 ;
        ZERO = 0 ;
        while ( SSS[h][j++] == 0 ) ZERO++;                   //定位小数点后无效的0
        for ( k = WWW[h] ; k > PPP[h] ; k-- ) printf( "%d" , SSS[h][k] );
        if ( PPP[h] > ZERO )
        {
            printf( "." ) ;
            for ( k = PPP[h] ; k > ZERO ; k-- ) printf( "%d" , SSS[h][k] );
        }
        printf( "\n" ) ;
    }

}