O(1)64位快速乘

原理十分简单，即$a\ast b\mathrm{m}\mathrm{o}\mathrm{d}p=a\ast b-\lfloor a\ast b/p\rfloor \ast p$。代码如下：

typedef long long ll;
inline ll mul(ll a, ll b, ll p) {
    a %= p, b %= p;
    ll c = (long double) a * b / p;
    cout << c << endl;
    ll ans = a * b - c * p;
    cout << ans << endl;
    if (ans < 0) ans += p;
    else if (ans >= p) ans -= p;
    return ans;
}

当$a,b<p$时，$\lfloor a\ast b/p\rfloor$肯定也小于$p$，我们可以先对$a,b$取模，使其满足这一条件。其实我不太明白(long double) a * b / p的具体表现，浮点数在上溢时会以类似科学计数法的形式舍弃低位，但为什么通过该运算不会导致结果的误差，希望能有大佬帮我解释一波OTZ