HDU 2035 人见人爱A^B（快速矩阵幂）

#include<bits/stdc++.h>
#define M 1000
using namespace std;
int del(int a,int b,int c)
{
    int ans=1;
    a=a%c;
    while(b>0)
    {
        if(b%2==1) ans=(ans*a)%c;
        b=b/2;
        a=(a*a)%c;
    }
    return ans;
}
int main()
{
    int a,b;
    while(cin>>a>>b)
    {
        if(a==0&&b==0)
            break;
        cout<<del(a,b,M)<<endl;
    }
    return 0;
}