很久以前做题的时候看到的。
inline long long multi(long long x,long long y,long long mod)
{
long long tmp=(x*y-(long long)((long double)x/mod*y+1.0e-8)*mod);
return tmp<0 ? tmp+mod : tmp;
}
利用了long double 在某些系统下大于64bit的特性。因此如果评测机的long
double是64bit时会WA。
由于浮点数的精度限制,mod不能太小,必须在1e18的级别。实测mod =
1e9+7时,x = 1e18,y = 1e18时会错。