光速乘

直接上代码

点击显/隐代码
1
2
3
4
5
6
7
8
typedef long long ll;
typedef long double lb;
const lb dx=1e-8;
inline ll multi(ll x,ll y,ll mod)
{
ll tmp=(x*y-(ll)((lb)x/mod*y+dx)*mod);
return tmp<0?tmp+mod:tmp;
}