Function reduce

src/compute-engine/numerics/bernoulli.ts:23–30

(num: bigint, den: bigint)

}

function reduce(num: bigint, den: bigint): [bigint, bigint] {
  if (den < 0n) {
    num = -num;
    den = -den;
  }
  const g = gcd(num, den);
  return g === 0n ? [0n, 1n] : [num / g, den / g];
}

function factorial(n: number): bigint {
  let f = 1n;