MOD = 1000000007
a = [1]
b = [1]
for i in xrange(1,2000001):
    a.append((i * a[i - 1]) % MOD)
    
for i in xrange(1,1000001):
    b.append(pow(a[i],MOD - 2,MOD))

def P(x,y):
    if x < y:
        return 0
    else:
        return a[x] * b[x - y] % MOD

def C(x,y):
    return P(x,y) * b[y] % MOD

def H(x,y):
    if x == 0 and y == 0:
        return 1
    else:
        return C(x + y - 1,y)

t = input()
for i in xrange(t):
    print eval(raw_input())