#include<bits/stdc++.h>
using namespace std;
int sol(){
	static long long f[655360],N;
	int A,B,C,g;
	cin >> N >> A >> B >> C;
	if(A>B)swap(A,B);
	if(A>C)swap(A,C);
	if(B>C)swap(B,C);
	if(N<2)return 0;
	g=__gcd(__gcd(A,B),C);
	A/=g;B/=g;C/=g;
	fill(f,f+B,1);
	for(int i=B;;++i){
		f[i]=f[i-A]+f[i-B];
		if(i>=C)f[i]+=f[i-C];
		if(f[i]>=N)return g*i;
	}
}
int main(){
	int T;
	for(cin >> T;T--;)cout << sol() << endl;	
}