題目的大意是在等腰三角形的高上堆圓,直到精度達到題目要求。並求出所有三角形的周長。
等腰三角形的內切圓半徑不難求,把裡面的三個三角形切開,即3個高相等的三角形,已知三角形總面積就可以求出三個三角形的高,也就是內切圓的半徑了。
然後要求堆砌內切圓,把三角形已經算過的部分切掉,根據相似就可以求出每個小三角形的內切圓,然後循環求到精度要求即可。總的周長就是總半徑*pi。
這裡pi=atan(1.0)*4.0
代碼:
#include <cstdio>
#include <cmath>
using namespace std;
const double pi = atan(1.0) * 4;
double b, h, l, th, r, k;
int main() {
int n;
scanf("%d", &n);
while (n--) {
scanf("%lf%lf", &b, &h);
l = sqrt(b * b / 4 + h * h);
r = b * h / (b + 2 * l);
th = 2 * r;
k = r / h;
l = h - r * 2;
while (l * k > 0.000001) {
l -= l * k * 2;
}
printf("%13.6lf\n", pi * (h - l));
if (n)
printf("\n");
}
return 0;