M=10**9+7
def inv(a):
    return pow(a,M-2,M)

T=int(input())
for t in range(T):
    n,ag,bg,ac,bc,ap,bp=map(int,input().split())
    prob=[ag*inv(bg),ac*inv(bc),ap*inv(bp)]
    anti=0
    for i in range(3):
        anti+=pow(1-prob[i],n,M)-pow(prob[(i+1)%3],n,M)-pow(prob[(i+2)%3],n,M)
    print((1-anti)%M)