bzoj1038

想法很简单

在每一条线段上三分横坐标

然后一定是对于左右的底层端点连到这个横坐标，然后求出高度

可以想到这个三分肯定是对的，因为在一条线段上，一定是个凸函数

但是精度被卡成傻子

在洛谷上eps开1e-9，80分，开1e-10,90分。。。

然后不用struct point 的向量写法，直接坐标运算就过了。

太毒瘤了。

#include<bits/stdc++.h>
#define maxl 310
#define eps 1e-12
using namespace std;

inline int sgn(double x)
{
	if(x>-eps && x<eps) return 0;
	if(x>0) return 1;
	else	return -1;
}

int n;
double x[maxl],y[maxl];
double ans=1e12;

inline void prework()
{
	int d;
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
		scanf("%d",&d),x[i]=d;
	for(int i=1;i<=n;i++)
		scanf("%d",&d),y[i]=d;
}

inline double jug(int i,double xi)
{
	double mx=0;
    double len = y[i]+(y[i+1]-y[i])/(x[i+1]-x[i])*(xi-x[i]);
    for(int j=1;j<i;j++)
	if(sgn(y[j+1]-y[j])>0)
    {
        double h = y[j]+(y[j+1]-y[j])/(x[j+1]-x[j])*(xi-x[j]);
        mx=max(mx,h-len);
    }
    for(int j=i+1;j<n;j++)
    if(sgn(y[j]-y[j+1])>0)
    {
    	double h=y[j+1]+(y[j]-y[j+1])/(x[j+1]-x[j])*(x[j+1]-xi);
    	mx=max(mx,h-len);
	}
    return mx;
}

inline void mainwork()
{
	double mxl,mxr;
	for(int i=1;i<n;i++)
	{
		double len,l=x[i],r=x[i+1],midl,midr;
		while(fabs(r-l)>1e-9)
		{
			len=(r-l)/3.0; 
			midl=l+len;
			midr=r-len;
			mxl=jug(i,midl);
			mxr=jug(i,midr);
			if(mxl>mxr)
				l=midl;
			else
				r=midr;
		}
		mxl=jug(i,l);
		ans=min(mxl,ans);
	}
}

inline void print()
{
	if(sgn(ans)<0)
		ans=eps;
	if(n==1)
		ans=eps;
	printf("%.3f",ans+eps);
}

int main()
{
	prework();
	mainwork();
	print();
	return 0;
}