#include<bits/stdc++.h>
using namespace std;
using ll=long long;
void solve()
{
    ll a,b,c,x; cin>>a>>b>>c>>x;
    ll lcm = (a*b)/__gcd(a,b);
    lcm = (lcm*c)/__gcd(lcm,c);
    
    ll gcm = (ceil((double)x/lcm)-1)*lcm;
    if(gcm==0) cout<<-1<<endl;
    else cout<<gcm<<endl;

}
int main(){int t; cin>>t; while(t--) solve();}