calculation of the normalising constant in the case of the initial distribution
somainit = 0
for(m in 1:2^K){
somainit = somainit + (alpha_init_k[m] + count_init_k[m] - 1)}