c - 为什么这个泊松分布实现在 msvc 上返回 0

Question

我很好奇为什么下面的实现使用 visual studio 编译器总是返回 0，但是用 gcc 和 clang 编译时它工作正常

int poissrand(double lambda){
  int k=0;
  double L=exp(-lambda), p=1;
  do {
    ++k;
    p *= rand()/(double)INT_MAX;
  } while (p > L);
  return --k;
}

示例 here (vc) 和 here (海合会)

Answer 1

在 MSVC 中，rand() 的范围不同。在所有版本中，它由 stdlib.h 中的 RAND_MAX 指定。

如果您将 p 乘以的值应该在 0..1 范围内，那么请尝试

p *= (double)rand() / RAND_MAX;