Codeforces 838E IndiaHacks 2nd Elimination 2017 mirror Convex Countour ：毛线DP啊

题意：顺时针给出一个严格凸包的n个点坐标。求最长的一条不交叉的哈密顿回路。

题解：每个哈密顿回路都有两个端点，其中一个作为起点，一个作为终点，就像一根毛线………………在这跟毛线没有结束的时候，他的终点是一个线头，那我们就拿这个线头来DP。dp[ i ][ k ]表是 连接i以及顺时针相邻 总共k个点的最短哈密顿路径。但是这个状态可以顺时针和逆时针拓展，所以考虑把二维变成三维dp[ i ][ k ][ 0 ]表示连接 i 以及顺时针相邻总共k个点，并且线头在逆时针那一端的哈密顿回路最短值。相应dp[ i ][ k ][ 1 ]表示连接这么多点，并且线头在顺时针那一边的哈密顿回路最短值。那么。。。比如说dp[ i ][ k ][ 0 ]就可以由dp[ i+1 ][ k-1 ][ 0 ]+dist( i , i+1)（线头在逆时针这边，继续往下边穿一个点）和dp[i+1 ][ k-1 ][ 1 ]+dist( i , i+1+k-1-1 )（线头在逆时针那边，横穿到左边来）来转移，相应dp[ i ][ k ][ 1 ]同理。

注意：哈密顿回路本身是NP问题，并没有什么好的办法能快速解决。但是如果是凸包且不交叉，那么就满足了DP的性质。因为此时按照顺时针或者逆时针来进行排序，就是这个DP的拓扑序（按照这个顺序进行DP可以保证不交叉），能够消除后效性。

Code：

#include<bits/stdc++.h>
#define sqr(_) (_)*(_)
using namespace std;
const int MAX = 2505;
struct Point{
	double x,y;
	double dist( Point p2){
//		printf("%d %.2lf %.2lf\n",sqr(p2.x-x)+sqr(p2.y-y),sqrt(sqr(p2.x-x)+sqr(p2.y-y)),sqrt(2));
		return sqrt((double)sqr(p2.x-x)+sqr(p2.y-y));
	}
}points[MAX];
int n;
double dp[MAX][MAX][2];
int main(){
	scanf("%d",&n);
	for (int i=0;i<n;i++){
		double x,y;
		scanf("%lf%lf",&x,&y);
		points[i].x=x,points[i].y=y;
	}
	for (int k=2;k<=n;k++){
		for (int i=0;i<n;i++){
			int j = (i+k-2)%n;
			int I = (i-1+n)%n;
			int J = (j+1)%n;
			Point IP = points[I];
			Point JP = points[J];
			dp[I][k][0] = max(dp[i][k-1][0]+IP.dist(points[i]),dp[i][k-1][1]+IP.dist(points[j]));
			dp[i][k][1] = max(dp[i][k-1][1]+JP.dist(points[j]),dp[i][k-1][0]+JP.dist(points[i]));
		}
	}
	double ans=-1;
	for (int i =0;i<n;i++){
		ans = max(ans,max(dp[i][n][0],dp[i][n][1]));
	}
	printf("%.12f\n",ans);
	return 0;
}