#!/usr/bin/env python3
# from typing import *

import sys
import io
import math
import collections
import decimal
import itertools
import bisect
import heapq


def input():
    return sys.stdin.readline()[:-1]


# sys.setrecursionlimit(1000000)

# _INPUT = """5
# 2 1 3 1 3 1 3
# 2 1 2 1 3 1 6
# 3 1 3 1 3 1 3
# 4 1 5 2 5 4 10
# 5 1 1 0 1 0 1

# """
# sys.stdin = io.StringIO(_INPUT)

MOD = 1000000007

def solve(N, A, B):
    P = 0
    bn = [pow(B[i], N, MOD) for i in range(3)]
    P = bn[0] * bn[1] * bn[2]
    P -= sum(pow(B[i]-A[i], N, MOD) * bn[(i+1)%3] * bn[(i+2)%3] for i in range(3))
    P += 2 * sum(pow(A[i], N, MOD) * bn[(i+1)%3] * bn[(i+2)%3] for i in range(3))
    P %= MOD

    Q = bn[0] * bn[1] * bn[2]
    Q %= MOD

    S = P * pow(Q, -1, MOD)
    return S % MOD


T0 = int(input())
A = [0, 0, 0]
B = [0, 0, 0]
for _ in range(T0):
    N, A[0], B[0], A[1], B[1], A[2], B[2] = map(int, input().split())
    print(solve(N, A, B))