HDU 5432 Pyramid Split(浮点数二分查找)

http://acm.hdu.edu.cn/showproblem.php?pid=5432

题意：
小明有n个金字塔，已知它们的高和宽，小明打算重铸，想将它重铸成两个体积一样的物体。
他现在把所有金字塔正放(即底面贴地放)在水平地面上，用屠龙刀切割一个平面，该平面与水平面平行，称之为割平面。
我们的任务是找到一个这样的割平面，使得这个割平面上方的体积之和等于下方的体积之和，该割平面称之为平均割平面。求平均割平面的高度。

解题思路：
1.论读题的重要性。题目所求的是将这几个金字塔放在水平桌上，然后找到某条线后计算该条线以上的体积是否和该线以下的体积相等，而不是对于每个金字塔都找到这条线在求这条线的平均值。
2.注意的是有的金字塔的总高度可能小于画的这条线，所以计算体积的时候要注意。
3.利用二分查找。

#include 
#include 
#include 
#include 
#include 
using namespace std;
#define EPS 1e-6
const int N = 10000+10;
int h[N],n,a[N];
double cal(int x) 
{
    double ans = 0;
    for(int i = 1;i<=n;i++)
        if(x<=h[i]) 
             ans += 1.0*a[i]*a[i]*(h[i]-x)*(h[i]-x)*(h[i]-x)/(1.0*h[i]*h[i]);
    return ans;
} 

double binsearch(double l,double r,double aim)
{
    while(fabs(l-r)>=EPS)
    {
        double mid = (l+r)/2;
        if(cal(mid)>=aim) l = mid;
        else r = mid;  
    }
    return l;
}

int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        int hight=0;
        double sum = 0;
        scanf("%d",&n);
        for(int i=1;i<=n;i++) scanf("%d",&h[i]),hight=max(hight,h[i]);
        for(int i=1;i<=n;i++) scanf("%d",&a[i]),sum+=a[i]*a[i]*h[i];
        double ans = binsearch(0,hight,sum/2.0);
        printf("%d\n",(int)ans);
    }
    return 0;
}