#include <iostream.h>
#include <stdio.h>
#include <math.h>

int main(void){
	double x, y, c, m, n1, n2, h1, h2, c1;
	while(cin >> x >> y >> c){
		c1 = 0;
		if(fabs(x-y)<0.00001){
			printf("%.3lf\n",(2*sqrt((x/2)*(x/2) - c*c)));
		}else if(c == 0){
			if(x <= y)
				printf("%.3lf\n",x);
			else 
				printf("%.3lf\n",y);
		}else{
			n1 = 0;
			n2 = 1;
			m = (n1 + n2)/2;
			while(1){
				h1 = sqrt(x*x - (c/m)*(c/m));
				h2 = sqrt(y*y - (c/(1-m))*(c/(1-m)));
				if(fabs(h1 - h2) < 0.00001) break;
				else{
					if(h1>h2){
						n2 = m;
						m = (m+n1)/2;
					}
					if(h1 <h2){
						n1 = m;
						m = (m+n2)/2;
					}
				}
			}
			printf("%.3lf\n",(h1 +h2)/2+0.0005);
		}
	}
	return 0;
}

• 有人救命吗？？ (7) smilezhen 2005-09-03 16:02:39
  • 列方程，二分法，求方程的解 (11) sunmoonstar_love 2005-09-03 16:21:07
  • AC了，但我原来怎么错了 (7) smilesmile 2005-09-09 18:13:44