r - 没有循环的增强矩阵

Question

我没有找到解决我的问题的函数。我有矩阵 A。

 A<-matrix(c(1,0,0,1,1,0,0,1,1),3,3)

我想要计算 A^13，不是元素与元素，而是行与列。A^2 的例子

A%*%A

我可以用循环来做，但我的问题是:有没有简单的方法或函数来计算这个？

我的 A^13 结果必须是:

[1 13 78]
[0  1 13]
[0  0  1]  

Answer 1

也许您正在寻找Reduce:

Reduce("%*%", replicate(13, A, FALSE))
#      [,1] [,2] [,3]
# [1,]    1   13   78
# [2,]    0    1   13
# [3,]    0    0    1