#include<bits/stdc++.h>
using namespace std;
const int N=2048;
// bitset<2048*2048>a[2048*2048];
int ch[2048*2048],n;int tt=0;
int pos[N][N];
bitset<N>F[N],b[N];
int G[N][N];
// bool solve()
// {
//     int n=tt;
//     for(int i=1;i<=n;++i)
//     {
//         // cerr<<i<<endl;
//         for(int j=i;j<=n;++j)
//             if(a[j][i]){swap(a[j],a[i]);break;}
//         if(!a[i][i])continue;
//         for(int j=1;j<=n;++j)if(i!=j&&a[j][i])a[j]^=a[i];
//     }
//     for(int i=1;i<=n;++i)if(a[i][tt+1])ch[i]=1;
//     return true;
// }
bitset<N>c[N];
vector<pair<int,int>>to;
int main()
{
    cin.tie(0)->sync_with_stdio(0);
    cin>>n;n=1<<n;
    for(int i=0;i<n;++i)
        for(int j=0;j<n;++j)
        {
            int w;cin>>w;
            b[i][j]=w;
        }
    F[1][1]=1;
    for(int u=4;u<=n;u<<=1)
    {
        // cerr<<"SOLVE: "<<u<<" "<<n<<endl;
        for(int i=0;i<u;++i)
            for(int j=0;j<u;++j)
                G[i][j]=0;
        for(int i=0;i<u/2;++i)
            for(int j=0;j<u/2;++j)
                G[i+u/4][j]^=F[i][j],
                G[i+u/4][j+u/4]^=F[i][j],
                G[i][j+u/4]^=F[i][j],
                G[i+u/2][j+u/4]^=F[i][j],
                G[i+u/4][j+u/2]^=F[i][j];
        int C=0;
        // for(int i=0;i<u;++i,cout<<endl)
        //     for(int j=0;j<u;++j)
        //         cout<<G[i][j],C+=G[i][j];cerr<<C<<endl;
        for(int i=0;i<u;++i)
            for(int j=0;j<u;++j)
                F[i][j]=G[i][j];
    }
    int C=0;
    for(int i=0;i<n;++i)
        for(int j=0;j<n;++j)
            if(F[i][j])to.emplace_back(i-n/2+n&n-1,j-n/2+n&n-1);
    for(int i=1;i<n-1;++i)
    {
        for(int j=0;j<n;++j)
            if(b[i-1][j])
            {
                c[i].flip(j),b[i].flip(j),b[i-1].flip(j),b[i].flip(j+n-1&n-1),b[i].flip(j+1&n-1),b[i+1].flip(j);

    // cerr<<"__\n";
    // for(int i=0;i<n;++i,cout<<endl)
    //     for(int j=0;j<n;++j)
    //         cout<<b[i][j];
            }
    }
    cerr<<clock()<<endl;
    for(auto k:to)
    {
        {
            const auto&A=b[n-2];
            auto&B=c[n-2+k.first&n-1];
            B^=(A<<k.second)|(A>>n-k.second);
        }
        {
            const auto&A=b[n-1];
            auto&B=c[n-1+k.first&n-1];
            B^=(A<<k.second)|(A>>n-k.second);
        }
        // for(int i=n-2;i<n;++i)
        //     for(int j=0;j<n;++j)if(b[i][j])
        //         c[i+k.first&n-1].flip(j+k.second&n-1);

    }
    cerr<<clock()<<endl;
    for(int i=0;i<n;++i)
        for(int j=0;j<n;++j)
            C+=c[i][j];
    cout<<C<<'\n';
    for(int i=0;i<n;++i)
        for(int j=0;j<n;++j)
            if(c[i][j])cout<<i<<" "<<j<<'\n';
    return 0;
}