ZOJ 2614 Bridge[三分求解]

题目链接：http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=2614

THE 30th ACM/ICPC ASIA REGIONAL 2005 HANGZHOU SITE

题目大意：给出桥的总长度B，线缆的总长度L，桥顶的高度H和桥东之间的最大距离D，

问线缆离桥面的最小距离：

解题方法：三分求解

涉及的知识：抛物线的面积与弧长公式；

面积 Area=2ab/3

　　弧长 Arc length ABC

=√(b^2+16a^2 )/2+b^2/8a ln((4a+√(b^2+16a^2 ))/b)

 

#include

#include

#include

#include

#include

const double eps=1e-9;

using namespace std;

double temp,len;

//求弧长

double cal(double a)

{

    double x=0;

    double y=sqrt(len*len+16*a*a);

    x=y/2+len*len/(8*a)*log((4*a+y)/len)-temp;

         return fabs(x);

}//三分求解

double div(double H)

{

   double left=eps,right=H;

   double ans1,ans2;

   double m,mm;

   do

   {

       m=(left+right)/2.0;

       mm=(right+m)/2.0;

       ans1=cal(m);

       ans2=cal(mm);

       if(ans1

           right=mm;

       else

           left=m;

   }while(fabs(ans1-ans2)>eps);

   return mm;

}

int main(){

    //freopen("C:\\Users\\acm\\Desktop\\001.in","r",stdin);

    int k=1,t;

    double D,H,B,L;

    scanf("%d",&t);

    while(t--){

        scanf("%lf%lf%lf%lf",&D,&H,&B,&L);

        len=ceil(1.0*B/D);

        temp=L/len;

        len=B/len;

        if(k!=1){

            printf("\n");

        }

        printf("Case %d:\n",k++);

        printf("%.2lf\n",H-div(H));

    }

    return 0;

}