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>
#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")
using namespace std;
using ll = long long;
using ull = unsigned long long;
using pii = pair<int, int>;
using pll = pair<long long, long long>;
constexpr char ln = '\n';
template <class T1, class T2>
inline bool chmax(T1 &a, T2 b) {
  if (a < b) {
    a = b;
    return true;
  }
  return false;
}
template <class T1, class T2>
inline bool chmin(T1 &a, T2 b) {
  if (a > b) {
    a = b;
    return true;
  }
  return false;
}
inline int topbit(int x) { return x == 0 ? -1 : 31 - __builtin_clz(x); }
inline int topbit(long long x) { return x == 0 ? -1 : 63 - __builtin_clzll(x); }
inline int botbit(int x) { return x == 0 ? 32 : __builtin_ctz(x); }
inline int botbit(long long x) { return x == 0 ? 64 : __builtin_ctzll(x); }
inline int popcount(int x) { return __builtin_popcount(x); }
inline int popcount(long long x) { return __builtin_popcountll(x); }
inline int kthbit(long long x, int k) { return (x >> k) & 1; }
inline void print() { cout << "\n"; }
template <class T>
inline void print(const vector<T> &v) {
  for (auto itr = v.begin(); itr != v.end(); ++itr) cout << *itr << " ";
  print();
}
template <class T, class... Args>
inline void print(const T &x, const Args &...args) {
  cout << x << " ";
  print(args...);
}
inline void debug() {}
template <class T>
inline void debug(const vector<T> &v) {}
template <class T, class... Args>
inline void debug(const T &x, const Args &...args) {}
const int dx[6] = {-2, -1, 0, 2, 1, 0};
const int dy[6] = {0, -1, -2, 0, 1, 2};
int main() {
  int N;
  cin >> N;
  vector<vector<int>> ans(N, vector<int>(N, -1));
  ans[0][0] = 1;
  ans[N - 1][N - 1] = 0;
  auto query = [](int x, int y, int nx, int ny) {
    cout << "? ";
    if (x <= nx and y <= ny) {
      cout << x + 1 << " " << y + 1 << " " << nx + 1 << " " << ny + 1 << endl;
    } else {
      cout << nx + 1 << " " << ny + 1 << " " << x + 1 << " " << y + 1 << endl;
    }
    int ret;
    cin >> ret;
    return ret;
  };
  queue<pii> que;
  que.emplace(0, 0);
  while (!que.empty()) {
    auto [x, y] = que.front();
    que.pop();
    for (int k = 0; k < (6); ++k) {
      int nx = x + dx[k];
      int ny = y + dy[k];
      if (nx < 0 or nx >= N or ny < 0 or ny >= N) continue;
      if (ans[nx][ny] != -1) continue;
      int ret = query(x, y, nx, ny);
      if (ret)
        ans[nx][ny] = ans[x][y];
      else
        ans[nx][ny] = 1 - ans[x][y];
      que.emplace(nx, ny);
    }
  }
  for (int i = 0; i < (N - 1); ++i) {
    for (int j = 0; j < (N - 1); ++j) {
      if (ans[i][j] != -1) continue;
      if (!que.empty()) continue;
      if (ans[i + 1][j] != ans[i][j + 1]) {
        if (j == 0) {
          int ret = query(i, j, i + 1, j + 2);
          if (ret) {
            ans[i][j] = ans[i + 1][j + 2];
            que.emplace(i, j);
          } else {
            ans[i][j] = 1 - ans[i + 1][j + 2];
            que.emplace(i, j);
          }
        } else {
          int ret = query(i, j - 1, i + 1, j + 1);
          if (ret) {
            ans[i + 1][j + 1] = ans[i][j - 1];
            que.emplace(i + 1, j + 1);
          } else {
            ans[i + 1][j + 1] = 1 - ans[i][j - 1];
            que.emplace(i + 1, j + 1);
          }
        }
      }
    }
  }
  for (int i = 0; i < (N - 2); ++i) {
    for (int j = 0; j < (N - 2); ++j) {
      if (ans[i][j] == -1) continue;
      if (!que.empty()) continue;
      if (ans[i][j] != ans[i + 2][j + 2]) {
        int ret1 = query(i, j, i + 2, j + 1);
        if (ret1) {
          ans[i + 2][j + 1] = ans[i][j];
          ans[i + 1][j] = ans[i + 1][j + 1];
          ans[i][j + 1] = ans[i + 1][j + 1];
          que.emplace(i + 2, j + 1);
          que.emplace(i + 1, j);
          que.emplace(i, j + 1);
          break;
        }
        int ret2 = query(i, j, i + 1, j + 2);
        if (ret2) {
          ans[i + 1][j + 2] = ans[i][j];
          ans[i + 1][j] = ans[i + 1][j + 1];
          ans[i][j + 1] = ans[i + 1][j + 1];
          que.emplace(i + 1, j + 2);
          que.emplace(i + 1, j);
          que.emplace(i, j + 1);
          break;
        }
        int ret3 = query(i + 2, j + 2, i, j + 1);
        if (ret3) {
          ans[i][j + 1] = ans[i + 2][j + 2];
          ans[i + 2][j + 1] = ans[i + 1][j + 1];
          ans[i + 1][j + 2] = ans[i + 1][j + 1];
          que.emplace(i, j + 1);
          que.emplace(i + 2, j + 1);
          que.emplace(i + 1, j + 2);
          break;
        }
        int ret4 = query(i + 2, j + 2, i + 1, j);
        if (ret4) {
          ans[i + 1][j] = ans[i + 2][j + 2];
          ans[i + 2][j + 1] = ans[i + 1][j + 1];
          ans[i + 1][j + 2] = ans[i + 1][j + 1];
          que.emplace(i + 1, j);
          que.emplace(i + 2, j + 1);
          que.emplace(i + 1, j + 2);
          break;
        }
        if (ans[i + 1][j + 1] == ans[i][j]) {
          ans[i + 1][j + 2] = ans[i + 2][j + 2];
          ans[i + 2][j + 1] = ans[i + 2][j + 2];
          que.emplace(i + 1, j + 2);
          que.emplace(i + 2, j + 1);
        } else {
          ans[i][j + 1] = ans[i + 1][j] = ans[i][j];
          que.emplace(i, j + 1);
          que.emplace(i + 1, j);
        }
      }
    }
  }
  while (!que.empty()) {
    auto [x, y] = que.front();
    que.pop();
    for (int k = 0; k < (6); ++k) {
      int nx = x + dx[k];
      int ny = y + dy[k];
      if (nx < 0 or nx >= N or ny < 0 or ny >= N) continue;
      if (ans[nx][ny] != -1) continue;
      int ret = query(x, y, nx, ny);
      if (ret)
        ans[nx][ny] = ans[x][y];
      else
        ans[nx][ny] = 1 - ans[x][y];
      que.emplace(nx, ny);
    }
  }
  cout << "!" << endl;
  for (int i = 0; i < (N); ++i) {
    for (int j = 0; j < (N); ++j) cout << ans[i][j];
    cout << endl;
  }
}