「HDU6158」 The Designer（圆的反演）

题目链接多校8-1009 HDU - 6158 The Designer

题意

T（<=1200）组，如图在半径R1、R2相内切的圆的差集位置依次绘制1,2,3,到n号圆，求面积之和（n<=1e7)。

题解

圆的反演：
（圆的反演就是半径为R，圆心O的圆为反演中心，点P的反演点就是在射线OP上满足\(|OP’|*|OP|=R^2\)的点P‘）
设切点为O，以O为圆心半径R的圆为反演点。将圆R1和R2反演得到两条直线，和两条直线相切的圆反演回去的圆就是1～n号圆的圆心。
那么它们的直径就是这些小圆的圆心和O的连线与小圆的交点反演回去的点的距离差。

再扔一次画图工具Desmos

比赛的时候想到这里就以为复杂度太高，不知道怎么预处理。其实到后面圆面积会收敛得很快。精度只要1e-5,就可以及时break掉。

代码

#include <bits/stdc++.h>
using namespace std;
const double pi = acos(-1);
const double R = 1;
int t,r1,r2,n;
double r0,d,a,b,r,s;
double ans;
int main() {
    scanf("%d",&t);
    while(t--){
        scanf("%d%d%d",&r1,&r2,&n);
        if(r2<r1)swap(r1,r2);
        d=R*(r1+r2)/r2/r1/4;
        r0=d-R/2/r1;
        r=r2-r1;
        ans=pi*r*r;
        for(int i=1;i<=n/2;++i){
            a=sqrt(d*d+i*r0*i*r0*4)-r0,b=a+r0*2;
            r=(R/a-R/b)/2;
            s=pi*r*r;
            ans+=s;
            if(i*2<n)ans+=s;
            if(s*(n-i*2)<1e-6){
                break;
            }
        }
        printf("%.5f\n",ans);
    }
    return 0;
}