这个题实际上求的是半个凸包的核,但只需要枚举直线两两交点和地上的点就可以很快跑过
根据线性规划原理答案应该在直线交点上,为什么还要看地上的点?
因为实际上可行域是有下限的,有可能下限与地线重合
码(注意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);
}