Problem: Peter had a cube with non-zero length of a side. He put the cube into three-dimensional space in such a way that its vertices lay at integer points (it is possible that the cube's sides are not parallel to the coordinate axes). Then he took a piece of paper and wrote down eight lines, each containing three integers — coordinates of cube's vertex (a single line contains coordinates of a single vertex, each vertex is written exactly once), put the paper on the table and left. While Peter was away, his little brother Nick decided to play with the numbers on the paper. In one operation Nick could swap some numbers inside a single line (Nick didn't swap numbers from distinct lines). Nick could have performed any number of such operations.

When Peter returned and found out about Nick's mischief, he started recollecting the original coordinates. Help Peter restore the original position of the points or else state that this is impossible and the numbers were initially recorded incorrectly.

Input

Each of the eight lines contains three space-separated integers — the numbers written on the piece of paper after Nick's mischief. All numbers do not exceed 106 in their absolute value.

Output

If there is a way to restore the cube, then print in the first line "YES". In each of the next eight lines print three integers — the restored coordinates of the points. The numbers in the i-th output line must be a permutation of the numbers in i-th input line. The numbers should represent the vertices of a cube with non-zero length of a side. If there are multiple possible ways, print any of them.

If there is no valid way, print "NO" (without the quotes) in the first line. Do not print anything else.

Examples

Input

0 0 0
0 0 1
0 0 1
0 0 1
0 1 1
0 1 1
0 1 1
1 1 1


Output

YES
0 0 0
0 0 1
0 1 0
1 0 0
0 1 1
1 0 1
1 1 0
1 1 1


Input

0 0 0
0 0 0
0 0 0
0 0 0
1 1 1
1 1 1
1 1 1
1 1 1


Output

NO

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
using namespace std;
using ll = long long;
using vvi = vector<vector<int>>;
using vi = vector<int>;
using vvll = vector<vector<long long>>;
using vll = vector<long long>;
using vd = vector<double>;
using vvd = vector<vector<double>>;
using pii = pair<int, int>;
using vpii = vector<pair<int, int>>;
void __print(int x) { cerr << x; }
void __print(long x) { cerr << x; }
void __print(long long x) { cerr << x; }
void __print(unsigned x) { cerr << x; }
void __print(unsigned long x) { cerr << x; }
void __print(unsigned long long x) { cerr << x; }
void __print(float x) { cerr << x; }
void __print(double x) { cerr << x; }
void __print(long double x) { cerr << x; }
void __print(char x) { cerr << '\'' << x << '\''; }
void __print(const char *x) { cerr << '\"' << x << '\"'; }
void __print(const string &x) { cerr << '\"' << x << '\"'; }
void __print(bool x) { cerr << (x ? "true" : "false"); }
template <typename T, typename V>
void __print(const pair<T, V> &x) {
  cerr << '{';
  __print(x.first);
  cerr << ',';
  __print(x.second);
  cerr << '}';
}
template <typename T>
void __print(const T &x) {
  int f = 0;
  cerr << '{';
  for (auto &i : x) cerr << (f++ ? "," : ""), __print(i);
  cerr << "}";
}
void _print() { cerr << "]\n"; }
template <typename T, typename... V>
void _print(T t, V... v) {
  __print(t);
  if (sizeof...(v)) cerr << ", ";
  _print(v...);
}
int SIX_TO_POWER_EIGHT = 1679616;
ll INF = 2e18 + 4;
struct P {
  int x, y, z;
  P() {}
  P(int _x, int _y, int _z) {
    x = _x;
    y = _y;
    z = _z;
  }
  bool operator==(const P &o) { return o.x == x && o.y == y && o.z == z; }
  P operator+(const P &o) { return P(x + o.x, y + o.y, z + o.z); }
  P operator-(const P &o) { return P(x - o.x, y - o.y, z - o.z); }
  P operator-() { return P(-x, -y, -z); }
  P rotate(int a) {
    switch (a) {
      case 0:
        return P(x, y, z);
      case 1:
        return P(x, z, y);
      case 2:
        return P(y, x, z);
      case 3:
        return P(y, z, x);
      case 4:
        return P(z, x, y);
      case 5:
        return P(z, y, x);
    }
  }
  ll square_length() { return (ll)x * x + (ll)y * y + (ll)z * z; }
  ll operator*(const P &o) { return (ll)x * o.x + (ll)y * o.y + (ll)z * o.z; }
};
void solve() {
  vector<P> points;
  for (int i = 0; i < 8; i++) {
    int x, y, z;
    cin >> x >> y >> z;
    points.push_back(P(x, y, z));
  }
  for (int mask = 0; mask < SIX_TO_POWER_EIGHT; mask++) {
    vector<P> np;
    int m = mask;
    for (int i = 0; i < 8; i++) {
      np.push_back(points[i].rotate(m % 6));
      m /= 6;
    }
    for (int i = 0; i < 8; i++) {
      ll md = INF;
      int mc = 0;
      vi min_idx;
      for (int j = 0; j < 8; j++) {
        if (i == j) continue;
        ll sd = (np[i] - np[j]).square_length();
        if (sd == md) {
          mc++;
          min_idx.push_back(j);
        } else if (sd < md) {
          md = sd;
          mc = 1;
          min_idx.clear();
          min_idx.push_back(j);
        }
      }
      if (mc != 3 || md == 0) goto next_mask;
      P a = np[i] - np[min_idx[1]];
      P b = np[i] - np[min_idx[0]];
      P c = np[i] - np[min_idx[2]];
      if (a * b != 0 || b * c != 0 || c * a != 0) goto next_mask;
    }
    cout << "YES" << endl;
    for (int i = 0; i < 8; i++) {
      cout << np[i].x << " " << np[i].y << " " << np[i].z << endl;
    }
    return;
  next_mask:;
  }
  cout << "NO" << endl;
}
int main() {
  ios_base::sync_with_stdio(false);
  cin.tie(NULL);
  solve();
  return 0;
}