Problem: The Trinitarian kingdom has exactly n = 3k cities. All of them are located on the shores of river Trissisipi, which flows through the whole kingdom. Some of the cities are located on one side of the river, and all the rest are on the other side.

Some cities are connected by bridges built between them. Each bridge connects two cities that are located on the opposite sides of the river. Between any two cities exists no more than one bridge.

The recently inaugurated King Tristan the Third is busy distributing his deputies among cities. In total there are k deputies and the king wants to commission each of them to control exactly three cities. However, no deputy can be entrusted to manage the cities, which are connected by a bridge — the deputy can set a too high fee for travelling over the bridge to benefit his pocket, which is bad for the reputation of the king.

Help King Tristan the Third distribute the deputies between the cities, if it is possible.

Input

The first line contains two integers n and m — the number of cities and bridges (3 ≤ n < 105, n = 3k, 0 ≤ m ≤ 105). Next m lines describe the bridges. The i-th line contains two integers ai and bi — the numbers of cities that are connected by the i-th bridge (1 ≤ ai, bi ≤ n, ai ≠ bi, 1 ≤ i ≤ m).

It is guaranteed that no bridge connects a city with itself and that any two cities are connected with no more than one bridge.

Output

If distributing the deputies in the required manner is impossible, print in a single line "NO" (without the quotes).

Otherwise, in the first line print "YES" (without the quotes), and in the second line print which deputy should be put in charge of each city. The i-th number should represent the number of the deputy (from 1 to k), who should be in charge of city numbered i-th in the input — overall there should be n numbers.

If there are multiple solutions, print any of them.

Examples

Input

6 6
1 2
4 1
3 5
6 5
2 6
4 6


Output

YES
1 2 1 2 2 1 

Input

3 1
1 2


Output

NO

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
using namespace std;
template <class T>
inline void read(T& x) {
  char c;
  bool f = 0;
  while ((c = getchar()) < '0' || '9' < c) f |= (c == '-');
  for (x = (c ^ 48); '0' <= (c = getchar()) && c <= '9';
       x = (x << 1) + (x << 3) + (c ^ 48))
    ;
  if (f) x = -x;
}
template <class T>
inline T read(const T sample) {
  T x = 0;
  char c;
  bool f = 0;
  while ((c = getchar()) < '0' || '9' < c) f |= (c == '-');
  for (x = (c ^ 48); '0' <= (c = getchar()) && c <= '9';
       x = (x << 1) + (x << 3) + (c ^ 48))
    ;
  return f ? -x : x;
}
template <class T>
void fwrit(const T x) {
  if (x < 0) return (void)(putchar('-'), fwrit(-x));
  if (x > 9) fwrit(x / 10);
  putchar(x % 10 ^ 48);
}
template <class T>
inline T Max(const T x, const T y) {
  return x > y ? x : y;
}
template <class T>
inline T Min(const T x, const T y) {
  return x < y ? x : y;
}
template <class T>
inline T fab(const T x) {
  return x > 0 ? x : -x;
}
inline int gcd(const int a, const int b) { return b ? gcd(b, a % b) : a; }
inline void getInv(int inv[], const int lim, const int MOD) {
  inv[0] = inv[1] = 1;
  for (int i = 2; i <= lim; ++i)
    inv[i] = 1ll * inv[MOD % i] * (MOD - MOD / i) % MOD;
}
inline long long mulMod(const long long a, const long long b,
                        const long long mod) {
  return ((a * b - (long long)((long double)a / mod * b + 1e-8) * mod) % mod +
          mod) %
         mod;
}
const int MAXN = 1e5;
const int MAXM = 1e5;
struct edge {
  int to, nxt;
  edge(const int T = 0, const int N = 0) : to(T), nxt(N) {}
} e[(MAXM << 1) + 5];
int tail[MAXN + 5], nsz[MAXN + 5], ecnt;
inline void add_edge(const int u, const int v) {
  e[++ecnt] = edge(v, tail[u]);
  tail[u] = ecnt;
  e[++ecnt] = edge(u, tail[v]);
  tail[v] = ecnt;
  ++nsz[u], ++nsz[v];
}
int n, m;
inline void Init() {
  n = read(1), m = read(1);
  for (signed i = (1), i_end_ = (m); i <= i_end_; ++i)
    add_edge(read(1), read(1));
}
vector<int> node[2];
int sz[2];
int col[MAXN + 5];
int g[MAXN + 5], gcnt;
void dfs_col(const int u, const int c) {
  col[u] = c;
  node[c].push_back(u);
  ++sz[c];
  for (signed i = tail[u], v = e[i].to; i; i = e[i].nxt, v = e[i].to)
    if (col[v] == -1) dfs_col(v, c ^ 1);
}
inline void Getcol() {
  memset(col, -1, sizeof col);
  for (signed i = (1), i_end_ = (n); i <= i_end_; ++i)
    if (col[i] == -1) dfs_col(i, 0);
}
inline void writAns() {
  for (signed i = (1), i_end_ = (n); i <= i_end_; ++i)
    fwrit(g[i]), putchar(' ');
  putchar('\n');
}
bool invalid[MAXN + 5];
bool vis[MAXN + 5];
inline bool Catch(const int c, const bool key) {
  int now = -1;
  for (signed i = (0), i_end_ = (sz[c] - 1); i <= i_end_; ++i)
    if (!invalid[node[c][i]] && nsz[node[c][i]] <= sz[c ^ 1] - 2) {
      now = node[c][i];
      break;
    }
  if (now == -1) {
    if (!key) return false;
    if (Catch(c ^ 1, 0) && Catch(c ^ 1, 0))
      return true;
    else
      return false;
  }
  memset(vis, 0, sizeof vis);
  for (int i = tail[now]; i; i = e[i].nxt) vis[e[i].to] = true;
  g[now] = ++gcnt;
  invalid[now] = true;
  int cnt = 0;
  for (signed i = (0), i_end_ = (sz[c ^ 1] - 1); i <= i_end_; ++i)
    if (!vis[node[c ^ 1][i]] && !invalid[node[c ^ 1][i]]) {
      invalid[node[c ^ 1][i]] = true;
      g[node[c ^ 1][i]] = gcnt, ++cnt;
      if (cnt == 2) break;
    }
  return true;
}
inline void Solve() {
  bool vis[MAXN + 5];
  memset(vis, false, sizeof vis);
  if (sz[1] % 3 == 1) {
    if (!Catch(1, 1)) {
      puts("NO");
      return;
    }
  } else if (sz[0] % 3 == 1) {
    if (!Catch(0, 1)) {
      puts("NO");
      return;
    }
  }
  ++gcnt;
  int nowsz = 0;
  for (signed i = (0), i_end_ = (sz[0] - 1); i <= i_end_; ++i)
    if (!g[node[0][i]]) {
      g[node[0][i]] = gcnt, ++nowsz;
      if (nowsz == 3) ++gcnt, nowsz = 0;
    }
  for (signed i = (0), i_end_ = (sz[1] - 1); i <= i_end_; ++i)
    if (!g[node[1][i]]) {
      g[node[1][i]] = gcnt, ++nowsz;
      if (nowsz == 3) ++gcnt, nowsz = 0;
    }
  puts("YES");
  writAns();
}
signed main() {
  Init();
  Getcol();
  Solve();
  return 0;
}