#include <bits/stdc++.h>
#define x first
#define y second
using namespace std;
using pii = pair<int, int>;
using ll = long long;

int n;
bool vis[3030];
pii p[3030];

int main(){
    cin.tie(nullptr)->sync_with_stdio(0);
    
    cin >> n;
    for(int i = 1; i <= n; i++) cin >> p[i].x >> p[i].y;

    vector<tuple<ll, int, int>> v;
    for(int i = 1; i <= n; i++) for(int j = i + 1; j <= n; j++){
        int dx = p[i].x - p[j].x;
        int dy = p[i].y - p[j].y;
        v.push_back({(ll)dx * dx + (ll)dy * dy, i, j});
    }
    sort(v.begin(), v.end());

    vector<array<double, 4>> ans;
    for(auto &[d, i, j] : v){
        if(vis[i] || vis[j]) continue;
        double mx = (double)(p[i].x + p[j].x) / 2;
        double my = (double)(p[i].y + p[j].y) / 2;
        double w = max(abs(p[i].x - p[j].x), abs(p[i].y - p[j].y));
        ans.push_back({mx - w / 2, my - w / 2, mx + w / 2, my + w / 2});
        vis[i] = 1;
        vis[j] = 1;
    }

    cout << "Yes\n";
    cout << fixed << setprecision(4);
    for(auto &[x1, y1, x2, y2] : ans){
        cout << x1 << ' ' << y1 << ' ' << x2 << ' ' << y2 << '\n';
    }
}