【jzoj 3510】 【NOIP2013模拟11.5B组】最短路径 {动态规划}

题目

Description
平面内给出 n 个点，记横坐标最小的点为 A，最大的点为 B，现在Zxd想要知道在每个点经过一次（A 点两次）的情况下从 A 走到 B，再回到 A 的最短路径。但他是个强迫症患者，他有许多奇奇怪怪的要求与限制条件：
1. 从 A 走到 B 时，只能由横坐标小的点走到大的点。
2. 由 B 回到 A 时，只能由横坐标大的点走到小的点。
3. 有两个特殊点 b1 和 b2， b1 在 0 到 n-1 的路上，b2 在 n-1 到 0 的路上。
请你帮他解决这个问题助他治疗吧！

Input
第一行三个整数 n,b1,b2,( 0 < b1，b2 < n-1 且 b1 <> b2)。n 表示点数，从 0 到 n-1 编号，b1 和 b2 为两个特殊点的编号。
以下 n 行，每行两个整数 x、y 表示该点的坐标(0 <= x，y <= 2000)，从 0 号点顺序
给出。Doctor Gao为了方便他的治疗，保证所有点横坐标不同，并且已经将给出的点按 x 增序排好了。

Output
仅一行，输出最短路径长度（精确到小数点后面 2 位）

---

解题思路

比赛时，想跑最短路（spfa），但限制太多，显然不是正解。正解是动态规划（需要注意一下题目的判断和细节）。

引：

(定义$dis[i][j]$表示点i到点j之间的直线距离)
把问题简化 后可以发现，问题其实可以看成这样一个问题：平面上n个点，确定一条连接各点的最短闭合旅程，这种旅程即为从最左点开始，严格地从左到右直至最右点，然后严格地从右到左直至出发点。
然后不妨把问题中从n到1路线看成是从1到n的路线，这样我们就有了两条从1到n的路线。
定义$f[i][j]$表示从点i到点j的路径，理解为：从点i开始，从右到左一直到点1，然后从左到右一直到点j。在这个路径上，会经过点1到点$max(i,j)$之间的所有点且只经过一次。
可以列出状态转移方程为：
定义k为$max(i,j)+1;$
$f[1][1]=0;$
$f[i][k]=min(f[i][k],f[i][j]+dis[j][k]);$
$f[k][j]=min(f[k][j],f[i][j]+dis[i][k]);$
还要考虑一些特殊情况，比如：
　　点b1,b2。
　　还有$k=max(i,j)+1;(i==n||j==n)$的情况。

---

代码

#include<cstdio> 
#include<algorithm>
#include<cmath>
#include<cstring>
#define fre(x) freopen(#x".in","r",stdin),freopen(#x".out","w",stdout);
using namespace std; 
int n,qq,ww,x[1001],y[1001]; double f[1001][1001],dis[1001][1001]; 
int cf(int x) { return x*x; }
int main()
{
//  fre(path); 
    memset(f,0x7f,sizeof(f)); 
    scanf("%d%d%d",&n,&qq,&ww); qq++; ww++; 
    for (int i=1;i<=n;i++) 
     scanf("%d%d",&x[i],&y[i]);
    for (int i=1;i<=n;i++)
     for (int j=1;j<=n;j++)
     dis[i][j]=sqrt(cf(x[i]-x[j])+cf(y[i]-y[j]));//求两点的直线距离
    f[1][1]=0; 
    for (int i=1;i<=n;i++)
     for (int j=1;j<=n;j++)
      if (i!=j||i==1)
      {
        int k=max(i,j)+1; 
        if (k-1==n) 
          {
            if (i!=n) f[n][n]=min(f[n][n],f[i][j]+dis[i][n]); 
            if (j!=n) f[n][n]=min(f[n][n],f[i][j]+dis[j][n]); 
            continue;
          }
        if (k!=qq) f[i][k]=min(f[i][k],f[i][j]+dis[j][k]); 
        if (k!=ww) f[k][j]=min(f[k][j],f[i][j]+dis[i][k]); 
      }
    printf("%.2lf",f[n][n]); 
}