相交圆的面积模板
纯数学知识,S1,S2是两个扇形的面积,S3是矩形的面积,
a1,a2是两个角,用来求扇形面积
q是周长的一半,用来求矩形面积
#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
using namespace std;
struct circle
{
double x,y,r;
};
int main()
{
circle a,b;
double S1,S2,S3,a1,a2,d,p,ans,mi;
while(cin>>a.x>>a.y>>a.r>>b.x>>b.y>>b.r)
{
d=(double)sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
mi=min(a.r,b.r);
if(d>a.r+b.r||d==a.r+b.r)
ans=0.0;
else if(d<abs(a.r-b.r)||d==abs(a.r-b.r))
ans=acos(-1.0)*mi*mi;
else
{
p=(a.r+b.r+d)/2.0;
a1=acos((a.r*a.r+d*d-b.r*b.r)/(2.0*a.r*d));
a2=acos((b.r*b.r+d*d-a.r*a.r)/(2.0*b.r*d));
S1=a1*a.r*a.r;
S2=a2*b.r*b.r;
S3=2*sqrt(p*(p-a.r)*(p-b.r)*(p-d));
ans=S1+S2-S3;
}
printf("%.3lf\n",ans);
}
return 0;
}
WA:输入double类型
相离相交内含的判断
pi 用acos(-1)表示