docs: public static double approximateBinomialCoefficient(int n, int k) {
      final int m = max(k, n - k);
      long temp = 1;
      for(int i = n, j = 1; i > m; i--, j++) {
        temp = temp * i / j;
      }
      return temp;
    }