timus 1874. Football Goal URAL 解題報告 計算幾何

題目大意:

想踢足球,沒有球門,找了一顆足夠高的樹,還有兩個木棍,組成一個球門,兩根門柱是大樹和一個木棍,橫樑是另一根木棍……  並且他們很有想象力,把球門組成了一個四邊形就好,沒必要是矩形或梯形等。這個從測試數據就可以看出來……   爲了的得分最大,想要球門面積最大!!  球門面積大得分就多嘛! 把中國隊放在上面再大的球門也白搭!



言歸正轉,這個題的解決不是我當時想出來的,我當時想這個圖形的面積可能和一個角度有關,然後二分這個角度利用面積和這個角度的關係二分,但是很遺憾沒有發現這個關係……

又來查找解題報告說是是枚舉邊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;
}


發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章