R 语言 | 计算可达矩阵

发布于

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

问题如下：

问题

计算过程：

计算过程

注意：是矩阵的乘法。

代码如下：

代码语言：javascript
复制
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

}

结果：

代码语言：javascript
复制
     [,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的使用。

#repeat