#include <boost/math/distributions/hypergeometric.hpp>
#include <boost/math/policies/policy.hpp>
#include <iostream>

using namespace std;
using namespace boost;

int main()
{
  unsigned N = 16086184;
  unsigned n = 256004;
  unsigned Q = 251138;
  math::hypergeometric_distribution<double> hyper(n, Q, N);
  cout <<  math::pdf<double>(hyper, 4000) << " " <<  math::pdf<double>(hyper, 4001) << " " <<  math::pdf<double>(hyper, 4002) << "\n";

  return 0;

}