zhxnb!!!
恩这道题是个数学题
以最下面和最左边两条边建系
可以用向量叉乘算出点(b,a)到矩形靠左c长为L的边的距离,本题实际上是当L边在坐标轴上动时求该距离的最小值
距离函数是一个单峰的(可以直接从图形中看出),用三分来解决
实型数据杀我……精度控制太难了errrrrrrrrrrrrr注意:三分的步长小于eps,否则永远跳不到停止那一步(永不停息!! )
郑tc讲课的时候说可以先*1e7变成longlong再来处理最后/1e7,精度更高
上代码~
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
double b,a,l,n;
double w,eps=0.000000001,d=0.000000000009;
double f(double m){
n = sqrt(l*l-m*m);
return ((a-n)*m+b*n)/l;
}
double getw(){
double L=-1.000000,R=1005.0000000;
while(fabs(R-L)>eps){
double mid = (L+R)/2;
if(f(mid-d)>f(mid+d)) L=mid-d;
else R=mid+d;
}
return f(L);
}
int main(){
freopen("b.in","r",stdin);
freopen("b.out","w",stdout);
scanf("%lf%lf%lf",&a,&b,&l);
w = getw();
if(l<=a) w = max(w,b);
if(l<=b) w = max(w,a);
if(w<0) {
printf("My poor head =(\n");
return 0;
}
if(w>l) w=l;
printf("%.7lf\n",w);
return 0;
}
依旧有精度问题
一时半会是A不了了
20191104
蒟蒻的csp倒计时11天