#include<cstdio>
#include<algorithm>
#include<bitset>
const int M=1<<12,mod=12289,g=11,invg=(mod*5+1)/g;typedef long long ll;
int qpow(int x,int y){int r=1;for(;y;x=x*x%mod,y>>=1)if(y&1)r=r*x%mod;return r;}
int L,invL,*rev,frv[M<<1],w[2][M<<1];
void initNTT()
{
	for(L=1;L<=M;L<<=1)
	{
		for(int i=0;i<L;++i)frv[L+i]=frv[L+(i>>1)]>>1|(i&1?L>>1:0);
		int gs=qpow(g,(mod-1)/L),ivgs=qpow(invg,(mod-1)/L);w[0][L]=w[1][L]=1;
		for(int i=1;i<L;++i)w[0][L+i]=w[0][L+i-1]*gs%mod;
		for(int i=1;i<L;++i)w[1][L+i]=w[1][L+i-1]*ivgs%mod;
	}
	
}
void preNTT(int n){L=(1<<n);rev=frv+L;invL=qpow(L,mod-2);}
void NTT(int*a,int g)
{
	for(int i=0;i<L;++i)if(i<rev[i])std::swap(a[i],a[rev[i]]);
	for(int i=1;i<L;i<<=1)
	{
		int*ws=w[g==invg]+i+i;
		for(int j=0;j<L;j+=i+i)
			for(int k=0;k<i;++k)
		{
			int x=a[j+k],y=(ll)a[i+j+k]*ws[k]%mod;
			a[j+k]=x+y;a[i+j+k]=x-y;
		}
	}
	for(int i=0;i<L;++i)a[i]=(a[i]%mod+mod)%mod;
	if(g==invg)for(int i=0;i<L;++i)a[i]=a[i]*invL%mod;
}
typedef std::bitset<(1<<11)> bts;
bts a[1<<11],ans[1<<11],ss[1<<11],dd[1<<11];
const int dx[]={0,-1,1,0,0},dy[]={0,0,0,1,-1};
inline void swp(bts*t,int i,int j){auto s=t[i][j];s=!s;}
int tma[2][M],tmb[1<<11][M],tmc[M];

void getdd(int n)
{
	if(n==1){dd[0].set(0);return;}
	getdd(n-1);int LS=(1<<(n-1))-1,AL=(1<<n)-1;for(int i=0;i<=LS;++i)ss[i]=dd[i],dd[i].reset();
	for(int i=LS;~i;--i)for(int j=LS;~j;--j)if(ss[i].test(j)){for(int d=0;d<5;++d)swp(dd,((i<<1)+dx[d])&AL,((j<<1)+dy[d])&AL);}
	
//	for(int i=0;i<=AL;++i,puts(""))
//		for(int j=0;j<=AL;++j)
//			printf("%d ",(int)dd[i][j]);
}
char s[1<<14];
#define SIZE (1<<20)
char out[SIZE],*p3=out;
#define flush() (fwrite(out,1,p3-out,stdout),p3=out)
#define putchar(ch) (p3==out+SIZE&&flush(),*p3++=(ch))
class Flush{public:~Flush(){flush();}}_;
void write(int x)
{
    static short Stack[50],top(0);
    do Stack[++top]=x%10,x/=10;while(x);
    while(top)putchar(Stack[top--]|48);
}
int main()
{
//	freopen("set.in","r",stdin);
//	freopen("set.out","w",stdout);
	int n;scanf("%d\n",&n);getdd(n);int AL=(1<<n)-1;
	for(int i=0;i<=AL;++i){fgets(s,1<<14,stdin);for(int j=0;j<(1<<n);++j)if(s[j<<1]&1)a[i].set(j);}
	for(int i=AL;i>1;--i){for(int j=0;j<(1<<n);++j)if(a[i].test(j)){swp(ans,i-1,j);for(int d=0;d<5;++d)swp(a,((i-1)+dx[d])&AL,(j+dy[d])&AL);}}
//	puts("");
//	for(int i=0;i<=AL;++i,puts(""))
//		for(int j=0;j<=AL;++j)
//			printf("%d ",(int)a[i][j]);
	initNTT();preNTT(n+1);
	for(int t=0;t<2;++t){for(int i=0;i<=AL;++i)tma[t][i]=a[t][i];NTT(tma[t],g);}
	for(int t=0;t<(1<<n);++t){for(int i=0;i<=AL;++i)tmb[t][i]=dd[t][i];NTT(tmb[t],g);}
	
	for(int ta=0;ta<2;++ta)for(int tb=0;tb<(1<<n);++tb)
	{
		for(int i=0;i<L;++i)tmc[i]=tma[ta][i]*tmb[tb][i]%mod;
		NTT(tmc,invg);for(int i=0;i<L;++i)if(tmc[i]&1)swp(ans,(ta+tb)&AL,i&AL);
	}
	
	int rescnt=0;for(int i=0;i<(1<<n);++i)rescnt+=ans[i].count();write(rescnt);putchar(10);
	for(int i=0;i<(1<<n);++i)for(int j=0;j<(1<<n);++j)if(ans[i].test(j))write(i),putchar(32),write(j),putchar(10);
	return 0;
}