那么 有位神仙已经说了O（1）的算法（当然不是我）
这是一种骚操作 直接放代码了啊
inline LL mul(LL a,LL b,LL Mod){ LL lf = a * ( b >> 25LL ) % Mod * ( 1LL << 25 ) % Mod ; LL rg = a * ( b & ( ( 1LL << 25 ) - 1 ) ) % Mod ; return ( lf + rg ) % Mod ; }