uva 10668

#include <iostream>
#include <cstdlib>
#include <cstdio>
#include <cmath>

using namespace std;

int main()
{
    double L,n,C,R;
    while (~scanf("%lf%lf%lf",&L,&n,&C)) {
        if (L <  && n <  && C < )
            break;

        double S = (+n*C)*L;
        double mid,l = 0.0,r = 0.5*L;
        while (r-l > 1e-) {
            mid = (l+r)/2.0;
            R = (*mid*mid+L*L)/8.0/mid;
            if (*R*asin(L/R/2.0) < S)
                l = mid;
            else r = mid;
        }

        printf("%.3lf\n",mid);
    }
    return ;
}

设升高h，则有：R^2 = （0.5L）^2 + （R-h）^2；

设弧长s，则有：s = 2*r*arcsin（L/（2R））；

然而我看见了显然  画一下图可以很容易推出公式，设圆弧扇形部弧度r，那么可以计算出铁棒长度为lr/sin(r)这个公式在[0,
pi/2]是单调递增的，所以可以用二分法去求解

#include <cstdio>
#include <cstring>
#include <cmath>

const double pi = acos(-1.0);

double l, n, c;

double cal(double r) {
    return l / sin(r) * r;
}

int main() {
    while (~scanf("%lf%lf%lf", &l, &n, &c)) {
    if (l < ) break;
    double x = , y = pi / , lx = ( + n * c) * l, m;
    for (int i = ; i < ; i++) {
        m = (x + y) / ;
        if (cal(m) < lx) x = m;
        else y = m;
    }
    printf("%.3lf\n", l /  / sin(m) * ( - cos(m)));
    }
    return ;
}