「模拟8.23」阴阳 DP

对于此题的性质我们考虑DP

分四种情况

黑色块在右侧单调降，单调升

还有在左侧

另外我们这样可能会记重，所以还要将重复记过的也就是边界线是横的和竖的

然后还要将全白全黑加上

#include<bits/stdc++.h>
#define MAXN 2100
#define int long long
using namespace std;
int f[5][MAXN][MAXN];
int sum[5][MAXN][MAXN];
int n,m;
const int mod=1e9+7;
char c[MAXN][MAXN];
int jud[2][MAXN][MAXN];
int jud_sum[2][MAXN][MAXN];
void work(int k,int cal)
{    
    f[k][0][m]=1;
    for(int i=1;i<=n;++i)//黑色单降靠左
    {
        sum[k][i-1][m]=f[k][i-1][m];
        for(int j=m-1;j>=0;--j)
        {
           sum[k][i-1][j]=(sum[k][i-1][j+1]+f[k][i-1][j])%mod;
           //printf("sum[%lld][%lld][%lld]=%lld
",k,i-1,j,sum[k][i-1][j]);
        }
        for(int j=0;j<=m;++j)
        {  
            //printf("i%lld j=%lld jud%lld %lld
",i,j,jud[cal][i][j],(jud[cal^1][i][m]-jud[cal^1][i][j]));
            if(j!=0)
                if(jud[cal][i][j]!=0||(jud[cal^1][i][m]-jud[cal^1][i][j])!=0){continue;}
            if(j==0)
                if(jud[cal][i][j]||jud[cal^1][i][m]){continue;}
                f[k][i][j]=sum[k][i-1][j]%mod;
            //printf("f[%lld][%lld][%lld]=%lld
",k,i,j,f[k][i][j]);
        }
    }
    f[k+1][0][0]=1;
    for(int i=1;i<=n;++i)//黑色单升靠左
    {
        sum[k+1][i-1][0]=f[k+1][i-1][0];
        for(int j=1;j<=m;++j)
        {
           sum[k+1][i-1][j]=(sum[k+1][i-1][j-1]+f[k+1][i-1][j])%mod;
        }
        for(int j=0;j<=m;++j)
        {  
            if(j)
                if(jud[cal][i][j]||(jud[cal^1][i][m]-jud[cal^1][i][j]))continue;
            if(j==0)
                if(jud[cal][i][j]||jud[cal^1][i][m])continue;
            f[k+1][i][j]=sum[k+1][i-1][j]%mod;
            //printf("f[%lld][%lld][%lld]=%lld
",k+1,i,j,f[k+1][i][j]);
        }
    }
}
int ans=0;
int get_sum(int k,int x1,int y1,int x2,int y2)
{
    if(x1>x2)return 0;if(y1>y2)return 0;
    return (jud_sum[k][x2][y2]-jud_sum[k][x1-1][y2]-jud_sum[k][x2][y1-1]+jud_sum[k][x1-1][y1-1]+mod)%mod;
}
void work2()
{
     for(int i=0;i<=n;++i)
     {
         if(get_sum(1,1,1,i,m)==0&&get_sum(0,i+1,1,n,m)==0)
         {
             ans--;
         }
     }
     for(int i=0;i<=n;++i)
     {
         if(get_sum(0,1,1,i,m)==0&&get_sum(1,i+1,1,n,m)==0)
         {
             ans--;
         }
     }
     for(int i=0;i<=m;++i)
     {
         if(get_sum(0,1,1,n,i)==0&&get_sum(1,1,i+1,n,m)==0)
         {
             ans--;
         }
     }
     for(int i=0;i<=m;++i)
     {
         if(get_sum(1,1,1,n,i)==0&&get_sum(0,1,i+1,n,m)==0)
         {
             ans--;
         }
     }     
     if(get_sum(1,1,1,n,m)==0)ans++;
     if(get_sum(0,1,1,n,m)==0)ans++;
}
signed main()
{
    scanf("%lld%lld",&n,&m);
    for(int i=1;i<=n;++i)
    {
        scanf("%s",c[i]+1);
        for(int j=1;j<=m;++j)
        {
            jud[0][i][j]+=jud[0][i][j-1];
            jud[1][i][j]+=jud[1][i][j-1];
            jud_sum[0][i][j]=(jud_sum[0][i-1][j]+jud_sum[0][i][j-1]-jud_sum[0][i-1][j-1])%mod;
            jud_sum[1][i][j]=(jud_sum[1][i-1][j]+jud_sum[1][i][j-1]-jud_sum[1][i-1][j-1])%mod;            
            if(c[i][j]=='B'){jud[0][i][j]++;jud_sum[0][i][j]++;}
            if(c[i][j]=='W'){jud[1][i][j]++;jud_sum[1][i][j]++;}
        }
    }
    work(1,1);
    work(3,0);
    for(int k=1;k<=4;++k)
    {
        for(int i=0;i<=m;++i)
        {
            ans=(ans+f[k][n][i])%mod;
            //printf("f[%lld][%lld][%lld]=%lld
",k,n,i,f[k][n][i]);
        }
    }
    //printf("初ans=%lld
",ans);
    work2();
    printf("%lld
",ans%mod);
}