在微博有位朋友问我可达矩阵的计算，于是发了点时间用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的使用。