题意:有两个内外径一样的环,给出两个环的圆心坐标,求重叠部分的面积。
思路:两圆相交求相交面积+容斥,
两圆相交求相交面积:这个可以用模版
容斥:就是大圆相交面积-2个大小圆相交面积+两个小圆相交面积
#include <iostream>
#include <cstring>
#include <cmath>
#include <cstdio>
#include <iomanip>
#define pi 4.0*atan(1.0)//圆周率
using namespace std;
double jiao(int r1,int r2,double d)//两圆相交模版
{
if(d>=(r1+r2)*(r1+r2))return 0;//相离和外切;
if(d <=(r1-r2)*(r1-r2))return r1<r2 ? pi*r1*r1 : pi*r2*r2;//内切;
d=sqrt(d);
double A1=acos((r1*r1+d*d-r2*r2)/(2.0*r1*d));//求半径是r1的圆所对扇形的圆心角;
double A2=acos((r2*r2+d*d-r1*r1)/(2.0*r2*d));//求半径是r2的圆所对扇形的圆心角;
double s;
s=r1*d*sin(A1);//多余三角形的面积;
double S;
S=r1*r1*A1+r2*r2*A2;//两扇形的面积和
S=S-s;
return S;//相交面积;
}
int main()
{
int t;
int f;
while(scanf("%d",&t)!=EOF)
{
f=0;
while(t--)
{
double r1,r2;
double x1,x2,y1,y2;
scanf("%lf%lf%lf%lf%lf%lf",&r1,&r2,&x1,&y1,&x2,&y2);
double ans,ans1,ans2,ans3,ans4;
double d;
d=(x1-x2)*(x1-x2)+1.0*(y1-y2)*(y1-y2);//两圆心的距离;
ans4=jiao(r2,r2,d);//大圆相交的面积;
ans1=jiao(r1,r1,d);//小圆相交的面积;
ans2=jiao(r1,r2,d);//一大一小两圆相交的面积;
ans=ans4-ans2*2+ans1;//一大一小两圆相交面积要减两次;
printf("Case #%d: ",++f);
printf("%.6f\n",ans);
}
}
return 0;
}