hdu 1063 Exponentiation（高精度幂）

poj这题数据很水很容易过，然而hdu的这题可能是因为case太多O(125*125*n)的复杂度TLE，最后用了快速幂的方法优化到O（125*125*logn）过了。

140MS。

/*
Author:Zhaofa Fang
Lang:C++
*/
#include <cstdio>
#include <cstdlib>
#include <sstream>
#include <iostream>
#include <cmath>
#include <cstring>
#include <algorithm>
#include <string>
#include <utility>
#include <vector>
#include <queue>
#include <stack>
#include <map>
#include <set>
using namespace std;

typedef long long ll;
#define DEBUG(x) cout<< #x << ':' << x << endl
#define REP(i,n) for(int i=0;i < (n);i++)
#define FOR(i,s,t) for(int i = (s);i <= (t);i++)
#define PII pair<int,int>
#define PB push_back
#define MP make_pair
#define ft first
#define sd second
#define lowbit(x) (x&(-x))
#define INF (1<<30)



void product(int a[],int b[],int c[])
{
    REP(j,125)c[j]=0;
    int e = 0;
    int N = 125-1;
    while(N>=0 && !b[N] )N--;
    REP(k,N+1)
    {
        for(int j=0;j+k<125;j++)
        {
            int tmp = a[j] * b[k];
            int tmp1 = tmp + c[j+k] + e;
            c[j+k] = (tmp1)%10;
            e = (tmp1)/10;
        }
    }
    REP(k,125)b[k] = c[k];
}
void print(int b[],int num)
{
    int k = 125 - 1,j=0;
    int end = 0;
    while(k>=0 && !b[k] )k--;
    while(end<num && !b[end])end++;
    k = max(num-1,k);
    //DEBUG(k);DEBUG(num);DEBUG(end);
    //puts("**********************************");
    for(int i=k;i>=end;i--)
    {
        j++;
        if(j == k-num+2)printf(".");
        printf("%d",b[i]);
    }
    if(end>k)printf("0");
    puts("");
}
void POW(int a[],int n,int num)
{
    int b[125],c[125];
    REP(i,125)b[i] = a[i];
    while(n)
    {
        if(n & 1)
        {
            product(a,b,c);
            REP(i,125)b[i] = c[i];
        }
        n >>= 1;
        product(a,a,c);
        REP(i,125)a[i] = c[i];
    }
    print(b,num);
}

int main()
{
    //freopen("in","r",stdin);
    char str[102];
    int a[125];
    int n;
    while(scanf("%s%d",str,&n)!=EOF)
    {
        int len = strlen(str);
        int num = 0;
        memset(a,0,sizeof(a));
        for(int i=len-1,j=0;i>=0;j++,i--)
        {
            if(str[i] == '.')
            {
                num = j;j--;
            }
            else
            {
                a[j] = str[i] - '0';
            }
        }
        num *= n;
        POW(a,n-1,num);

    }
    return 0;
}