// g++-13 1.cpp -std=c++17 -O2 -I .
#include <bits/stdc++.h>
using namespace std;

#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")

#include<ext/pb_ds/assoc_container.hpp>
#include<ext/pb_ds/tree_policy.hpp>
#include<ext/pb_ds/tag_and_trait.hpp>
using namespace __gnu_pbds;

using ll = long long;
using ld = long double;
 
using vi = vector<int>;
using vvi = vector<vi>;
using vll = vector<ll>;
using vvll = vector<vll>;
using vld = vector<ld>;
using vvld = vector<vld>;
using vst = vector<string>;
using vvst = vector<vst>;
 
#define fi first
#define se second
#define pb push_back
#define eb emplace_back
#define pq_big(T) priority_queue<T,vector<T>,less<T>>
#define pq_small(T) priority_queue<T,vector<T>,greater<T>>
#define all(a) a.begin(),a.end()
#define rep(i,start,end) for(ll i=start;i<(ll)(end);i++)
#define per(i,start,end) for(ll i=start;i>=(ll)(end);i--)
#define uniq(a) sort(all(a));a.erase(unique(all(a)),a.end())

random_device seed;
mt19937_64 randint(seed());

ll grr(ll mi, ll ma) { // [mi, ma)
    return mi + randint() % (ma - mi);
}

ll mx = 1e9;

using P = pair<int,int>;
using T = tuple<P,P,P>;

P mid(P &a,P &b){
  auto [ax,ay]=a;
  auto [bx,by]=b;
  return {(ll)(ax+bx)/2,(ll)(ay+by)/2};
}

vector<T> plus3(T &tr){
  auto [i,j,k]=tr;
  auto ni=mid(j,k),nj=mid(k,i),nk=mid(i,j);
  vector<T> res;
  res.emplace_back(i,nj,nk);
  res.emplace_back(ni,j,nk);
  res.emplace_back(ni,nj,k);
  res.emplace_back(ni,nj,nk);
  return res;
}

void sub_output(P &p){
  cout<<p.first<<" "<<p.second<<" ";
}

ll len(P &a,P &b){
  auto [ax,ay]=a;
  auto [bx,by]=b;
  return (ll)(ax-bx)*(ax-bx)+(ll)(ay-by)*(ay-by);
}

void is_ok(T &tr){
  auto [i,j,k]=tr;
  vll l={len(i,j),len(j,k),len(k,i)};
  sort(all(l));
  if(l[0]+l[1]<=l[2]){
    cerr<<"FAILED!"<<endl;
    cerr<<l[0]<<" "<<l[1]<<" "<<l[2]<<endl;
  }
  return;
}

void output(T &tr){
  is_ok(tr);
  auto [i,j,k]=tr;
  sub_output(i);
  sub_output(j);
  sub_output(k);
  cout<<endl;
}

int main(){
  ios::sync_with_stdio(false);
  cin.tie(nullptr);

  // https://kennji72019.blog.fc2.com/blog-entry-9326.html
  // https://ics.uci.edu/~eppstein/junkyard/acute-square/

  int k;cin>>k;
  if(k<=7){
    cout<<"No"<<endl;
    return 0;
  }
  cout<<"Yes"<<endl;

  int md=k%3;
  vector<T> ans;

  auto push=[&](int ax,int ay,int bx,int by,int cx,int cy,int l=1){
    ax*=l;ay*=l;bx*=l;by*=l;cx*=l;cy*=l;
    ans.push_back({{ax,ay},{bx,by},{cx,cy}});
  };

  if(md==0){
    // k = 9
    int l=mx/100;
    push(0,0,39,0,30,27,l);
    push(0,0,30,27,0,40,l);
    push(30,27,0,40,32,44,l);
    push(0,40,0,100,32,44,l);
    push(0,100,100,100,40,30,l);
    push(100,100,40,30,100,0,l);
    push(39,0,100,0,40,30,l);
    push(39,0,30,27,40,30,l);
    push(30,27,40,30,32,44,l);
  }
  else if(md==1){
    // k = 10
    int l=mx/100;
    push(0,100,30,30,100,100,l);
    push(100,0,30,30,100,100,l);
    push(0,24,0,100,30,30,l);
    push(24,0,100,0,30,30,l);
    push(0,0,15,15,0,24,l);
    push(0,0,15,15,24,0,l);
    push(0,24,15,15,20,28,l);
    push(24,0,15,15,28,20,l);
    push(30,30,28,20,20,28,l);
    push(15,15,20,28,28,20,l);
  }
  else{
    // k = 8
    int l=mx/100;
    push(0,0,50,0,49,90,l);
    push(100,0,50,0,51,90,l);
    push(49,90,50,0,51,90,l);
    push(0,100,50,100,49,90,l);
    push(100,100,50,100,51,90,l);
    push(49,90,50,100,51,90,l);
    push(0,0,0,100,49,90,l);
    push(100,0,100,100,51,90,l);
  }

  while(ans.size()<k){
    int n=ans.size();
    vector<T> nans;
    rep(i,1,n)nans.emplace_back(ans[i]);
    auto add=plus3(ans[0]);
    rep(i,0,4)nans.emplace_back(add[i]);
    ans=nans;
  }

  rep(i,0,k){
    output(ans[i]);
  }
}