https://123dok.net/document/yr31no4v-fast-algorithm-computing-characteristic-polynomial-p-curvature.html