AM_prior_K_NegBin function

computes the prior number of clusters