def factorielle(n): assert n==int(n) and n>=0 p = 1 for k in range(1,n+1): p *= k return p