快速乘法取模算法

原理：

32+16+4=52

 1 LL qmul(LL x, LL y, LL mod) { // 乘法防止溢出， 如果p * p不爆LL的话可以直接乘； O(1)乘法或者转化成二进制加法
 2 //快速乘法取模算法
 3
 4     LL ret = 0;
 5     while(y) {
 6         if(y & 1)
 7             ret = (ret + x) % mod;
 8         x = x * 2 % mod;
 9         y >>= 1;
10     }
11     return ret;
12 }

原文地址：https://www.cnblogs.com/Fy1999/p/8908522.html