poj 1001高精浮点数乘法

http://poj.org/problem?id=1001

据说在10天之内做完poj前十题的就可以迈向神之路……

作为打头阵的一道题自然是不会防水

本题乍眼一看是高精度浮点数次幂的计算，如果直接计算自然是比较难受，所以我们的数学基础在这里又体现了它非常重要的作用

对于两个浮点数相乘的计算我们会记录这两个数的小数点后分别有几个数，之后将这两部分相加作为最终答案小数点的位置（即最终答案小数点后数字的数量为计算前两个数小数点后数字数量之和）

那么对于一个浮点数次幂的计算我们也很容易就得到最终答案小数点后的数字的数量为最初的浮点数中小数点后数字的数量乘上要进行计算的次幂

这样一来，在进行高精度计算的时候便可以忽略小数点的存在，只是在最后输出的时候注意一下在某个位置输出这个小数点即可。

另外此题还有几个比较恶心的情况:

1.输入的时候可能会出现类似1.000000的情况，要注意不要计算上后面无效的0

2.输出的时候要求去掉前导0，即输出.001而非0.001

3.本题依然有不存在小数点的情况，处理的时候要注意

这道题思想上没什么难度，但很考验代码功底（据说java做这种题特别方便……c++已哭晕在厕所…………

---

#include<iostream>
#include<cstdio>
#include<queue>
#include<cstring>
#include<algorithm>
using namespace std;
char in[2333];
struct bign {
    int len;
    int num[2333];
};
int cnt;
int q;
bool flag;
bool flag_wa;
bign mul(bign a,bign b)
{
    bign c;
    memset(c.num,0,sizeof(c.num));
    int len;
    for(int i=0;i<a.len;i++)
    {
        for(int j=0;j<b.len;j++)
        {
            c.num[i+j]+=a.num[i]*b.num[j];
            if(c.num[i+j]>9)
            {
                c.num[i+j+1]+=c.num[i+j]/10;
                c.num[i+j]%=10;
            }
        }
    } 
    len=a.len+b.len-1;
    while(!c.num[len]&&len>1)
    len--;
    while(c.num[len])
    len++;
    c.len=len;
    return c;
}
bign ksm(bign a,int q)
{
    bign c;
    c.len=1;
    c.num[0]=1;
    while(q)
    {
        if(q&1)
        {
            c=mul(c,a);
        }
        a=mul(a,a);
        q>>=1;
    }
    return c;
} 
void init(bign &a,char *in)
{
    flag=0;
    flag_wa=0;
    int len=strlen(in)-1;
    for(int i=0;i<=len;i++)
    {
        if(in[i]=='.')
        {
            flag=1;
            for(int j=i+1;j<=len;j++)
            {
                if(in[j]>='1')
                flag_wa=1;
            }
        }

    }
    if(flag)
    a.len=len;
    else
    a.len=len+1;
    int i=0;
    while(len>=0)
    {
        if(in[len]=='.')
        {
            cnt=i;
            //cout<<cnt<<endl;
            len--;
            continue;
        }
        a.num[i]=in[len]-'0';
        i++;
        len--;
    }
}
void print(bign b)
{
    int dleta=b.len-cnt*q;
    if(dleta<0)
    {
        cout<<".";
        for(int i=1;i<=-dleta;i++)
        cout<<"0";
    }
    //cout<<endl<<dleta<<endl;
    for(int i=b.len-1;i>=0;i--)
    {
        cout<<b.num[i];
        if(i==cnt*q&&flag&&flag_wa)
        cout<<".";
    }
    //cout<<endl;
    //cout<<cnt;
}
int main()
{
    while(scanf("%s%d",in,&q)!=EOF)
    {
        cnt=0;
        int y=strlen(in);
        for(int i=y-1;i>=0;i--)
        {
            if(in[i]>='1'||in[i]=='.')
            break;
            else
            in[i]-='0';
        }
        bign a;
        init(a,in);
        bign b=ksm(a,q);
        print(b);
        cout<<endl;
    }
    return 0;
}