BZOJ3707 圈地

我们可以枚举每条线，用离这个线段最近的两个点更新答案，考虑如何快速找到最近的两个点，我们把所有线段按斜率排序，把所有点按到第一个点的距离排序（如果点在线段左边那么距离算负的）然后按斜率枚举每个线段，那么从一个线段l1转到另一个线段l2时，和l2的距离的相对顺序可能改变的点只有l1的两个端点，因为如果有别的点对相对顺序改变，那么这个点对的斜率一定介于l1和l2之间

然后就是n^2的了

#include<iostream>
#include<cstdlib>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<ctime>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<stack>
#include<queue>
#include<map>
#include<set>
#include<bitset>
using namespace std;
#define MAXN 1010
#define MAXM 1010
#define ll long long
#define INF 1000000000
#define MOD 1000000007
#define eps 1e-8
struct pt{
	double x;
	double y;
	int num;
	pt(){
		
	}
	pt(double _x,double _y){
		x=_x;
		y=_y;
	}
};
struct line{
	pt x;
	pt y;
	double k;
	double b;
	line(){
		
	}
	line(pt x,pt y);
	friend bool operator <(line x,line y){
		return x.k<y.k;
	}
};
pt p[MAXN];
line l[MAXN*MAXN];
int pos[MAXN];
int n;
int tot;
double xl(pt x,pt y){
	if(x.x==y.x){
		return 1e60;
	}
	return (x.y-y.y)/(x.x-y.x);
}
line::line(pt _x,pt _y){
	x=_x;
	y=_y;
	k=xl(x,y);
	b=x.y-k*x.x;
}
double dis(pt x,pt y){
	return sqrt((x.x-y.x)*(x.x-y.x)+(x.y-y.y)*(x.y-y.y));
}
double dis(pt x,line y){
	double K=-1.0/y.k;
	double B=x.y-K*x.x;
	double X=(B-y.b)/(y.k-K);
	double Y=K*X+B;
	if(x.x*y.k+y.b<=x.y){
		return dis(pt(X,Y),x);
	}else{
		return -dis(pt(X,Y),x);
	}
}
int now;
bool cmp(pt x,pt y){
	return dis(x,l[1])<dis(y,l[now]);
}
int main(){
	int i,j;
	scanf("%d",&n);
	for(i=1;i<=n;i++){
		scanf("%lf%lf",&p[i].x,&p[i].y);
		p[i].num=i;
		for(j=1;j<i;j++){
			l[++tot]=line(p[i],p[j]);
		}
	}
	sort(l+1,l+tot+1);
	now=1;
	sort(p+1,p+n+1,cmp);
	for(i=1;i<=n;i++){
		pos[p[i].num]=i;
	}
	double ans=1e60;
	for(i=1;i<=tot;i++){
		int x=min(pos[l[i].x.num],pos[l[i].y.num]);
		int y=max(pos[l[i].x.num],pos[l[i].y.num]);
		if(x!=1){
			ans=min(ans,fabs(dis(p[x-1],l[i]))*dis(l[i].x,l[i].y)/2);
		}
		if(y!=n){
			ans=min(ans,fabs(dis(p[y+1],l[i]))*dis(l[i].x,l[i].y)/2);
		}
		if(i!=tot){
			if(dis(p[x],l[i+1])>dis(p[y],l[i+1])){
				swap(p[x],p[y]);
				swap(pos[l[i].x.num],pos[l[i].y.num]);
			}
		}
	}
	printf("%.2lf\n",ans);
	return 0;
}

/*
3
0 0
0 1
1 0

*/