【五校联考2day1】池塘

Description

从前有两个青蛙王国，两个王国商业都非常繁荣。但是一块池塘阻碍了两国的商业往来。一次，两只青蛙在池塘的两岸，他们都希望到对岸去。我们可以将池塘看做一个n×m的矩形，在每个格子里，可能会有荷叶。青蛙必须踩在荷叶上，不能跳进水里。如图青蛙可以向他前方的5个有荷叶的地方跳去。

由于有的地方荷叶比较小，当一个青蛙从该荷叶上跳走之后，荷叶会沉入水底，两个青蛙也不能同时跳上这种荷叶。两个青蛙想知道有多少种方式使他们都到达对岸。第一个青蛙可以从第一行任何一个有荷叶的格子出发。第二个青蛙可以从最后一行任何一个有荷叶的格子出发。当第一个青蛙到达最后一行任何一个有荷叶的格子时，他就算到达了对岸。当第二个青蛙到达第一行任何一个有荷叶的格子时，他也算到达了对岸。请你帮助青蛙们计算有多少种方案可以让他们都到达对岸。
注：第一个青蛙只能向下跳，第二个青蛙只能向上跳。青蛙并不能跳出矩形区域。

Input

输入的第一行包含两个整数n，m。
第2至n+1行包含m个整数。若该数为0，表示该格子上没有荷叶，青蛙不能通过。若该数为1，表示该格子上的荷叶只允许一个青蛙通过。若该数为2，表示该格子上的荷叶可以允许两个青蛙都通过。

Output

输出的第一行包含一个整数，表示两个青蛙都到达对岸的方案数。由于结果可能非常大，输出答案模1000000007的结果。

分析

首先，这是一道很明显的dp题，但是关键却是如何转移状态。
当然我们还是先把两只青蛙不同方向地走改成两只青蛙从同向去走，这样方便我们判断路径是否合法。
明显的：f[i][j][k][l]表示一只青蛙走到(i,j)，另一只青蛙走到(k,l)的可行的方案数。
那么难道我们的转移就直接每只青蛙都转移一次吗？不！这样会重复！
一只青蛙走的很快时，你就不知道另一只慢的青蛙往哪里走是合法的。所以我们可以限制它们的abs(i-k)<=3，然后再让走的慢的那只先走就可以了。

这样我们便可以解决这个问题了，但我还是很傻地调了很久，原因是我在转移状态时一直在纠结两只青蛙的坐标放到哪里。改着改着我还甚至硬性规定了左边坐标大于右边最后在乘2。。但其实没有那么复杂，比如（i，j）转到（i’，j’），（x，y）转到（x’，y’），那么转移时就直接f[i’,j’,x’,y’]….f[i,j,x,y]当然我想就只有我这种蒟蒻才会犯这种错误吧。。

这题的初始化我也纠结了很久，有个地方一直打错了，是个小细节。这再次暴露了我的粗心，浮躁。要打一步想一步。

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#define ll long long
using namespace std;
const int N=55;
const int mo=1000000007;
ll n,m,w[N],a[N][N],f[N][N][N][N],x,y,ans,t[5][2],i,j,k,l;
void make(int i,int j,int k,int l,int z,int op){
    for(int r=0;r<=4;r++){
        x=i+t[r][0];y=j+t[r][1];
        if (x<1||y<1||x>n||y>m||a[x][y]==0) continue;
        if (x==k&&y==l&&a[x][y]<2) continue;
        if (op)f[k][l][x][y]+=z,f[k][l][x][y]%=mo;
        else f[x][y][k][l]+=z,f[x][y][k][l]%=mo;
    }
}
int main(){
    scanf("%d %d",&n,&m);
    t[0][0]=1;t[0][1]=-2;t[1][0]=2;t[1][1]=-1;t[2][0]=3;t[3][0]=2;t[3][1]=1;t[4][0]=1;t[4][1]=2;
    for(i=1;i<=n;i++)
        for(j=1;j<=m;j++)
            scanf("%lld",&a[i][j]);

    for(i=1;i<=m;i++)
        for(j=1;j<=m;j++)
        if (i==j&&a[1][i]>1) f[1][i][1][i]=1;
        else if (i!=j&&a[1][i]&&a[1][j])f[1][i][1][j]=1;
    for(i=1;i<=n;i++)
       for(j=1;j<=m;j++)
           for(k=max(1*1LL,i-3);k<=min(n,i+3);k++)
               for(l=1;l<=m;l++){
               if (f[i][j][k][l]>0){
                    ll we=f[i][j][k][l];
                   if (i<k) make(i,j,k,l,we,0);
                    else make(k,l,i,j,we,1);
               }
            }
    for(i=1;i<=m;i++)
       for(j=1;j<=m;j++) 
       ans+=f[n][i][n][j],ans%=mo;
    printf("%lld",ans);
}