百度之星 1002 （动态规划，双调旅行商问题）

                                                                                                          Disk Schedule 

 Problem Description
 有很多从磁盘读取数据的需求，包括顺序读取、随机读取。为了提高效率，需要人为安排磁盘读取。然而，在现实中，这种做法很复杂。我们考虑一个相对简单的场景。
磁盘有许多轨道，每个轨道有许多扇区，用于存储数据。当我们想在特定扇区来读取数据时，磁头需要跳转到特定的轨道、具体扇区进行读取操作。为了简单，我们假设磁头可以在某个轨道顺时针或逆时针匀速旋转，旋转一周的时间是360个单位时间。磁头也可以随意移动到某个轨道进行读取，每跳转到一个相邻轨道的时间为400个单位时间，跳转前后磁头所在扇区位置不变。一次读取数据的时间为10个单位时间，读取前后磁头所在的扇区位置不变。磁头同时只能做一件事：跳转轨道，旋转或读取。
现在，需要在磁盘读取一组数据，假设每个轨道至多有一个读取请求，这个读取的扇区是轨道上分布在 0到359内的一个整数点扇区，即轨道的某个360等分点。磁头的起始点在0轨道0扇区，此时没有数据读取。在完成所有读取后，磁头需要回到0轨道0扇区的始点位置。请问完成给定的读取所需的最小时间。
Input
输入的第一行包含一个整数M（0<M<=100），表示测试数据的组数。
对于每组测试数据，第一行包含一个整数N（0<N<=1000），表示要读取的数据的数量。之后每行包含两个整数T和S（0<T<=1000，0<= S<360），表示每个数据的磁道和扇区，磁道是按升序排列，并且没有重复。
 Output
对于每组测试数据，输出一个整数，表示完成全部读取所需的时间。
 Sample Input
3
1
1 10
3
1 20
3 30
5 10
2
1 10
2 11
 Sample Output
830
4090
1642

解题思路：

dp

双调欧几里得旅行商问题。等价成两个人从起点开始走。

dp[i][j]表示一个人走到第i个位置，一个人走到第j个位置时的最小花费。保证i>j 且前i个位置都走遍了。

dp[i][j]=min(dp[i][j],dp[i-1][j]+dis(i-1,i)) ;要么给快人从i-1走到了i

dp[i][i-1]=min(dp[i][i-1],dp[i-1][j]+dis(j,i); 要么给慢人从j走到了i 

（此图来自网络，感谢作者）

源代码（后续，欢迎大家来帖）：

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cmath>
#define INF 0x3f3f3f3f
const int maxn=1000+10;

using namespace std;

int dp[maxn][maxn],save[maxn];
int dis[maxn][maxn];
int n;

int cal(int a,int b)
{
    return min(abs(save[a]-save[b]),360-abs(save[a]-save[b]));
}



int bitonicTSP()
{
    int i,j;
    //for(i=1;i<=n;i++)
    dp[1][0]=dis[0][1];  //init

    for(i=1;i<n;i++)
    {
        dp[i+1][i]=INF;
        for(j=0;j<i;j++)
        {
            //dp[i+1][j]=dp[i][j]+dis[i][i+1];
            dp[i+1][j]=dp[i][j]+dis[i][i+1];
            dp[i+1][i]=min(dp[i+1][i],dp[i][j]+dis[j][i+1]);
        }
    }
    return dp[n][n-1]+dis[n-1][n];

}



int main()
{
   int T,Max;
   scanf("%d",&T);
   while(T--)
   {
       scanf("%d",&n);
       for(int i=1;i<=n;i++) //离散化
       {
           int a,b;
           scanf("%d%d",&a,&b);
           save[i]=b;
           Max=a;
       }
       save[0]=0;
       for(int i=0;i<=n;i++) //预处理出距离
       {
           for(int j=i+1;j<=n;j++)
               dis[i][j]=cal(i,j);

       }
       cout<<bitonicTSP()+Max*400*2+10*n<<endl;

   }
   return 0;


}

（终于可以贴上自己的代码了，2014/5/20 11：42，表示终于看懂了，TSP）