题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=6158
题意:
给你两个内切的圆的半径,一大一小,在两个圆之间放n个小圆(如图),求n个小圆的总面积。
思路:
设给定的两个圆,大圆的半径为r1,小圆的半径为r2,那么r1与r2之间圆1的半径为r2-r1,
由于每个小圆的半径求法类似,只要能求出圆2的半径,就可以求出圆3,4,5,...,n的半径。
现在问题转化为了,四个圆两两相切,已知三个圆的半径,求第四个圆的半径。
赛后学习了一下,这个可以用圆的反演来求,下面介绍一下圆的反演:
首先有一个圆O作为反演圆,取圆O的圆心作为反演中心,圆O的半径R作为反演半径,现在有一点P,P经过反演中心O反演之后变成P2,那么,,如图,O、P、P2三点共线且|OP| * |OP2| = R^ 2
不经过反演中心的圆反演之后还是圆:
经过反演中心的圆反演后是一条直线:
对于这个题,取大圆r1与小圆r2的切点作为反演中心,任取一个反演半径,分别对大圆r1,r2,圆1,圆2做反演变换:
设直线BL与圆2的交点为P1, Q1,直线 BL与圆2的反演圆交点为P2,Q2,
BM,ML的长度根据已知条件可以求出,所以可以求出线段BL,BP2, BQ2的长度,
由于P2,Q2是P1,Q1的反演点,所以可以求出BP1,BQ1的长度,这样圆2的半径r = (BQ1 - BP1) / 2,
其他圆半径的求法类似。。。
知道半径如何求就可以求面积了。。。这个题给的n特别大,直接枚举求面积会超时,容易知道,当n特别大时,
小圆的半径会特别小,它的面积可以忽略不计,所以设置一个eps来控制精度(这个地方特别坑,需要选取合适的精度,精度调的太低了会WA,调的太高了会TLE。。。),当面积小到一定程度时,不再计算,直接退出。
代码:
# include <iostream>
# include <algorithm>
# include <cstdio>
# include <cstring>
# include <cmath>
using namespace std;
typedef long long ll;
const double eps = 1e-13;
const double pi = acos(-1.0);
int r1, r2, n;
double sqr(double x) {
return x * x;
}
int main(void)
{
int T; scanf("%d", &T);
while (T-- && scanf("%d %d %d", &r1, &r2, &n)) {
if (r1 < r2) swap(r1, r2);
double r3 = r1 - r2;
double ans = sqr(r3);
--n;
double R = 2 * r1;
double x1 = R, x2 = sqr(R) / (2 * r2);
double r = (x2 - x1) / 2;
double L = R + r;
for (int i = 0, j = 1; i < n; i += 2, ++j) {
double t = sqrt(sqr(L) + sqr(2 * j * r));
double t1 = t - r, t2 = t + r;
double s1 = sqr(R) / t1, s2 = sqr(R) / t2;
double tmpr = (s1 - s2) / 2;
double tmp = sqr(tmpr);
if (tmp < eps) break;
if (i + 1 == n) ans += tmp;
else ans += 2 * tmp;
}
printf("%.5lf\n", ans * pi);
}
return 0;
}