Problem: I'm the Map, I'm the Map! I'm the MAP!!!

Map

In anticipation of new adventures Boots wanted to do a good deed. After discussion with the Map and Backpack, they decided to gift Dora a connected graph. After a long search, Boots chose t graph's variants, which Dora might like. However fox Swiper wants to spoil his plan.

The Swiper knows, that Dora now is only able to count up to 3, so he has came up with a following idea. He wants to steal some non-empty set of vertices, so that the Dora won't notice the loss. He has decided to steal some non-empty set of vertices, so that after deletion of the stolen vertices and edges adjacent to them, every remaining vertex wouldn't change it's degree modulo 3. The degree of a vertex is the number of edges it is adjacent to. It would've been suspicious to steal all the vertices, so Swiper needs another plan.

Boots are sure, that the crime can not be allowed. However they are afraid, that they won't be able to handle this alone. So Boots decided to ask for your help. Please determine for every graph's variant whether the Swiper can perform the theft or not.

Input

The first line contains a single integer t (1 ≤ t ≤ 100 000) — the number of graph variants.

The first line of each variant contains integers n, m (1 ≤ n ≤ 500 000, 0 ≤ m ≤ 500 000), the number of vertexes and edges in the graph.

Then m lines follow, each containing integers a_i, b_i (1 ≤ a_i, b_i ≤ n), the indices of the vertices connected with a corresponding edge.

It's guaranteed, that the graph is connected and doesn't contain multiple edges or self-loops.

It's guaranteed, that the sum of n over all variants is at most 500 000 and that the sum of m over all variants is at most 500 000.

Descriptions of graph's variants are separated with an empty line.

Output

For each variant:

  * In case the answer exists, print "Yes" and then the answer itself.

The first line should contain an integer c (1 < c < n), the number of vertices the Crook can steal, without Dora noticing the loss. On the next line print c distinct integers, the indices of the graph's vertices in arbitrary order.

  * Otherwise print "No". 



In case there are several correct ways to steal the vertices, print any of them.

Please note, that it's not required to maximize the number of stolen vertices.

Example

Input


3
3 3
1 2
2 3
3 1

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

8 12
1 2
1 3
2 3
1 4
4 5
5 1
3 6
3 7
3 8
6 1
7 1
8 1


Output


No
Yes
3
4 5 6
Yes
3
6 7 8

Note

The picture below shows the third variant from the example test. The set of the vertices the Crook can steal is denoted with bold. 

<image>

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
using namespace std;
const int N = 2e6 + 5;
int Case;
namespace G {
int n, m;
int hed[N], to[N], nxt[N], cnt;
int deg[N], u[N], v[N], d[N], vis[N], par[N], del[N];
bool vis2[N];
inline void adde(int u, int v) {
  del[cnt + 1] = 0;
  ++cnt;
  to[cnt] = v, nxt[cnt] = hed[u];
  hed[u] = cnt;
}
inline void read() {
  cnt = 1;
  cin >> n >> m;
  for (int i = 1; i <= m; i++) {
    scanf("%d%d", &u[i], &v[i]);
    deg[u[i]]++, deg[v[i]]++;
  }
}
inline void output(vector<int> del) {
  vector<bool> ans(n + 1, 0);
  if ((int)del.size() == n) {
    puts("No");
  } else {
    puts("Yes");
    printf("%d\n", n - (int)del.size());
    for (size_t i = 0; i < del.size(); i++) ans[del[i]] = 1;
    for (int i = 1; i <= n; i++) {
      if (!ans[i]) printf("%d ", i);
    }
    puts("");
  }
}
vector<int> ret;
int rt, Rt[N], dep[N];
inline void dfs(int x, int pre) {
  dep[x] = dep[pre] + 1;
  Rt[x] = rt;
  par[x] = pre;
  vis[x] = 1;
  int T = 0;
  for (int i = hed[x]; i; i = nxt[i]) {
    int v = to[i];
    if (v == pre) continue;
    if (vis[v]) {
      if (dep[v] < dep[x]) {
        if (!ret.size()) {
          if (!T)
            T = v;
          else if (dep[v] > dep[T])
            T = v;
        }
      }
    }
  }
  if (T) {
    int cur = x;
    ret.push_back(cur);
    while (cur != T) {
      cur = par[cur];
      ret.push_back(cur);
    }
  }
  for (int i = hed[x]; i; i = nxt[i]) {
    int v = to[i];
    if (!vis[v]) {
      dfs(v, x);
    }
  }
}
vector<int> getpath(int x, int y) {
  bool tag = 0;
  if (dep[x] < dep[y]) swap(x, y), tag = 1;
  vector<int> a, b;
  while (dep[x] > dep[y]) {
    a.push_back(x);
    x = par[x];
  }
  while (x != y) {
    a.push_back(x), b.push_back(y);
    x = par[x], y = par[y];
  }
  a.push_back(x);
  reverse(b.begin(), b.end());
  for (size_t i = 0; i < b.size(); i++) a.push_back(b[i]);
  if (tag) reverse(a.begin(), a.end());
  return a;
}
bool sx[N], sy[N];
void Answer(int x, vector<int> path, int y) {
  for (int i = 1; i <= n; i++) sx[i] = sy[i] = 0;
  for (int i = 1; i <= m; i++) {
    if (u[i] == x) sx[v[i]] = 1;
    if (u[i] == y) sy[v[i]] = 1;
  }
  for (int i = 0; i < (int)path.size(); i++)
    if (sy[path[i]]) {
      path.resize(i + 1);
      break;
    }
  vector<int> s;
  for (int i = (int)path.size() - 1; ~i; i--) {
    s.push_back(path[i]);
    if (sx[path[i]]) break;
  }
  s.push_back(x), s.push_back(y);
  if (s.size() > n) cout << "A" << endl;
  output(s);
}
vector<pair<int, int> > s[N];
vector<int> RRR;
int delnode[N];
inline int DF(int x, int pre) {
  if (sx[x]) {
    RRR.push_back(x);
    return x;
  }
  int T = 0;
  for (int i = hed[x]; i; i = nxt[i])
    if (to[i] != pre && !delnode[to[i]]) {
      int ret = DF(to[i], x);
      if (ret > 0) {
        T = 1;
        break;
      }
    }
  if (T) RRR.push_back(x);
  return T;
}
void Sol(int nd1, int cnt0) {
  for (int i = 1; i <= n; i++) sx[i] = 0;
  for (int i = 1; i <= m; i++)
    if (u[i] == cnt0) sx[v[i]] = 1;
  RRR.push_back(nd1);
  int R = 2;
  if (sx[nd1]) R--;
  for (int i = hed[nd1]; i; i = nxt[i]) {
    int v = to[i];
    if (DF(v, nd1)) R--;
    if (!R) break;
  }
}
int Sz[N], _rt;
void _ddd(int x, int pre) {
  Sz[x] = sx[x];
  Rt[x] = _rt;
  for (int i = hed[x]; i; i = nxt[i]) {
    int v = to[i];
    if (v == pre) continue;
    _ddd(v, x);
    Sz[x] += Sz[v];
  }
}
void Del(int x, int &b) {
  delnode[x] = 1;
  Rt[x] = 0;
  for (int i = hed[x]; i; i = nxt[i]) {
    int v = to[i];
    _rt = 0;
    _ddd(v, x);
    if (Sz[v] > 1) {
      b++, _rt = v;
      _ddd(v, x);
    }
  }
}
inline void solve() {
  int BLOCKS = 0;
  for (int i = 1; i <= n; i++) {
    if (deg[i] % 3 == 0) {
      output(vector<int>(1, i));
      return;
    }
  }
  for (int i = 1; i <= m; i++) {
    if (deg[u[i]] % 3 == 1 && deg[v[i]] % 3 == 1) {
      vector<int> r(2, 0);
      r[0] = u[i], r[1] = v[i];
      output(r);
      return;
    }
    if (deg[u[i]] % 3 == 2 && deg[v[i]] % 3 == 2) {
      adde(u[i], v[i]), adde(v[i], u[i]);
      d[u[i]]++, d[v[i]]++;
    }
  }
  ret.clear();
  for (int i = 1; i <= n; i++)
    if (deg[i] % 3 == 2 && !vis[i]) {
      BLOCKS++;
      rt = i;
      dfs(i, 0);
    }
  if (ret.size()) {
    if (ret.size() > n) cout << "RET";
    output(ret);
    return;
  }
  vector<int> cnt1;
  for (int i = 1; i <= n; i++)
    if (deg[i] % 3 == 1) cnt1.push_back(i);
  if ((int)cnt1.size() > 1) {
    for (int i = 1; i <= m; i++) {
      if (deg[u[i]] % 3 == 2 && deg[v[i]] % 3 == 1) swap(u[i], v[i]);
      if (deg[u[i]] % 3 == 1 && deg[v[i]] % 3 == 2) {
        int r = Rt[v[i]];
        s[r].push_back(pair<int, int>(u[i], v[i]));
      }
    }
    for (int i = 1; i <= n; i++)
      if (s[i].size()) {
        sort(s[i].begin(), s[i].end());
        for (int j = 1; j < (int)s[i].size(); j++) {
          if (s[i][j].first != s[i][j - 1].first) {
            vector<int> ret = getpath(s[i][j].second, s[i][j - 1].second);
            Answer(s[i][j].first, ret, s[i][j - 1].first);
            return;
          }
        }
      }
  } else {
    RRR.clear();
    for (int i = 1; i <= m; i++)
      if (deg[u[i]] % 3 == 2 && deg[v[i]] % 3 == 1) swap(u[i], v[i]);
    for (int i = 1; i <= n; i++) sx[i] = 0, delnode[i] = 0, vis2[i] = 0;
    for (int i = 1; i <= m; i++)
      if (u[i] == cnt1[0]) sx[v[i]] = 1;
    if (BLOCKS == 1) {
      for (int i = 1; i <= n; i++)
        if (deg[i] % 3 == 2) {
          if (deg[i] >= 5) {
            Del(i, BLOCKS);
            break;
          }
        }
    }
    if (BLOCKS > 1) {
      RRR.push_back(cnt1[0]);
      for (int i = 1; i <= n; i++) vis2[i] = 0;
      int T = 2;
      for (int i = 1; i <= n; i++)
        if (deg[i] % 3 == 2) {
          if (!vis2[Rt[i]]) {
            vis2[Rt[i]] = 1, Sol(i, cnt1[0]), T--;
          }
          if (T == 0) {
            break;
          }
        }
      if (RRR.size() > n) cout << "RRR" << endl;
      output(RRR);
      return;
    }
    return;
  }
}
void clear() {
  for (int i = 1; i <= n; i++)
    dep[i] = Rt[i] = par[i] = vis[i] = hed[i] = 0, d[i] = deg[i] = 0;
  cnt = 1;
  for (int i = 1; i <= n; i++) s[i].clear();
}
}  // namespace G
int main() {
  int T;
  cin >> T;
  while (T--) {
    Case++;
    G::read();
    G::solve();
    G::clear();
  }
}