被数论怒虐了一天 心力憔悴啊 感觉脑细胞已经快消耗殆尽了>_< 但是今天还是会了很多之前觉得特别神的东西

比如BSGS 之前听了两遍 好像都因为听得睡着了没听懂-。- 今天终于硬着头皮学会了~ 做个总结吧 免得又忘记- –

BSGS：

BSGS就是求 A^x=B(mod C) 0<=x<C的解(C为素数)

做一个转换 设m*i+j=x (m=trunc(sqrt(C)))

将A^i(0<=i<m) 存入hash表中(i,A^i)

这样我们就能O(1)求出A^x=B 对应的x的值

设AA=(A^m^i)mod C 枚举i的值(0<=i<m)

原方程就变为 AA*(A^j)=B(mod C)

因为C为素数 所以(AA,C)=1

故 可以用拓展欧几里德求出A^j值为多少 再在hash中找A^j对应的j 答案get

其实不难发现这是个神奇的分块思想 将n分为√n块

预处理使得计算每块的复杂度降为O(1) 到达加速

拓展BSGS：

上题中有一个限制条件 C为素数 如果没有这个条件要怎么做呢

以下是AC大神的做法 由AK大神跟我讲解（orz两个神犇）

A^x=B(mod C)

先将这个方程转换为 A^x+Cy=B

设 A’=A/(A,C),C’=C/((A,C)^z) (C’与(A,C)^z互质)

将方程左右都除以(A,C)^z （如果B不能整除(A,C)^z 则无解）

会得到：

(A’^z)*A^(x-z)+C’y=B/((A,C)^z)

这时 将A’^z视为系数 A与C’互质就能用上面的方法求解了