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的使用。