[JOISC2016]サンドイッチ

题目大意：
　　一个$n\times m(n,m\leq400)$的网格图中，每个格子上放了两个三明治，摆放的方式分为'N'和'Z'两种。一个三明治可以被拿走当且仅当与该三明治的两条直角边相邻的三明治均被拿走或与该三明治斜边相邻的三明治被拿走。问对于每个格子，想要拿走这个格子中的所有三明治至少需要先拿走多少三明治。

思路：
　　对于同一个格子，不难发现这一个格子中两个三明治接连被拿走一定是最优的。
　　于是这题就和每个单独的三明治取走顺序没什么关系了，只和每个方格取走顺序及三明治的摆放方式有关。
　　$O(n^2)$枚举每个格子$(x,y)$，假设它是因为$(x-1,y)$和$(x,y-1)$被取走后才被取走，我们可以$O(n^2)$DFS出最优情况下，取走每个格子之前一定要取走哪些格子。时间复杂度$O(n^4)$，bitset优化为$O(\frac{n^4}{\omega})$。
　　不难发现，若$(x,y)$是因为$(x-1,y)$被取走才被取走的，$(x-1,y)$不可能因为$(x,y)$被取走才被取走。因此对于同一行的格子，我们可以让后面的DFS重复利用前面DFS出的信息。DFS是$O(n^2)$的，每一行要重新DFS，时间复杂度是$O(n^3)$。
　　具体实现上，可以用$0\sim3$来表示不同的方向。若摆放方式为'N'的格子，一个直角边的方向为$d$，则另一个直角边的方向为$d\oplus3$；若摆放方式为'Z'的格子，一个直角边的方向为$d$，则另一个直角边的方向为$d\oplus1$。搜索时的一系列分类讨论可以通过简单的位运算实现。

#include<cstdio>
#include<cctype>
#include<algorithm>
inline int getint() {
    register char ch;
    while(!isdigit(ch=getchar()));
    register int x=ch^'0';
    while(isdigit(ch=getchar())) x=(((x<<2)+x)<<1)+(ch^'0');
    return x;
}
inline bool getblock() {
    register char ch;
    while(!isalpha(ch=getchar()));
    return ch=='N';
}
const int N=400,inf=0x7fffffff;
const int dx[]={0,-1,0,1},dy[]={-1,0,1,0};
bool a[N][N];
int n,m,f[N][N],ans[N][N],tmp,flag;
inline bool check(const int &x,const int &y) {
    return x>=0&&x<n&&y>=0&&y<m;
}
void dfs(const int &x,const int &y,const int &d) {
    if(!~f[x][y]) {
        flag=true;
        return;
    }
    if(f[x][y]) return;
    tmp+=2;
    f[x][y]=-1;
    const int p=a[x][y]?3:1;
    if(check(x-dx[d],y-dy[d])) dfs(x-dx[d],y-dy[d],d);
    if(check(x-dx[d^p],y-dy[d^p])) dfs(x-dx[d^p],y-dy[d^p],d^p);
    f[x][y]=1;
}
int main() {
    for(register int T=getint();T;T--) {
        n=getint(),m=getint();
        for(register int i=0;i<n;i++) {
            for(register int j=0;j<m;j++) {
                a[i][j]=getblock();
            }
        }
        for(register int i=0;i<n;i++) {
            flag=tmp=0;
            for(register int i=0;i<n;i++) {
                for(register int j=0;j<m;j++) {
                    f[i][j]=0;
                }
            }
            for(register int j=0;j<m;j++) {
                if(!flag) dfs(i,j,2);
                ans[i][j]=flag?inf:tmp;
            }
            flag=tmp=0;
            for(register int i=0;i<n;i++) {
                for(register int j=0;j<m;j++) {
                    f[i][j]=0;
                }
            }
            for(register int j=m-1;~j;j--) {
                if(!flag) dfs(i,j,0);
                ans[i][j]=std::min(ans[i][j],flag?inf:tmp);
            }
        }
        for(register int i=0;i<n;i++) {
            for(register int j=0;j<m;j++) {
                printf("%d%c",ans[i][j]!=inf?ans[i][j]:-1," \n"[j==m-1]);
            }
        }
    }
    return 0;
}