// what is matter? never mind. 
#pragma GCC optimize("Ofast")
#pragma GCC optimize("unroll-loops")
//#pragma GCC target("sse,sse2,sse3,sse4,popcnt,abm,mmx,avx,avx2")
#include<bits/stdc++.h>
#define For(i,a,b) for(int i=(a);i<=(b);++i)
#define Rep(i,a,b) for(int i=(a);i>=(b);--i)
#define ll long long
#define ull unsigned long long
//#define int long long
#define SZ(x) ((int)((x).size()))
#define ALL(x) (x).begin(),(x).end()
using namespace std;
inline int read()
{
	char c=getchar();int x=0;bool f=0;
	for(;!isdigit(c);c=getchar())f^=!(c^45);
	for(;isdigit(c);c=getchar())x=(x<<1)+(x<<3)+(c^48);
	return f?-x:x;
}

#define fi first
#define se second
#define pb push_back
#define mkp make_pair
typedef pair<int,int>pii;
typedef vector<int>vi;

#define maxn 100005
#define inf 0x3f3f3f3f

int n;
struct P{
    int x,y;
    P(int xx=0,int yy=0){x=xx,y=yy;}
};
P operator +(P a,P b){
    return P((a.x+b.x)/2,(a.y+b.y)/2);
}
int V=1e7;

int now=0;
bool isok(P a,P b){
    return ((a.x+b.x)%2==0) && ((a.y+b.y)%2==0);
}
struct tri{
    P a,b,c;
    bool iso(){
        return isok(a,b) && isok(b,c) && isok(c,a);
    }
    void out(){
        cout<<a.x<<" "<<a.y<<" "<<b.x<<" "<<b.y<<" "<<c.x<<" "<<c.y<<"\n";
    }
};
vector<tri>res;
bool mul;
void add(P a,P b,P c){
    if(mul){
        a.x*=V,a.y*=V;
        b.x*=V,b.y*=V;
        c.x*=V,c.y*=V;
    }
    res.pb((tri){a,b,c});
}
void add(int a,int b,int c,int d,int e,int f){
    return add(P(a,b),P(c,d),P(e,f));
}

void op8(){
    mul=1;

    add(0,0,0,100,48,90);
    add(100,100,100,0,52,90);

//OP
    add(0,100,50,100,48,90);
    add(50,100,48,90,52,90);
    add(52,90,100,100,50,100);

    add(0,0,48,90,50,0);
    add(48,90,52,90,50,0);
    add(52,90,50,0,100,0);

    mul=0;
    now=8;
}

void op10(){
    mul=1;

    add(0,0,0,100,48,90);
    add(100,100,100,0,52,90);

    P a=P(48,90),b=P(52,90);
    P c=P(49,100),d=P(51,100),e=P(50,90);
    
    add(a,c,P(0,100));

    add(a,c,e);
    add(c,d,e);
    add(d,e,b);

    add(d,b,P(100,100));

    add(0,0,48,90,50,0);
    add(48,90,52,90,50,0);
    add(52,90,50,0,100,0);

    mul=0;
    now=10;
}

void op12(){
    mul=1;

    add(0,0,0,100,48,90);
    add(100,100,100,0,52,90);

    P a=P(48,90),b=P(52,90);
    P c=P(49,100),d=P(51,100),e=P(50,90);
    
    add(a,c,P(0,100));

    add(a,c,e);
    add(c,d,e);
    add(d,e,b);

    add(d,b,P(100,100));

    c.y=d.y=0;

    add(a,c,P(0,0));

    add(a,c,e);
    add(c,d,e);
    add(d,e,b);

    add(d,b,P(100,0));

    mul=0;
    now=10;
}

void op9()
{
	//angry
	mul=1;
	add(40,80,100,0,100,100);
    add(40,80,0,0,100,0);
    add(40,80,39,100,100,100);
    add(39,100,40,80,20,80);
    add(39,100,0,100,20,80);
    add(30,60,20,80,40,80);
    add(0,0,0,65,30,60);
    add(0,65,30,60,20,80);
    add(0,65,20,80,0,100);
    mul=0;
}

void add3(){
    vector<tri>nw;
    For(i,0,res.size()-1){
        if(res[i].iso()){
            P a=res[i].a,b=res[i].b,c=res[i].c;
            For(j,0,res.size()-1) if(i!=j) nw.pb(res[j]);
            res=nw;
            P ab=(a+b),bc=(b+c),ca=(a+c);
            add(ab,bc,ca);
            add(a,ab,ca);
            add(b,ab,bc);
            add(c,ca,bc);
            return; 
        }
    }
    assert(0);
}

signed main(){
    cin>>n;
    if(n<=7) {
        puts("No"); exit(0);
    }
    if(n%3==2) op8();
    else if(n%3==1) op10();
    else if(n%3==0) op9();
    while(res.size()<n) add3();
    assert(res.size()==n);
    puts("Yes");
    for(auto it:res)it.out();
    return 0;
}
/*
6
1000000000
1
3
7
10
20221218
*/