POJ 1001 高精度实数相乘与快速幂

题目链接：http://poj.org/problem?id=1001

题意：给一个100以内的六位有效数字的正实数x，再给一个指数n，求x的n次幂。

分析：先把小数点去掉，就是大数相乘加上快速幂，要注意的就是输出的时候加小数点，和去前导0。答案的小数位数就是x的小数位数乘以n。大数相乘有点不熟悉了，用这个来复习一下，注意点就是先char型数组读入，再用int型数组倒着存数据，这样后面的计算比较方便。快速幂是我第一次敲，感觉好神奇。

代码：

#include<cstdio>
#include<stack>
#include<cstring>
#include<iostream>
#include<cmath>
using namespace std;
const int N=150;
char num[N];

struct bignum//用结构体存大数，这样调用函数传参时比较方便
{
    int val[N*2];
    bignum()//初始化
    {
        memset(val,0,sizeof(val));
    }
    bignum& operator = (bignum a)//重载等号，后来发现没什么用，C++没学清楚。。。
    {
        for(int i=0;i<N;i++)
            val[i]=a.val[i];
        return *this;//要有返回值，codeblocks这里没返回值居然不会报错，坑~
    }
};

bignum mul(bignum a,bignum b)
{
    bignum temp;
    int jin=0,i,j;
    for(i=0;i<N;i++)//全部都乘以遍，懒得优化了，感觉也优化不了多少
    {
        for(j=0;j<N;j++)
        {
            int temp0=temp.val[i+j]+a.val[i]*b.val[j]+jin;//大数相乘的精髓，这个地方自己坑了好久
            temp.val[i+j]=temp0%10;
            jin=temp0/10;
        }
    }
    return temp;
}

bignum fastmul(bignum a,int times)
{
    bignum temp;
    temp.val[0]=1;
    while(times)
    {
        if(times&1) temp=mul(a,temp);//快速幂的精髓。。。
        a=mul(a,a);
        times/=2;
    }
    return temp;
}

int main ()
{
    //freopen("in.txt","r",stdin);
    while(cin>>num)
    {

        int i=0,j=0,n=N-1,len_point=0;
        bool point=true,first=true;
        bignum num2;
        for(i=5;i>=0;i--)
        {
            if(num[i]=='0'&&first)//处理后导0
                continue;
            if(num[i]!='0')
                first=false;
            if(num[i]=='.')
                point=false;
            else
                num2.val[j++]=num[i]-'0';//倒着存
            if(point)//记录小数位数
                len_point++;
        }
        int time;
        cin>>time;


        bignum ans;
        ans=fastmul(num2,time);

        for(i=N-1;i!=-1&&i>len_point*time-1;i--)//处理小数点前的前导0,
            if(ans.val[i])
                break;
        if(i==-1)
            printf("0");
        else
        {
            for(;i>=0;i--)
            {
                if(i==len_point*time-1)
                    cout<<'.';
                cout<<ans.val[i];
            }
            cout<<endl;
        }
    }
    return 0;
}