Problem: Now that Heidi has made sure her Zombie Contamination level checker works, it's time to strike! This time, the zombie lair is a strictly convex polygon on the lattice. Each vertex of the polygon occupies a point on the lattice. For each cell of the lattice, Heidi knows the level of Zombie Contamination – the number of corners of the cell that are inside or on the border of the lair.

Given this information, Heidi wants to know the exact shape of the lair to rain destruction on the zombies. Help her!

<image>

Input

The input contains multiple test cases.

The first line of each test case contains one integer N, the size of the lattice grid (5 ≤ N ≤ 500). The next N lines each contain N characters, describing the level of Zombie Contamination of each cell in the lattice. Every character of every line is a digit between 0 and 4. 

Cells are given in the same order as they are shown in the picture above: rows go in the decreasing value of y coordinate, and in one row cells go in the order of increasing x coordinate. This means that the first row corresponds to cells with coordinates (1, N), ..., (N, N) and the last row corresponds to cells with coordinates (1, 1), ..., (N, 1).

The last line of the file contains a zero. This line should not be treated as a test case. The sum of the N values for all tests in one file will not exceed 5000.

Output

For each test case, give the following output:

The first line of the output should contain one integer V, the number of vertices of the polygon that is the secret lair. The next V lines each should contain two integers, denoting the vertices of the polygon in the clockwise order, starting from the lexicographically smallest vertex.

Examples

Input

8
00000000
00000110
00012210
01234200
02444200
01223200
00001100
00000000
5
00000
01210
02420
01210
00000
7
0000000
0122100
0134200
0013200
0002200
0001100
0000000
0


Output

4
2 3
2 4
6 6
5 2
4
2 2
2 3
3 3
3 2
3
2 5
4 5
4 2

Note

It is guaranteed that the solution always exists and is unique. It is guaranteed that in the correct solution the coordinates of the polygon vertices are between 2 and N - 2. A vertex (x1, y1) is lexicographically smaller than vertex (x2, y2) if x1 < x2 or <image>.

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
using namespace std;
static const int INF = 0x3f3f3f3f;
static const long long INFL = 0x3f3f3f3f3f3f3f3fLL;
template <typename T, typename U>
static void amin(T& x, U y) {
  if (y < x) x = y;
}
template <typename T, typename U>
static void amax(T& x, U y) {
  if (x < y) x = y;
}
struct P {
  typedef long long T2;
  int x, y;
  P(int x_, int y_) : x(x_), y(y_) {}
  P() : x(0), y(0) {}
};
inline bool operator==(const P& a, const P& b) {
  return a.x == b.x && a.y == b.y;
}
inline bool operator<(const P& a, const P& b) {
  return a.x < b.x || (a.x == b.x && a.y < b.y);
}
inline P operator+(const P& a, const P& b) { return P(a.x + b.x, a.y + b.y); }
inline P operator-(const P& a, const P& b) { return P(a.x - b.x, a.y - b.y); }
inline P operator-=(P& a, const P& b) {
  a.x -= b.x, a.y -= b.y;
  return a;
}
inline P::T2 cross(const P& a, const P& b) {
  return (P::T2)a.x * b.y - (P::T2)a.y * b.x;
}
inline P::T2 dot(const P& a, const P& b) {
  return (P::T2)a.x * b.x + (P::T2)a.y * b.y;
}
inline P::T2 norm(const P& a) { return (P::T2)a.x * a.x + (P::T2)a.y * a.y; }
ostream& operator<<(ostream& out, const P& x) {
  out << "(" << x.x << ", " << x.y << ")";
  return out;
}
inline int ccw(const P& a, const P& b, const P& c) {
  int ax = b.x - a.x, ay = b.y - a.y, bx = c.x - a.x, by = c.y - a.y;
  P::T2 t = (P::T2)ax * by - (P::T2)ay * bx;
  if (t > 0) return 1;
  if (t < 0) return -1;
  if ((P::T2)ax * bx + (P::T2)ay * by < 0) return +2;
  if ((P::T2)ax * ax + (P::T2)ay * ay < (P::T2)bx * bx + (P::T2)by * by)
    return -2;
  return 0;
}
vector<P> convex_hull(vector<P> ps) {
  if ((int)ps.size() <= 1) return ps;
  int n = ps.size(), k = 0;
  sort(ps.begin(), ps.end());
  vector<P> ch(2 * n);
  for (int i = 0; i < n; ch[k++] = ps[i++])
    while (k >= 2 && ccw(ch[k - 2], ch[k - 1], ps[i]) <= 0) --k;
  for (int i = n - 2, t = k + 1; i >= 0; ch[k++] = ps[i--])
    while (k >= t && ccw(ch[k - 2], ch[k - 1], ps[i]) <= 0) --k;
  ch.resize(k - 1);
  return ch;
}
int main() {
  char v[500][501];
  int N;
  while (~scanf("%d", &N) && N != 0) {
    for (int(i) = 0; (i) < (int)(N); ++(i)) scanf("%s", v[i]);
    vector<vector<int> > cnt(N + 1, vector<int>(N + 1));
    for (int(i) = 0; (i) < (int)(N); ++(i))
      for (int(j) = 0; (j) < (int)(N); ++(j)) cnt[i][j] = v[i][j] - '0';
    vector<vector<int> > sum1(N + 1, vector<int>(N + 1));
    for (int(i) = 0; (i) < (int)(N + 1); ++(i))
      for (int j = N - 1; j >= 0; --j) sum1[i][j] = cnt[i][j] - sum1[i][j + 1];
    vector<vector<int> > point(N + 1, vector<int>(N + 1));
    for (int i = N - 1; i >= 0; --i)
      for (int(j) = 0; (j) < (int)(N + 1); ++(j))
        point[i][j] = sum1[i][j] - point[i + 1][j];
    vector<P> points;
    for (int(i) = 0; (i) < (int)(N + 1); ++(i))
      for (int(j) = 0; (j) < (int)(N + 1); ++(j))
        if (point[i][j] > 0) points.emplace_back(j, N - i);
    vector<P> ch = convex_hull(points);
    reverse(ch.begin(), ch.end());
    rotate(ch.begin(), min_element(ch.begin(), ch.end()), ch.end());
    printf("%d\n", (int)ch.size());
    for (int(i) = 0; (i) < (int)(ch.size()); ++(i))
      printf("%d %d\n", ch[i].x, ch[i].y);
  }
  return 0;
}