快速乘法取模算法

原理：

32+16+4=52

 

LL qmul(LL x, LL y, LL mod) { // 乘法防止溢出， 如果p * p不爆LL的話可以直接乘； O(1)乘法或者轉化成二進制加法
//快速乘法取模算法

    LL ret = 0;
    while(y) {
        if(y & 1)
            ret = (ret + x) % mod;
        x = x * 2 % mod;
        y >>= 1;
    }
    return ret;
}