import sys

# sys.setrecursionlimit(200005)
int1 = lambda x: int(x)-1
pDB = lambda *x: print(*x, end="\n", file=sys.stderr)
p2D = lambda x: print(*x, sep="\n", end="\n\n", file=sys.stderr)
def II(): return int(sys.stdin.readline())
def LI(): return list(map(int, sys.stdin.readline().split()))
def LLI(rows_number): return [LI() for _ in range(rows_number)]
def LI1(): return list(map(int1, sys.stdin.readline().split()))
def LLI1(rows_number): return [LI1() for _ in range(rows_number)]
def SI(): return sys.stdin.readline().rstrip()

# dij = [(0, 1), (-1, 0), (0, -1), (1, 0)]
dij = [(0, 1), (-1, 0), (0, -1), (1, 0), (1, 1), (1, -1), (-1, 1), (-1, -1)]
inf = (1 << 63)-1
# inf = (1 << 31)-1
md = 10**9+7
# md = 998244353

#   f(x)=[a0, a1, a2, a3...]
# pxf(x)=[0 ,pa0,pa1,pa2,pa3...]
# x2f(x)=[0,   0, a0, a1, a2,a3...]
# (x2+px-1)f(x)=[-a0,-a1+pa0]
# (x4+2px3+(p^2-2)x2-2px+1)f^2(x)=[....]

p = II()
aa = [0, 1]
for _ in range(2):
    aa.append((aa[-2]+aa[-1]*p)%md)
ans = [0]*4
for i in range(4):
    for s in range(i+1):
        ans[i] += aa[s]*aa[i-s]
    ans[i] %= md

a, b = 2*p%md, (pow(p, 2, md)-2)%md
for _ in range(2000000):
    ans.append((-ans[-4]-a*ans[-3]-b*ans[-2]+a*ans[-1])%md)

for _ in range(II()):
    q = II()
    print(ans[q-2])