hdu-6158-圆的反演

题目链接：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;
}