注意到对于矩阵，通常意义上的乘法分配律成立，也即 A(B+C)=AB+AC 那么就有

（A+…+A^n)

=A(E+…+A^(n-1))

=A*(A^n-E)*(A-E)^(-1)

这样不会快很多吗