Disk Schedule HDU - 4824 (双调欧几里得旅行商问题)

有很多从磁盘读取数据的需求，包括顺序读取、随机读取。为了提高效率，需要人为安排磁盘读取。然而，在现实中，这种做法很复杂。我们考虑一个相对简单的场景。 
磁盘有许多轨道，每个轨道有许多扇区，用于存储数据。当我们想在特定扇区来读取数据时，磁头需要跳转到特定的轨道、具体扇区进行读取操作。为了简单，我们假设磁头可以在某个轨道顺时针或逆时针匀速旋转，旋转一周的时间是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

很明显双调欧几里得旅行商问题，其实只要理解了原理，这类题也就是个模板题，其中原理紫书269页讲的很清楚，这里就不复述了，知识上代码看一下模板和具体实现方法。

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cmath>
using namespace std;
int dp[1010][1010];
struct node{
	int x;
	int y;
};
node op[1010];
//int dis(op[a],op[b]);
int main()
{
	int T;
	scanf("%d",&T);
	while(T--)
	{
		int n;
		scanf("%d", &n);
		op[0].x = op[0].y = 0;
		for(int i=1; i<=n; i++)
		scanf("%d %d", &op[i].x,&op[i].y);
		dp[1][0] = dp[0][1] = op[1].x*400 + min(op[1].y,360-op[1].y);
		for(int i=2; i<=n; i++)
		{
			dp[i][i-1] = 100000000;
			for(int j=0; j<i-1; j++)
			{
				dp[i][i-1] = min(dp[i][i-1], dp[i-1][j] + (op[i].x-op[j].x)*400 + min(abs(op[j].y-op[i].y),360-abs(op[j].y-op[i].y)));
				dp[i][j] = dp[i-1][j] + (op[i].x-op[i-1].x)*400 + min(abs(op[i-1].y-op[i].y),360-abs(op[i-1].y-op[i].y));
			}
		}
		int cc = 100000000;
		for(int i=0; i<=n-1; i++)
		{
			cc = min(cc, dp[n][i] + (op[n].x-op[i].x)*400 + min(abs(op[i].y-op[n].y),360-abs(op[n].y-op[i].y)) );
		}
		printf("%d\n",cc+10*n);
	}
	return 0;
}

水波.