[CODE FESTIVAL 2017 Final H]Poor Penguin

题意：在一个$n\times m$的网格上，每个格子是薄冰或冰山（网格外什么都没有），有一片薄冰上站着一只企鹅，对于薄冰$(i,j)$，如果不满足（$(i-1,j),(i 1,j)$都有东西或$(i,j-1),(i,j 1)$都有东西），那么它会消失，并且会发生连锁反应，现在你可以把一些冰山削成薄冰，问最少多少次操作可以使得企鹅掉入水中

先考虑什么时候企鹅所在的薄冰会消失（以下的图片全部来自官方题解）

如果一个格子的右下角没有冰山，那么它最终会消失，对其他方向也是这样

如果能把整个网格用十字分开，使得某两个相对区域中都没有冰山，那么另外两个区域可以被分开考虑，且之后互相独立，这种分割可以递归地进行

所以对于一个包含企鹅的矩形，我们DP出让它独立于其他格子所需的最小操作次数，再枚举删掉企鹅的四个方向的冰山来更新答案即可

设$f_{i,j,k,l}$表示让$(i,j),(k,l)$这个矩形独立的最小操作次数，枚举它里面的一个点$(x,y)$，以它为中心画十字分开原矩形来转移即可

总时间复杂度$O((nm)^3)$，感觉Atcoder评测机挺快的？

#include<stdio.h>#include<string.h>#include<algorithm>using namespace std;void fmin(int&a,int b){	if(b<a)a=b;}int s[41][41],f[41][41][41][41];char str[41];int get(int i,int j,int k,int l){	if(i>k||j>l)return 0;	return s[k][l]-s[i-1][l]-s[k][j-1] s[i-1][j-1];}int main(){	int n,m,i,j,k,l,x,y,sx,sy,ans;	scanf("%d%d",&n,&m);	for(i=1;i<=n;i  ){		scanf("%s",str 1);		for(j=1;j<=m;j  ){			if(str[j]=='P'){				sx=i;				sy=j;			}			s[i][j]=s[i-1][j] s[i][j-1]-s[i-1][j-1] (str[j]=='#');		}	}	memset(f,63,sizeof(f));	ans=f[0][0][0][0];	f[1][1][n][m]=0;	for(i=1;i<=sx;i  ){		for(j=1;j<=sy;j  ){			for(k=n;k>=sx;k--){				for(l=m;l>=sy;l--){					fmin(ans,f[i][j][k][l] min(min(get(i,j,sx,sy),get(i,sy,sx,l)),min(get(sx,j,k,sy),get(sx,sy,k,l))));					for(x=i;x<=k;x  ){						for(y=j;y<=l;y  ){							if(sx<=x&&sy<=y)fmin(f[i][j][x][y],f[i][j][k][l] get(i,y 1,x,l) get(x 1,j,k,y));							if(sx<=x&&y<=sy)fmin(f[i][y][x][l],f[i][j][k][l] get(i,j,x,y-1) get(x 1,y,k,l));							if(x<=sx&&y<=sy)fmin(f[x][y][k][l],f[i][j][k][l] get(i,y,x-1,l) get(x,j,k,y-1));							if(x<=sx&&sy<=y)fmin(f[x][j][k][y],f[i][j][k][l] get(i,j,x-1,y) get(x,y 1,k,l));						}					}				}			}		}