fork download
  1. from scipy.special import*
  2. f=lambda a,b,p:(binom(a*p,b*p)-binom(a,b))/p**3
  3.  
  4. for i in [[6,2,5], [3,1,13], [7,3,13]]:
  5. print('%2s %2s %2s -> %s'%(i[0],i[1],i[2],f(*i)))
Success #stdin #stdout 0.27s 31928KB
stdin
Standard input is empty
stdout
 6  2  5 -> 240360.0
 3  1 13 -> 3697053.0
 7  3 13 -> 3.74036217417e+22