size_t binomial(size_t n, size_t k)
{
    size_t c = 1;

    for (size_t i = 0; i < k; i++)
    {
        c *= n - i;
        c /= i + 1;
    }

    return c;
}