BinomialCoefficients.cc

Go to the documentation of this file.

#include <utl/BinomialCoefficients.h>
#include <utl/ErrorLogger.h>

using namespace utl;


void
BinomialCoefficients::Generate(const int n)
{
  if (n > 500)
    WARNING("Are you sure you know what you're doing?");

  fCoefficients.resize((n+1)*(n+2)/2);
  for (int i = fMaxN+1; i <= n; ++i) {
    GetUnchecked(i, 0) = 1;
    for (int j = 1; j < i; ++j)
      GetUnchecked(i, j) = GetUnchecked(i-1, j-1) + GetUnchecked(i-1, j);
    GetUnchecked(i, i) = 1;
  }
  fMaxN = n;
}