先日の計算機代数勉強会からだいぶ経ってしまったが、整数係数多項式の因数分解をする Berlekamp-Zassenhaus アルゴリズムをようやく実装できた!

有限体上での因数分解をするBerlekampアルゴリズムは実装できてたのだけれど、pを法とする因子からp^kを法とする因子に持ち上げる Hensel lifting と、p^k を法とする因子から最終的に整数上での因子を求める部分の理解にだいぶ苦戦してしまっていた。

Kronecker's method で数分とかかかっていた処理が一瞬で処理できるようになって、ちょっと感動。