Problem: Bob has a simple undirected connected graph (without self-loops and multiple edges). He wants to learn whether his graph is bipartite (that is, you can paint all vertices of the graph into two colors so that there is no edge connecting two vertices of the same color) or not. As he is not very good at programming, he asked Alice for help. He does not want to disclose his graph to Alice, but he agreed that Alice can ask him some questions about the graph.

The only question that Alice can ask is the following: she sends s — a subset of vertices of the original graph. Bob answers with the number of edges that have both endpoints in s. Since he doesn't want Alice to learn too much about the graph, he allows her to ask no more than 20000 questions. Furthermore, he suspects that Alice might introduce false messages to their communication channel, so when Alice finally tells him whether the graph is bipartite or not, she also needs to provide a proof — either the partitions themselves or a cycle of odd length.

Your task is to help Alice to construct the queries, find whether the graph is bipartite. 

Input

The first line contains a single integer n (1 ≤ n ≤ 600) — the number of vertices in Bob's graph.

Interaction

First, read an integer n (1≤ n≤ 600) — the number of vertices in Bob's graph.

To make a query, print two lines. First of which should be in the format "? k" (1 ≤ k ≤ n), where k is the size of the set to be queried. The second line should contain k space separated distinct integers s_1, s_2, ..., s_k (1 ≤ s_i ≤ n) — the vertices of the queried set.

After each query read a single integer m (0 ≤ m ≤ (n(n-1))/(2)) — the number of edges between the vertices of the set \\{s_i\}.

You are not allowed to ask more than 20000 queries.

If m = -1, it means that you asked more queries than allowed, or asked an invalid query. Your program should immediately terminate (for example, by calling exit(0)). You will receive Wrong Answer; it means that you asked more queries than allowed, or asked an invalid query. If you ignore this, you can get other verdicts since your program will continue to read from a closed stream.

After printing a query do not forget to print end of line and flush the output. Otherwise, you will get Idleness limit exceeded. To do this, use:

  * fflush(stdout) or cout.flush() in C++;
  * System.out.flush() in Java;
  * flush(output) in Pascal;
  * stdout.flush() in Python;
  * see documentation for other languages.



When you know the answer, you need to print it.

The format of the answer depends on whether the graph is bipartite or not.

If the graph is bipartite, print two lines. The first should contain the letter "Y" (short for "YES") followed by a space, and then a single integer s (0 ≤ s ≤ n) — the number of vertices in one of the partitions. Second line should contain s integers a_1, a_2, ..., a_s — vertices belonging to the first partition. All a_i must be distinct, and all edges in the main graph must have exactly one endpoint in the set \\{a_i\}.

If the graph is not bipartite, print two lines. The first should contain the letter "N" (short for "NO") followed by a space, and then a single integer l (3 ≤ l ≤ n) — the length of one simple cycle of odd length. Second line should contain l integers c_1, c_2, ..., c_l — the vertices along the cycle. It must hold that for all 1 ≤ i ≤ l, there is an edge \\{c_i, c_{(i mod l)+1}\} in the main graph, and all c_i are distinct.

If there are multiple possible answers, you may print any of them.

Hacks format For hacks, use the following format:

The first line contains two integers n and m~(1 ≤ n ≤ 600, 0 ≤ m ≤ (n(n-1))/(2)) — the number of vertices and edges of the graph, respectively.

Each of the next m lines contains two integers u_i and v_i~(1 ≤ u_i, v_i ≤ n) mean that there is an edge between u_i and v_i. There must not be any multiple edges, no loops, and the graph must be connected.

For example, you need to use this test to get the first sample:
    
    
      
    4 4  
    4 1  
    1 3  
    3 2  
    2 4  
    

Examples

Input

4
4
0
1
1
1
0

Output

? 4 
1 2 3 4
? 2
1 2
? 2
1 3
? 2
1 4
? 2
2 4
? 2
3 4
Y 2
1 2

Input

4
4
3

Output

? 4
1 4 2 3
? 3
1 2 4
N 3
2 1 4

Note

In the first case, Alice learns that there are 4 edges in the whole graph. Over the course of the next three queries, she learns that vertex 1 has two neighbors: 3 and 4. She then learns that while vertex 2 is adjacent to 4, the vertex 3 isn't adjacent to 4. There is only one option for the remaining edge, and that is (2, 3). This means that the graph is a cycle on four vertices, with (1, 2) being one partition and (3, 4) being the second. Here, it would be also valid to output "3 4" on the second line.

In the second case, we also have a graph on four vertices and four edges. In the second query, Alice learns that there are three edges among vertices (1, 2, 4). The only way this could possibly happen is that those form a triangle. As the triangle is not bipartite, Alice can report it as a proof. Notice that she does not learn where the fourth edge is, but she is able to answer Bob correctly anyway.

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
using namespace std;
template <class T, size_t D>
struct vec : vector<vec<T, D - 1>> {
  template <class... Args>
  vec(size_t n = 0, Args... args)
      : vector<vec<T, D - 1>>(n, vec<T, D - 1>(args...)) {}
};
template <class T>
struct vec<T, 1> : vector<T> {
  template <class... Args>
  vec(Args... args) : vector<T>(args...) {}
};
template <class T>
inline bool Minimize(T& a, const T& b) {
  return a > b ? a = b, true : false;
}
template <class T>
inline bool Maximize(T& a, const T& b) {
  return a < b ? a = b, true : false;
}
inline int Next(int i, int n) { return i == n - 1 ? 0 : i + 1; }
inline int Prev(int i, int n) { return !i ? n - 1 : i - 1; }
mt19937 rng(static_cast<uint32_t>(
    chrono::steady_clock::now().time_since_epoch().count()));
int Ask(const vec<int, 1>& vertices) {
  cout << "? " << static_cast<int>((vertices).size()) << endl;
  for (auto& u : vertices) {
    cout << u + 1 << ' ';
  }
  cout << endl;
  int ret;
  cin >> ret;
  return ret;
}
struct DSU {
  int n_sets;
  vec<int, 1> parents, cardinalities;
  DSU(int t_n_sets)
      : n_sets(t_n_sets), parents(n_sets), cardinalities(n_sets, 1) {
    iota((parents).begin(), (parents).end(), 0);
  }
  int Parent(int i) {
    return i == parents[i] ? i : parents[i] = Parent(parents[i]);
  }
  bool SameSet(int i, int j) { return Parent(i) == Parent(j); }
  bool Unite(int i, int j) {
    i = Parent(i);
    j = Parent(j);
    if (i == j) {
      return false;
    }
    if (cardinalities[i] > cardinalities[j]) {
      swap(i, j);
    }
    --n_sets;
    parents[i] = j;
    cardinalities[j] += cardinalities[i];
    return true;
  }
};
int main() {
  ios_base::sync_with_stdio(0);
  cin.tie(0);
  int n_vertices;
  cin >> n_vertices;
  vec<int, 2> neighbors_vertices(n_vertices);
  vec<int, 1> n_edges_prefixes(n_vertices);
  DSU dsu(n_vertices);
  vec<int, 2> elements_sets(n_vertices);
  elements_sets[0].emplace_back(0);
  for (int u = 1; u < n_vertices; ++u) {
    vec<int, 1> sets;
    for (int v = 0; v < u; ++v) {
      if (dsu.Parent(v) == v) {
        sets.emplace_back(v);
      }
    }
    vec<int, 1> prefixes(static_cast<int>((sets).size()));
    for (int i = 0; i < static_cast<int>((sets).size()); ++i) {
      prefixes[i] = n_edges_prefixes[elements_sets[sets[i]].back()] +
                    (i ? prefixes[i - 1] : 0);
    }
    vec<int, 1> connected_sets;
    for (int i = 0; i < static_cast<int>((sets).size()); ++i) {
      int left_index = i, right_index = static_cast<int>((sets).size());
      while (left_index < right_index) {
        int middle_index = (left_index + right_index) >> 1;
        vec<int, 1> vertices;
        vertices.emplace_back(u);
        for (int j = left_index; j <= middle_index; ++j) {
          vertices.insert(vertices.end(), (elements_sets[sets[j]]).begin(),
                          (elements_sets[sets[j]]).end());
        }
        if (Ask(vertices) != prefixes[middle_index] -
                                 (left_index ? prefixes[left_index - 1] : 0)) {
          right_index = middle_index;
        } else {
          left_index = middle_index + 1;
        }
      }
      i = left_index;
      if (i >= static_cast<int>((sets).size())) {
        continue;
      }
      connected_sets.emplace_back(sets[i]);
      left_index = 0,
      right_index = static_cast<int>((elements_sets[sets[i]]).size()) - 1;
      while (left_index < right_index) {
        int middle_index = (left_index + right_index) >> 1;
        vec<int, 1> vertices;
        vertices.emplace_back(u);
        for (int j = 0; j <= middle_index; ++j) {
          vertices.emplace_back(elements_sets[sets[i]][j]);
        }
        if (Ask(vertices) !=
            n_edges_prefixes[elements_sets[sets[i]][middle_index]]) {
          right_index = middle_index;
        } else {
          left_index = middle_index + 1;
        }
      }
      neighbors_vertices[u].emplace_back(elements_sets[sets[i]][left_index]);
      neighbors_vertices[elements_sets[sets[i]][left_index]].emplace_back(u);
    }
    vec<int, 1> new_elements;
    for (auto& connected_set : connected_sets) {
      vec<int, 1> vertices((elements_sets[connected_set]).begin(),
                           (elements_sets[connected_set]).end());
      vertices.emplace_back(u);
      n_edges_prefixes[u] += Ask(vertices);
      dsu.Unite(connected_set, u);
      if (static_cast<int>((new_elements).size())) {
        for (int i = 0;
             i < static_cast<int>((elements_sets[connected_set]).size()); ++i) {
          n_edges_prefixes[elements_sets[connected_set][i]] +=
              n_edges_prefixes[new_elements.back()];
        }
      }
      new_elements.insert(new_elements.end(),
                          (elements_sets[connected_set]).begin(),
                          (elements_sets[connected_set]).end());
    }
    new_elements.emplace_back(u);
    elements_sets[dsu.Parent(u)] = new_elements;
  }
  array<vec<int, 1>, 2> vertices_parities;
  function<void(int, int, bool)> DFS = [&](int u, int parent, bool parity) {
    vertices_parities[parity].emplace_back(u);
    for (auto& v : neighbors_vertices[u]) {
      if (v == parent) {
        continue;
      }
      DFS(v, u, parity ^ 1);
    }
  };
  DFS(0, -1, 0);
  for (int parity = 0; parity < 2; ++parity) {
    if (static_cast<int>((vertices_parities[parity]).size()) &&
        Ask(vertices_parities[parity])) {
      int left_index = 1,
          right_index =
              static_cast<int>((vertices_parities[parity]).size()) - 1;
      while (left_index < right_index) {
        int middle_index = (left_index + right_index) >> 1;
        vec<int, 1> vertices;
        for (int i = 0; i <= middle_index; ++i) {
          vertices.emplace_back(vertices_parities[parity][i]);
        }
        if (Ask(vertices)) {
          right_index = middle_index;
        } else {
          left_index = middle_index + 1;
        }
      }
      int index_1 = left_index;
      left_index = 0, right_index = index_1 - 1;
      while (left_index < right_index) {
        int middle_index = (left_index + right_index + 1) >> 1;
        vec<int, 1> vertices;
        for (int i = middle_index; i <= index_1; ++i) {
          vertices.emplace_back(vertices_parities[parity][i]);
        }
        if (Ask(vertices)) {
          left_index = middle_index;
        } else {
          right_index = middle_index - 1;
        }
      }
      int index_2 = left_index;
      vec<int, 1> parents(n_vertices, -1);
      queue<int> q;
      q.emplace(vertices_parities[parity][index_1]);
      while (static_cast<int>((q).size())) {
        int u = q.front();
        q.pop();
        for (auto& v : neighbors_vertices[u]) {
          if (v == parents[u]) {
            continue;
          }
          parents[v] = u;
          q.emplace(v);
        }
      }
      vec<int, 1> answer;
      for (int u = vertices_parities[parity][index_2]; ~u; u = parents[u]) {
        answer.emplace_back(u);
      }
      cout << "N " << static_cast<int>((answer).size()) << endl;
      for (auto& u : answer) {
        cout << u + 1 << ' ';
      }
      cout << endl;
      return 0;
    }
  }
  cout << "Y " << static_cast<int>((vertices_parities[0]).size()) << endl;
  for (auto& u : vertices_parities[0]) {
    cout << u + 1 << ' ';
  }
  cout << endl;
  return 0;
}