sum(binomial(n-1, r-1-k) * binomial(r, k)) for k = 0 to (r-1)