桂电第一届程序设计 D 双线程DP

Problem 1198 # 情迷意乱
ediszhao添加于2017-12-22 15:01:51
时间限制 : 1000MS 内存限制 : 65536KB
×提示 : 你还没有解决此题目!
问题描述
小明和小红是同班同学，上课了，小明心里乱慌慌的，因为刚才下课时小红盯着小明看了一会，小明想急切知道小红刚才为什么盯着他。他们坐在一个m行n列的教室里，而小明和小红被安排在矩阵对角线的两端，因此小明需要通过传纸条来问小红缘由。纸条要经由许多同学传到对方手里，小明坐在教室的左上角，坐标为(1,1)，小红坐在教室的右下角，坐标为(m,n)。从小明传到小红的纸条只可以向下或者向右传递，从小红传给小明的纸条只可以向上或者向左传递。班里每个同学都可以帮他们传递，但只会帮他们一次，也就是说如果此人在小明递给小红纸条的时候帮忙，那么在小红递给小明的时候就不会再帮忙。反之亦然。还有需要注意的是，全班每个同学愿意帮忙的心情指数有高有低，可以用一个0到100的自然数来表示，数越大表示心情越好。小明和小红的心情指数皆为0。小明希望尽可能找心情指数高的同学来帮忙传纸条让小红早点回复他。因此他想找出心情指数和达到最大的来回传递路径，并想知道它的的值是多少。现在，请你帮帮心猿意乱的小明。

输入
第一行输入两个整数m和n，表示班里有m行n列（1<=m,n<=50）。

接下来的m行是一个m*n的式子，式子中第i行j列的整数表示坐在第i行j列的学生的心情指数。

输出
输出来回传递途径的心情指数的和达到最大的值。

输入范例
3 3

0 3 9

2 8 5

5 7 0

输出范例
34

首先这肯定是DP题目。可以看错两条线。
用四维数组存放嘛。然后就把DP方程列错了。
这样的每次只走一步。就不是线程了好吗。。尴尬。放上来警示自己吧。

                        dp[i][j][k][t] = max(dp[i][j][k][t],dp[i-1][j][k][t]+M[i][j]);
                        dp[i][j][k][t] = max(dp[i][j][k][t],dp[i][j-1][k][t]+M[i][j]);
                        dp[i][j][k][t] = max(dp[i][j][k][t],dp[i][j][k-1][t]+M[k][t]);
                        dp[i][j][k][t] = max(dp[i][j][k][t],dp[i][j][k][t-1]+M[k][t]);

正确的DP方程应该是
dp[i][j][x][y] = max( max(dp[i-1][j][x-1][y] , dp[i-1][j][x][y-1] ), max( dp[i][j-1][x-1][j] , dp[i][j-1][x][j-1]))

但是如果用4个循环怎么都感觉要超时，所以肯定要优化。优化就看看有没有没办法优化方程或者说减少循环了。至于说减少循环我们可以通过观察两条线的位置发现。由于走的步数相同所以在与当前知道的一格子的斜线上。所以就可以设k=i+1,做另一线程端的x，然后i+j-k就是y了。所以用三个循环就表示了。
代码如下：

#include<stdio.h>  
#include<string.h>  
#define max 50+5  
int dp[max][max][max][max];  
int a[max][max];  
int MAX(int x,int y)  
{  
    return x>y?x:y;  
}  
int main()  
{  
    int n,m,sum;  
    int i,j,k,l;  

    while(scanf("%d%d",&n,&m))  
    {  
        for(i=1;i<=n;i++)  
        {  
            for(j=1;j<=m;j++)  
            {  
                scanf("%d",&a[i][j]);  
            }  
        }  
        for(i=1;i<=n;i++)  
        {  
            for(j=1;j<=m;j++)  
            {  
                for(k=i+1;k<=n;k++)  
                {  
                    l=i+j-k;  
                    if(l<=0)  
                    continue;  
                    dp[i][j][k][l]=MAX(MAX(dp[i-1][j][k-1][l],dp[i][j-1][k][l-1]),MAX(dp[i-1][j][k][l-1],dp[i][j-1][k-1][l]))+a[i][j]+a[k][l];  
                }  
            }  
        }  
        sum=MAX(MAX(dp[n-1][m][n-1][m],dp[n][m-1][n][m-1]),MAX(dp[n-1][m][n][m-1],dp[n][m-1][n-1][m]))+a[n][m];  
        printf("%d\n",sum);  
    }  
    return 0;  
} 

ps：上述代码在桂电的OJ上通过不了。因为是while（）输入。他们OJ判断不了？？？会超时？？？所以我一直说他们OJ是真的。。咳咳。。