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;
int n, m;
vector<int> adj[100050];
bool side[100050];
int i1, i2, sz1, sz2;
vector<int> vset[2];
int visited[100050] = {0}, visid = 0;
int visited2[100050] = {0}, visid2 = 0;
int dep[100050];
void dfs(int v, bool c) {
  visited[v] = visid;
  vset[c].push_back(v);
  side[v] = c;
  for (int i = 0; i < adj[v].size(); i++) {
    int u = adj[v][i];
    if (visited[u] == visid) continue;
    dfs(u, !c);
  }
}
inline void pre() {
  ++visid;
  for (int i = 1; i <= n; i++)
    if (visited[i] != visid) dfs(i, 0);
}
inline void output_no() { puts("NO"); }
inline void output_yes() {
  puts("YES");
  for (int i = 1; i <= n; i++) printf("%d ", dep[i]);
  puts("");
}
inline void solve() {
  pre();
  if (vset[0].size() % 3 == 0) {
    int c = 0;
    for (int i = 0; i < vset[0].size(); i++) dep[vset[0][i]] = (c++) / 3 + 1;
    for (int i = 0; i < vset[1].size(); i++) dep[vset[1][i]] = (c++) / 3 + 1;
    output_yes();
    return;
  } else {
    if (vset[0].size() % 3 == 1) {
      i1 = 0;
      i2 = 1;
    } else {
      i1 = 1;
      i2 = 0;
    }
    sz1 = vset[i1].size();
    sz2 = vset[i2].size();
    for (int i = 0; i < sz1; i++) {
      int v = vset[i1][i];
      if ((int)adj[v].size() <= sz2 - 2) {
        int j, u, w;
        ++visid;
        for (j = 0; j < adj[v].size(); j++) visited[adj[v][j]] = visid;
        for (j = 0; j < sz2; j++)
          if (visited[vset[i2][j]] < visid) break;
        u = vset[i2][j];
        for (j = j + 1; j < sz2; j++)
          if (visited[vset[i2][j]] < visid) break;
        w = vset[i2][j];
        dep[v] = dep[u] = dep[w] = n / 3;
        int c = 0;
        for (int j = 0; j < vset[i1].size(); j++) {
          int x = vset[i1][j];
          if (x == v || x == u || x == w) continue;
          dep[x] = (c++) / 3 + 1;
        }
        for (int j = 0; j < vset[i2].size(); j++) {
          int x = vset[i2][j];
          if (x == v || x == u || x == w) continue;
          dep[x] = (c++) / 3 + 1;
        }
        output_yes();
        return;
      }
    }
    int vplus[2], cnt2plus = 0;
    for (int i = 0; i < sz2; i++) {
      int v = vset[i2][i];
      if ((int)adj[v].size() <= sz1 - 2) {
        vplus[cnt2plus++] = v;
      }
      if (cnt2plus >= 2) break;
    }
    if (cnt2plus >= 2) {
      ++visid;
      int c = 0;
      for (int t = 0; t < 2; t++) {
        int v = vplus[t];
        ++visid2;
        for (int i = 0; i < adj[v].size(); i++) visited2[adj[v][i]] = visid2;
        int j, u, w;
        for (j = 0; j < sz1; j++)
          if (visited2[vset[i1][j]] < visid2) break;
        u = vset[i1][j];
        for (j = j + 1; j < sz1; j++)
          if (visited2[vset[i1][j]] < visid2) break;
        w = vset[i1][j];
        dep[v] = dep[u] = dep[w] = c / 3 + 1;
        visited[v] = visited[u] = visited[w] = visid;
        c += 3;
      }
      for (int j = 0; j < vset[i1].size(); j++) {
        int x = vset[i1][j];
        if (visited[x] == visid) continue;
        dep[x] = (c++) / 3 + 1;
      }
      for (int j = 0; j < vset[i2].size(); j++) {
        int x = vset[i2][j];
        if (visited[x] == visid) continue;
        dep[x] = (c++) / 3 + 1;
      }
      output_yes();
      return;
    }
    output_no();
  }
}
int main(void) {
  scanf("%d %d", &n, &m);
  for (int i = 0; i < m; i++) {
    int v, u;
    scanf("%d %d", &v, &u);
    adj[v].push_back(u);
    adj[u].push_back(v);
  }
  solve();
  return 0;
}