Modular exponentiation

# compute B^E mod M
def prodmodul(B, E, M):
    buf=1
    while E > 0:
        buf=(buf * (( B**(E%10) ) % M) ) % M
        B=(B**10) % M
        E//=10L
    return buf