Problem: Developer Petr thinks that he invented a perpetual motion machine. Namely, he has a lot of elements, which work in the following way.

Each element has one controller that can be set to any non-negative real value. If a controller is set on some value x, then the controller consumes x2 energy units per second. At the same time, any two elements connected by a wire produce y·z energy units per second, where y and z are the values set on their controllers.

Petr has only a limited number of wires, so he has already built some scheme of elements and wires, and is now interested if it's possible to set the controllers in such a way that the system produces at least as much power as it consumes, and at least one controller is set on the value different from 0. Help him check this, and if it's possible, find the required integer values that should be set.

It is guaranteed that if there exist controllers' settings satisfying the above conditions, then there exist required integer values not greater than 106.

Input

There are several (at least one) test cases in the input. The first line contains single integer — the number of test cases.

There is an empty line before each test case. The first line of test case contains two integers n and m (1 ≤ n ≤ 105, 0 ≤ m ≤ 105) — the number of elements in the scheme and the number of wires.

After that, m lines follow, each of them contains two integers a and b (1 ≤ a, b ≤ n) — two elements connected by a wire. No element is connected with itself, no two elements are connected by more than one wire.

It is guaranteed that the sum of n and the sum of m over all test cases do not exceed 105.

For hacks you can only use tests with one test case.

Output

Print answer for each test case.

For each test case print "YES" if it's possible to set the controllers in such a way that the consumed power is not greater than the power produced, and the required values on the next line. The settings should be integers from 0 to 106, inclusive, and at least one value should be different from 0. If there are multiple answers, print any of them.

If it's not possible to set the controllers in the required way, print one line "NO".

Example

Input

4
 
4 4
1 2
2 3
3 4
4 2
 
3 2
2 3
3 1
 
4 6
1 2
3 4
4 2
1 4
1 3
3 2
 
10 9
2 1
3 2
5 2
6 2
2 7
2 8
2 9
2 10
4 2


Output

YES
1 2 2 1
NO
YES
1 1 1 1
YES
1 5 1 1 1 1 1 1 1 1

Note

In the first example it's possible to set the controllers in the required way, for example, in the following way: set 1 on the first element, set 2 on the second and on the third, set 1 on the fourth. The consumed power is then equal to 12 + 22 + 22 + 12 = 10 energy units per second, the produced power is equal to 1·2 + 2·2 + 2·1 + 2·1 = 10 energy units per second. Thus the answer is "YES".

In the second test case it's not possible to set the controllers in the required way. For example, if we set all controllers to 0.5, then the consumed powers equals 0.75 energy units per second, while produced power equals 0.5 energy units per second.

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
using namespace std;
template <typename T>
inline T read(register T& t) {
  register T f = 1;
  register char ch = getchar();
  t = 0;
  while (ch < '0' || ch > '9') {
    if (ch == '-') f = -f;
    ch = getchar();
  }
  while (ch >= '0' && ch <= '9') t = t * 10 + ch - '0', ch = getchar();
  t *= f;
  return t;
}
template <typename T, typename... Args>
inline void read(T& t, Args&... args) {
  read(t);
  read(args...);
}
const long long p = 998244353;
inline long long power(register long long x, register long long k = p - 2) {
  register long long re = 1;
  for (; k; k >>= 1, x = (x * x) % p)
    if (k & 1) re = (re * x) % p;
  return re;
}
int n, m;
int head[100005], ne[200005], to[200005], total = 1;
int deg[100005];
inline void add(int a, int b) {
  deg[a]++;
  total++;
  to[total] = b;
  ne[total] = head[a];
  head[a] = total;
}
int fa[100005];
vector<int> cir;
vector<int> thd;
vector<int> mth;
void find_cir(int now) {
  if (cir.size()) return;
  if (deg[now] > 3) mth.push_back(now);
  if (deg[now] == 3) thd.push_back(now);
  for (int i = head[now]; i; i = ne[i])
    if (to[i] != fa[now]) {
      if (cir.size()) return;
      if (fa[to[i]]) {
        int qwq = now;
        while (qwq != to[i]) cir.push_back(qwq), qwq = fa[qwq];
        cir.push_back(to[i]);
        return;
      }
      fa[to[i]] = now;
      find_cir(to[i]);
    }
}
void reroot(int now) {
  for (int i = head[now]; i; i = ne[i])
    if (to[i] != fa[now]) fa[to[i]] = now, reroot(to[i]);
}
vector<int> lin[3];
void find_lin(int now, int fa, vector<int>& lin) {
  lin.push_back(now);
  for (int i = head[now]; i; i = ne[i])
    if (to[i] != fa) find_lin(to[i], now, lin);
}
void work() {
  memset(deg, 0, sizeof(deg));
  total = 1;
  memset(head, 0, sizeof(head));
  read(n, m);
  for (int i = 1, x, y; i <= m; i++) read(x, y), add(x, y), add(y, x);
  cir.clear();
  memset(fa, 0, sizeof(fa));
  static long long ans[100005];
  memset(ans, 0, sizeof(ans));
  for (int qwq = 1; qwq <= n; qwq++)
    if (!fa[qwq]) {
      thd.clear();
      mth.clear();
      fa[qwq] = qwq;
      find_cir(qwq);
      if (cir.size()) {
        for (auto i : cir) ans[i] = 1;
        printf("YES\n");
        for (int i = 1; i <= n; i++) printf("%lld ", ans[i]);
        printf("\n");
        return;
      }
      if (mth.size()) {
        int now = mth[0];
        ans[now] = 2;
        for (int i = head[now]; i; i = ne[i]) ans[to[i]] = 1;
        printf("YES\n");
        for (int i = 1; i <= n; i++) printf("%lld ", ans[i]);
        printf("\n");
        return;
      }
      if (thd.size() >= 2) {
        int a = thd[0], b = thd[1];
        memset(fa, 0, sizeof(fa));
        reroot(b);
        int aa = a;
        while (aa != b) ans[aa] = 2, aa = fa[aa];
        ans[b] = 2;
        for (int i = head[a]; i; i = ne[i])
          if (ans[to[i]] == 0) ans[to[i]] = 1;
        for (int i = head[b]; i; i = ne[i])
          if (ans[to[i]] == 0) ans[to[i]] = 1;
        printf("YES\n");
        for (int i = 1; i <= n; i++) printf("%lld ", ans[i]);
        printf("\n");
        return;
      }
      if (thd.size()) {
        int now = thd[0];
        lin[0].clear();
        lin[1].clear();
        lin[2].clear();
        int tnt = -1;
        for (int i = head[now]; i; i = ne[i]) find_lin(to[i], now, lin[++tnt]);
        sort(lin, lin + 3, [](vector<int>& a, vector<int>& b) {
          return a.size() < b.size();
        });
        long long a = lin[0].size() + 1, b = lin[1].size() + 1,
                  c = lin[2].size() + 1;
        if (a * c + b * c + a * b > a * b * c) continue;
        if (a * b * c > 1e6) {
          if (lin[0].size() >= 1 && lin[1].size() >= 2 && lin[2].size() >= 5) {
            while (lin[0].size() > 1) lin[0].pop_back();
            while (lin[1].size() > 2) lin[1].pop_back();
            while (lin[2].size() > 5) lin[2].pop_back();
          }
          if (lin[0].size() >= 1 && lin[1].size() >= 3 && lin[2].size() >= 3) {
            while (lin[0].size() > 1) lin[0].pop_back();
            while (lin[1].size() > 3) lin[1].pop_back();
            while (lin[2].size() > 3) lin[2].pop_back();
          }
          if (lin[0].size() >= 2 && lin[1].size() >= 2 && lin[2].size() >= 2) {
            while (lin[0].size() > 2) lin[0].pop_back();
            while (lin[1].size() > 2) lin[1].pop_back();
            while (lin[2].size() > 2) lin[2].pop_back();
          }
        }
        a = lin[0].size() + 1, b = lin[1].size() + 1, c = lin[2].size() + 1;
        ans[now] = a * b * c;
        long long tp;
        tp = ans[now];
        for (auto i : lin[0]) tp -= b * c, ans[i] = tp;
        tp = ans[now];
        for (auto i : lin[1]) tp -= a * c, ans[i] = tp;
        tp = ans[now];
        for (auto i : lin[2]) tp -= a * b, ans[i] = tp;
        printf("YES\n");
        for (int i = 1; i <= n; i++) printf("%lld ", ans[i]);
        printf("\n");
        return;
      }
    }
  printf("NO\n");
  return;
}
int main() {
  int T;
  read(T);
  while (T--) work();
  return 0;
}