[SPOJ PUTNIK]区间DP_putnik 区间dp-CSDN博客

本文链接：https://blog.csdn.net/acmore_xiong/article/details/54003047

[SPOJ PUTNIK]区间DP

知识点：dynamic programing uncommon TSP

1. 题目链接

2. 题意描述

旅行商问题TSP( Travelling Salesman Problem)概都清楚，经典的NP难问题之一。这个题目是一个变式的TSP问题。
经典的TSP问题是指：城市与城市之间有边权，求经过所有城市有且仅有一次，最小边权之和。
这个题目增加了一个约束条件：对于标号为 $K$ 的城市，所有标号小于 $K$ 的城市，要么在 $K$ 之前被访问，要么都在 $K$ 之后被访问。
城市数目 $2\le n\le 1500$ ，从 $i$ 到 $j$ 的边权 $cost[i][j] \in[0, 1000]$ .

3. 解题思路

题目关键之所在，在于分析这个题目的特殊性。对于标号为 $K$ 的城市，所有标号小于 $K$ 的城市，要么在 $K$ 之前被访问，要么都在 $K$ 之后被访问。对于整个访问序列，可以这么理解，先放好标号为 $1$ 的城市，接下来放编号为 $k+1$ 的城市，它一定是只能放在当前序列的最左端或者最右端，整个序列可以看成从 $1$ 开始向两端拓展【也可以像其中一端拓展】。
然后，用 $dp[i][j]$ 表示起点为 $i$ 终点为 $j$ 的最小访问边权之和。从 $2$ 到 $n$ 开始枚举顶点。有如下转移方程：

d p [u] [i] = min {d p [u - 1] [i] + c o s t [u] [u - 1]} d p [i] [u] = min {d p [i] [u - 1] + c o s t [u] [u - 1]}} i \in [1, u - 1] d p [u] [u - 1] = min {d p [i] [u - 1] + c o s t [u] [i]} d p [u - 1] [u] = min {d p [u - 1] [i] + c o s t [i] [u]}} i \in [1, u - 1]

$\left. \begin{array}{l} dp[u][i]=\min\{dp[u-1][i]+cost[u][u-1]\}\\ dp[i][u]=\min\{dp[i][u-1]+cost[u][u-1]\} \end{array} \right\} i\in [1,u-1]\\ \left. \begin{array}{l} dp[u][u-1]=\min\{dp[i][u-1]+cost[u][i]\}\\ dp[u-1][u]=\min\{dp[u-1][i]+cost[i][u]\} \end{array} \right\} i\in [1,u-1]$

4. 实现代码

#include <bits/stdc++.h>
using namespace std;

typedef long long LL;

const int MAXN = 1500 + 5;
const int INF = 0x3f3f3f3f;

int n, cost[MAXN][MAXN];
int dp[MAXN][MAXN];

int main() {
#ifdef ___LOCAL_WONZY___
    freopen("input.txt", "r", stdin);
#endif // ___LOCAL_WONZY___
    while(~scanf("%d", &n)) {
        for(int i = 1; i <= n; ++i) {
            for(int j = 1; j <= n; ++j) {
                scanf("%d", &cost[i][j]);
                dp[i][j] = INF;
            }
        }
        dp[1][1] = dp[1][1] = 0;
        dp[1][2] = dp[2][1] = cost[1][2];
        for(int u = 3; u <= n; ++u) {
            for(int i = 1; i < u - 1; i++) {
                dp[u][u - 1] = min(dp[u][u - 1], dp[i][u - 1] + cost[u][i]);
                dp[u - 1][u] = min(dp[u - 1][u], dp[u - 1][i] + cost[i][u]);
                dp[u][i] = min(dp[u][i], dp[u - 1][i] + cost[u][u - 1]);
                dp[i][u] = min(dp[i][u], dp[i][u - 1] + cost[u - 1][u]);
            }
        }
        int ans = INF;
        for(int i = 1; i < n; i++) {
            ans = min(ans, dp[i][n]);
            ans = min(ans, dp[n][i]);
        }
        printf("%d\n", ans);
    }
    return 0;
}