Problem: This is an interactive task.

Scientists are about to invent a new optimization for the Floyd-Warshall algorithm, which will allow it to work in linear time. There is only one part of the optimization still unfinished.

It is well known that the Floyd-Warshall algorithm takes a graph with n nodes and exactly one edge between each pair of nodes. The scientists have this graph, what is more, they have directed each edge in one of the two possible directions.

To optimize the algorithm, exactly m edges are colored in pink color and all the rest are colored in green. You know the direction of all m pink edges, but the direction of green edges is unknown to you. In one query you can ask the scientists about the direction of exactly one green edge, however, you can perform at most 2 ⋅ n such queries.

Your task is to find the node from which every other node can be reached by a path consisting of edges of same color. Be aware that the scientists may have lied that they had fixed the direction of all edges beforehand, so their answers may depend on your queries.

Input

The first line contains two integers n and m (1 ≤ n ≤ 100 000, 0 ≤ m ≤ 100 000) — the number of nodes and the number of pink edges.

The next m lines describe the pink edges, the i-th of these lines contains two integers u_i, v_i (1 ≤ u_i, v_i ≤ n, u_i ≠ v_i) — the start and the end of the i-th pink edge. It is guaranteed, that all unordered pairs (u_i, v_i) are distinct.

Output

When you found the answer, print "!" and the number of the node from which every other node can be reached by a single-colored path.

Interaction

To ask the direction of the green edge between nodes a and b, print "?", a and b in single line, separated by the space characters. 

In answer to this read a single integer, which will be 1 if the edge is directed from a to b and 0 if the edge is directed from b to a.

You can ask at most 2 ⋅ n queries, otherwise you will get Wrong Answer.

After printing a query do not forget to output 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. 



Answer -1 instead of 0 or 1 means that you made an invalid query or exceeded the query limit. Exit immediately after receiving -1 and you will see Wrong answer verdict. Otherwise you can get an arbitrary verdict because your solution will continue to read from a closed stream.

Hacks

Hacks should be formatted as follows.

The first line should contain two integers n and m (1 ≤ n ≤ 300, 0 ≤ m ≤ n ⋅ (n - 1) / 2) — the number of nodes and number of pink edges.

The next m lines should describe the pink edges, the i-th line should contain 2 integers a_i, b_i (1 ≤ a_i, b_i ≤ n, a_i ≠ b_i), which means that there is a pink edge from a_i to b_i. All unordered pairs (a_i, b_i) should be distinct.

The next (n ⋅ (n - 1) / 2 - m) lines should describe the green edges, the i-th line should contain two integers a_i, b_i (1 ≤ a_i, b_i ≤ n, a_i ≠ b_i)), which means that there is a green edge from a_i to b_i. All unordered pairs of (a_i, b_i) should be distinct and should also be different from the pairs for the pink edges.

Example

Input


4 2
1 2
3 4
0
1
1

Output


? 1 3
? 4 2
? 3 2
! 3

Note

In the example above the answer for the query "? 1 3" is 0, so the edge is directed from 3 to 1. The answer for the query "? 4 2" is 1, so the edge is directed from 4 to 2. The answer for the query "? 3 2" is 1, so the edge is directed from 3 to 2. So there are green paths from node 3 to nodes 1 and 2 and there is a pink path from node 3 to node 4.

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
using namespace std;
const int inf = (int)1.01e9;
const long long infll = (long long)1.01e18;
const long double eps = 1e-9;
const long double pi = acos((long double)-1);
mt19937 mrand(chrono::steady_clock::now().time_since_epoch().count());
int rnd(int x) { return mrand() % x; }
void precalc() {}
const int maxn = (int)1e5 + 5;
int n, m;
vector<int> g[maxn], gr[maxn];
bool read() {
  if (scanf("%d%d", &n, &m) < 2) {
    return false;
  }
  for (int i = 0; i < n; i++) {
    g[i].clear();
    gr[i].clear();
  }
  for (int i = 0; i < m; i++) {
    int v, u;
    scanf("%d%d", &v, &u);
    v--;
    u--;
    g[v].push_back(u);
    gr[u].push_back(v);
  }
  return true;
}
int used[maxn];
vector<int> p;
void dfs0(int v) {
  used[v] = true;
  for (int i = 0; i < ((int)(g[v]).size()); i++) {
    int u = g[v][i];
    if (!used[u]) {
      dfs0(u);
    }
  }
  p.push_back(v);
}
int col[maxn];
vector<int> c[maxn];
int deg[maxn];
int pos[maxn];
int cs;
void dfs1(int v, int cc) {
  used[v] = true;
  col[v] = cc;
  c[cc].push_back(v);
  for (int i = 0; i < ((int)(gr[v]).size()); i++) {
    int u = gr[v][i];
    if (!used[u]) {
      dfs1(u, cc);
    }
  }
}
int ask(int v, int u) {
  printf("? %d %d\n", v + 1, u + 1);
  fflush(stdout);
  int res;
  scanf("%d", &res);
  return res;
}
vector<int> q;
void del(int cc) {
  for (int i = 0; i < ((int)(c[cc]).size()); i++) {
    int v = c[cc][i];
    for (int j = 0; j < ((int)(g[v]).size()); j++) {
      int u = g[v][j];
      int cu = col[u];
      if (cu == cc) {
        continue;
      }
      deg[cu]--;
      if (!deg[cu]) {
        q.push_back(cu);
      }
    }
  }
}
int solveComps(int a, int b) {
  while (pos[a] < ((int)(c[a]).size()) && pos[b] < ((int)(c[b]).size())) {
    if (ask(c[a][pos[a]], c[b][pos[b]])) {
      pos[b]++;
    } else {
      pos[a]++;
    }
  }
  if (pos[a] >= ((int)(c[a]).size())) {
    swap(a, b);
  }
  del(b);
  return a;
}
void solve() {
  for (int i = 0; i < n; i++) {
    used[i] = false;
  }
  p.clear();
  for (int v = 0; v < n; v++) {
    if (!used[v]) {
      dfs0(v);
    }
  }
  for (int i = 0; i < n; i++) {
    used[i] = false;
  }
  cs = 0;
  for (int i = n - 1; i >= 0; i--) {
    int v = p[i];
    if (!used[v]) {
      c[cs].clear();
      deg[cs] = 0;
      pos[cs] = 0;
      dfs1(v, cs++);
    }
  }
  for (int v = 0; v < n; v++) {
    int cv = col[v];
    for (int i = 0; i < ((int)(g[v]).size()); i++) {
      int u = g[v][i];
      int cu = col[u];
      if (cv == cu) {
        continue;
      }
      deg[cu]++;
    }
  }
  for (int i = 0; i < cs; i++) {
    if (!deg[i]) {
      q.push_back(i);
    }
  }
  assert(!q.empty());
  int curc = q[0];
  for (int l = 1; l < ((int)(q).size()); l++) {
    int cc = q[l];
    curc = solveComps(curc, cc);
  }
  int res = c[curc][pos[curc]];
  printf("! %d\n", res + 1);
  fflush(stdout);
}
int main() {
  precalc();
  read();
  solve();
  return 0;
}