#include<bits/stdc++.h>
#define PII pair<int,int>
#define lim 1e7
using namespace std;
inline int read(){int s=0,f=1;char c=getchar();while(c<'0'||c>'9'){if(c=='-')f*=-1;c=getchar();}while(c>='0'&&c<='9')s=(s<<1)+(s<<3)+(c^48),c=getchar();return s*f;}
const int N=2010;
int n,l,r;
struct Point{
    int x,y;
    friend Point operator + (Point a,Point b){return {a.x+b.x,a.y+b.y};}
    friend Point operator * (Point a,int k){return {a.x*k,a.y*k};}
    friend Point operator / (Point a,int k){return {a.x/k,a.y/k};}
};
struct Tri{
    Point a,b,c;
}q[N];
int main(){
    n=read(),l=1,r=0;
    if(n<8){printf("No\n");return 0;}
    printf("Yes\n");
    if(!(n%3)){/*构造三角形*/
        n-=9;
        q[++r]={{0,100},{42,70},{55,100}};
		q[++r]={{0,100},{70,50},{0,0}};
		q[++r]={{0,0},{70,50},{100,0}};
		q[++r]={{100,0},{70,50},{100,60}};
		q[++r]={{100,60},{70,75},{100,100}};
		q[++r]={{100,100},{55,100},{70,75}};
		q[++r]={{55,100},{70,75},{42,70}};
		q[++r]={{42,70},{70,50},{70,75}};
		q[++r]={{70,50},{100,60},{70,75}};
    }else if(n%3==1){
        n-=10;
        q[++r]={{0,0},{40,0},{30,18}};
		q[++r]={{30,18},{50,14},{40,0}};
		q[++r]={{40,0},{50,14},{60,0}};
		q[++r]={{60,0},{50,14},{70,18}};
		q[++r]={{60,0},{70,18},{100,0}};
		q[++r]={{30,18},{50,14},{50,30}};
		q[++r]={{50,30},{70,18},{50,14}};
		q[++r]={{0,0},{50,30},{0,100}};
		q[++r]={{0,100},{100,100},{50,30}};
		q[++r]={{100,100},{50,30},{100,0}};
    }else{
        n-=8;
        q[++r]={{0,0},{45,20},{50,0}};
		q[++r]={{45,20},{50,0},{55,20}};
		q[++r]={{50,0},{55,20},{100,0}};
		q[++r]={{0,0},{0,100},{45,20}};
		q[++r]={{0,100},{45,20},{50,100}};
		q[++r]={{50,100},{45,20},{55,20}};
		q[++r]={{50,100},{100,100},{55,20}};
		q[++r]={{100,100},{55,20},{100,0}};
    }
    for(int i=l;i<=r;++i)q[i].a=q[i].a*lim,q[i].b=q[i].b*lim,q[i].c=q[i].c*lim;
    while(n){
        n-=3;
        Tri hd=q[l];++l;
        Point ma=(hd.a+hd.b)/2,mb=(hd.b+hd.c)/2,mc=(hd.a+hd.c)/2;
        q[++r]={ma,mb,mc},q[++r]={hd.a,ma,mc},q[++r]={hd.b,mb,ma},q[++r]={hd.c,mb,mc};
    }
    for(int i=l;i<=r;++i)printf("%d %d %d %d %d %d\n",q[i].a.x,q[i].a.y,q[i].b.x,q[i].b.y,q[i].c.x,q[i].c.y);
    return 0;
}