棋盘覆盖问题 分治

NKOJ 2307 棋盘覆盖问题

问题描述

在一个2k x 2k(1<=k<=10)个方格组成的棋盘中，恰有一个方格是黑色的，其他方格都是白色。你的任务是要用图示的由三个方格构成的4种不同形态的L型骨牌覆盖所有的白色方格，且任何一个白色方格不能同时被两次或多个次覆盖。

输入格式

第一行一个整数k,表示棋盘的大小为2k x 2k
第二行是两个整数，代表特殊方格所在行号和列号。

输出格式

一个以空格为间隔的2k x 2k数字矩阵，表示一种可行的覆盖方案
同一块骨牌用相同的数字表示，用数字0表示黑色方格
按矩阵从左到右，从上到下的顺序将数字由小到大输出，详情见样例。

样例输入

样例输入1：

1
1 1

样例输入2：

2
2 1

样例输入3：

3
3 2

样例输出

样例输出1：

0 1
1 1

样例输出2：

1 1 2 2
0 1 3 2
4 3 3 5
4 4 5 5

样例输出3：

1 1 2 2 3 3 4 4
1 5 5 2 3 6 6 4
7 0 5 8 9 9 6 10
7 7 8 8 11 9 10 10
12 12 13 11 11 14 15 15
12 16 13 13 14 14 17 15
18 16 16 19 20 17 17 21
18 18 19 19 20 20 21 21

---

直接上正解：

容易想到从小规模的情况出发寻找做法。

对于最小的规模，即k=1时，显然就是用一个L形骨牌把方形补齐。

k=2呢？很容易得到，对于黑色方格所在的2$\times$2的方格，可以用一个L形补齐。对于剩下的三个2$\times$2的方格，可以先在中心处放一个L形骨牌，如果把这个L形骨牌看成黑色的，那么问题就转换为了前面的形式。

对于k=3有类似的操作：

这是原来的样子：

这是我们的操作：

这时可以注意到，这样做实际上是用分治的思想把原来的问题分成四个规模更小的问题，这样做一定有解。所以我们得到这道题的核心操作方法：

采用分治的思想，当问题规模为2$\times$2的时候，用L形骨牌填满方格即可；

当问题规模更大时，找到当前格子里的黑色方格，根据它的位置在中心处放一个L形骨牌，并把这个骨牌的颜色视为黑色，之后对四个更小的方格进行相同的操作步骤，直至问题规模为2$\times$2.

分治完之后，重新编号后输出即可。

分治操作过程中，记下当前矩阵的上下左右边界实现比较方便。编号从0~Len-1似乎比1~Len更好，这似乎在很多题，特别是有循环的题中更不容易错，而且写起来比较方便。比如说NOIP2016D1T1，另一道分治题分形之城。

---

代码：

#include<stdio.h>

int K,A,B,Len,Map[1050][1050],Ans[1050][1050],dx[8]={1,-1,0,0,-1,-1,1,1},dy[8]={0,0,1,-1,-1,1,-1,1};

int id,tot;

void FZ(int cur,int x,int y,int u,int d,int l,int r)//cur 阶数  x y 空白格所在方位 u d l r 矩阵四边 
{
    tot++;
    int i,j,t,cmid,rmid,X[4],Y[4],pos=0;

    if(cur==1)
    {

        for(i=0;i<2;i++)
        for(j=0;j<2;j++)if((u+i!=x)||(l+j!=y))Map[u+i][l+j]=tot;
        return;
    }

    cmid=u+d>>1;rmid=l+r>>1;

    if(x>cmid)pos+=2;
    if(y>rmid)pos++;

    X[0]=cmid;Y[0]=rmid;if(pos!=0)Map[cmid][rmid]=tot;
    X[1]=cmid;Y[1]=rmid+1;if(pos!=1)Map[cmid][rmid+1]=tot;
    X[2]=cmid+1;Y[2]=rmid;if(pos!=2)Map[cmid+1][rmid]=tot;
    X[3]=cmid+1;Y[3]=rmid+1;if(pos!=3)Map[cmid+1][rmid+1]=tot;
    X[pos]=x;Y[pos]=y;
    //填中心L骨牌，找到四个小正方形里黑方格的坐标

    FZ(cur-1,X[0],Y[0],u,cmid,l,rmid);//左上角
    FZ(cur-1,X[1],Y[1],u,cmid,rmid+1,r);//右上角
    FZ(cur-1,X[2],Y[2],cmid+1,d,l,rmid);//左下角
    FZ(cur-1,X[3],Y[3],cmid+1,d,rmid+1,r);//右下角
}
/*
    u
l       r
    d
*/ 

/*
pos:
0 1
2 3
*/

int main()
{
    int i,j,k,tx,ty;

    scanf("%d%d%d",&K,&A,&B);

    Len=1<<K;
    FZ(K,A-1,B-1,0,Len-1,0,Len-1);

    for(i=0;i<Len;i++)
    for(j=0;j<Len;j++)
    {
        if((i==A-1&&j==B-1)||Ans[i][j])continue;
        Ans[i][j]=++id;
        for(k=0;k<8;k++)
        {
            tx=dx[k]+i;
            ty=dy[k]+j;
            if(tx>=0&&ty>=0&&tx<Len&&ty<Len&&Map[i][j]==Map[tx][ty])Ans[tx][ty]=Ans[i][j];
        }
    }
    //重编号

    for(i=0;i<Len;i++,putchar('\n'))
    for(j=0;j<Len;j++)printf("%d ",Ans[i][j]);
}