Description
给出大小圆半径,按此顺序一个个放 n 个小圆,问所有小圆面积之和
Input
第一行一整数
Output
输出小圆面积之和,保留五位小数
Sample Input
2
5 4
1
4 5
1
Sample Output
3.14159
3.14159
Solution
首先给出圆的曲率 k=±1r ,其中 r 为圆的半径,如果该圆外切于其他圆则取正,内切于其他圆则取负
然后给出
对于该题,由于加入的小圆除第一个外两两对称,只求一边即可,取 k1=−1R1 为外侧大圆的曲率, k2=1R2 为内侧大圆的曲率,加入的第一个圆半径为 R1−R2 ,故其曲率为 k3=1R1−R2 ,根据 Descartes 定理求出第二三个圆的曲率为 k4=k1+k2+k3+k1k2+k1k3+k2k3−−−−−−−−−−−−−−−√=k1+k2+k3 ,求出 k4 后,为了求第四五个圆的曲率,相当于把 k3 换成 k4 继续该过程,注意到
(k1+k2+k3+k4)2=2(k21+k22+k23+k24)
(k1+k2+k4+k5)2=2(k21+k22+k24+k25)
两式做差得 k3+k5=2(k1+k2+k4) ,在求其他小圆过程中 k1,k2 不变, k3,k4 为上两步求出的圆的曲率,故每次可以线性求出新加入的圆的曲率
Code
#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
#define eps 1e-13
#define PI acos(-1.0)
int main()
{
int T,r1,r2,n;
scanf("%d",&T);
while(T--)
{
scanf("%d%d%d",&r1,&r2,&n);
if(r1<r2)swap(r1,r2);
double k1=-1.0/r1,k2=1.0/r2,k3=1.0/(r1-r2),k4=k1+k2+k3;
double ans=(r1-r2)*(r1-r2);
for(int i=2;i<=n;i+=2)
{
double r=(1.0/k4)*(1.0/k4);
if(r<eps)break;
ans+=r;
if(i+1<=n)ans+=r;
double k5=2.0*(k1+k2+k4)-k3;
k3=k4,k4=k5;
}
printf("%.5f\n",PI*ans);
}
return 0;
}