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;
const int maxn = 200005;
inline pair<int, int> operator-(const pair<int, int> &p,
                                const pair<int, int> &q) {
  return pair<int, int>(p.first - q.first, p.second - q.second);
}
int pmy[4 * maxn], pMy[4 * maxn], smy[4 * maxn], sMy[4 * maxn];
pair<int, int> p[maxn];
int sz;
vector<pair<int, int> > lo, up;
vector<pair<int, int> > qx;
inline long long cross(const pair<int, int> &p, const pair<int, int> &q) {
  return 1LL * p.first * q.second - 1LL * p.second * q.first;
}
inline int hull(vector<pair<int, int> > &h) {
  if (sz <= 1) return sz;
  h.clear();
  h.push_back(p[0]);
  for (int i = 1; i < sz; i++) {
    while (h.size() > 1 &&
           cross(h.back() - h[(int)h.size() - 2], p[i] - h.back()) <= 0)
      h.pop_back();
    h.push_back(p[i]);
  }
  return h.size();
}
inline bool fnd(const pair<int, int> &p, const vector<pair<int, int> > &v) {
  int i = lower_bound((v).begin(), (v).end(), p) - v.begin();
  return i < int(v.size()) && v[i] == p;
}
inline int convex_hull() {
  lo.clear();
  up.clear();
  if (sz <= 1) return sz;
  sort(p, p + sz);
  sz = unique(p, p + sz) - p;
  int L = hull(lo);
  reverse(p, p + sz);
  int H = hull(up);
  sz = 0;
  for (int i = 0; i < int(lo.size()); i++) p[sz++] = lo[i];
  for (int i = 1; i < int(up.size()) - 1; i++) p[sz++] = up[i];
  return L + H - 2;
}
bool in(pair<int, int> q) {
  int it =
      lower_bound(qx.begin(), qx.end(), make_pair(q.first, 0)) - qx.begin();
  if (it == qx.size() || (pMy[it] <= q.second || pmy[it] >= q.second))
    return false;
  it = lower_bound(qx.begin(), qx.end(), make_pair(q.first + 1, 0)) -
       qx.begin() - 1;
  if (it < 0 || (sMy[it] <= q.second || smy[it] >= q.second)) return false;
  return true;
}
int main() {
  ios::sync_with_stdio(false);
  int n, m;
  while (cin >> n >> m && n) {
    map<pair<int, int>, int> cnt;
    int y, x;
    qx.clear();
    for (int i = 0; i < m; i++) {
      cin >> x >> y;
      x--, y--;
      cnt[make_pair(x, y)]++;
      cnt[make_pair(x + 1, y)]++;
      cnt[make_pair(x, y + 1)]++;
      cnt[make_pair(x + 1, y + 1)]++;
      qx.push_back(make_pair(x, y));
      qx.push_back(make_pair(x + 1, y));
      qx.push_back(make_pair(x, y + 1));
      qx.push_back(make_pair(x + 1, y + 1));
    }
    sort(qx.begin(), qx.end());
    qx.resize(unique(qx.begin(), qx.end()) - qx.begin());
    pmy[0] = n, pMy[0] = 0;
    for (int i = 1; i <= qx.size(); i++) {
      pmy[i] = min(pmy[i - 1], qx[i - 1].second);
      pMy[i] = max(pMy[i - 1], qx[i - 1].second);
    }
    smy[qx.size() - 1] = n, sMy[qx.size() - 1] = 0;
    for (int i = int(qx.size()) - 2; i >= 0; i--) {
      smy[i] = min(smy[i + 1], qx[i + 1].second);
      sMy[i] = max(sMy[i + 1], qx[i + 1].second);
    }
    sz = 0;
    for (auto q : cnt) {
      if (q.second == 4 || (q.second == 3 && in(q.first))) {
        p[sz++] = q.first;
      }
    }
    convex_hull();
    cout << sz << endl;
    int st = 0;
    for (int i = 0; i < sz; i++)
      if (p[i] < p[st]) st = i;
    int j = st;
    do {
      cout << p[j].first << ' ' << p[j].second << '\n';
      j--;
      if (j < 0) j = sz - 1;
    } while (j != st);
  }
  return 0;
}