POJ2677 DP tour 双调欧几里得旅行商问题

Tour

Time Limit: 1000MS		Memory Limit: 65536K
Total Submissions: 2699		Accepted: 1193

Description

John Doe, a skilled pilot, enjoys traveling. While on vacation, he rents a small plane and starts visiting beautiful places. To save money, John must determine the shortest closed tour that connects his destinations. Each destination is represented by a point in the plane pi = < xi,yi >. John uses the following strategy: he starts from the leftmost point, then he goes strictly left to right to the rightmost point, and then he goes strictly right back to the starting point. It is known that the points have distinct x-coordinates. Write a program that, given a set of n points in the plane, computes the shortest closed tour that connects the points according to John's strategy.

Input

The program input is from a text file. Each data set in the file stands for a particular set of points. For each set of points the data set contains the number of points, and the point coordinates in ascending order of the x coordinate. White spaces can occur freely in input. The input data are correct.

Output

For each set of data, your program should print the result to the standard output from the beginning of a line. The tour length, a floating-point number with two fractional digits, represents the result. An input/output sample is in the table below. Here there are two data sets. The first one contains 3 points specified by their x and y coordinates. The second point, for example, has the x coordinate 2, and the y coordinate 3. The result for each data set is the tour length, (6.47 for the first data set in the given example).

Sample Input

3
1 1
2 3
3 1
4
1 1
2 3
3 1
4 2

Sample Output

6.47
7.89

Source

Southeastern Europe 2005

思路【转】：

  欧几里得旅行商问题是对平面上给定的n个点确定一条连接各点的最短闭合旅程的问题。如图（a）给出了一个7个点问题的解。这个问题的一般形式是NP完全的，故其解需要多于多项式的时间。

J.L. Bentley 建议通过只考虑双调旅程(bitonic tour)来简化问题,这种旅程即为从最左点开始，严格地从左到右直至最右点，然后严格地从右到左直至出发点。下图(b)显示了同样的7个点的最短双调路线。在这种情况下，多项式的算法是可能的。事实上，存在确定的最优双调路线的O(n*n)时间的算法。

 图a            图b

注:在一个单位栅格上显示的平面上的七个点。 a)最短闭合路线,长度大约是24.89。这个路线不是双调的。b)相同点的集合上的最短双调闭合路线。长度大约是25.58。

这是一个算导上的思考题15-1。

首先将给出的点排序，关键字x，重新编号，从左至右1，2，3，…，n。

定义p[i][j]，表示结点i到结点j之间的距离。

定义d[i][j]，表示从i连到1，再从1连到j，（注意，i>j，且并没有相连。）

对于任意一个点i来说，有两种连接方法，一种是如图（a）所示，i与i-1相连，另一种呢是如图（b），i与i-1不相连。

根据双调旅程，我们知道结点n一定与n相连，那么，如果我们求的d[n][n-1]，只需将其加上p[n-1][n]就是最短双调闭合路线。

根据上图，很容易写出方程式：

dp[i][j]=dp[i-1][j]+dist[i][i-1];

dp[i][i-1]=min(dp[i][i-1],dp[i-1][j]+dist[j][i]);

View Code

 1 #include 
 2 #include 
 3 #include 
 4 #include 
 5 #include <string.h>
 6 using namespace std;
 7 struct Point
 8 {
 9     int x,y;
10     void input()
11     {
12         scanf("%d%d",&x,&y);
13     }
14 }p[101];
15 double Dis(Point p1,Point p2)
16 {
17     return sqrt((double)(p2.x-p1.x)*(p2.x-p1.x)+(double)(p2.y-p1.y)*(p2.y-p1.y));
18 }
19 bool cmp(Point m,Point n)
20 {
21     return m.x<n.x;
22 }
23 int main()
24 {
25     int n;
26     double a[101][101],dp[101][101];
27     while(scanf("%d",&n)!=EOF)
28     {
29         p[0].input();
30         if(n==1) {printf("0.00\n");continue;}
31         for(int i=1;i)
32         {
33             p[i].input();
34             if(p[0].x>p[i].x) swap(p[0],p[i]);
35         }
36         sort(p,p+n,cmp);
37         for(int i=0;i)
38         {
39             a[i][i]=0;
40             for(int j=i+1;j)
41             {
42                a[i][j]=a[j][i]=Dis(p[i],p[j]);
43             }
44         }
45         dp[0][0]=0;
46         for(int i=1;i)
47         dp[i][0]=a[i][0];
48         for(int i=1;i1;i++)
49         {
50             dp[i+1][i]=1000000000;
51             for(int j=0;j<=i-1;j++)
52             {
53                 dp[i+1][j]=dp[i][j]+a[i][i+1];
54                 dp[i+1][i]=min(dp[i+1][i],dp[i][j]+a[j][i+1]);
55             }
56         }
57         printf("%.2lf\n",dp[n-1][n-2]+a[n-1][n-2]);
58     }
59     return 0;
60 }

 

转载于:https://www.cnblogs.com/-sunshine/archive/2012/07/23/2605251.html