HDU 5954 Do not pour out (

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

参考博客：http://blog.csdn.net/danliwoo/article/details/53002695

#include <bits/stdc++.h>
using namespace std;
const double eps = 1e-10;   //1e-8就会挂....
const double PI = acos(-1.0);
bool equ(double x, double y) {
	if(fabs(x - y) < eps) return 1;
	return 0;
}
double getalp(double the) {
	return acos(2*tan(the)-1);
}
double f(double the) {
	double alp = getalp(the);
	double cosalp = cos(alp);
	double sinalp = sin(alp);
	return (PI*cosalp - alp*cosalp + sinalp - sinalp*sinalp*sinalp/3) / tan(the);
}
int main() {
	int t;
	scanf("%d", &t);
	double d, ans;
	while(t--) {
		scanf("%lf", &d);
		if(equ(d, 0)) {
			puts("0.00000");
			continue;
		}
		if(equ(d, 1)) ans = sqrt(2.0)*PI;
		else if(d > 1) ans = sqrt(4*d*d-16*d+20)*PI/2;
		else {
			double l = 0, r = PI / 4;
			double mid, res;
			while(!equ(l, r)) {
				mid = (l + r) / 2;
				res = f(mid);
				if(equ(res, PI*d)) break;
				if(res > PI*d) r = mid;
				else l = mid;
			}
			double alp = getalp(mid);
			double s1 = PI - alp + (2.0*tan(mid) - 1)*sin(alp);
			ans = s1 / sin(mid);
		}
		printf("%.5lf\n", ans);
	}
	return 0;
}