Comment by asplake

Oops yes! And unfortunately it still doesn't work for n=0

Edit: prod(range(1, n+1)) does work