#include<bits/stdc++.h>
namespace {
#pragma GCC diagnostic ignored "-Wunused-function"
#include<atcoder/all>
#pragma GCC diagnostic warning "-Wunused-function"
using namespace std;
using namespace atcoder;
#define rep(i,n) for(int i = 0; i < (int)(n); i++)
#define rrep(i,n) for(int i = (int)(n) - 1; i >= 0; i--)
#define all(x) begin(x), end(x)
#define rall(x) rbegin(x), rend(x)
template<class T> bool chmax(T& a, const T& b) { if (a < b) { a = b; return true; } else return false; }
template<class T> bool chmin(T& a, const T& b) { if (b < a) { a = b; return true; } else return false; }
using ll = long long;
using P = pair<int,int>;
using VI = vector<int>;
using VVI = vector<VI>;
using VL = vector<ll>;
using VVL = vector<VL>;

} int main() {
  ios::sync_with_stdio(false);
  cin.tie(0);
  vector<P> ans;
  auto ask = [](int xl, int xr, int yl, int yr) {
    cout << "? " << xl << ' ' << xr - 1 << ' ' << yl << ' ' << yr - 1 << endl;
    int res;
    cin >> res;
    assert(res >= 0);
    return res;
  };
  auto dfs = [&](auto&& self, int xl, int xr, int yl, int yr) -> void {
    if (!ask(xl, xr, yl, yr)) return;
    if (xl + 1 == xr && yl + 1 == yr) {
      ans.emplace_back(xl, yl);
      return;
    }
    if (xr - xl > yr - yl) {
      int xc = (xl + xr) / 2;
      self(self, xl, xc, yl, yr);
      self(self, xc, xr, yl, yr);
    } else {
      int yc = (yl + yr) / 2;
      self(self, xl, xr, yl, yc);
      self(self, xl, xr, yc, yr);
    }
  };
  dfs(dfs, 0, 1000001, 0, 1000001);
  cout << "! " << ssize(ans) << '\n';
  for (auto [x, y] : ans) cout << x << ' ' << y << '\n';
}