ll fast_mult(ll a,ll b,ll mod) {
    return (a*b - (ll)((long double)a/mod*b)*mod+mod)%mod;
}

参考：https://blog.csdn.net/m0_38013346/article/details/81435369

注意事项：

  1. double可能会挂，最好long double。
  2. u,v>=p可能会挂，必要时先%p。
  3. 用浮点数算出u*v/p的值时事实上允许了±1的误差，因此可能出现负数，所以必须+p再%p。因此理论上不需要+eps。

还有一种模板：

inline ll fast_mult(ll x,ll y,ll mod) {
	ll tmp=(x*y-(ll)((long double)x/mod*y+1.0e-8)*mod);
	return tmp<0 ? tmp+mod : tmp;
}

另： 

ll n,k,p;
ll kre(ll a,ll b){
    a=a%p;
    b=b%p;
    return ((a*b-(ll)(((double)a*b+0.5)/p)*p)+p)%p;
}
ll kru(ll a,ll b){
    ll ret=1;
    while(b){
        if(b&1){
            ret=kre(ret,a)%p;
            b--;
        }
        b/=2;
        a=kre(a,a)%p;
    }
    return ret;
}

参考：https://blog.csdn.net/axuhongbo/article/details/82011464