不怕溢出の乘法取模【转】

$\Theta (1)$模取乘

分析

我们把两个 $10^{18}$ 的数乘起来会发生什么？会炸掉？会循环，$9223372036854775807+x$ 会变成 $-9223372036854775808+x-1$ （因为 $-9223372036854775808-1=9923372036854775807$ ）。

怎么规避负数然后取模呢？

这样：$-9223372036854775808+ab-1-(-9223372036854775808+\lfloor \frac{ab}{c}\rfloor \times c-1)=ab-\lfloor \frac{ab}{c}\rfloor \times c=ab\mathrm{m}\mathrm{o}\mathrm{d}c$ 。

如果 $longdouble$ 也救不了你的话就算了。

代码

typedef long long LL;
LL mul(LL a,LL b,LL c)
{
	return ((a*b-(LL)((long double)a/c*b+1e-8)*c)%c+c)%c; 
}

---

版权声明：本文为CSDN博主「C20182030Epic」的原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接及本声明。

---

备注

上面的说法好像并不完全正确？

因为 $ab$ 可能“爆很多圈”，也就是说如果 $ab=10^{36}$ 的话……？

所以作者认为，原理应该是：

• 第一点， 0 ≤ A − B < C ⇒ ( A   m o d   C ) − ( B   m o d   C ) ∈ { A − B + k C ∣ k ∈ { 0 , 1 } } 0\le A-B<C\Rightarrow(A\bmod C)-(B\bmod C)\in\{A-B+kC|k\in\{0,1\}\} 0≤A−B<C⇒(AmodC)−(BmodC)∈{A−B+kC∣k∈{0,1}} 。

很简单的道理：

• 设 $B=kC+m(0\le m<C)$
• 则 $A=kC+m+n(0\le n<C)$
• 故 A   m o d   C ∈ { m + n , m + n − C } A\bmod C\in\{m+n,m+n-C\} AmodC∈{m+n,m+n−C}
• 而 $B\mathrm{m}\mathrm{o}\mathrm{d}C=m$

所以得证。

• 第二点，“溢出”的本质是 $x$ 变为 $x\mathrm{m}\mathrm{o}\mathrm{d}{2}^{64}$ 。

• 第三点，无论如何，$ab\mathrm{m}\mathrm{o}\mathrm{d}c=ab-c\lfloor \frac{ab}{c}\rfloor$ 。

将第三点中的 $ab$ 当做 $A$ ，$c\lfloor \frac{ab}{c}\rfloor$ 当做 $B$ ，$2^{64}$ 当做 $C$ ，代入第一点，即可证明。