FZU 2234 牧场物语【Dp】好题！好题！

 Problem 2234 牧场物语

Accept: 172    Submit: 876
Time Limit: 1000 mSec    Memory Limit : 32768 KB

 Problem Description

小茗同学正在玩牧场物语。该游戏的地图可看成一个边长为n的正方形。

小茗同学突然心血来潮要去砍树，然而，斧头在小茗的右下方。

小茗是个讲究效率的人，所以他会以最短路程走到右下角，然后再返回到左上角。并且在路上都会捡到/踩到一些物品，比如说花朵，钱和大便等。

物品只能被取最多一次。位于某个格子时，如果格子上还有物品，就一定要取走。起点和终点上也可能有物品。

每种物品我们将为其定义一个价值，当然往返之后我们取得的物品的价值和越大越好。但是小茗同学正在认真地玩游戏，请你计算出最大的价值和。

 Input

多组数据(<=10)，处理到EOF。

第一行输入正整数N(N≤100)，表示正方形的大小。

接下来共N行，每行N个整数Ai,j(|Ai,j|≤10^9)，表示相应对应位置上物品的价值。值为0表示没有物品。

 Output

每组数据输出一个整数，表示最大价值和。

 Sample Input

2

11 14

16 12

 Sample Output

53

 Source

福州大学第十三届程序设计竞赛

思路：

FZU校赛的质量确实很高。

好多ACM强校的校赛都可以拿来当省赛打了= =

1、一去一回，我们可以看做两个去。

所以这类题可以直接跑费用流。但是一定要考虑是否存在负环。

这个题存在，所以直接跑费用流是不行的。

2、那么考虑Dp.类似这种Dp直接跑多线程Dp就行，即设定dp【i】【j】【k】【l】表示第一个人走到了（i，j）第二个人走到（k，l）的最优方案，那么其每个人都有两种走法，所以四个状态转移很轻松搞定问题。

然而时间复杂度和空间复杂度都不允许，虽然我们可以对数组进行滚动，但是毕竟时间也是一大障碍。

我们要Dp的就是位子，又已知每次行走都是最短路，那么我们不妨设定dp【i】【j】【k】表示第一个人走到第i行，第二个人走到第j行，一共两个人一起走了k步的方案。

那么我们可以通过k和i，j之间的关系计算出两个人现在的位子，同理，对四种走法Dp一下转移即可。

思路很宽广的一道题。

Ac代码：

#include<stdio.h>
#include<iostream>
#include<string.h>
using namespace std;
#define ll __int64
ll dp[105][105][220];
ll a[105][105];
int main()
{
    int n;
    while(~scanf("%d",&n))
    {
        for(int i=0;i<=n;i++)
        {
            for(int j=0;j<=n;j++)
            {
                for(int k=0;k<=210;k++)
                {
                    dp[i][j][k]=-1000000000000000000ll;
                }
            }
        }
        for(int i=1;i<=n;i++)
        {
            for(int j=1;j<=n;j++)
            {
                scanf("%I64d",&a[i][j]);
            }
        }
        dp[1][1][0]=a[1][1];
        for(int k=1;k<=(n-1)*2;k++)
        {
            for(int i=1;i<=n;i++)
            {
                int c1=2+k-i;
                if(c1>=1&&c1<=n)
                for(int j=1;j<=n;j++)
                {
                    int c2=2+k-j;
                    if(c2>=1&&c2<=n)
                    {
                        ll add=a[i][c1]+a[j][c2];
                        if(i==j&&c1==c2)add/=2;
                        dp[i][j][k]=max(dp[i][j][k],dp[i][j][k-1]+add);
                        dp[i][j][k]=max(dp[i][j][k],dp[i][j-1][k-1]+add);
                        dp[i][j][k]=max(dp[i][j][k],dp[i-1][j][k-1]+add);
                        dp[i][j][k]=max(dp[i][j][k],dp[i-1][j-1][k-1]+add);
                    }
                }
            }
        }
        printf("%I64d\n",dp[n][n][(n-1)*2]);
    }
}