R语言-快速提取所有的组合

以 1， 2，3为例，要提取他们的所有组合的可能性，那么一共有

$$C^1_3+C^2_3 + C^3_3$$ 种情况，计算结果为：

considerNum <- sapply(1:3, function(x) choose(3,x))
#  3 3 1
sum(considerNum)
# 7

为快速提取他们所有组合，用矩阵的思考方法，利用Bit矩阵进行提取。原理以1，2，3为例子。需用矩阵方法进行提取1，2，3的所有组合。那么需生成一个这样的矩阵：

1	2	3
0	0	1
0	1	0
0	1	1
1	0	0
1	0	1
1	1	0
1	1	1

它其实是一个$1：（2^{n}-1）$的一个二进制表达

BitMatrix <- function(n){
  # 作用：返还Bit矩阵
  # Args:n：数据长度

  set <- 0:(2^n-1)
  rst <- matrix(0,ncol = n,nrow = 2^n)
  for (i in 1:n){
    rst[, i] = ifelse((set-rowSums(rst*rep(c(2^((n-1):0)), each=2^n)))/(2^(n-i))>=1, 1, 0)
  }
  rst
}

k <- c(1,2,3)
selectMatrix <- BitMatrix(length(k))
apply(selectMatrix, 1, function(x){ k[which(x==1)] })
#[[1]]
#[1] 3
#[[2]]
#[1] 2
#[[3]]
#[1] 2 3
#[[4]]
#[1] 1
#[[5]]
#[1] 1 3
#[[6]]
#[1] 1 2
#[[7]]
#[1] 1 2 3