Problem: An undirected graph is called k-regular, if the degrees of all its vertices are equal k. An edge of a connected graph is called a bridge, if after removing it the graph is being split into two connected components.

Build a connected undirected k-regular graph containing at least one bridge, or else state that such graph doesn't exist.

Input

The single line of the input contains integer k (1 ≤ k ≤ 100) — the required degree of the vertices of the regular graph.

Output

Print "NO" (without quotes), if such graph doesn't exist. 

Otherwise, print "YES" in the first line and the description of any suitable graph in the next lines.

The description of the made graph must start with numbers n and m — the number of vertices and edges respectively. 

Each of the next m lines must contain two integers, a and b (1 ≤ a, b ≤ n, a ≠ b), that mean that there is an edge connecting the vertices a and b. A graph shouldn't contain multiple edges and edges that lead from a vertex to itself. A graph must be connected, the degrees of all vertices of the graph must be equal k. At least one edge of the graph must be a bridge. You can print the edges of the graph in any order. You can print the ends of each edge in any order.

The constructed graph must contain at most 106 vertices and 106 edges (it is guaranteed that if at least one graph that meets the requirements exists, then there also exists the graph with at most 106 vertices and at most 106 edges). 

Examples

Input

1


Output

YES
2 1
1 2

Note

In the sample from the statement there is a suitable graph consisting of two vertices, connected by a single edge.

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
struct custom_hash {
  static uint64_t splitmix64(uint64_t x) {
    x += 0x9e3779b97f4a7c15;
    x = (x ^ (x >> 30)) * 0xbf58476d1ce4e5b9;
    x = (x ^ (x >> 27)) * 0x94d049bb133111eb;
    return x ^ (x >> 31);
  }
  size_t operator()(uint64_t x) const {
    static const uint64_t FIXED_RANDOM =
        std::chrono::steady_clock::now().time_since_epoch().count();
    return splitmix64(x + FIXED_RANDOM);
  }
};
using namespace std;
const long long INF = 9223372036854775807;
const long long N = 600005;
const long long mod = 1e9 + 7;
const long long M = 30;
vector<pair<long long, long long> > v;
vector<long long> vn;
void linear(long long k, vector<long long> vk) {
  queue<pair<long long, long long> > q;
  for (long long i = 0; i < vk.size(); i++) {
    for (long long j = 0; j < k - 1; j++) {
      v.push_back({vk[i], vn.back() + 1});
      q.push({vn.back() + 1, 1});
      vn.push_back(vn.back() + 1);
    }
  }
  while (!q.empty()) {
    long long x = q.front().first;
    long long y = q.front().second;
    q.pop();
    if (y == k) continue;
    while (y != k) {
      long long a = q.front().first;
      long long b = q.front().second;
      q.pop();
      v.push_back({x, a});
      b++;
      q.push({a, b});
      y++;
    }
  }
}
int main() {
  ios_base::sync_with_stdio(false);
  cin.tie(NULL);
  long long k;
  cin >> k;
  if (k % 2 == 0)
    cout << "NO" << endl;
  else if (k == 1) {
    cout << "YES" << endl;
    cout << "2 1" << endl;
    cout << "1 2" << endl;
  } else {
    cout << "YES" << endl;
    vn.push_back(1);
    vn.push_back(2);
    v.push_back({1, 2});
    vector<long long> vk;
    for (long long i = 1; i <= k - 1; i++) {
      v.push_back({1, vn.back() + 1});
      vk.push_back(vn.back() + 1);
      vn.push_back(vn.back() + 1);
    }
    linear(k, vk);
    vk.clear();
    for (long long i = 1; i <= k - 1; i++) {
      v.push_back({2, vn.back() + 1});
      vk.push_back(vn.back() + 1);
      vn.push_back(vn.back() + 1);
    }
    linear(k, vk);
    cout << vn.size() << " " << v.size() << endl;
    for (long long i = 0; i < v.size(); i++) {
      printf("%lld %lld\n", v[i].first, v[i].second);
    }
  }
}