題目描述

做法

首先加入所有詢問座標，求出答案。
初始答案很好求，你可以設f[i,j]表示以(i,j)為左下角的最大正方形邊長。
那麼f[i,j]=min(f[i−1,j],f[i,j+1])(+1)後面是否+1要判斷右上角。
當然左下角是障礙f[i,j]=0。
接下來來考慮時間倒流。
首先我們維護left和right表示每個位置可以往左往右延伸多長。
時間倒流後，相當於刪除障礙，答案一定是遞增的。
刪除一個障礙只需要暴力更改這一行的left和right。
接下來我們不斷判斷是否存在ans+1（ans是當前答案）邊長的正方形，存在將ans+1，繼續做。
容易得知這個正方形一定跨過障礙所在列。
於是可以單調佇列，得到障礙所在列每一個長度為ans+1的區間，向左向右能延伸多長，如果加起來達到ans+1，則判定成功。
複雜度是n^2。

#include<cstdio>
#include<algorithm>
#define fo(i,a,b) for(i=a;i<=b;i++)
#define fd(i,a,b) for(i=a;i>=b;i--)
using namespace std;
const int maxn=2000+10;
int f[maxn][maxn],left[maxn][maxn],right[maxn][maxn];
int ans[maxn],b[maxn],dl[maxn],ask[maxn][2];
int pd[maxn][maxn];
int i,j,k,l,t,n,m,q,head,tail,now,x,y;
int read(){
    int x=0,f=1;
 

    char ch=getchar();
    while (ch<'0'||ch>'9'){
        if (ch=='-') f=-1;
        ch=getchar();
    }
    while (ch>='0'&&ch<='9'){
        x=x*10+ch-'0';
        ch=getchar();
    }
    return x*f;
}
char get(){
    char ch=getchar();
    while (ch!='.'&&ch!='X') ch=getchar();
 

    return ch;
}
bool check(int ans){
    int i,j,k;
    head=1;tail=0;
    fo(i,1,n){
        while (head<=tail&&left[dl[tail]][y]>=left[i][y]) tail--;
        dl[++tail]=i;
        while (dl[head]<=i-ans) head++;
        b[i]=left[dl[head]][y];
    }
    head=1;tail=0;
    fo(i,1,n){
        while (head<=tail&&right[dl[tail]][y]>=right[i][y]) tail--;
        dl[++tail]=i;
        while (dl[head]<=i-ans) head++;
        b[i]+=right[dl[head]][y]-1;
    }
    fo(i,1,n)
        if (i>=ans&&b[i]>=ans) return 1;
    return 0;
}
int main(){
    scanf("%d%d%d",&n,&m,&q);
    fo(i,1,n)
        fo(j,1,m)
            if (get()=='X') pd[i][j]++;
    fo(i,1,q){
        ask[i][0]=read();ask[i][1]=read();
        pd[ask[i][0]][ask[i][1]]++;
    }
    fo(i,1,n)
        fd(j,m,1){
            if (pd[i][j]) continue;
            if (i==1||j==m){
                f[i][j]=1;
                now=max(now,f[i][j]);
                continue;
            }
            f[i][j]=min(f[i-1][j],f[i][j+1]);
            if (!pd[i-f[i][j]][j+f[i][j]]) f[i][j]++;
            now=max(now,f[i][j]);
        }
    fo(i,1,n)
        fo(j,1,m)
            if (pd[i][j]) left[i][j]=0;else left[i][j]=left[i][j-1]+1;
    fo(i,1,n)
        fd(j,m,1)
            if (pd[i][j]) right[i][j]=0;else right[i][j]=right[i][j+1]+1;
    fd(i,q-1,1){
        ans[i+1]=now;
        x=ask[i+1][0];y=ask[i+1][1];
        pd[x][y]--;
        fo(j,1,m)
            if (pd[x][j]) left[x][j]=0;else left[x][j]=left[x][j-1]+1;
        fd(j,m,1)
            if (pd[x][j]) right[x][j]=0;else right[x][j]=right[x][j+1]+1;
        while (check(now+1)) now++;
    }
    ans[1]=now;
    fo(i,1,q) printf("%d\n",ans[i]);
}