// 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 1000005
#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
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);
}
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);
}
int 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
*/