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>
#pragma GCC optimize("O3")
using namespace std;
unsigned seed = chrono::system_clock::now().time_since_epoch().count();
mt19937 eng(seed);
long long random2() { return (1ll << 31ll) * eng() + eng(); }
long long n, m, k, q, T;
const long long big = 1000000007;
const long long big2 = 1000000009;
const long long mod = 998244353;
const int MAXN = 300001;
vector<vector<int> > C(MAXN, vector<int>());
long long A[MAXN] = {0};
int CC[MAXN] = {0};
int cc = 0;
bool comp(long long i, long long j) {
  return (int)(C[i]).size() > (int)(C[j]).size();
}
bool mark1[MAXN] = {0};
int par[MAXN] = {0};
vector<int> I;
int cyc1 = -1;
int cyc2 = -1;
void dfs1(int i) {
  CC[i] = cc;
  if (mark1[i]) return;
  mark1[i] = 1;
  I.push_back(i);
  for (auto& y : C[i]) {
    if (y == par[i]) continue;
    if (mark1[y]) {
      if (cyc1 == -1) {
        cyc1 = i;
        cyc2 = y;
      }
      continue;
    }
    par[y] = i;
    dfs1(y);
  }
}
bool mark2[MAXN] = {0};
void dfs2(int i) {
  if (mark2[i]) return;
  mark2[i] = 1;
  for (auto& y : C[i]) {
    if (!mark2[y]) {
      par[y] = i;
      dfs2(y);
    }
  }
}
void fill_path(int from, int to, int with) {
  while (1) {
    A[from] = with;
    if (from == to) break;
    from = par[from];
  }
}
void fill_around(int i, int x, int with) {
  for (auto& y : C[i]) {
    if (A[y] == 0 && x > 0) {
      A[y] = with;
      x--;
    }
  }
}
map<long long, pair<int, int> > M;
pair<int, int> dp(long long i, long long prev, long long par) {
  if (prev <= 1 && (int)(C[i]).size() != 3) return {0, 0};
  long long h = i * n * 1234 + prev * n + par;
  if (M.find(h) != M.end()) return M[h];
  long long lim = prev;
  if ((int)(C[i]).size() == 3) lim = 11;
  pair<int, int> best = {-big, 0};
  for (int val = 0; val < (lim); ++val) {
    if (val == 0 && (int)(C[i]).size() == 3) continue;
    int num = val * prev - val * val;
    for (auto& y : C[i]) {
      if (y != par) {
        num += dp(y, val, i).first;
      }
    }
    if (num > best.first) best = {num, val};
  }
  M[h] = best;
  return best;
}
void get_ans(long long i, long long prev, long long par) {
  long long h = i * n * 1234 + prev * n + par;
  if (M.find(h) == M.end()) return;
  pair<int, int> nxt = M[h];
  A[i] = nxt.second;
  for (auto& y : C[i]) {
    if (y != par) {
      get_ans(y, A[i], i);
    }
  }
}
void print_ans() {
  cout << "YES\n";
  for (int c1 = 0; c1 < (n); ++c1) {
    if (CC[c1] != cc) A[c1] = 0;
    cout << A[c1] << " ";
  }
  cout << "\n";
}
int main() {
  ios_base::sync_with_stdio(0);
  cin.tie(0);
  long long a, b, c, d, e;
  cin >> T;
  for (int _ = 0; _ < (T); ++_) {
    cin >> n >> m;
    for (int c1 = 0; c1 < (n); ++c1) {
      C[c1].clear();
      mark1[c1] = 0;
      mark2[c1] = 0;
      A[c1] = 0;
      par[c1] = -1;
      CC[c1] = cc;
    }
    M.clear();
    for (int c1 = 0; c1 < (m); ++c1) {
      cin >> a >> b;
      a--;
      b--;
      C[a].push_back(b);
      C[b].push_back(a);
    }
    bool solved = 0;
    for (int c1 = 0; c1 < (n); ++c1) {
      if (!mark1[c1]) {
        I.clear();
        cyc1 = -1;
        cyc2 = -1;
        dfs1(c1);
        if (cyc1 != -1) {
          fill_path(cyc1, cyc2, 1);
          print_ans();
          solved = 1;
          break;
        }
        sort(I.begin(), I.end(), comp);
        int i = I[0];
        if ((int)(I).size() <= 1) continue;
        int j = I[1];
        if ((int)(C[i]).size() <= 2) {
          continue;
        }
        if ((int)(C[i]).size() > 3) {
          A[i] = 2;
          fill_around(i, 4, 1);
          print_ans();
          solved = 1;
          break;
        }
        if ((int)(C[j]).size() < 3) {
          pair<int, int> thing = dp(i, 0, -1);
          if (thing.first >= 0) {
            get_ans(i, 0, -1);
            print_ans();
            solved = 1;
            break;
          }
          for (auto& y : I) {
            A[y] = 0;
          }
          continue;
        }
        dfs2(i);
        fill_path(j, i, 2);
        fill_around(i, 2, 1);
        fill_around(j, 2, 1);
        print_ans();
        solved = 1;
        break;
      }
    }
    if (!solved) {
      cout << "NO\n";
    }
  }
  return 0;
}