def rademacher_penalty_bound(n, delta, growth_function): return sqrt(2 * log(2 * n * growth_function(n)) / n) + sqrt(2/n * log(1/delta)) + 1/n