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

const int N = 1e5 + 10;
const double eps = 1e-6;

const int dA[] = {1, 0, 1, 1};
const int dB[] = {0, 1, 1, -1};

struct Point{
	double x, y;
}P[N];
int n;

struct Line{double A, B, C;};

bool check_on_line(Point p, Line l){return abs(l.A * p.x + l.B * p.y + l.C) <= eps;}
bool legal_point(Point p){return abs(p.x) <= 1e9 && abs(p.y) <= 1e9 && (p.x - floor(p.x)) <= eps && (p.y - floor(p.y)) <= eps;}
Line calc_line(Point p, int fx){return (Line){dA[fx], dB[fx], - dA[fx] * p.x - dB[fx] * p.y};}
Point cross_line(Line l1, Line l2){
	if(l1.A == 0) swap(l1, l2);
	if(l2.A){l2.A -= l1.A, l2.B -= l1.B, l2.C -= l1.C;}
	double y = -l2.C / l2.B;
	double x = (-l1.B * y - l1.C) / l1.A;
	return (Point){x, y};
}


map<double, int>ids;

bool check(Point p){
	for (int i = 1; i <= n; ++ i){
		int flag = 0;
		for (int fx = 0; fx < 4; ++ fx) flag |= check_on_line(p, (Line){dA[fx], dB[fx], - dA[fx] * P[i].x - dB[fx] * P[i].y});
		if(flag == 0) return 0;
	}
	return 1;
}

pair<int, Point> two_point_check(Point p1, Point p2, int disable_fx){
	for (int fx1 = 0; fx1 < 4; ++ fx1)
		for (int fx2 = 0; fx2 < 4; ++ fx2){
			if(fx1 == disable_fx || fx2 == disable_fx || fx1 == fx2) continue;
			Point p = cross_line((Line){dA[fx1], dB[fx1], - dA[fx1] * p1.x - dB[fx1] * p1.y}, (Line){dA[fx2], dB[fx2], - dA[fx2] * p2.x - dB[fx2] * p2.y});
			if(legal_point(p) && check(p)) return make_pair(1, p);
		}
	return make_pair(-1, (Point){0, 0});
}

pair<int, Point> find(int A, int B, int disable_fx){
	ids.clear();
	
	for (int i2 = 1; i2 <= n; ++ i2){
		if(ids[A * P[i2].x + B * P[i2].y]){
			int i1 = ids[A * P[i2].x + B * P[i2].y];
			Line l = (Line){A, B, - A * P[i2].x - B * P[i2].y};
			for (int i = 1; i <= n; ++ i)
				if (!check_on_line(P[i], l)){
					for (int fx = 0; fx < 4; ++ fx){
						Point p = cross_line(l, (Line){dA[fx], dB[fx], - dA[fx] * P[i].x - dB[fx] * P[i].y});
						if(legal_point(p) && check(p)) return make_pair(1, p);
					}
					break;
				}
			pair<int, Point> re = two_point_check(P[i1], P[i2], disable_fx);
			return re;
		}
		ids[A * P[i2].x + B * P[i2].y] = i2;
	}
	return make_pair(0, (Point){0, 0});
}

void print(pair<int, Point> ans){
	cout << (ans.first == 1 ? "Yes\n" : "No\n");
	if(ans.first == 1) cout << fixed << setprecision(0) << ans.second.x << " " << ans.second.y << "\n";
}

int main(){
	ios::sync_with_stdio(0); cin.tie(0), cout.tie(0);
	int T; cin >> T;
	while(T --){
		cin >> n;
		
		for (int i = 1; i <= n; ++ i) cin >> P[i].x >> P[i].y;
		
		if(n == 1) {print(make_pair(1, P[1])); continue;}
		
		
		int flag = 0;
		
		for (int i = 0; i < 4; ++ i){
			pair<int, Point> ans = find(dA[i], dB[i], i);
			if(ans.first != 0){print(ans); flag = 1; break;}
		}
		
		if(flag) continue;
		
		pair<int, Point> ans = two_point_check(P[1], P[2], -1);
		
		print(ans);
		
	}
}