HDU 2289 Cup（可以二分法，但是除了它的一半？）

这道题目。运营商做数学题？算上两个子主题做？顶多算一个水主要议题...

首先，没有实际的二分法，但是，我们发现了一个新的解决方案，以取代二分法。

若果按照i从0，每次添加0.00000001我一直枚举h如果，绝对逼暂停。枚举太大

然而，它可以分为两个步骤：

#include<cstdio>

#include<cmath>

#define pai acos(-1.0)

double r1,r2,h,v;

double get_v(double temp)

{

    double rr=r1+(r2-r1)*temp/h;

    return pai*temp*(r1*r1+rr*rr+r1*rr)/3;

}

int main()

{

    int t;

    scanf("%d",&t);

    while(t--)

    {

        scanf("%lf%lf%lf%lf",&r1,&r2,&h,&v);

        double i;

        for(i=0;i<100;i+=0.0001)

        {

            if(get_v(i)>v)

                break;

        }

        double j;

        for(j=i-0.0001;j<=i;j+=0.00000001)

        {

            if(get_v(j)>v)

                break;

        }

        printf("%.6lf\n",j);

    }

    return 0;

}

尽管时间久了一点，单这道题的数据是能够过得。

自己想到方法解决这个问题，非常开心。

以下是迪神一下子就想到的二分解法，说来迪神也真是牛逼……在我们这样的菜鸟云集的地方。应该算的上是鹤立鸡群了。。。

#include<iostream>

#include<algorithm>

#include<cstring>

#include<cstdio>

#include<cmath>

using namespace std;

#define PI acos(-1.0)

double r,R,H,v,m,x,y;

int L()

{

    if((PI * m * (r*r + r*((R-r)/H*m+r) + ((R-r)/H*m+r)*((R-r)/H*m+r))/3-v)>1e-7) return 1;

    else return 0;

}

int main()

{

    int T;

    cin>>T;

    while(T--)

    {

        cin>>r>>R>>H>>v;

        x=0,y=H;

        while((y-x)>1e-7)

        {

            m=x+(y-x)/2;

            if(L()) y=m;

            else x=m;

        }

        printf("%.6lf\n",x);

    }

}

就是这个样子，打电话久了会困……它被证明是真实的。。。