T
def v(k, n):
x=1
for i in range(k):
x*=(n-i)/n
return x
for i in range(1,5000,1):
print(i, 1-v(i, 16**6))
x=1
for i in range(k):
x*=(n-i)/n
return x
for i in range(1,5000,1):
print(i, 1-v(i, 16**6))
Ну вот тут по сути написано
n!/k!/n**k
n!/k!/n**k
