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;
const double PI = 3.141592653589793;
const long long I_MAX = 1LL << 60;
long long mod = 1000000007;
const long long maxN = 2E5;
vector<long long> inv(maxN + 5, 1);
vector<long long> fac(maxN + 5, 1);
void ADD(long long &x, long long y) {
  x += y;
  if (x >= mod) x -= mod;
  if (x < 0) x += mod;
}
void MUL(long long &x, long long y) {
  x *= y;
  x = (x % mod + mod) % mod;
}
long long mypow(long long b, long long e) {
  long long ans = 1;
  while (e) {
    if (e & 1) MUL(ans, b);
    MUL(b, b);
    e >>= 1;
  }
  return ans;
}
long long modinv(long long n) { return mypow(n, mod - 2); }
void calcInv(long long n) {
  for (int i = 2; i <= n; ++i) {
    inv[i] = (mod - mod / i) * inv[mod % i] % mod;
  }
}
void calcFac(long long n) {
  for (int i = 2; i <= n; ++i) {
    fac[i] = fac[i - 1] * i % mod;
  }
}
long long cm(long long a, long long b) {
  long long ans = 1;
  if (a < b) return 0;
  MUL(ans, inv[fac[a - b]]);
  MUL(ans, inv[fac[b]]);
  MUL(ans, fac[a]);
  return ans;
}
long long Lucas(long long n, long long m) {
  if (m == 0) return 1;
  return cm(n % mod, m % mod) * Lucas(n / mod, m / mod) % mod;
}
long long doLucas(long long n, long long m) {
  calcInv(mod);
  calcFac(mod);
  return Lucas(n, m);
}
const long long maxn = 1E5;
void precal() { return; }
void solve() {
  long long n;
  scanf("%I64d", &n);
  int pln[n + 5][n + 5];
  memset(pln, -1, sizeof(pln));
  for (int i = 1; i <= n; ++i)
    for (int j = 1; j <= n; ++j) {
      if ((i & 1) and (j & 1))
        pln[i][j] = 0;
      else if ((i & 1))
        pln[i][j] = 2;
      else if ((j & 1))
        pln[i][j] = 4;
      else
        pln[i][j] = 6;
    }
  pln[1][1] = 1;
  pln[n][n] = 0;
  int rep;
  for (int i = 1; i <= n; i += 2) {
    for (int j = 1; j <= n - 2 - 2 * (i == n); j += 2) {
      printf("? %d %d %d %d\n", i, j, i, j + 2);
      fflush(stdout);
      scanf("%d", &rep);
      if (rep == 1) {
        pln[i][j + 2] = pln[i][j];
      } else {
        pln[i][j + 2] = pln[i][j] ^ 1;
      }
    }
    if (i + 2 <= n) {
      printf("? %d %d %d %d\n", i, 1, i + 2, 1);
      fflush(stdout);
      scanf("%d", &rep);
      if (rep == 1) {
        pln[i + 2][1] = pln[i][1];
      } else {
        pln[i + 2][1] = pln[i][1] ^ 1;
      }
    }
  }
  for (int i = 2; i <= n; i += 2) {
    for (int j = 1; j <= n - 2; j += 2) {
      printf("? %d %d %d %d\n", i, j, i, j + 2);
      fflush(stdout);
      scanf("%d", &rep);
      if (rep == 1) {
        pln[i][j + 2] = pln[i][j];
      } else {
        pln[i][j + 2] = pln[i][j] ^ 1;
      }
    }
    if (i + 2 <= n) {
      printf("? %d %d %d %d\n", i, 1, i + 2, 1);
      fflush(stdout);
      scanf("%d", &rep);
      if (rep == 1) {
        pln[i + 2][1] = pln[i][1];
      } else {
        pln[i + 2][1] = pln[i][1] ^ 1;
      }
    }
  }
  for (int i = 1; i <= n; i += 2) {
    for (int j = 2; j <= n - 2; j += 2) {
      printf("? %d %d %d %d\n", i, j, i, j + 2);
      fflush(stdout);
      scanf("%d", &rep);
      if (rep == 1) {
        pln[i][j + 2] = pln[i][j];
      } else {
        pln[i][j + 2] = pln[i][j] ^ 1;
      }
    }
    if (i + 2 <= n) {
      printf("? %d %d %d %d\n", i, 2, i + 2, 2);
      fflush(stdout);
      scanf("%d", &rep);
      if (rep == 1) {
        pln[i + 2][2] = pln[i][2];
      } else {
        pln[i + 2][2] = pln[i][2] ^ 1;
      }
    }
  }
  for (int i = 2; i <= n; i += 2) {
    for (int j = 2; j <= n - 2; j += 2) {
      printf("? %d %d %d %d\n", i, j, i, j + 2);
      fflush(stdout);
      scanf("%d", &rep);
      if (rep == 1) {
        pln[i][j + 2] = pln[i][j];
      } else {
        pln[i][j + 2] = pln[i][j] ^ 1;
      }
    }
    if (i + 2 <= n) {
      printf("? %d %d %d %d\n", i, 2, i + 2, 2);
      fflush(stdout);
      scanf("%d", &rep);
      if (rep == 1) {
        pln[i + 2][2] = pln[i][2];
      } else {
        pln[i + 2][2] = pln[i][2] ^ 1;
      }
    }
  }
  printf("? %d %d %d %d\n", 1, 1, 2, 2);
  fflush(stdout);
  scanf("%d", &rep);
  if (rep == 1) {
    for (int i = 2; i <= n; i += 2)
      for (int j = 2; j <= n; j += 2) {
        pln[i][j] ^= 7;
      }
  } else {
    for (int i = 2; i <= n; i += 2)
      for (int j = 2; j <= n; j += 2) {
        pln[i][j] ^= 6;
      }
  }
  printf("? %d %d %d %d\n", 1, 2, 2, 3);
  fflush(stdout);
  scanf("%d", &rep);
  if (rep == 1) {
    int x = pln[2][3] ^ 2;
    for (int i = 2; i <= n; i += 2)
      for (int j = 1; j <= n; j += 2) {
        pln[i][j] ^= x;
      }
  } else {
    int x = pln[2][3] ^ 3;
    for (int i = 2; i <= n; i += 2)
      for (int j = 1; j <= n; j += 2) {
        pln[i][j] ^= x;
      }
  }
  int px = n, py = n;
  for (int i = 2; i <= n; ++i) {
    for (int j = 3; j <= n; ++j) {
      if (pln[i - 1][j - 2] == pln[i][j - 1] and pln[i][j - 2] == pln[i][j]) {
        px = i;
        py = j;
      }
      if (pln[i - 1][j - 2] != pln[i][j - 1] and pln[i][j - 2] != pln[i][j]) {
        px = i;
        py = j;
      }
      if (pln[i - 1][j - 2] == pln[i][j - 1] and
          pln[i - 1][j - 1] == pln[i][j]) {
        px = i;
        py = j;
      }
      if (pln[i - 1][j - 2] != pln[i][j - 1] and
          pln[i - 1][j - 1] != pln[i][j]) {
        px = i;
        py = j;
      }
      if (pln[i - 1][j - 2] == pln[i - 1][j] and pln[i][j - 2] == pln[i][j]) {
        px = i;
        py = j;
      }
      if (pln[i - 1][j - 2] != pln[i - 1][j] and pln[i][j - 2] != pln[i][j]) {
        px = i;
        py = j;
      }
    }
  }
  printf("? %d %d %d %d\n", px - 1, py - 2, px, py);
  fflush(stdout);
  scanf("%d", &rep);
  if (px % 2 != py % 2) {
    if (rep == 1) {
      int x = pln[px][py] ^ pln[px - 1][py - 2];
      for (int i = 2; i <= n; i += 2)
        for (int j = 1; j <= n; j += 2) {
          pln[i][j] ^= x;
        }
      for (int i = 1; i <= n; i += 2)
        for (int j = 2; j <= n; j += 2) {
          pln[i][j] ^= x;
        }
    } else {
      int x = pln[px][py] ^ pln[px - 1][py - 2] ^ 1;
      for (int i = 2; i <= n; i += 2)
        for (int j = 1; j <= n; j += 2) {
          pln[i][j] ^= x;
        }
      for (int i = 1; i <= n; i += 2)
        for (int j = 2; j <= n; j += 2) {
          pln[i][j] ^= x;
        }
    }
  } else {
    if (rep == 1) {
      int x = pln[px][py] ^ pln[px - 1][py - 2];
      for (int i = 2; i <= n; i += 2)
        for (int j = 1; j <= n; j += 2) {
          pln[i][j] ^= x;
        }
      for (int i = 1; i <= n; i += 2)
        for (int j = 2; j <= n; j += 2) {
          pln[i][j] ^= x;
        }
    } else {
      int x = pln[px][py] ^ pln[px - 1][py - 2] ^ 1;
      for (int i = 2; i <= n; i += 2)
        for (int j = 1; j <= n; j += 2) {
          pln[i][j] ^= x;
        }
      for (int i = 1; i <= n; i += 2)
        for (int j = 2; j <= n; j += 2) {
          pln[i][j] ^= x;
        }
    }
  }
  printf("!\n");
  for (int i = 1; i <= n; ++i) {
    for (int j = 1; j <= n; ++j) {
      printf("%d", pln[i][j]);
    }
    printf("\n");
  }
  fflush(stdout);
}
int main() {
  precal();
  int t = 1;
  for (int i = 1; i <= t; ++i) {
    solve();
  }
  return 0;
}