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 <class T>
void splitstr(const string &s, vector<T> &out) {
  istringstream in(s);
  out.clear();
  copy(istream_iterator<T>(in), istream_iterator<T>(), back_inserter(out));
}
template <class T>
T gcd(T a, T b) {
  return b ? gcd(b, a % b) : a;
}
static void redirect(int argc, const char **argv) {
  ios::sync_with_stdio(false);
  cin.tie(NULL);
  if (argc > 1) {
    static filebuf f;
    f.open(argv[1], ios::in);
    cin.rdbuf(&f);
    if (!cin) {
      cerr << "Failed to open '" << argv[1] << "'" << endl;
      exit(1);
    }
  }
  if (argc > 2) {
    static filebuf f;
    f.open(argv[2], ios::out | ios::trunc);
    cout.rdbuf(&f);
    if (!cout) {
      cerr << "Failed to open '" << argv[2] << "'" << endl;
    }
  }
}
template <typename Derived>
class field_mixin {
 public:
  Derived operator-(const Derived &other) const {
    return static_cast<const Derived &>(*this) + -other;
  }
  Derived operator/(const Derived &other) const {
    return static_cast<const Derived &>(*this) * other.inverse();
  }
  Derived &operator+=(const Derived &other) {
    Derived &self = static_cast<Derived &>(*this);
    return self = self + other;
  }
  Derived &operator-=(const Derived &other) {
    Derived &self = static_cast<Derived &>(*this);
    return self = self - other;
  }
  Derived &operator*=(const Derived &other) {
    Derived &self = static_cast<Derived &>(*this);
    return self = self * other;
  }
  Derived &operator/=(const Derived &other) {
    Derived &self = static_cast<Derived &>(*this);
    return self = self / other;
  }
};
template <typename Derived>
class ordered_mixin {
 public:
  bool operator!=(const Derived &other) const {
    const Derived &self = static_cast<const Derived &>(*this);
    return !(self == other);
  }
  bool operator<=(const Derived &other) const {
    const Derived &self = static_cast<const Derived &>(*this);
    return !(other < self);
  }
  bool operator>=(const Derived &other) const {
    const Derived &self = static_cast<const Derived &>(*this);
    return !(self < other);
  }
  bool operator>(const Derived &other) const {
    const Derived &self = static_cast<const Derived &>(*this);
    return other < self;
  }
};
template <typename T>
class rational : public field_mixin<rational<T>>,
                 public ordered_mixin<rational<T>> {
  T num;
  T den;

 public:
  using field_mixin<rational<T>>::operator-;
  rational(T num = 0, T den = 1) {
    assert(den != 0);
    if (den < T(0)) {
      num = -num;
      den = -den;
    }
    while (num % 2 == 0 && den % 2 == 0) {
      num >>= 1;
      den >>= 1;
    }
    this->num = num;
    this->den = den;
  }
  rational operator-() const { return rational(-num, den); }
  rational inverse() const { return rational(den, num); }
  rational operator+(const rational &other) const {
    return rational(num * other.den + den * other.num, den * other.den);
  }
  rational operator*(const rational &other) const {
    return rational(num * other.num, den * other.den);
  }
  bool operator==(const rational &other) const {
    return num * other.den == other.num * den;
  }
  bool operator<(const rational &other) const {
    return num * other.den < other.num * den;
  }
  T numerator() const { return num; }
  T denominator() const { return den; }
  explicit operator double() const { return double(num) / double(den); }
};
template <typename T>
static std::ostream &operator<<(std::ostream &o, const rational<T> &r) {
  return o << r.numerator() << '/' << r.denominator();
}
struct node {
  bool seen = false;
  int parent = -1;
  rational<long long> lambda = 1;
  vector<int> ch;
};
static pair<long long, long long> recurse(vector<node> &nodes,
                                          const vector<vector<int>> &edges,
                                          int cur, int parent) {
  node &n = nodes[cur];
  n.seen = true;
  n.parent = parent;
  rational<long long> r = 1;
  for (int y : edges[cur]) {
    if (y == parent) continue;
    if (nodes[y].seen) return pair<long long, long long>(cur, y);
    auto sub = recurse(nodes, edges, y, cur);
    if (sub.first != -1) return sub;
    r -= rational<long long>(1, 4) / nodes[y].lambda;
    n.ch.push_back(y);
  }
  if (r <= 0)
    return pair<long long, long long>(cur, -1);
  else if (r > rational<long long>(1, 2) && r < rational<long long>(4, 7)) {
    r = nodes[n.ch[0]].lambda;
    nodes[n.ch[0]].lambda = rational<long long>(1, 2);
  }
  n.lambda = r;
  return {-1, -1};
}
static void pushdown(const vector<node> &nodes, vector<rational<long long>> &c,
                     int cur, rational<long long> value) {
  const node &n = nodes[cur];
  c[cur] = value;
  for (int y : n.ch) {
    rational<long long> x = rational<long long>(1, 2) / nodes[y].lambda;
    pushdown(nodes, c, y, value * x);
  }
}
int main(int argc, const char **argv) {
  redirect(argc, argv);
  int T = 0;
  cin >> T;
  for (int cas = 0; cas < T; cas++) {
    int N, M;
    cin >> N >> M;
    vector<vector<int>> edges(N);
    for (int i = 0; i < M; i++) {
      int a, b;
      cin >> a >> b;
      a--;
      b--;
      edges[a].push_back(b);
      edges[b].push_back(a);
    }
    vector<node> nodes(N);
    pair<long long, long long> base{-1, -1};
    for (int i = 0; i < N; i++) {
      if (!nodes[i].seen) {
        base = recurse(nodes, edges, i, -1);
        if (base.first != -1) break;
      }
    }
    bool win = false;
    vector<long long> c(N);
    if (base.second != -1) {
      win = true;
      int cur = base.first;
      while (cur != base.second) {
        c[cur] = 1;
        cur = nodes[cur].parent;
      }
      c[cur] = 1;
    } else if (base.first != -1) {
      win = true;
      vector<rational<long long>> cr(N);
      pushdown(nodes, cr, base.first, 1);
      long long lcm = 1;
      for (int i = 0; i < N; i++) {
        long long d = cr[i].denominator();
        if (d != 1) lcm = lcm * (d / gcd(lcm, d));
      }
      for (int i = 0; i < N; i++)
        c[i] = cr[i].numerator() * (lcm / cr[i].denominator());
    }
    if (win) {
      cout << "YES\n";
      long long produce = 0;
      for (int i = 0; i < N; i++) {
        cout << c[i] << ' ';
        produce -= c[i] * c[i];
        for (int trg : edges[i])
          if (trg > i) produce += c[i] * c[trg];
      }
      cout << '\n';
    } else
      cout << "NO\n";
  }
  return 0;
}