基于另一个矩阵修改矩阵值而不使用循环或应用函数

Question

我有两个矩阵S 和A，我想根据S 的值修改A 的值。矩阵S 包含我想修改的矩阵A 的索引。比如

S <- matrix(c(2,3,3,1,1,2),nrow=3, byrow=TRUE)

   [,1] [,2]
[1,]    2    3
[2,]    3    1
[3,]    1    2

生成的矩阵A 将是

     [,1] [,2] [,3]
[1,]    0    1    1
[2,]    1    0    1
[3,]    1    1    0

换句话说，在矩阵S中，每一行代表一个图顶点。每行中的值（比如ith）是连接到ith 顶点的顶点索引。矩阵A是图的邻接矩阵。

问题：如何解决问题，即获取矩阵A，而不使用for循环或apply/sapply/mapply等函数？最快的方法是什么？

编辑：一个例子

S <- matrix(c(2,3,3,1,1,2,1,2),nrow=4, byrow=TRUE)

    [,1] [,2]
[1,]    2    3
[2,]    3    1
[3,]    1    2
[4,]    1    2

结果应该是：

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

Answer 1

我们可以使用'S'作为行/列索引来改变A中的值

A <- matrix(0, 4, 4)
A[cbind(c(row(S)), c(S))] <- 1
+(A|t(A))
#     [,1] [,2] [,3] [,4]
#[1,]    0    1    1    1
#[2,]    1    0    1    1
#[3,]    1    1    0    0
#[4,]    1    1    0    0