既然标题为偷懒。。。
那么肯定不用一大堆神奇的定理啦
如果令
a
表示被逆元的数,
a×a′≡1modQ
令 d 满足
假如我们一开始知道 d 的逆元
d×d′≡1modQ
(Q−a⌊Qa⌋)×d′≡1modQ
−a⌊Qa⌋d′≡1modQ
(Q−⌊Qa⌋)d′×a≡1modQ
于是 a 的逆元就是
于是程序就很简单啦。
int inverse(int a) {
return a == 1 ? 1 : (Q - Q / a) * inverse(Q % a) % Q;
}
由于递归的特性使问题规模不断减小,因此还可以用来预处理1~N的逆元。
inv[1] = 1;
for (int a = 2; a < Q; ++a) inv[i] = (Q - Q / a) * inv[Q % a] % Q;
与Q是否为质数应该无关。
所以懒得打extgcd或phi+qpow的。。。。
实际上就长一点。。
不会证复杂度。。
求教。。
猜
O(lgn)
,逃
感觉Q和a特别搞一个奇怪的数字就可以卡了。。