利用了快速幂方法。

#include <iostream>
#include <string>
#include <cstdio>

using namespace std;

int fun(int a, int b, int mod){
    int ans = 1;
    while(b != 0){
        if(b % 2 != 0){
            ans *= a;
            ans %= mod;
        }
        b /= 2;
         a *= a;
         a %= mod;
    }
    return ans;
}

int main(){
    int a, b;
    int ans; 
    while (cin >> a >> b) {
        if(a == 0 || b == 0) return 0;
        ans = fun(a, b, 1000);
        cout << ans << endl;
    }
    return 0; 
}