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

int A[8][6]={
    0,0,9,4,0,20,
    0,20,9,4,10,20,
    10,20,9,4,11,4,
    10,20,11,4,20,20,
    20,0,20,20,11,4,
    0,0,10,0,9,4,
    10,0,11,4,9,4,
    10,0,20,0,11,4
};
int B[9][6]={
    0,0,20,0,16,12,
    0,0,16,12,0,20,
    0,20,12,14,13,20,
    20,0,20,13,16,12,
    16,16,16,12,20,13,
    16,16,20,13,20,20,
    16,16,20,20,13,20,
    16,16,13,20,12,14,
    16,16,12,14,16,12
};
int C[10][6]={
    0,0,10,8,0,20,
    20,0,10,8,20,20,
    0,20,10,8,20,20,
    0,0,8,0,5,4,
    10,3,5,4,8,0,
    10,3,10,8,5,4,
    10,3,8,0,12,0,
    10,3,12,0,15,4,
    10,3,15,4,10,8,
    20,0,15,4,12,0
};
void gen(auto T[][6],int t){
    queue<vector<int>>q;
    for(int i=0;i<t;i++) {
        vector<int>x;
        for(int j=0;j<6;j++) x.push_back(T[i][j]*((int)5e7));
        q.push(x);
    }
    while(q.size()<k){
        auto x=q.front();
        q.pop();
        int x1=x[0],y1=x[1],x2=x[2],y2=x[3],x3=x[4],y3=x[5];
        int x4=(x1+x2)/2,y4=(y1+y2)/2;
        int x5=(x2+x3)/2,y5=(y2+y3)/2;
        int x6=(x3+x1)/2,y6=(y3+y1)/2;
        q.push({x1,y1,x4,y4,x6,y6});
        q.push({x2,y2,x5,y5,x4,y4});
        q.push({x3,y3,x6,y6,x5,y5});
        q.push({x4,y4,x5,y5,x6,y6});
    }
    for(;!q.empty();){
        auto x=q.front();
        q.pop();
        printf("%d %d %d %d %d %d\n", x[0],x[1],x[2],x[3],x[4],x[5]);
    }
}

int main(){
    cin>>k;
    if(k<8) {cout<<"No"<<endl;}
    else{
        cout<<"Yes"<<endl;
        if(k%3==2) gen(A,8);
        if(k%3==0) gen(B,9);
        if(k%3==1) gen(C,10);
    }
    return 0;
}