BZOJ3808 : Neerc2012 Labyrinth of the Minotaur

左上角和右下角不四连通等价于左下角和右上角八连通

枚举正方形的左上角，先二分出最大的边长，使得里面不含障碍物

然后再二分出最小的边长，使得两部分连通，用前缀和判断

 

这题WA了好久…一直对拍都没问题…于是去看原题，发现有SPJ…

然后改了个枚举姿势就过了…

 

#include<cstdio>
#define N 1510
int n,m,i,j,dx[8]={-1,-1,-1,0,0,1,1,1},dy[8]={-1,0,1,-1,1,-1,0,1},x[N*N],y[N*N],h,t,sum[3][N][N],ans,ansx,ansy,l,r,mid;
char s[N];
bool a[N][N],v[N][N];
inline int min(int&a,int b){if(a>b)a=b;}
inline int max(int&a,int b){if(a<b)a=b;}
inline int ask(int p,int x,int y,int k){
  int x1=x-1,y1=y-1,x2=x+k-1,y2=y+k-1;
  max(x1,0),max(y1,0),min(x2,n),min(y2,m);
  return sum[p][x2][y2]-sum[p][x1][y2]-sum[p][x2][y1]+sum[p][x1][y1];
}
inline void add(int nx,int ny){
  if(nx<1||nx>n||ny<1||ny>m||!a[nx][ny]||v[nx][ny])return;
  v[nx][ny]=1,x[++t]=nx,y[t]=ny;
}
inline void bfs(int nx,int ny){for(h=1,t=0,add(nx,ny);h<=t;h++)for(j=0;j<8;j++)add(x[h]+dx[j],y[h]+dy[j]);}
int main(){
  scanf("%d%d",&m,&n);gets(s);
  for(i=1;i<=n;i++)for(gets(s+1),j=1;j<=m;j++)a[i][j]=s[j]=='#';
  a[1][1]=a[n][m]=1;
  for(i=1;i<=n;i++)for(j=1;j<=m;j++)sum[0][i][j]=sum[0][i-1][j]+sum[0][i][j-1]-sum[0][i-1][j-1]+a[i][j];
  a[1][1]=a[n][m]=0;
  for(i=1;i<=n;i++)bfs(i,1);
  for(i=1;i<=m;i++)bfs(n,i);
  for(i=1;i<=n;i++)for(j=1;j<=m;j++)sum[1][i][j]=sum[1][i-1][j]+sum[1][i][j-1]-sum[1][i-1][j-1]+v[i][j],v[i][j]=0;
  for(i=1;i<=n;i++)bfs(i,m);
  for(i=1;i<=m;i++)bfs(1,i);
  for(i=1;i<=n;i++)for(j=1;j<=m;j++)sum[2][i][j]=sum[2][i-1][j]+sum[2][i][j-1]-sum[2][i-1][j-1]+v[i][j];
  ans=N;
  for(i=1;i<=n;i++)for(j=1;j<=m;j++){
    l=1,r=n-i+1,min(r,m-j+1),min(r,ans-1),t=0;
    while(l<=r){
      mid=(l+r)>>1;
      if(!ask(0,i,j,mid))l=(t=mid)+1;else r=mid-1;
    }
    l=1,r=t;
    while(l<=r){
      mid=(l+r)>>1;
      if((ask(1,i-1,j-1,mid+2)||j==1||i+mid-1==n)&&(ask(2,i-1,j-1,mid+2)||i==1||j+mid-1==m))ans=mid,ansx=i,ansy=j,r=mid-1;else l=mid+1;
    }
  }
  if(ans==N)puts("Impossible");else printf("%d %d %d",ans,ansy,ansx);
  return 0;
}