传纸条 NOIP2008（线性DP）

题干：
小渊和小轩是好朋友也是同班同学，他们在一起总有谈不完的话题。一次素质拓展活动中，班上同学安排做成一个m行n列的矩阵，而小渊和小轩被安排在矩阵对角线的两端，因此，他们就无法直接交谈了。幸运的是，他们可以通过传纸条来进行交流。纸条要经由许多同学传到对方手里，小渊坐在矩阵的左上角，坐标(1,1)，小轩坐在矩阵的右下角，坐标(m,n)。从小渊传到小轩的纸条只可以向下或者向右传递，从小轩传给小渊的纸条只可以向上或者向左传递。
在活动进行中，小渊希望给小轩传递一张纸条，同时希望小轩给他回复。班里每个同学都可以帮他们传递，但只会帮他们一次，也就是说如果此人在小渊递给小轩纸条的时候帮忙，那么在小轩递给小渊的时候就不会再帮忙。反之亦然。
还有一件事情需要注意，全班每个同学愿意帮忙的好感度有高有低（注意：小渊和小轩的好心程度没有定义，输入时用0表示），可以用一个0-100的自然数来表示，数越大表示越好心。小渊和小轩希望尽可能找好心程度高的同学来帮忙传纸条，即找到来回两条传递路径，使得这两条路径上同学的好心程度值和最大。现在，请你帮助小渊和小轩找到这样的两条路径。

数据范围及提示 Data Size & Hint
30%的数据满足：1<=m,n<=10
100%的数据满足：1<=m,n<=50

输入格式：
的第一行有2个用空格隔开的整数m和n，表示班里有m行n列（1<=m,n<=50）。
接下来的m行是一个m*n的矩阵，矩阵中第i行j列的整数表示坐在第i行j列的学生的好心程度。每行的n个整数之间用空格隔开。
输出格式：
共一行，包含一个整数，表示来回两条路上参与传递纸条的学生的好心程度之和的最大值。

输入样例：
3 3
0 3 9
2 8 5
5 7 0

输出样例：
34

思路：
题目大意：给定一个N*M的矩阵A，每个格子中有一个数，现在需要找到两条从左上角（1,1）到右下角（N，M）的路径，路径上的每一步只能向右向下走，路径经过的格子中的数会被取走，若两条路径同时经过一个格子，只算一次。求取得的数之和最大是多少。N，M<=50。

根据Mobile Service TYVJ1061（线性DP）这题我们类比一下，把每走的步数i当成阶段，扩展到下一步还需要知道上一步两条路径的坐标。设上一次路径长度为i，第一条路径末尾位于（x1，y1），第二条路径末尾位于（x2，y2）。我们需要找到能够覆盖这个状态空间的“最小维度”。

我们检查上面五个值有何种关系，能否互相能导出。容易发现有下面等式成立：
x1+y1=x2+y2=i+2

过i，x1，x2三个维度就可以描述一个状态。
设F[i，x1，x2]表示路径长度为i，第一条路径末尾在第x1行，第二条路径末尾在x2行时，已经取得的数的和的最大值。
此时可以推出y1=i+2-x1，y2=i+2-x2。
每条路径有向右、向下两种扩展方法，故共有2*2=4种转移。
两条路径的组合为右右，下下，右下，下右。
以右右为例:
如果x1=x2，并且y1+1=y2+1, 那么两条路径进入同一格子，只累加一次。
F[i+1,x1,x2]=max（F[i+1,x1,x2],F[i,x1,x2]+A[x1,y1+1]）
否则，两条路径走不同的格子，两个格子中的数都进行累加：

F[i+1,x1,y1]=max(F[i+1,x1,y1],F[i,x1,x2]+A[x1,y1+1]+A[x2,y2+1])
初值为F[0,1，1]=A[1,1].
目标为F[N+M-2，N，N]。

#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
int f[105][55][55],a[55][55],n,m,i,x1,x2,j,ans=0;
int ff(int k,int x1,int x2)
{
    if (k<0||x1<0||x2<0) 
       return 0;
    return f[k][x1][x2];
}
int main()
{   scanf("%d%d",&m,&n);
     for (i=1;i<=m;i++)
     for (j=1;j<=n;j++)
         scanf("%d",&a[i][j]);
     memset(f,0,sizeof(f));
     for (i=1;i<=n+m-2;i++)
         for (x1=0;x1<=min(i,n);x1++)
             for (x2=0;x2<=min(i,n);x2++)
                 {
				   if (x1==x2) continue;
                   f[i][x1][x2]=max(max(ff(i-1,x1,x2),ff(i-1,x1-1,x2)),max(ff(i-1,x1,x2-1),ff(i-1,x1-1,x2-1)))+a[i-x1+1][1+x1]+a[i-x2+1][1+x2];
                   if (f[i][x1][x2]>ans) ans=f[i][x1][x2];
                 }
cout<<ans;
}
/*
3 3
0 3 9
2 8 5
5 7 0
*/ 

为了加深对上面的理解，把O(n⁴)复杂度的代码放出来：

设dp[i][j][k][l]是第一张纸条到达(i,j)，第二张到达(k,l)时好心值最大

那么方程就是dp[i][j][k][l] = map[i][j] + map[k][l] + max(dp[i - 1][j][k - 1][l],dp[i -1][j][k][l - 1],dp[i][j - 1][k -1][l],dp[i][j -1][k][]l -1);

还有一点注意的是，如果i == k && j == l，也就是第二张走了第一张的路径，那么就要减去，dp[i][j][k][l] -= map[i][j];

接下来枚举每个i,j,k,l，最后输出dp[m][n][m][n]就行了，但这样做时间复杂度是O(n⁴)

#include<cstdio>
#include<iostream>
using namespace std;

int dp[55][55][55][55];
int n, m;
int map[55][55];

int main()
{
    scanf("%d%d", &n, &m);
    for(int i = 1; i <= n; i++)
    for(int j = 1; j <= m; j++)
    scanf("%d", &map[i][j]);
    for(int i = 1; i <= n; i++)
    {
        for(int j = 1; j <= m; j++)
        {
            for(int k = 1; k <= n; k++)
            {
                int l;
                if(i + j - k > 0)l = i + j - k;
                else continue;
                dp[i][j][k][l]=max(max(dp[i-1][j][k-1][l],dp[i-1][j][k][l-1]),max(dp[i][j-1][k-1][l],dp[i][j-1][k][l-1]))+map[i][j]+map[k][l];
                if(i == k && j == l)dp[i][j][k][l] -= map[i][j];
            }
        }
    }
    printf("%d\n",dp[n][m][n][m]);
    return 0;
}
/*
3 3
0 3 9
2 8 5
5 7 0
*/