zoukankan      html  css  js  c++  java
  • HDU 2224 The shortest path

    The shortest path

    Time Limit: 1000ms
    Memory Limit: 32768KB
    This problem will be judged on HDU. Original ID: 2224
    64-bit integer IO format: %I64d      Java class name: Main
     
    There are n points on the plane, Pi(xi, yi)(1 <= i <= n), and xi < xj (i<j). You begin at P1 and visit all points then back to P1. But there is a constraint: 
    Before you reach the rightmost point Pn, you can only visit the points those have the bigger x-coordinate value. For example, you are at Pi now, then you can only visit Pj(j > i). When you reach Pn, the rule is changed, from now on you can only visit the points those have the smaller x-coordinate value than the point you are in now, for example, you are at Pi now, then you can only visit Pj(j < i). And in the end you back to P1 and the tour is over.
    You should visit all points in this tour and you can visit every point only once.
     

    Input

    The input consists of multiple test cases. Each case begins with a line containing a positive integer n(2 <= n <= 200), means the number of points. Then following n lines each containing two positive integers Pi(xi, yi), indicating the coordinate of the i-th point in the plane.
     

    Output

    For each test case, output one line containing the shortest path to visit all the points with the rule mentioned above.The answer should accurate up to 2 decimal places.
     

    Sample Input

    3
    1 1
    2 3
    3 1

    Sample Output

    6.47
    
    Hint: The way 1 - 3 - 2 - 1 makes the shortest path.

    解题:双调TSP

     1 #include <bits/stdc++.h>
     2 using namespace std;
     3 const int INF = 0x3f3f3f3f;
     4 const int maxn = 500;
     5 typedef long long LL;
     6 double dp[maxn][maxn],ds[maxn][maxn];
     7 struct Point {
     8     int x,y;
     9     bool operator<(const Point &t) const {
    10         return x < t.x;
    11     }
    12 } p[maxn];
    13 double calc(int a,int b) {
    14     LL x = p[a].x - p[b].x;
    15     LL y = p[a].y - p[b].y;
    16     double tmp = x*x + y*y;
    17     return sqrt(tmp);
    18 }
    19 
    20 int main() {
    21     int n;
    22     while(~scanf("%d",&n)) {
    23         for(int i = 1; i <= n; ++i)
    24             scanf("%d%d",&p[i].x,&p[i].y);
    25         sort(p,p+n);
    26         for(int i = 1; i <= n; ++i)
    27             for(int j = 1; j <= n; ++j)
    28                 ds[i][j] = calc(i,j);
    29         dp[2][1] = ds[1][2];
    30         for(int i = 3; i <= n; ++i) {
    31             dp[i][i-1] = INF;
    32             for(int j = 1; j < i - 1; j++) {
    33                 dp[i][j] = dp[i-1][j] + ds[i][i-1];
    34                 dp[i][i-1] = min(dp[i][i-1],dp[i-1][j] + ds[j][i]);
    35             }
    36         }
    37         printf("%.2f
    ",dp[n][n-1]+ds[n-1][n]);
    38     }
    39     return 0;
    40 }
    View Code
  • 相关阅读:
    摩托罗拉挪动诉TiVo数字录像手艺专利侵权
    动静称Verizon三月末推出WP7手机
    2月25日中国观点股全线下跌 优酷网涨7.06%
    Verizon CEO称C版iPhone销售微弱
    外来往戏出海案例:热酷日本掘金之路
    摩根大通旗下基金或以4.5亿美元入股Twitter
    TCL通讯2010年净利润6.11亿元 同比增200%
    马克·扎克伯格和Facebook星球
    2月25日中国观点股评级:维持网易买入评级
    美媒评全球十家增速最快IT办事公司 当当网居首
  • 原文地址:https://www.cnblogs.com/crackpotisback/p/4621827.html
Copyright © 2011-2022 走看看