Description
Input
Output
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
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]);
![](https://i-blog.csdnimg.cn/blog_migrate/81178cc93a2a3bb5048d90d76e7ec935.gif)
1 #include <algorithm> 2 #include <iostream> 3 #include <cmath> 4 #include <stdio.h> 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<n;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<n;i++) 38 { 39 a[i][i]=0; 40 for(int j=i+1;j<n;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<n;i++) 47 dp[i][0]=a[i][0]; 48 for(int i=1;i<n-1;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 }