#include <bits/stdc++.h>
using namespace std;

vector<int> Xs;
vector<pair<int, int>> ans;

const int MX = 1000000;

bool ask(int xl, int xr, int yl, int yr)
{
    cout << "? " << xl << " " << xr << " " << yl << " " << yr << endl;
    int ret;
    cin >> ret;
    assert(ret == 0 || ret == 1);
    return ret == 1;
}

void dfs(int xl, int xr, int yl, int yr)
{
    if(xl + 1 == xr && yl + 1 == yr)
    {
        ans.emplace_back(xl, yl);
        return;
    }
    if(xr - xl >= yr - yl)
    {
        int m = (xl + xr)/2;
        if(ask(xl, m - 1, yl, yr - 1)) dfs(xl, m, yl, yr);
        if(ask(m, xr - 1, yl, yr - 1)) dfs(m, xr, yl, yr);
    }
    else
    {
        int m = (yl + yr)/2;
        if(ask(xl, xr - 1, yl, m - 1)) dfs(xl, xr, yl, m);
        if(ask(xl, xr - 1, m, yr - 1)) dfs(xl, xr, m, yr);
    }
    return;
}

int main()
{
    dfs(0, MX + 1, 0, MX + 1);
    int N = ans.size();
    cout << "! " << N << endl;
    for(auto [x, y] : ans)
    {
        cout << x << " " << y << endl;
    }
}