c - Cにおけるポアソン分布の計算

Question

k の値が 720 までのポアソン分布を計算する C 関数が必要です。非常に効率的なソリューションが必要です。

Answer 1

GSLを試してみてください:gsl_ran_poisson_pdf

Answer 2

ポアソン乱数発生器

int poissonRandom(double expectedValue) {
  int n = 0; //counter of iteration
  double limit; 
  double x;  //pseudo random number
  limit = exp(-expectedValue);
  x = rand() / INT_MAX; 
  while (x > limit) {
    n++;
    x *= rand() / INT_MAX;
  }
  return n;
}

私はあなたの緊急の要求にかなり遅れていると思います.

Answer 3

ライブラリを使用せずに自分で計算する場合は、式を使用して計算できます.. e^k*e^(-lambda)/k!
log(n!) = log(n)+log(n-1!) と動的計画法を使用できます