bzoj 3235: [Ahoi2013]好方的蛇

题意：

给一个n*n的网格，有黑白两色，问不重合的黑色矩阵对有多少。

题解：

不难的题，然而想了很久……
算出以某个点为左上角，右上角，左下角，右下角的矩形个数。
这个暴力枚举要$n^3$，但是用单调栈可以$n^2$，挺容易想的就不累述了。
然后枚举点(i,j)，那么就可以用右下角的f乘上右边和下边的左上角f。

就是红色这部分的和。
但是有重复算的。

这样位置关系的就会算两次。
所以就有左下角乘左下方右上角的就可以了。
code：

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#define LL long long
using namespace std;
const LL mod=10007;
LL top,sta[1010];
LL n,map[1010][1010],MAP[1010][1010],g[1010][1010],f[1010][1010][4];
LL sum[1010][1010],num=0;
char s[1010];
//0右下角 1左上角 2右上角 3左下角
void pre(LL op)
{
    for(LL i=1;i<=n;i++)
        for(LL j=1;j<=n;j++)
            g[i][j]=(MAP[i][j]==1)?g[i-1][j]+1:0;
    for(LL i=1;i<=n;i++)
    {
        top=0;LL sum=0;
        for(LL j=1;j<=n;j++)
        {
            while(top!=0&&g[i][sta[top]]>=g[i][j]) sum-=(sta[top]-sta[top-1])*g[i][sta[top]],top--;
            sum+=(j-sta[top])*g[i][j];sta[++top]=j;
            f[i][j][op]=sum;
            if(f[i][j][op]!=0) f[i][j][op]--;
        }
    }
}
int main()
{
    scanf("%lld",&n);
    for(LL i=1;i<=n;i++)
    {
        scanf("%s",s+1);
        for(LL j=1;j<=n;j++) map[i][j]=(s[j]=='B'),num+=(s[j]=='B');
    }
    for(LL i=1;i<=n;i++) for(LL j=1;j<=n;j++) MAP[i][j]=map[i][j];pre(0);
    for(LL i=1;i<=n;i++) for(LL j=1;j<=n;j++) MAP[n-i+1][n-j+1]=map[i][j];pre(1);
    for(LL i=1;i<=n;i++) for(LL j=1;j<=n;j++) MAP[j][n-i+1]=map[i][j];pre(2);
    for(LL i=1;i<=n;i++) for(LL j=1;j<=n;j++) MAP[n-j+1][i]=map[i][j];pre(3);
    for(LL i=n;i>=1;i--)
        for(LL j=n;j>=1;j--) sum[i][j]=(sum[i][j+1]+f[n-i+1][n-j+1][1])%mod;
    for(LL i=n;i>=1;i--)
        for(LL j=n;j>=1;j--) (sum[i][j]+=sum[i+1][j])%=mod;
    LL ans=0;
    for(LL i=1;i<=n;i++)
        for(LL j=1;j<=n;j++)
            ans=(ans+f[i][j][0]*(sum[1][j+1]+sum[i+1][1]-sum[i+1][j+1])%mod)%mod;
    memset(sum,0,sizeof(sum));
    for(LL i=n;i>=1;i--)
        for(LL j=1;j<=n;j++) sum[i][j]=(sum[i][j-1]+f[j][n-i+1][2])%mod;
    for(LL i=n;i>=1;i--)
        for(LL j=1;j<=n;j++) (sum[i][j]+=sum[i+1][j])%=mod;
    for(LL i=1;i<=n;i++)
        for(LL j=1;j<=n;j++)
            ans-=(sum[i+1][j-1]*f[n-j+1][i][3])%mod,ans%=mod;
    printf("%lld",(ans+mod)%mod);
}