Problem: This is an interactive problem

You are given a grid n× n, where n is odd. Rows are enumerated from 1 to n from up to down, columns are enumerated from 1 to n from left to right. Cell, standing on the intersection of row x and column y, is denoted by (x, y).

Every cell contains 0 or 1. It is known that the top-left cell contains 1, and the bottom-right cell contains 0.

We want to know numbers in all cells of the grid. To do so we can ask the following questions: 

"? x_1 y_1 x_2 y_2", where 1 ≤ x_1 ≤ x_2 ≤ n, 1 ≤ y_1 ≤ y_2 ≤ n, and x_1 + y_1 + 2 ≤ x_2 + y_2. In other words, we output two different cells (x_1, y_1), (x_2, y_2) of the grid such that we can get from the first to the second by moving only to the right and down, and they aren't adjacent.

As a response to such question you will be told if there exists a path between (x_1, y_1) and (x_2, y_2), going only to the right or down, numbers in cells of which form a palindrome.

For example, paths, shown in green, are palindromic, so answer for "? 1 1 2 3" and "? 1 2 3 3" would be that there exists such path. However, there is no palindromic path between (1, 1) and (3, 1).

<image>

Determine all cells of the grid by asking not more than n^2 questions. It can be shown that the answer always exists.

Input

The first line contains odd integer (3 ≤ n < 50) — the side of the grid.

Interaction

You begin the interaction by reading n.

To ask a question about cells (x_1, y_1), (x_2, y_2), in a separate line output "? x_1 y_1 x_2 y_2".

Numbers in the query have to satisfy 1 ≤ x_1 ≤ x_2 ≤ n, 1 ≤ y_1 ≤ y_2 ≤ n, and x_1 + y_1 + 2 ≤ x_2 + y_2. Don't forget to 'flush', to get the answer.

In response, you will receive 1, if there exists a path going from (x_1, y_1) to (x_2, y_2) only to the right or down, numbers in cells of which form a palindrome, and 0 otherwise.

In case your query is invalid or you asked more than n^2 queries, program will print -1 and will finish interaction. You will receive Wrong answer verdict. Make sure to exit immediately to avoid getting other verdicts.

When you determine numbers in all cells, output "!".

Then output n lines, the i-th of which is a string of length n, corresponding to numbers in the i-th row of the grid.

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.



Hack Format

To hack, use the following format.

The first line should contain a single odd integer n (side of your grid).

The i-th of n following lines should contain a string of length n corresponding to the i-th row of the grid. Top left element of the grid has to be equal to 1, bottom right has to be equal to 0.

Example

Input


3
0
1
0
1
1
1
1

Output


? 1 1 1 3
? 1 1 2 3
? 2 1 2 3
? 3 1 3 3
? 2 2 3 3
? 1 2 3 2
? 1 2 3 3
!
100
001
000

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
using namespace std;
int ask(int x, int y, int xx, int yy) {
  if ((x + y) > (xx + yy)) {
    swap(x, xx);
    swap(y, yy);
  }
  cout << "? " << (x + 1) << " " << (y + 1) << " " << (xx + 1) << " "
       << (yy + 1) << "\n";
  cout.flush();
  int res;
  cin >> res;
  return res;
}
void print(vector<vector<int>> &ta) {
  cout << "!"
       << "\n";
  int n = ta.size();
  for (int i = 0; i < n; i++) {
    for (int j = 0; j < n; j++) {
      cout << ta[i][j];
    }
    cout << "\n";
  }
  cout.flush();
}
void make_table(vector<vector<int>> &ta2, vector<vector<bool>> &dp2) {
  int n = ta2.size();
  for (int i = 0; i < n * n; i++) {
    dp2[i][i] = true;
  }
  for (int i = 0; i < n; i++) {
    for (int j = 0; j < n; j++) {
      int nx = i + 1, ny = j;
      if (nx < n && ta2[i][j] == ta2[nx][ny]) {
        dp2[n * i + j][nx * n + ny] = true;
      }
      nx = i, ny = j + 1;
      if (ny < n && ta2[i][j] == ta2[nx][ny]) {
        dp2[n * i + j][nx * n + ny] = true;
      }
    }
  }
  for (int len = 3; len < 2 * n; len += 2) {
    for (int i = 0; i < n; i++) {
      for (int j = 0; j < n; j++) {
        for (int k = 0; k <= len; k++) {
          int h = len - k;
          int nx = i + k;
          int ny = j + h;
          if (nx < 0 || ny < 0 || nx >= n || ny >= n) continue;
          if (ta2[i][j] != ta2[nx][ny]) {
            dp2[i * n + j][nx * n + ny] = false;
            continue;
          }
          if ((i + 1) < n && (nx - 1) >= 0) {
            if (dp2[(i + 1) * n + j][(nx - 1) * n + ny])
              dp2[i * n + j][nx * n + ny] = true;
          }
          if ((i + 1) < n && (ny - 1) >= 0) {
            if (dp2[(i + 1) * n + j][nx * n + (ny - 1)])
              dp2[i * n + j][nx * n + ny] = true;
          }
          if ((j + 1) < n && (nx - 1) >= 0) {
            if (dp2[i * n + (j + 1)][(nx - 1) * n + ny])
              dp2[i * n + j][nx * n + ny] = true;
          }
          if ((j + 1) < n && (ny - 1) >= 0) {
            if (dp2[i * n + (j + 1)][nx * n + (ny - 1)])
              dp2[i * n + j][nx * n + ny] = true;
          }
        }
      }
    }
  }
}
int sol(vector<vector<int>> &ta2, vector<vector<int>> &ta3) {
  int n = ta2.size();
  vector<vector<bool>> dp2(n * n, vector<bool>(n * n, false)), dp3 = dp2;
  make_table(ta2, dp2);
  make_table(ta3, dp3);
  for (int i = 0; i < n; i++) {
    for (int j = 0; j < n; j++) {
      for (int nx = 0; nx < n; nx++) {
        for (int ny = 0; ny < n; ny++) {
          if ((abs(i - nx) + abs(j - ny)) < 2) continue;
          if (dp2[i * n + j][nx * n + ny] != dp3[i * n + j][nx * n + ny]) {
            int tt = ask(i, j, nx, ny);
            bool bo = (tt == 1) ? true : false;
            if (bo == dp2[i * n + j][nx * n + ny])
              return 0;
            else
              return 1;
          }
        }
      }
    }
  }
  return 0;
}
int main() {
  int n;
  cin >> n;
  vector<vector<int>> ta(n, vector<int>(n, -1));
  ta[0][0] = 1;
  ta[n - 1][n - 1] = 0;
  stack<pair<int, int>> st;
  st.push(make_pair(0, 0));
  vector<vector<int>> hou = {{0, 2}, {1, 1}, {2, 0}};
  while (st.size()) {
    pair<int, int> p0 = st.top();
    st.pop();
    int cx = p0.first, cy = p0.second;
    for (vector<int> v0 : hou) {
      int nx = cx + v0[0];
      int ny = cy + v0[1];
      if (nx >= n || ny >= n || nx < 0 || ny < 0) continue;
      if (ta[nx][ny] != -1) continue;
      int tt = ask(cx, cy, nx, ny);
      if (tt) {
        ta[nx][ny] = ta[cx][cy];
      } else {
        ta[nx][ny] = ta[cx][cy] ^ 1;
      }
      st.push(make_pair(nx, ny));
    }
  }
  vector<vector<int>> ta2 = ta;
  vector<vector<int>> hou2 = {{0, 2},  {1, 1},   {2, 0},
                              {0, -2}, {-1, -1}, {-2, 0}};
  ta2[1][2] = 0;
  st.push(make_pair(1, 2));
  while (st.size()) {
    pair<int, int> p0 = st.top();
    st.pop();
    int cx = p0.first, cy = p0.second;
    for (vector<int> v0 : hou2) {
      int nx = cx + v0[0];
      int ny = cy + v0[1];
      if (nx >= n || ny >= n || nx < 0 || ny < 0) continue;
      if (ta2[nx][ny] != -1) continue;
      int tt = ask(cx, cy, nx, ny);
      if (tt) {
        ta2[nx][ny] = ta2[cx][cy];
      } else {
        ta2[nx][ny] = ta2[cx][cy] ^ 1;
      }
      st.push(make_pair(nx, ny));
    }
  }
  vector<vector<int>> ta3 = ta2;
  for (int i = 0; i < n; i++) {
    for (int j = 0; j < n; j++) {
      if ((i + j) & 1) ta3[i][j] ^= 1;
    }
  }
  for (int i = 0; i < n; i++) {
    for (int j = 0; j < n; j++) {
    }
  }
  if (sol(ta2, ta3)) {
    print(ta3);
  } else
    print(ta2);
  return 0;
}