题意:
给一个n*n的网格,有黑白两色,问不重合的黑色矩阵对有多少。
题解:
不难的题,然而想了很久……
算出以某个点为左上角,右上角,左下角,右下角的矩形个数。
这个暴力枚举要
n3
,但是用单调栈可以
n2
,挺容易想的就不累述了。
然后枚举点(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);
}