bzoj1038 [ZJOI2008]瞭望塔 枚举

这个题实际上求的是半个凸包的核，但只需要枚举直线两两交点和地上的点就可以很快跑过

根据线性规划原理答案应该在直线交点上，为什么还要看地上的点？

因为实际上可行域是有下限的，有可能下限与地线重合

码（注意eps和运算中的括号）：

#include<iostream>
#include<cstdio>
using namespace std;
long double ans=99999999999999,X,Y,K1,K2,B1,B2,lin;
#define eps 0.000000001
int n,i,j,k;
bool ky;
double lliinn;
struct dian
{
long	double x,y;
}D[400];
int main()
{
	scanf("%d",&n);
	for(i=1;i<=n;i++)scanf("%lf",&lliinn),D[i].x=lliinn;
	for(i=1;i<=n;i++)scanf("%lf",&lliinn),D[i].y=lliinn;
	for(i=1;i<=n;i++)
	{
		lin=0;
	 for(k=1;k<n;k++)
	 {
   	 K1=(D[k].y-D[k+1].y)/(D[k].x-D[k+1].x);
	 B1=D[k].y-K1*D[k].x;  
	 lin=max(lin,K1*D[i].x+B1-D[i].y);	 	
	 if(lin+eps>=ans)break;	
	 }		
	ans=min(ans,lin);	
	}	
	for(i=1;i<n;i++)
	for(j=i+1;j<n;j++)//都是和后一个点 
	{
	 K1=(D[i].y-D[i+1].y)/(D[i].x-D[i+1].x);
	 B1=D[i].y-K1*D[i].x;
	 K2=(D[j].y-D[j+1].y)/(D[j].x-D[j+1].x);
	 B2=D[j].y-K2*D[j].x;
	 if(K1+eps>K2&&K1-eps<K2)continue;	
	 X=(B2-B1)/(K1-K2);	
	 Y=K2*X+B2;
	 lin=999999999999;
	 for(k=1;k<n;k++)
	 {
	 	if(D[k].x<=X&&X<=D[k+1].x)
	 	{
	K1=(D[k].y-D[k+1].y)/(D[k].x-D[k+1].x);
	 B1=D[k].y-K1*D[k].x;
	lin=Y-(K1*X+B1);
		}	
	 }
	 if(lin>=ans||lin<0)continue;
	 ky=1;
	 for(k=1;k<n;k++)
	 {
   	 K1=(D[k].y-D[k+1].y)/(D[k].x-D[k+1].x);
	 B1=D[k].y-K1*D[k].x;  	 	
	 if(K1*X+B1>Y+eps){ky=0;break;}	
	 }	
		if(ky)
		{
			ans=min(ans,lin);			
		}		
	}
	lliinn=ans;
	printf("%.3lf",lliinn);
}