题目链接: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;
}