2016ACM/ICPC亚洲区沈阳站 HDU5954 Do not pour out (积分+二分)

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=5954

题意：有一个圆柱形杯子，地面为半径为1的圆，高为2.现在给定一个实数d为杯子中水面的高度，问将杯子倾斜到最大角度且杯中水不溢出水面时的水面表面面积

思路：

代码：

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <vector>
#include <string>
#include <queue>
#include <map>
#include <algorithm>
using namespace std;

typedef long long ll;
const int maxn = 1e5 + 5;
const double eps = 1e-100;
const double pi = acos(-1);

int cmp(double x, double y)
{
	if (fabs(x - y) <= eps)return 0;
	if (x < y)return -1;
	return 1;
}

double cal(double x)
{
	return -x * cos(x) + sin(x) - sin(x)*sin(x)*sin(x) / 3.0;
}

int main()
{
	int t;
	scanf("%d", &t);
	while (t--)
	{
		double d;
		scanf("%lf", &d);
		if (cmp(d, 1) >= 0)
		{
			printf("%.5lf\n", pi*sqrt((2 - d)*(2 - d) + 1));
			continue;
		}
		if (cmp(d, 0) == 0)
		{
			printf("0.00000\n");
			continue;
		}
		double l = 0, r = pi / 4.0;
		double res;
		int cnt = 2e3;
		while (cmp(l, r) <= 0)
		{
			cnt--;
			double mid = (l + r) / 2.0;
			double a = acos(1.0-2.0*tan(mid));
			double s = (cal(a) - cal(0))/tan(mid);
			if (cmp(s, pi*d) == 0 || cnt<=0)
			{
				res = mid;
				break;
			}
			else if (cmp(s, pi*d) >= 0)
			{
				r = mid;
			}
			else
			{
				l = mid;
			}
		}
		double a = acos(1.0-2.0*tan(res));
		double ans = (a - sin(2*a)/2.0) /(sin(res));
		printf("%.5lf\n", ans);
	}
}