UVa 116 - Unidirectional TSP(动态规划)

紫书上的原题，第270页。

书上已经讲解的很详细了。也没啥好说的。

首先dp(i, j) 表示从格子(i, j)出发到最后一列的开销。 这样我们就确定了状态了。

然后根据逆推，从后往前倒着走，得到动态规划方程 dp(i, j) =min { dp (i,j+1) , dp(i, j) , dp(i, j-1) } + 格子(i, j) 的值。

然后这道题就没有什么难点了。

也提示我思考问题不要太那啥了。

上代码：

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

const int maxn = 100 + 5;
const int INF = 1000000000;
int m, n, a[maxn][maxn], d[maxn][maxn], next[maxn][maxn]; // next 用于保存路径 

int main() {
  while(scanf("%d%d", &m, &n) == 2 && m) {
    for(int i = 0; i < m; i++)
      for(int j = 0; j < n; j++)
        scanf("%d", &a[i][j]);

    int ans = INF, first = 0;
    for(int j = n-1; j >= 0; j--) {    // 从最后一列往前倒着走 
      for(int i = 0; i < m; i++) {
        if(j == n-1) d[i][j] = a[i][j];  // 最后一列的初始开销就是自己的值 
        else {
          int rows[3] = {i, i-1, i+1};
          if(i == 0) rows[1] = m-1;
          if(i == m-1) rows[2] = 0; // 根据题目要求 
          sort(rows, rows+3);   // 重新排序，确保字典序最小 
          d[i][j] = INF;
          for(int k = 0; k < 3; k++) {
            int v = d[rows[k]][j+1] + a[i][j];
            if(v < d[i][j]) { d[i][j] = v; next[i][j] = rows[k]; }
          }
        }
        if(j == 0 && d[i][j] < ans) { ans = d[i][j]; first = i; }   
      }
    }
    printf("%d", first+1);
    for(int i = next[first][0], j = 1; j < n; i = next[i][j], j++) printf(" %d", i+1);
    printf("\n%d\n", ans);
  }
  return 0;
}

水波。