UVa 375 Inscribed Circles and Isosceles Triangles

【转载出处】

http://blog.csdn.net/frankiller/article/details/7729784

【解题思路】

求等腰三角形的内切圆周长，内切圆不止一个，切完一个继续内切直到r<=0.000001；

设初始的内切圆半径为R；

勾股定理 斜边L=sqrt（R*R+B*B\4）；

接下来求R有2种方法：

1继续用勾股定理：（L-B/2）^2+R^2=(H-R）^2    R=sqrt(H*H+B*B/4);   R=R-B/2;     R=(R*R-H*H)/(-2*H)；

2.面积法：R*L+R*B/2=B*H/2；                                   R =B*H/(2*L+B)；

【代码】

#include<stdio.h>
#include<math.h>
#define pi asin(1.0)*2

void main() {
	double B,H,r,h,R;
	int t;

	scanf("%d",&t);
	while (t--) {
		scanf("%lf%lf",&B,&H);
		R = sqrt(H*H+B*B/4);
		R = B*H/(2*R+B);
		h=H;  r=R;
		while (r >= 0.000001) {
			h = h-2*r;
			r = R*h/H;
		}
		printf("%13lf\n",(H-h)*pi);

	if (t) 
		printf("\n");
	}
}