当前位置: 代码迷 >> 综合 >> 二分 hdu5432 Pyramid Split
  详细解决方案

二分 hdu5432 Pyramid Split

热度:75   发布时间:2023-12-14 03:49:42.0

传送门:点击打开链接

题意:给你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;
}


  相关解决方案