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;
const double pi = acos(0.0) * 2.0;
const long double eps = 1e-10;
const int step[8][2] = {{-1, 0}, {0, 1}, {1, 0},  {0, -1},
                        {-1, 1}, {1, 1}, {1, -1}, {-1, -1}};
template <class T>
inline T abs1(T a) {
  return a < 0 ? -a : a;
}
template <typename t, typename t1>
t min1(t a, t1 b) {
  return a < b ? a : b;
}
template <typename t, typename... arg>
t min1(t a, arg... arr) {
  return min1(a, min1(arr...));
}
template <typename t, typename t1>
t max1(t a, t1 b) {
  return a > b ? a : b;
}
template <typename t, typename... arg>
t max1(t a, arg... arr) {
  return max1(a, max1(arr...));
}
inline int jud(double a, double b) {
  if (abs(a) < eps && abs(b) < eps)
    return 0;
  else if (abs1(a - b) / max(abs1(a), abs1(b)) < eps)
    return 0;
  if (a < b) return -1;
  return 1;
}
template <typename t>
inline int jud(t a, t b) {
  if (a < b) return -1;
  if (a == b) return 0;
  return 1;
}
template <typename it, typename t1>
inline int find(t1 val, it a, int na, bool f_small = 1, bool f_lb = 1) {
  if (na == 0) return 0;
  int be = 0, en = na - 1;
  if (*a <= *(a + na - 1)) {
    if (f_lb == 0)
      while (be < en) {
        int mid = (be + en + 1) / 2;
        if (jud(*(a + mid), val) != 1)
          be = mid;
        else
          en = mid - 1;
      }
    else
      while (be < en) {
        int mid = (be + en) / 2;
        if (jud(*(a + mid), val) != -1)
          en = mid;
        else
          be = mid + 1;
      }
    if (f_small && jud(*(a + be), val) == 1) be--;
    if (!f_small && jud(*(a + be), val) == -1) be++;
  } else {
    if (f_lb)
      while (be < en) {
        int mid = (be + en + 1) / 2;
        if (jud(*(a + mid), val) != -1)
          be = mid;
        else
          en = mid - 1;
      }
    else
      while (be < en) {
        int mid = (be + en) / 2;
        if (jud(*(a + mid), val) != 1)
          en = mid;
        else
          be = mid + 1;
      }
    if (!f_small && jud(*(a + be), val) == -1) be--;
    if (f_small && jud(*(a + be), val) == 1) be++;
  }
  return be;
}
template <class T>
inline T lowb(T num) {
  return num & (-num);
}
inline int bitnum(unsigned int nValue) { return __builtin_popcount(nValue); }
inline int bitnum(int nValue) { return __builtin_popcount(nValue); }
inline int bitnum(unsigned long long nValue) {
  return __builtin_popcount(nValue) + __builtin_popcount(nValue >> 32);
}
inline int bitnum(long long nValue) {
  return __builtin_popcount(nValue) + __builtin_popcount(nValue >> 32);
}
inline int bitmaxl(unsigned int a) {
  if (a == 0) return 0;
  return 32 - __builtin_clz(a);
}
inline int bitmaxl(int a) {
  if (a == 0) return 0;
  return 32 - __builtin_clz(a);
}
inline int bitmaxl(unsigned long long a) {
  int temp = a >> 32;
  if (temp) return 32 - __builtin_clz(temp) + 32;
  return bitmaxl(int(a));
}
inline int bitmaxl(long long a) {
  int temp = a >> 32;
  if (temp) return 32 - __builtin_clz(temp) + 32;
  return bitmaxl(int(a));
}
long long pow(long long n, long long m, long long mod = 0) {
  if (m < 0) return 0;
  long long ans = 1;
  long long k = n;
  while (m) {
    if (m & 1) {
      ans *= k;
      if (mod) ans %= mod;
    }
    k *= k;
    if (mod) k %= mod;
    m >>= 1;
  }
  return ans;
}
template <class t1, class t2>
inline void add(t1 &a, t2 b, int mod = -1) {
  if (mod == -1) mod = 1000000007;
  a += b;
  while (a >= mod) a -= mod;
  while (a < 0) a += mod;
}
template <class t>
void output1(t arr) {
  for (int i = 0; i < (int)arr.size(); i++) cerr << arr[i] << ' ';
  cerr << endl;
}
template <class t>
void output2(t arr) {
  for (int i = 0; i < (int)arr.size(); i++) output1(arr[i]);
}
const int maxn = 610;
int fa[maxn], cate[maxn];
int CheckOneNode(int no, vector<int> others) {
  printf("? %d\n", 1 + (int)others.size());
  for (int i = 0; i < (int)others.size(); i++) {
    printf("%d ", others[i]);
  }
  printf("%d\n", no);
  fflush(stdout);
  int cnt;
  scanf("%d", &cnt);
  int rcnt = 0;
  if (others.size() > 1 && cnt > 0) {
    printf("? %d\n", (int)others.size());
    for (int i = 0; i < (int)others.size(); i++) {
      printf("%d ", others[i]);
    }
    fflush(stdout);
    scanf("%d", &rcnt);
  }
  return cnt - rcnt;
}
int FindOneEdge(int no, vector<int> others) {
  if (others.size() == 0) {
    return -1;
  }
  if (CheckOneNode(no, vector<int>(others.begin(),
                                   others.begin() + (others.size() + 1) / 2))) {
    if (((int)others.size() + 1) / 2 == 1) {
      return others[0];
    }
    return FindOneEdge(
        no, vector<int>(others.begin(),
                        others.begin() + ((int)others.size() + 1) / 2));
  } else {
    return FindOneEdge(
        no, vector<int>(others.begin() + ((int)others.size() + 1) / 2,
                        others.end()));
  }
}
void CheckCircle(vector<int> &all) {
  int no = *all.rbegin();
  int rno =
      CheckOneNode(no, vector<int>(all.begin(), all.begin() + all.size() - 1));
  if (rno > 0) {
    rno =
        FindOneEdge(no, vector<int>(all.begin(), all.begin() + all.size() - 1));
    for (int j = 0; j < (int)all.size(); j++)
      if (rno == all[j]) {
        vector<int> fa1, fa2;
        while (no != -1) {
          fa1.push_back(no);
          no = fa[no];
        }
        while (rno != -1) {
          fa2.push_back(rno);
          rno = fa[rno];
        }
        for (int i = 0; i < (int)fa1.size(); i++)
          for (int j = 0; j < (int)fa2.size(); j++)
            if (fa1[i] == fa2[j]) {
              printf("N %d\n", (int)i + j + 1);
              for (int k = 0; k < i; k++) {
                printf("%d ", fa1[k]);
              }
              for (int k = j; k >= 0; k--) {
                printf("%d ", fa2[k]);
              }
              puts("");
              fflush(stdout);
              exit(0);
            }
      }
  }
}
int main() {
  ios_base::sync_with_stdio(0);
  int n;
  scanf("%d", &n);
  vector<int> all;
  for (int i = 1; i <= n; i++) all.push_back(i);
  memset(fa, -1, sizeof(fa));
  vector<int> ans[2];
  while (all.size()) {
    ans[0].push_back(*all.rbegin());
    cate[*all.rbegin()] = 0;
    int q[maxn], lq = 0;
    q[lq++] = *all.rbegin();
    all.pop_back();
    for (int i = 0; i < lq; i++) {
      int no = q[i], c = cate[no];
      int nxt;
      while ((nxt = FindOneEdge(no, all)) != -1) {
        ans[1 - c].push_back(nxt);
        cate[nxt] = 1 - c;
        fa[nxt] = no;
        q[lq++] = nxt;
        for (int i = 0; i < (int)all.size(); i++)
          if (all[i] == nxt) {
            all.erase(all.begin() + i);
            break;
          }
        CheckCircle(ans[1 - c]);
      }
    }
  }
  printf("Y %d\n", (int)ans[0].size());
  for (auto no : ans[0]) {
    printf("%d ", no);
  }
  puts("");
  fflush(stdout);
  return 0;
}