我有一个 7000x7000 的稀疏矩阵（scipy），我想对其取幂。我尝试过使用 scipy.sparse.linalg.expm，它适用于较小的矩阵（1000x1000 矩阵需要几秒钟），但计算相关矩阵需要很长时间。有没有办法优化这个或解决这个算法的时间复杂度（显然是 O(n^2)）？对于背景，我已经尝试了与 R 包 {expm} 相同的结果，因此它显然更多地是关于问题的算法方法而不是使用的特定软件。你有什么想法？