[USACO Dec08]巨大的围栏Largest Fence解题报告

题目

http://cogs.pro/cogs/problem/problem.php?pid=2008

分析

这是道挺有意思的计算几何题……

算法就是：枚举凸包的最左下点B，然后将所有合法的点按相对于B的极角序排序，这样凸包上的点一定是顺序选的。然后记忆化搜索：设状态f(B,p2,p3)代表：最左下点为B，凸包上选的最后一个点为p2，正在测试选不选p3.那么有两种转移：①不选p3，将p3改为p2->p3这条射线沿逆时针旋转卡到的第一个点。②选p3，将p2改为p3，p3改为以p3为起点，方向为p2->p3的射线沿逆时针旋转卡到的第一个点。

这个过程很像是Graham和卷包裹法的某种综合。

“卡到的第一个点”暴力计算是O(N^3)的，那么如果在每次枚举B时都计算一次会是O(N^4)的。不过实际上我们只需要开始计算一下，记忆化搜索中跳过那些不合法的点即可。

ps：我怎么也开始写意识流题解了……身败名裂……

pps：我擦，意识流题解真省事……

代码

<span style="font-size:14px;">#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
using namespace std;
const int SIZEN=310,INF=0x7fffffff/2;
class Point{
public:
	int x,y;
	Point(int _x=0,int _y=0){
		x=_x;
		y=_y;
	}
};
void print(const Point &p){
	cout<<"("<<p.x<<" "<<p.y<<")";
}
bool cmp_coor(const Point &a,const Point &b){
	if(a.y==b.y) return a.x<b.x;
	return a.y<b.y;
}
Point operator - (const Point &a,const Point &b){
	return Point(a.x-b.x,a.y-b.y);
}
int cross(const Point &a,const Point &b){
	return a.x*b.y-b.x*a.y;
}
int N;
Point P[SIZEN];
int from_scan[SIZEN][SIZEN]={0},to_scan[SIZEN][SIZEN]={0};
//直观地想象"一个引出扫描射线的转动枢轴......"例如机械式雷达的波束
void prepare(void){
	for(int i=1;i<=N;i++){
		for(int j=1;j<=N;j++){
			if(i==j) continue;
			//如果没有"三点不共线"的条件,这段可能有潜在问题
			int t1=0,t2=0;
			for(int k=1;k<=N;k++){
				if(k==i||k==j) continue;
				if(cross(P[j]-P[i],P[k]-P[i])>0){
					if(t1==0||cross(P[k]-P[i],P[t1]-P[i])>0) t1=k;
				}
				if(cross(P[j]-P[i],P[k]-P[j])>0){
					if(t2==0||cross(P[k]-P[j],P[t2]-P[j])>0) t2=k;
				}
			}
			from_scan[i][j]=t1;
			to_scan[i][j]=t2;
		}
	}
}
bool valid[SIZEN]={0};
int dp[SIZEN][SIZEN];
int del_step(const int &p1,int p2,int p3){//不选p3
	int x=from_scan[p2][p3];
	while(x&&!(valid[x]&&cross(P[p2]-P[p1],P[x]-P[p1])>0)){
		if(x==p3||x==p1) return 0;
		x=from_scan[p2][x];
	}
	if(!x||x==p3||x==p1) return 0;
	if(cross(P[p3]-P[p2],P[x]-P[p2])>0) return x;
	else return 0;
}
int sel_step(const int &p1,int p2,int p3){//选p3
	int x=to_scan[p2][p3];
	while(x&&!(valid[x]&&cross(P[p3]-P[p1],P[x]-P[p1])>0)){
		if(x==p3||x==p1) return 0;
		x=from_scan[p3][x];
	}
	if(!x||x==p3||x==p1) return 0;
	if(cross(P[p3]-P[p2],P[x]-P[p3])) return x;
	else return 0;
}
int DFS(const int &p1,int p2,int p3){//凸包起点p1,目前钦定了p2,待探查p3
	if(!p2||!p3) return -INF;
	//p1它其实是确定的......
	if(cross(P[p1]-P[p2],P[p3]-P[p2])>0) return -INF;
	int &ans=dp[p2][p3];
	if(ans!=-1) return ans;
	ans=1;//直接结束凸包
	ans=max(ans,DFS(p1,p3,sel_step(p1,p2,p3))+1);
	ans=max(ans,DFS(p1,p2,del_step(p1,p2,p3)));
	return ans;
}
int calc(int d){
	memset(dp,-1,sizeof(dp));
	for(int i=1;i<=N;i++) valid[i]=cmp_coor(P[d],P[i]);
	valid[0]=false;valid[d]=true;
	int ans=0;
	for(int i=1;i<=N;i++){
		for(int j=1;j<=N;j++){
			if(valid[i]&&valid[j]&&cross(P[i]-P[d],P[j]-P[d])>0) ans=max(ans,DFS(d,i,j)+2);
		}
	}
	return ans;
}
void work(void){
	int ans=0;
	for(int i=1;i<=N;i++) ans=max(ans,calc(i));
	printf("%d\n",ans);
}
void read(void){
	scanf("%d",&N);
	for(int i=1;i<=N;i++) scanf("%d%d",&P[i].x,&P[i].y);
}
int main(){
	freopen("fence1.in","r",stdin);
	freopen("fence1.out","w",stdout);
	read();
	prepare();
	work();
	return 0;
}
</span>