R 语言 | 计算可达矩阵
R 语言 | 计算可达矩阵
努力在北京混出人样
发布于 2019-02-18 16:50:20
发布于 2019-02-18 16:50:20
代码可运行
运行总次数：0
代码可运行
在微博有位朋友问我可达矩阵的计算，于是发了点时间用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的使用。
本文参与 腾讯云自媒体同步曝光计划，分享自作者个人站点/博客。
原始发表：2017年04月30日，如有侵权请联系 cloudcommunity@tencent.com 删除
评论
登录后参与评论
推荐阅读
腾讯云开发者
