#include <bits/stdc++.h>
using namespace std;
#define pb push_back

vector<pair<vector<int>,int>> Cand(int x,int A,int SQ){
	vector<pair<vector<int>,int>> ans;
	for(int mod=2;mod<=min(SQ,A);mod++){
		vector<int> tmp;
		int Z=x;
		while(Z>0){
			tmp.pb(Z%mod);
			Z/=mod;
		}
		ans.pb({tmp,mod});
	}
	sort(ans.begin(),ans.end());
	return ans;
}

bool Solve(int x,int y,int A,int B,int&M1,int&M2){
	if(x==y){
		M1=2;
		M2=2;
		return true;
	}

	int SQ=sqrt(x)+3;
	vector<pair<vector<int>,int>> candA=Cand(x,A,SQ);
	vector<pair<vector<int>,int>> candB=Cand(y,B,SQ);
	for(int i=0,j=0;i<candA.size();i++){
		while(j<candB.size() && candB[j].first<candA[i].first){
			j++;
		}
		if(j<candB.size() && candA[i].first==candB[j].first){
			M1=candA[i].second;
			M2=candB[j].second;
			return true;
		}
	}

	for(int a1=1;a1<=SQ;a1++){
		int L1=x/(a1+1)+1;
		int R1=x/a1;

		L1=max(L1,2);
		R1=min(R1,A);

		int L2=y/(a1+1)+1;
		int R2=y/a1;

		L2=max(L2,2);
		R2=min(R2,B);

		int diff=y-x;
		if(diff%a1==0){
			int gap=diff/a1;
			int L3=max(L2,L1+gap);
			int R3=min(R2,R1+gap);
			if(L3<=R3 && R3-gap>a1){
				M1=R3-gap;
				M2=R3;
				return true;
			}
		}
	}

	return false;
}

int main(){
	int t;
	scanf("%i",&t);
	while(t--){
		int x,y,A,B;
		scanf("%i %i %i %i",&x,&y,&A,&B);
		int M1,M2;
		bool ok;
		if(x<=y){
			ok=Solve(x,y,A,B,M1,M2);
		}else{
			ok=Solve(y,x,B,A,M2,M1);
		}
		if(ok){
			printf("YES\n");
			printf("%i %i\n",M1,M2);
		}else{
			printf("NO\n");
		}
	}
	return 0;
}