设船最宽行列的交点为船的重心,那么只要预处理出重心在每个位置是否可行,以及在边界上走出边界所需的最小值之后,进行一遍BFS即可。
枚举每个点$(x,y)$,求出它上下最近的障碍物的距离。考虑重心在第$x$行左右移动的时候,不能走的部分是一个区间,可以$O(1)$计算。
对于左右障碍物也同理,因此可以通过维护差分前缀和可以做到$O(n^2)$预处理出所有可行位置。
对于边界,首先一定可以通过一直往一个方向走得到一个解;其次对于四个角落上的点,可以通过$O(n)$枚举它卡住船的哪个角落来得到走出地图的最优解。
总时间复杂度为$O(n^2)$。
#include<cstdio>
const int N=2010,inf=100000000;
int n,i,j,x,y,z,xl=inf,xr,yl=inf,yr,xm,ym,d,mu,md,l,r;
int f[N][N],h,t,q[N*N][2],ans=inf;
int ul[N],ur[N],dl[N],dr[N],fu[N],fd[N],f11,f1n,fn1,fnn;
char a[N][N],ban[N][N];
inline void umin(int&a,int b){if(a>b)a=b;}
inline void umax(int&a,int b){if(a<b)a=b;}
inline void add(int x,int y,int z){
if(x<1||x>n||y<1||y>n)return;
if(ban[x][y]||~f[x][y])return;
if(x==1)umin(ans,z+xr-xm+1);
if(x==n)umin(ans,z+xm-xl+1);
if(y==1)umin(ans,z+yr-ym+1);
if(y==n)umin(ans,z+ym-yl+1);
if(x==1&&y==1)umin(ans,z+f11);
if(x==1&&y==n)umin(ans,z+f1n);
if(x==n&&y==1)umin(ans,z+fn1);
if(x==n&&y==n)umin(ans,z+fnn);
q[++t][0]=x,q[t][1]=y,f[x][y]=z;
}
int main(){
scanf("%d",&n);
for(i=1;i<=n;i++)for(scanf("%s",a[i]+1),j=1;j<=n;j++)if(a[i][j]=='r'){
umin(xl,i);
umax(xr,i);
umin(yl,j);
umax(yr,j);
}
for(i=xl;i<=xr;i++){
for(j=yl;j<=yr;j++)if(a[i][j]!='r')break;
if(j>yr)xm=i;
}
for(i=yl;i<=yr;i++){
for(j=xl;j<=xr;j++)if(a[j][i]!='r')break;
if(j>xr)ym=i;
}
for(i=xl;i<=xm;i++){
d=xm-i;
ul[d]=inf;ur[d]=0;
for(j=yl;j<=yr;j++)if(a[i][j]=='r'){
umin(ul[d],j);
umax(ur[d],j);
}
ul[d]=ym-ul[d];
ur[d]=ur[d]-ym;
}
for(i=xm;i<=xr;i++){
d=i-xm;
dl[d]=inf;dr[d]=0;
for(j=yl;j<=yr;j++)if(a[i][j]=='r'){
umin(dl[d],j);
umax(dr[d],j);
}
dl[d]=ym-dl[d];
dr[d]=dr[d]-ym;
}
mu=xm-xl;
md=xr-xm;
for(i=1;i<=n;i++)for(j=1;j<=n;j++)f[i][j]=0;
for(i=1;i<=n;i++){
fu[0]=fd[n+1]=inf;
for(j=1;j<=n;j++)if(a[j][i]=='X')fu[j]=0;else fu[j]=fu[j-1]+1;
for(j=n;j;j--)if(a[j][i]=='X')fd[j]=0;else fd[j]=fd[j+1]+1;
for(j=1;j<=n;j++){
l=n;
if(fu[j]<=mu)umin(l,i-ur[fu[j]]);
if(fd[j]<=md)umin(l,i-dr[fd[j]]);
r=1;
if(fu[j]<=mu)umax(r,i+ul[fu[j]]);
if(fd[j]<=md)umax(r,i+dl[fd[j]]);
umax(l,1),umin(r,n);
if(l<=r)f[j][l]++,f[j][r+1]--;
}
}
for(i=1;i<=n;i++)for(j=1;j<=n;j++){
f[i][j]+=f[i][j-1];
if(f[i][j])ban[i][j]=1;
}
f11=f1n=fn1=fnn=inf;
for(i=1;i<=md;i++){
umin(f11,i+dr[i]+1);
umin(f1n,i+dl[i]+1);
umin(fn1,i+ur[i]+1);
umin(fnn,i+ul[i]+1);
}
for(i=yl;i<=ym;i++){
d=ym-i;
ul[d]=inf;ur[d]=0;
for(j=xl;j<=xr;j++)if(a[j][i]=='r'){
umin(ul[d],j);
umax(ur[d],j);
}
ul[d]=xm-ul[d];
ur[d]=ur[d]-xm;
}
for(i=ym;i<=yr;i++){
d=i-ym;
dl[d]=inf;dr[d]=0;
for(j=xl;j<=xr;j++)if(a[j][i]=='r'){
umin(dl[d],j);
umax(dr[d],j);
}
dl[d]=xm-dl[d];
dr[d]=dr[d]-xm;
}
mu=ym-yl;
md=yr-ym;
for(i=1;i<=n;i++)for(j=1;j<=n;j++)f[i][j]=0;
for(i=1;i<=n;i++){
fu[0]=fd[n+1]=inf;
for(j=1;j<=n;j++)if(a[i][j]=='X')fu[j]=0;else fu[j]=fu[j-1]+1;
for(j=n;j;j--)if(a[i][j]=='X')fd[j]=0;else fd[j]=fd[j+1]+1;
for(j=1;j<=n;j++){
l=n;
if(fu[j]<=mu)umin(l,i-ur[fu[j]]);
if(fd[j]<=md)umin(l,i-dr[fd[j]]);
r=1;
if(fu[j]<=mu)umax(r,i+ul[fu[j]]);
if(fd[j]<=md)umax(r,i+dl[fd[j]]);
umax(l,1),umin(r,n);
if(l<=r)f[l][j]++,f[r+1][j]--;
}
}
for(i=1;i<=n;i++)for(j=1;j<=n;j++){
f[i][j]+=f[i-1][j];
if(f[i][j])ban[i][j]=1;
}
for(h=i=1;i<=n;i++)for(j=1;j<=n;j++)f[i][j]=-1;
add(xm,ym,0);
while(h<=t){
x=q[h][0],y=q[h++][1],z=f[x][y]+1;
add(x+1,y,z);
add(x-1,y,z);
add(x,y+1,z);
add(x,y-1,z);
}
if(ans<inf)printf("%d",ans);else puts("NIE");
return 0;
}