快速幂取模

//a^n的结果对MOD取余
const int MOD=10003;
int PowMOD(int a,int n){
    int ret=1;
    a%=MOD;
    while(n){
        if(n&1){
            ret =(ret*a)%MOD;
        }
        n>>=1;
        a=(a*a)%MOD;
    }
    return ret;
}