#include<bits/stdc++.h>
using namespace std; 
typedef long long ll;
typedef array<int,2> S;
int n;
S a[3005];
void Solve(){
	scanf("%d",&n);
	for(int i=1;i<=n;i++)scanf("%d%d",&a[i][0],&a[i][1]);
	sort(a+1,a+n+1);
	puts("Yes");
	for(int i=n/2;i>=1;i--){
		int x=2*i,y=2*i-1;
		if(a[x][0]==a[y][0]){
			cout<<a[x][0]<<' '<<a[x][1]<<' ';
			cout<<a[x][0]+(a[x][1]-a[y][1])<<' '<<a[y][1]<<'\n';
			continue;
		}
		if(a[x][1]>=a[y][1]){
			cout<<a[y][0]<<' '<<a[y][1]<<' ';
			int w=max(a[x][0]-a[y][0],a[x][1]-a[y][1]);
			cout<<a[y][0]+w<<' '<<a[y][1]+w<<'\n';
			continue;
		}
		int z=2*i-2;
		if(z==0||a[z][0]!=a[y][0]||a[z][1]<a[x][1]){
			cout<<a[y][0]<<' '<<a[x][1]<<' ';
			int w=max(a[x][0]-a[y][0],a[y][1]-a[x][1]);
			cout<<a[y][0]+w<<' '<<a[x][1]+w<<'\n';
			continue;
		}
		if(a[z][1]>a[x][1]){
			cout<<a[z][0]<<' '<<a[z][1]<<' ';
			cout<<a[z][0]+(a[y][1]-a[z][1])<<' '<<a[y][1]<<'\n';
			swap(a[x],a[z]);
			continue;
		}
		if(a[y][1]-a[z][1]<=a[x][0]-a[z][0]){
			printf("%.1lf %d ",a[y][0]-0.5,a[y][1]);
			printf("%.1lf %d\n",a[y][0]-0.5+(a[y][1]-a[z][1]),a[z][1]);
			swap(a[x],a[z]);
		}
		else {
			cout<<a[z][0]<<' '<<a[z][1]<<' ';
			cout<<a[x][0]<<' '<<a[z][1]+(a[x][0])<<'\n';
			swap(a[y],a[z]);
		}
	}
}
int main(){
	int t=1;
	//scanf("%d",&t);
	while(t--)Solve();
}