$total *= lower($sizes->[$_], $sampled[$_]) / fact($sampled[$_])

##</code><code>##

  $total *= lower($sizes->[$_], $sampled[$_]) / fact($sampled[$_]) for 0..$n-1;