Problem: Zombies have found out about the Zombie Contamination level checker and managed to damage it! Now detecting the shape of their main compound will be a real challenge for Heidi. As before, a lair can be represented as a strictly convex polygon on a lattice. Each vertex of the polygon occupies a point on the lattice. However, the damaged Zombie Contamination level checker can only tell, for each cell, whether the level of Zombie Contamination for that cell is in the set {1, 2, 3}. In other words, Heidi knows all the cells of the lattice for which the Contamination level is not 0 and not 4.

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

Input

The input contains multiple test cases.

The first line of each test case contains two space-separated integers N and M, where N is the size of the lattice grid (5 ≤ N ≤ 100000) and M is the number of lattice points for which the Zombie Contamination level is 1, 2, or 3 (8 ≤ M ≤ 200000).

The second line of each test case contains M pairs of integers x1, y1, ..., xM, yM – coordinates of the cells with Zombie Contamination level not equal to 0 nor 4. It is guaranteed that 1 ≤ xi, yi ≤ N. All pairs xi, yi are different.

Cells are enumerated based on the coordinates of their upper right corner. This means that the bottommost leftmost cell that touches the origin has coordinates (1, 1), and the uppermost leftmost cell is identified as (1, N).

The last line of the file contains two zeroes. This line should not be treated as a test case. The sum of the M values for all tests in one file will not exceed 200000.

Output

For each test case, the following output is expected:

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.

Example

Input

8 19
2 3 2 4 2 5 3 3 3 5 4 3 4 5 4 6 5 2 5 3 5 6 6 2 6 3 6 4 6 5 6 6 6 7 7 6 7 7
5 8
2 2 2 3 2 4 3 2 3 4 4 2 4 3 4 4
0 0


Output

4
2 3
2 4
6 6
5 2
4
2 2
2 3
3 3
3 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 1 and N - 1. 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;
template <class T>
int chkmax(T& a, T b) {
  if (b > a) {
    a = b;
    return 1;
  }
  return 0;
}
template <class T>
int chkmin(T& a, T b) {
  if (b < a) {
    a = b;
    return 1;
  }
  return 0;
}
template <class iterator>
void output(iterator begin, iterator end, ostream& out = cerr) {
  while (begin != end) {
    out << (*begin) << " ";
    begin++;
  }
  out << endl;
}
template <class T>
void output(T x, ostream& out = cerr) {
  output(x.begin(), x.end(), out);
}
void fast_io() {
  ios_base::sync_with_stdio(false);
  cin.tie(nullptr);
  cout.tie(nullptr);
}
struct pt {
  long long x, y;
  pt() : x(0), y(0) {}
  pt(long long _x, long long _y) : x(_x), y(_y) {}
};
pt operator+(pt a, pt b) { return pt(a.x + b.x, a.y + b.y); }
pt operator-(pt a, pt b) { return pt(a.x - b.x, a.y - b.y); }
long long operator*(pt a, pt b) { return a.x * b.x + a.y * b.y; }
long long operator/(pt a, pt b) { return a.x * b.y - a.y * b.x; }
istream& operator>>(istream& is, pt& p) {
  is >> p.x >> p.y;
  return is;
}
ostream& operator<<(ostream& os, const pt& p) {
  os << p.x << " " << p.y;
  return os;
}
bool operator<(pt a, pt b) { return make_pair(a.y, a.x) < make_pair(b.y, b.x); }
bool operator==(pt a, pt b) { return !(a < b) && !(b < a); }
int half_plane(pt a) { return !(a.y > 0 || (a.y == 0 && a.x > 0)); }
bool polar_cmp(pt a, pt b) {
  int ha = half_plane(a), hb = half_plane(b);
  if (ha < hb) {
    return true;
  }
  if (ha > hb) {
    return false;
  }
  return (a / b > 0);
}
pt base;
long long len(pt a) { return a * a; }
bool base_polar_cmp(pt a, pt b) {
  long long cross = (a - base) / (b - base);
  return (cross > 0 || (cross == 0 && len(a - base) < len(b - base)));
}
vector<pt> convex_hull(vector<pt>& a) {
  sort(a.begin(), a.end());
  a.erase(unique(a.begin(), a.end()), a.end());
  base = a[0];
  sort(a.begin() + 1, a.end(), base_polar_cmp);
  vector<pt> st;
  int n = a.size();
  for (int i = 0; i < n; ++i) {
    while (st.size() >= 2) {
      pt last = st.back() - st[st.size() - 2];
      pt next_vect = a[i] - st.back();
      if (last / next_vect <= 0) {
        st.pop_back();
      } else {
        break;
      }
    }
    st.push_back(a[i]);
  }
  return st;
}
const vector<pt> dir = {{0, 0}, {-1, 0}, {0, -1}, {-1, -1}};
const vector<pt> dir2 = {{0, 0}, {1, 0}, {0, 1}, {1, 1}};
int main() {
  fast_io();
  while (true) {
    int n, m;
    cin >> n >> m;
    if (n == 0) return 0;
    vector<pt> a(m);
    for (int i = 0; i < m; i++) cin >> a[i];
    map<int, int> mn, mx;
    for (int i = 0; i < m; i++) {
      int x = a[i].x, y = a[i].y;
      if (mn.find(x) == mn.end())
        mn[x] = y;
      else
        mn[x] = min(mn[x], y);
      if (mx.find(x) == mx.end())
        mx[x] = y;
      else
        mx[x] = max(mx[x], y);
    }
    vector<pt> pts;
    for (const pt& p : a) {
      for (const pt& delta : dir) {
        pt node = p + delta;
        bool is_ok = true;
        for (const pt& delta2 : dir2) {
          pt cell = node + delta2;
          if (mn.find(cell.x) == mn.end()) {
            is_ok = false;
            break;
          }
          if (mn[cell.x] > cell.y || mx[cell.x] < cell.y) {
            is_ok = false;
            break;
          }
        }
        if (is_ok) {
          pts.push_back(node);
        }
      }
    }
    auto ans = convex_hull(pts);
    reverse(ans.begin(), ans.end());
    int pos = min_element(ans.begin(), ans.end(),
                          [](const pt& a, const pt& b) {
                            return make_pair(a.x, a.y) < make_pair(b.x, b.y);
                          }) -
              ans.begin();
    rotate(ans.begin(), ans.begin() + pos, ans.end());
    cout << ans.size() << endl;
    for (const pt& p : ans) {
      cout << p << '\n';
    }
  }
}