R 语言 | 计算可达矩阵

在微博有位朋友问我可达矩阵的计算，于是发了点时间用R语言写出来了。

问题如下：

问题

计算过程：

计算过程

注意：是矩阵的乘法。

代码如下：

A <- matrix(c(0,0,0,0,0,0,0,1,0,0,1,1,0,0,0,0,0,0,0,0,0,0,1,1,0),nrow=5)
A1 <- A+diag(5)

repeat{
  B <- A1 * A1;
  if(identical(B,A1)) {print(B);break} else A1 <- B
}

结果：

     [,1] [,2] [,3] [,4] [,5]
[1,]    1    0    1    0    0
[2,]    0    1    1    0    0
[3,]    0    1    1    0    1
[4,]    0    0    0    1    1
[5,]    0    0    0    0    1

重点在repeat的使用。