HDU 2686 (双线程) Matrix

这也是当初卡了很久的一道题

题意：从左上角的格子出发选一条路径到右上角然后再回到左上角，而且两条路径除了起点和终点不能有重合的点。问所经过的格子中的最大和是多少

状态设计：我们可以认为是从左上角出发了两条路径，然后同时到达右下角。容易看出，第k个阶段所有可能到达的格子构成一条斜线而且满足x1 + y1 = x2 + y2 = k + 1

dp[k][x1][x2]表示第k阶段两条路分别到达(x1, y1) (x2, y2)所能取得的最大值(y1 y2根据上面的等量关系来计算)，如果x1 = x2表示两条路径汇合于一点了

状态转移方程：

在求解第k阶段最优解时，首先枚举k-1阶段时的状态(x1‘, x2‘)，然后向四个方向延伸出(x1, x2)

dp[k][x1][x2] = max{dp[k-1][x1‘][x2‘] + num[x1][y1] | (x1=x2),  dp[k-1][x1‘][x2‘] + num[x1][y1] + num[x2][y2] | (x1≠x2)}

在循环的过程中可能会出现从一个点延伸的下一个点的情况，即两条路径可能重合，不过不要紧，因为格子里面都是正数，所以得到的最优解一定是两条不重合的路径

优化：

因为两条路径具有任意性，所以不妨规定第一条路在第二条路的下方(或同一水平线)，时间从92MS优化到62MS

第k阶段状态的转移只依赖第k-1阶段的状态，可以用滚动数组(由于原来的空间本不大，所以代码中没有实现滚动数组)

 1 //#define LOCAL
 2 #include <iostream>
 3 #include <cstdio>
 4 #include <cstring>
 5 using namespace std;
 6
 7 int dp[70][35][35], num[35][35], n;
 8
 9 bool islegal(int x, int y)
10 {
11     return (x>=1 && x<=n && y>=1 && y<=n);
12 }
13
14 int main(void)
15 {
16     #ifdef LOCAL
17         freopen("2686in.txt", "r", stdin);
18     #endif
19
20     int i, j, k, d1, d2, x1, x2, y1, y2;
21     while(scanf("%d", &n) == 1)
22     {
23         memset(dp, 0, sizeof(dp));
24         for(i = 1; i <= n; ++i)
25             for(j = 1; j <= n; ++j)
26                 scanf("%d", &num[i][j]);
27
28         dp[1][1][1] = num[1][1];
29         for(k = 2; k <= 2*n-1; ++k)
30         {
31             for(i = 1; i <= k - 1; ++i)
32                 for(j = 1; j <= i; ++j)
33                 {//枚举两条路径在k-1步时的状态
34                     for(d1 = 0; d1 <= 1; ++d1)
35                         for(d2 = 0; d2 <= 1; ++d2)
36                         {//共四个延伸方向
37                             x1 = i + d1, y1 = k + 1 - x1;
38                             x2 = j + d2, y2 = k + 1 - x2;
39                             if(islegal(x1, y1) && islegal(x2, y2))
40                             {
41                                 if(x1 == x2)
42                                     dp[k][x1][x2] = max(dp[k][x1][x2], dp[k-1][i][j] + num[x1][y1]);
43                                 else
44                                     dp[k][x1][x2] = max(dp[k][x1][x2], dp[k-1][i][j] + num[x1][y1] + num[x2][y2]);
45                             }
46                         }
47                 }
48         }
49         printf("%d\n", dp[2*n-1][n][n]);
50     }
51     return 0;
52 }

代码君