实现gibbs 等离散分布的sampling的code

实现gibbs sampling的code

方法一：

double x = rand() * 1.0 / (1.0 + RAND_MAX);

        while (true)
        {
          x -= topicScores[newTopic];
          if (x < 0)
            break;
          newTopic++;
        }