【CF838E】 Convex Countour

【CF838E】 Convex Countour

首先观察题目的性质

由于是凸包，因此不自交路径中的一条边\((x, y)\)的两端点只能向与\(x\)或\(y\)相邻的结点连边。

举个栗子，若选取了一条边\((x, y)\)，且假设编号从\(x\)到\(y\)结点已经在一条不自交路径中（不考虑特殊情况），那么向外扩展路径只能连向相邻的点，即只能连边\((x+1, y)\)或\((x, x+1)\)或\((x, y-1)\)或\((y-1, y)\)

很容易用反证法证明。假设连边\((x-2, y)\)，那么点\(x-1\)则无法通过一条不与\((x, y)\)或\((x-2, y)\)相交的路径与其他点连通。而此题路径要覆盖所有点，即所有点之间连通，则矛盾。因此上述结论成立。

由于选取的路径每次只能向外扩展一个点，那么此题就变成了区间动态规划问题。

设\(f_{l, r, 0/1}\)表示区间\([l, r]\)的最长路径长度，\(0\)表示路径终点在\(l\)， \(1\)表示路径终点在\(r\)。

那么可以得到

\[f_{l, r, 0}=max\{f_{l+1, r, 0}+dis(l, l+1), f_{l+1, r, 1}+dis(l, r)\}\\f_{l, r, 1}=max\{f_{l, r-1, 0}+dis(r, l), f_{l, r-1, 1}+dis(r, r-1)\} \]

且易知\(f_{x, x, 0}=f_{x, x, 1}=0\)

此题卡空间，不能开两倍大小，将下标取模后再dp即可

代码如下

#include <cstdio>
#include <cmath>
#include <algorithm>
using namespace std;

const int N=2510;
struct Point {
	double x, y;
	Point(int x=0, int y=0):x(x), y(y){}
} p[N];
double dis(Point a, Point b) {
	return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
double f[N][N][2];
int n;

int main() {
	scanf("%d", &n);
	for (int i=0; i<n; i++) {
		int x, y; scanf("%d%d", &x, &y);
		p[i]=Point(x, y);
	}
	for (int len=2; len<=n; len++)
		for (int l=0; l<n; l++) {
			int r=(l+len-1)%n;
			f[l][r][0]=max(f[(l+1)%n][r][0]+dis(p[l], p[(l+1)%n]), f[(l+1)%n][r][1]+dis(p[l], p[r]));
			f[l][r][1]=max(f[l][(r-1+n)%n][0]+dis(p[r], p[l]), f[l][(r-1+n)%n][1]+dis(p[r], p[(r-1+n)%n]));
		}
	double ans=0;
	for (int i=0; i<n; i++) ans=max(ans, max(f[i][(i+n-1)%n][0], f[i][(i+n-1)%n][1]));
	printf("%.10lf", ans);
	return 0;
}