# 1 ... g^n + c^n + p^n
# 3 ... 1 - ((g + c)^n + (c + p)^n + (p + g)^n) + g^n + c^n + p^n
MOD = 10 ** 9 + 7
for _ in range(int(input())):
    n, ag, bg, ac, bc, ap, bp = map(int, input().split())
    g = ag * pow(bg, MOD-2, MOD) % MOD
    c = ac * pow(bc, MOD-2, MOD) % MOD
    p = ap * pow(bp, MOD-2, MOD) % MOD
    t1 = sum(pow(x, n, MOD) for x in [g, c, p]) % MOD
    t2 = (1 - sum(pow(x + y, n, MOD) for x, y in [(g, c), (c, p), (p, g)]) + t1) % MOD
    ans = (t1 + t2) % MOD
    print(ans)