import math as m
t=int(input())
for _ in range(t):
    a,b,c,x = map(int,input().split())
    lcmm = m.lcm(a,b)
    lcmm = m.lcm(lcmm,c)
    gcm = (m.ceil(x/lcmm)-1)*lcmm
    if(gcm):
        print(gcm)
    else:
        print(-1)