二分 hdu5432 Pyramid Split

传送门：点击打开链接

题意：给你n个四棱锥，一起放在一个平面上，然后在高度为h的位置，用一个平行平面去切所有的四棱锥，得到上下两部分，使上下两部分体积相等，求h的大小(取整数部分输出)

思路：很裸的二分，推出高度为h时四棱锥上部分体积和下部分体积，然后利用二分就可以了，注意可以直接循环一定的次数这样就可以避免精度问题了。还有一个要注意的地方，当h>A[i]时，上部分的体积应该是0，而不是负数，话说如果不考虑这个情况竟然可以过样例，但是小数据会wa，如果小数据里没有这些数据，估计还会挂一片→_→

#include<map>
#include<set>
#include<cmath>
#include<stack>
#include<queue>
#include<cstdio>
#include<cctype>
#include<string>
#include<vector>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<functional>
#define fuck printf("fuck")
#define FIN freopen("input.txt","r",stdin)
#define FOUT freopen("output.txt","w+",stdout)
using namespace std;
typedef long long LL;
typedef pair<int, int> PII;

const int MX = 1e5 + 5;

double A[MX], B[MX];

int main() {
    int T, n; //FIN;
    scanf("%d", &T);
    while(T--) {
        scanf("%d", &n);
        for(int i = 1; i <= n; i++) {
            scanf("%lf", &A[i]);
        }
        for(int i = 1; i <= n; i++) {
            scanf("%lf", &B[i]);
        }

        double L = 0, R = 1000, m;
        for(int i = 1; i <= 40; i++) {
            m = (L + R) / 2;
            double s, sum1 = 0, sum2 = 0;
            for(int j = 1; j <= n; j++) {
                if(A[j] <= m) s = 0;
                else s = B[j] * B[j] * (A[j] - m) * (A[j] - m) * (A[j] - m) / (3.0 * A[j] * A[j]);
                sum1 += s;
                sum2 += B[j] * B[j] * A[j] / 3.0 - s;
            }

            if(sum1 < sum2) R = m;
            else L = m;
        }
        printf("%d\n", (int)m);
    }
    return 0;
}