#include <bits/stdc++.h>
#define endl "\n"
typedef long long i64;
typedef double lb;

signed main(){
    std::cin.tie(nullptr) -> sync_with_stdio(false);
    int n;
    std::cin >> n;
    std::cout << "Yes" << endl;
    std::map<int, std::vector<int> > mp;
    for(int i = 0, x, y;i < n;i++){
        std::cin >> x >> y;
        mp[x].emplace_back(y);
    }

    std::cout << std::fixed << std::setprecision(3);
    int tx = -1, ty = -1;
    for(auto &[x, f] : mp){
        sort(f.begin(), f.end());
        int i = 0, m = f.size();
        if(tx == -1){
            for(;i + 1 < m;i += 2){
                lb l = f[i], r = f[i + 1];
                lb y = x - (r - l);
                std::cout << y << ' ' << l << ' ' << x << ' ' << r << endl;
            }
            if(i == m - 1) tx = x, ty = f.back();
        }
        else{
            for(;i + 1 < m && f[i + 1] < ty;i += 2){
                lb l = f[i], r = f[i + 1];
                lb y = x - (r - l);
                std::cout << y << ' ' << l << ' ' << x << ' ' << r << endl;
            }

            if(i < m && f[i] <= ty){
                int delx = x - tx;
                int dely = ty - f[i];
                if(i == m - 1 || f[i + 1] > ty){
                    int len = std::max(delx, dely);
                    std::cout << x - len << ' ' << ty - len << ' ' << x << ' ' << ty << endl;
                    tx = ty = -1;
                    i ++;
                }
                else{
                    if(delx > dely){
                        int l = f[i], r = f[i + 1];
                        int y = x - (r - l);
                        std::cout << y << ' ' << l << ' ' << x << ' ' << r << endl;
                        i += 2;
                    }
                    else{
                        int len = x - tx; 
                        std::cout << x - len << ' ' << ty - len << ' ' << x << ' ' << ty << endl;
                        tx = ty = -1;
                        i ++;
                    }
                }
            }
            if(tx >= 0 && i < m){
                int len = std::max(x - tx, f[i] - ty);
                std::cout << x - len << ' ' << f[i] - len << ' ' << x << ' ' << f[i] << endl;
                tx = ty = -1;
                i ++;
            }

            for(;i + 1 < m;i += 2){
                lb l = f[i], r = f[i + 1];
                lb y = x - (r - l);
                std::cout << y << ' ' << l << ' ' << x << ' ' << r << endl;
            }
            if(i == m - 1) tx = x, ty = f.back();
        }
    }

    return 0;
}