要证明它的正确性,我们只需要稍微了解 long double 的精度以及 long long 值域超出后的行为。
支持数值大小:1e18*1e18
#define int long long
inline int Mul(int x, int y, int p) {
int tmp = x * y - (int)((long double)x * y / p) * p;
tmp = (tmp < 0 ? tmp + p : tmp);
tmp = (tmp >= p ? tmp - p : tmp);
return tmp;
}