import numpy as np

mod=10**9+7

d=np.zeros(100100,dtype=np.int64)

d[0]=1

for i in range(1,10**5+1):

  d[i]=d[i-1]*i%mod

print((d[int(eval(input()))]))
d,mod=[1]+[0]*100100,10**9+7

for i in range(1,10**5+1):d[i]=d[i-1]*i%mod

print((d[int(eval(input()))]))
import numpy as np

mod = 10**9 + 7
d = np.zeros(100100, dtype=np.int64)
d[0] = 1
for i in range(1, 10**5 + 1):
    d[i] = d[i - 1] * i % mod
print((d[int(eval(input()))]))

d, mod = [1] + [0] * 100100, 10**9 + 7
for i in range(1, 10**5 + 1):
    d[i] = d[i - 1] * i % mod
print((d[int(eval(input()))]))

-import numpy as np - -mod = 10**9 + 7 -d = np.zeros(100100, dtype=np.int64) -d[0] = 1 +d, mod = [1] + [0] * 100100, 10**9 + 7