xjb——洛谷P1191 矩形

https://www.luogu.org/problem/show?pid=1191
来一个用xjb算法的题解；
xjb是什么算法？
这种算法基本上可以解决所有题目；
也许你不会lct；
也许你不会插头dp；
也许你不会莫比乌斯反演；
但是，只要你会xjb算法，noip甚至省选noi都不成问题！

我们对于这道题
要有一个初步的想法；
对于(i,j)这个点；
如果可以求出以这个点为右下角的所有矩阵；
那么这道题就好办了；
但是这个类似与dp的思想，却好像无法实现；
那怎么办呢？
我们对于(i,j)
维护x表示这个点向左最多连续几个W
维护y表示这个点向上最多连续几个W
显然(i,j)这个点最大贡献就是x*y
但是我们咋么求这个点的实际贡献呢？
我们从x*y这个矩阵的左下角开始向右上角搜索；
勾出一条边界线；
边界线左下放全部是W
这样我们一遍勾出边界，一边计算答案就好了；
时间复杂度
150*150*300
满

#include<bits/stdc++.h>
using namespace std;
const int N=155;
int a[N][N],x[N][N],y[N][N],s[N][N];
int ans,n,m,xx,yy;
char c;
bool check(int x,int y,int xx,int yy){
    return (s[xx][yy]-s[xx][y-1]-s[x-1][yy]+s[x-1][y-1])==(xx-x+1)*(yy-y+1);
}
int dfs(int a1,int a2,int xx,int yy){
    int ans=0;
    int x=xx;
    int y=yy-a1+1;
    while(a1&&a2){
        ans+=a1;
        while(!check(x-1,y,xx,yy))y+=1,a1--;
        x-=1;
    }
    return ans;
}
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++)for(int j=1;j<=n;j++){cin>>c;if(c=='W')a[i][j]=1;s[i][j]=a[i][j];}
    for(int i=1;i<=n;i++)for(int j=1;j<=n;j++)s[i][j]+=s[i][j-1];
    for(int i=1;i<=n;i++)for(int j=1;j<=n;j++)s[i][j]+=s[i-1][j];
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
            if(a[i][j]){
                x[i][j]=x[i][j-1]+1;
                y[i][j]=y[i-1][j]+1;
                ans+=dfs(x[i][j],y[i][j],i,j);
            }
    printf("%d",ans);
}