poj 1001 求高精度幂

#include <iostream>
#include <string.h>
using namespace std;
int n,a[1000000],b[1000000],c[1000000];
int dot,l,r,lena,lenb,lens;
char s[10];

void mul(){ //高精度乘法
int i,j;
memset(c,0,sizeof(c));
for(j=0;j<lenb;j++){
for(i=0;i<lena;i++){
c[i+j]+=a[i]*b[j];
c[i+j+1]+=c[i+j]/10; ///
c[i+j]%=10;
}
}
lena+=lenb;
for(i=lena-1;i>=0;i--){
if(c[i]){
lena=i+1;
for(;i>=0;i--){a[i]=c[i];}
break;
}
}
}

int main(){
int i,j;
while(scanf("%s %d",s,&n)!=EOF){
dot=-1;
lens=lena=lenb=strlen(s);
for(i=lens-1,j=0;i>=0;i--){ //string->int数组
if(s[i]!='.'){
a[j]=b[j]=s[i]-'0';
j++;
}
else dot=j;
}
if(dot!=-1){ //整数
lena=lenb=lens-1;
dot*=n; //小数部分长度
}
for(i=0;i<n-1;i++) mul(); //n次相乘

if(dot==-1){
for(i=lena-1;i>=0;i--) printf("%d",a[i]);
}
else{
l=0;
for(i=0;i<lena;i++){
if(a[i]!=0){
r=i;
break;
}
}
if(dot>=lena){
printf(".");
for(i=dot-lena;i>0;i--)
printf("0");
for(i=lena-1;i>=r;i--) printf("%d",a[i]);
}
else{
for(i=lena-1;i>dot-1;i--) printf("%d",a[i]);
if(i>=r)printf(".");
for(;i>=r;i--) printf("%d",a[i]);
}
}
putchar(10);
}
return 0;
}