What you need is a multiset permutation.
> library(iterpc)
> I <- iterpc(c(16, 14), labels=c(0,1), ordered=TRUE)
> getlength(I)
[1] 145422675
The matrix is 145422675 by 30, too large to be stored in a single matrix. Use getnext
to get the next 10000000 permutations.
> getnext(I, 10000000)
You could repeat the above 15 times in order to get all the permutations.
与恶龙缠斗过久,自身亦成为恶龙;凝视深渊过久,深渊将回以凝视…