51nod 1503 猪和回文串(动态规划)

1503 猪和回文

题目来源： CodeForces

基准时间限制：2 秒 空间限制：131072 KB 

 

一只猪走进了一个森林。很凑巧的是，这个森林的形状是长方形的，有n行，m列组成。我们把这个长方形的行从上到下标记为1到n，列从左到右标记为1到m。处于第r行第c列的格子用(r,c)表示。

刚开始的时候猪站在(1,1),他的目标是走到(n,m)。由于猪回家心切，他在(r,c)的时候，只会往(r+1,c)或(r,c+1)走。他不能走出这个森林。

这只猪所在的森林是一个非同寻常的森林。有一些格子看起来非常相似，而有一些相差非常巨大。猪在行走的过程中喜欢拍下他经过的每一个格子的照片。一条路径被认为是漂亮的当且仅当拍下来的照片序列顺着看和反着看是一样的。也就是说，猪经过的路径要构成一个回文。

数一数从(1,1)到(n,m)有多少条漂亮路径。答案可能非常巨大，请输出对 109+7 取余后的结果。

样例解释：有三种可能

  

Input

单组测试数据。
第一行有两个整数 n,m (1≤n,m≤500)，表示森林的长和宽。
接下来有n行，每行有m个小写字母，表示每一个格子的类型。同一种类型用同一个字母表示，不同的类型用不同的字母表示。

Output

输出答案占一行。

Input示例

3 4
aaab
baaa
abba

Output示例

3

dp[i][j][k]代表从起点走到(i,j)和从终点走到(k,t);由于已知i,j,k可以推导出t,所以空间上可以优化

　　dp[i][j][k]=dp[i][j][k]+dp[i-1][j][k]
　　dp[i][j][k]=dp[i][j][k]+dp[i-1][j][k+1]
　　dp[i][j][k]=dp[i][j][k]+dp[i][j-1][k]
　　dp[i][j][k]=dp[i][j][k]+dp[i][j-1][k+1]

共四种情况,可以发现dp[i][][]由dp[i-1][][]得出，所以空间上还可以优化,滚动数组,记得滚动的时候不忘清0

奇偶可以统一讨论，不要想的太复杂

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<string>
#include<algorithm>
#include<iostream>
#include<queue>
#include<map>
#include<set>
#define ll long long
#define max(x,y) (x)>(y)?(x):(y)
#define min(x,y) (x)>(y)?(y):(x)
#define cls(name,x) memset(name,x,sizeof(name))
using namespace std;
const int inf=1<<28;
const int maxn=510;
const int maxm=1e2+10;
const int mod=1e9+7;
int dp[2][maxn][maxn];
char mapp[maxn][maxn];
bool judge(int x1,int y1,int x2,int y2)
{
    if(x1==x2&&y1==y2)
        return true;
    if(x1+1==x2&&y1==y2)
        return true;
    if(x1==x2&&y1+1==y2)
        return true;
    return false;
}
int main()
{
    int n,m;
    //freopen("in.txt","r",stdin);
    scanf("%d%d",&n,&m);
    {
        int ans=0;
        for(int i=1;i<=n;i++)
        {
            scanf("%s",mapp[i]);
            for(int j=m;j>=1;j--)
                mapp[i][j]=mapp[i][j-1];
        }
        if(mapp[1][1]==mapp[n][m])
            dp[1][1][n]=1;
        for(int i=1;i<=n;i++)
        {
            for(int j=1;(i+j-1)<=(n-i+m-j+1);j++)//始终在上三角，包括对角线
                for(int k=n;(i+j-1)>=n-k+1&&k>=1;k--)
                    if(mapp[i][j]==mapp[k][n-k+m-i-j+2])
                    {
                        dp[i%2][j][k]=(dp[i%2][j][k]+dp[(i-1)%2][j][k])%mod;
                        dp[i%2][j][k]=(dp[i%2][j][k]+dp[(i-1)%2][j][k+1])%mod;
                        dp[i%2][j][k]=(dp[i%2][j][k]+dp[i%2][j-1][k])%mod;
                        dp[i%2][j][k]=(dp[i%2][j][k]+dp[i%2][j-1][k+1])%mod;if(judge(i,j,k,n-k+m-i-j+2))
                                ans=(ans+dp[i%2][j][k])%mod;
                    }
            cls(dp[(i-1)%2],0);
        }
        printf("%d\n",ans);
    }
    return 0;
}