the m1 indicates that m is followed by 1
p_m0 <- matrix(0, dd, 2^K)
p_m1 <- matrix(0, dd, 2^K)