題目大意:
想踢足球,沒有球門,找了一顆足夠高的樹,還有兩個木棍,組成一個球門,兩根門柱是大樹和一個木棍,橫樑是另一根木棍…… 並且他們很有想象力,把球門組成了一個四邊形就好,沒必要是矩形或梯形等。這個從測試數據就可以看出來…… 爲了的得分最大,想要球門面積最大!! 球門面積大得分就多嘛! 把中國隊放在上面再大的球門也白搭!
言歸正轉,這個題的解決不是我當時想出來的,我當時想這個圖形的面積可能和一個角度有關,然後二分這個角度利用面積和這個角度的關係二分,但是很遺憾沒有發現這個關係……
又來查找解題報告說是是枚舉邊AB,然後AB確定之後,三角形ABD面積就確定了,用海倫_秦九韶公式去運算; 至於三角形ABO,由於其是個直角三角形,AB邊長度確定之後,AO,BO可以變化,再加上勾股定理,求其面積的時候利用消元得到一個關於AB的一元二次方程,其最值是AO=BO的時候得到的,所以三角巷AOB的面積最大值時z*z/4;
三角形ABD的面積用海倫公式知道是關於Z的一元二次方程,具體很難化簡到一般形式,所以用三分求面積,因爲他是個凸型函數!
z是變長,範圍是[0,x+y].
參考博客:http://blog.csdn.net/lwbaptx/article/details/7175161
#include <stdio.h>
#include <stdlib.h>
#include<math.h>
#include<iostream>
#define MAXSIZE 100
#define EPS 1e-9
const double pi=acos(-1.0);
using namespace std;
double x,y;
double gets(double z)
{
double p=(x+y+z)/2.0;
return z*z/4.0+sqrt(p*(p-x)*(p-y)*(p-z));
}
int main()
{
while(cin>>x>>y)
{
double l=0,r=(x+y);
double ans=0;
while(fabs(r-l)>EPS)
{
double mid1=(l+r)/2.0;
double mid2=(r+mid1)/2.0;
double s1=gets(mid1);
double s2=gets(mid2);
if(s1>s2)
{
r=mid2;
}else l=mid1;
ans=s1;
}
printf("%.9f\n",ans+EPS);
}
return 0;
}