3039.sagews

n = 160
N = 212480000
r = 7000
d = 52

p = (1-(falling_factorial(N-r, n) / falling_factorial(N, n))^d) * 1.0
p

0.239744716797366

1/p - 1

3.17110338596202