f <- function(n){
  if(n==1) return(matrix(1))
  else if(n==2) return(matrix(c(1,4,2,3),nrow=2))
  else{
    a <- (4*n-4):(3*n-1):-1
    b <- (n+1):(2*n-2)
    c <- 1:n
    d <- (3*n-2):(2*n-1):-1
    return(rbind(c,cbind(a,f(n-2)+4*n-4,b),d))
  }
}

上楼是代码

你运行下
1：5：-1
就知道原因了，R中没有这样的表达

顺便问一下怎么把这个矩阵的标签去掉呢？ 我想了很久names(x) <- NULL在上面递归函数里面好像加不上去