链接1
链接2
通过消元，转化为与之互为相似矩阵的上海森堡矩阵，完成内求解矩阵的特征多项式。
由此，可以优化向量×矩阵的p次幂的问题，估计时间复杂度为，没问题。