#include<bits/stdc++.h>
#define int long long
using namespace std;
int T,x,y,A,B,W[100],cntw;
bool flag,NO;
bool num_is_ok(int num,int lim){
	if(num>=2&&num<=lim) return 1;
	return 0;
}
bool check_w0(int a,int b,int w1){
	if(x-a*w1!=y-b*w1) return 0;
	int w0=x-a*w1;
	if(w1<1||w1>=min(a,b)) return 0;
	if(w0<1||w0>=min(a,b)) return 0;
	if(a*w1+w0==x&&b*w1+w0==y&&num_is_ok(a,A)&&num_is_ok(b,B)){
		if(flag) swap(a,b);
		printf("YES\n%lld %lld\n",a,b),NO=0;
		return 1;
	}
	return 0;
}
bool cal(int w1){
	int a=x/w1,b=y/w1;
	if(check_w0(a,b,w1)) return 1;
	if(check_w0(a+1,b+1,w1)) return 1;
	if(check_w0(a-1,b-1,w1)) return 1;
	return 0;
}
void DI(int Num,int bas){
	cntw=0;
	while(Num){
		W[++cntw]=Num%bas;
		Num/=bas;
	}
}
int num_base_on(int di){
	int sum=0,mi=1;
	for(int i=1;i<=cntw;i++){
		sum+=W[i]*mi;
		mi*=di;
	}
	return sum;
}
int bin_sea(){
	int l=2,r=A;
	while(l<=r){
		int mid=(l+r)>>1;
		int NUM_BASE_ON=num_base_on(mid);
		if(NUM_BASE_ON==x) return mid; 
		else if(NUM_BASE_ON<x) l=mid+1;
		else r=mid-1;
	} 
	return -1;
}
bool work_on_less_than_sqrt(int b){
	DI(y,b);
	int a=bin_sea();
	if(a!=-1){
		if(flag) swap(a,b);
		printf("YES\n%lld %lld\n",a,b),NO=0;
		return 1;
	}
	return 0;
}
signed main(){
	scanf("%lld",&T);
	while(T--){
		flag=0;
		scanf("%lld%lld%lld%lld",&x,&y,&A,&B);
		if(x<y) swap(x,y),swap(A,B),flag=1;
		NO=1;
		for(int i=1;i*i<=(x-y);i++)
			if((x-y)%i==0){
				if(cal((x-y)/i)) break;
				if(cal(i)) break;
			}
		for(int b=2;b<=B&&b*b<=y;b++)
			if(work_on_less_than_sqrt(b)) break;
		if(NO&&x==y&&A>=x&&B>=y){
			if(x==1&&num_is_ok(x+1,A)&&num_is_ok(x+1,B)) printf("YES\n%lld %lld\n",x+1,x+1),NO=0;
			else if(num_is_ok(x,A)&&num_is_ok(y,B)) printf("YES\n%lld %lld\n",x,x),NO=0;
		}
		if(NO) printf("NO\n");
	}
	return 0;
}