http://icpc.upc.edu.cn/problem.php?cid=1418&pid=9
#include<cstdio>
using namespace std;
int k,x,y,j,kk,i,a[2000][2000];
void dfs(int x0,int y0,int x,int y,int t) {
if (t==0) return;
int k=t;
if (x0+k<=x&&y0+k<=y) {
a[k+x0-1][y0+k-1]=a[k+x0][k+y0-1]=a[k+x0-1][k+y0]=1;
dfs(x0,y0,k+x0-1,k+y0-1,k/2);
dfs(x0,y0+k,k+x0-1,k+y0,k/2);
dfs(x0+k,y0,k+x0,k+y0-1,k/2);
dfs(x0+k,y0+k,x,y,k/2);
} else if (x>=k+x0&&y<k+y0) {
a[k+x0-1][y0+k-1]=a[k+x0][k+y0]=a[k+x0-1][k+y0]=2;
dfs(x0,y0,k+x0-1,k+y0-1,k/2);
dfs(x0,y0+k,k+x0-1,k+y0,k/2);
dfs(x0+k,y0+k,k+x0,k+y0,k/2);
dfs(x0+k,y0,x,y,k/2);
} else if (x<k+x0&&y>=k+y0) {
a[k+x0-1][y0+k-1]=a[k+x0][k+y0]=a[k+x0][k+y0-1]=3;
dfs(x0,y0,k+x0-1,k+y0-1,k/2);
dfs(x0+k,y0,k+x0,k+y0-1,k/2);
dfs(x0+k,y0+k,k+x0,k+y0,k/2);
dfs(x0,y0+k,x,y,k/2);
} else {
a[k+x0][y0+k]=a[k+x0-1][k+y0]=a[k+x0][k+y0-1]=4;
dfs(x0,y0,x,y,k/2);
dfs(x0+k,y0,k+x0,k+y0-1,k/2);
dfs(x0,y0+k,k+x0-1,k+y0,k/2);
dfs(x0+k,y0+k,k+x0,y0+k,k/2);
}
}
int main() {
scanf("%d%d%d",&k,&x,&y);
a[x][y]=7;
kk=1;
for (i=1; i<=k; i++) {
kk=kk*2;
}
dfs(1,1,x,y,kk);
for (i=1; i<=kk; i++) {
for (j=1; j<=kk; j++) {
if (j==kk) printf("%d \n",a[i][j]);
else printf("%d ",a[i][j]);
}
}
}