HDU 6158 圆的反演

上一篇介绍了院的反演，这一篇就是一道网络赛圆的反演的题目。其实感觉这道题目有点…. 感觉它卡的东西不是太好 反正计算几何你一定要注意精度问题。emmmm…… 这个题目就是圆的反演，你只需要找对里面的几何关系就好了

注意由于反演 中具有一一对应的性质，只有原来相交的点到后面才会可能相交，剩下的就是圆的反演了 （这个东西没有模版）.

还有我们用反演变换的时候一般关注两个东西，一个是过反演中心的圆 （因为这种圆的反形是一条直线） 另一种情况我们关注不过反演中心的直线，因为这个的反形是圆。也一般在这种情况下我们才使用反演变换 多个圆的时候有时也会去考虑这种变换。

#include <iostream>#include <algorithm>#include <cstring>#include <cstdio>#include <cmath>using namespace std;const double eps = 1e-7;const double PI = 3.141592653589793;int main () {    int T;    scanf ("%d",&T);    while (T --) {        double r1,r2;        int n;        scanf ("%lf%lf%d",&r1,&r2,&n);        if (r1 < r2) swap (r1,r2);        double r3 = r1 - r2;        double ans = r3 * r3;        double R = 4 * r1;        double l = R / (2 * r1),r = R / (2 * r2);        double rr = (r - l) / 2;        int u = (n - 1)/ 2;        int v = 0;        if (n % 2 == 0) v = 1;        double add = 2 * rr;        double len = rr + l;        for (int i = 1;i <= u + v; ++ i) {            double dis = sqrt (add * add + len * len);            add += 2 * rr;            double l1 = dis - rr;            double r1 = dis + rr;            double ansr = (R / l1 - R / r1) / 2;            ans += ansr * ansr * 2;            if (i == u + v && v) {                ans -= ansr * ansr;            }            if (ansr <= eps) break;        }        printf ("%.5f\n",ans*PI);    }    return 0;}