题目大意:给出2圆的圆心坐标及半径,求相交部分的面积。
解题思路:分情况讨论相离,包含,及相交关系。相交时用四边形面积-2扇形面积之和
#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
#include <algorithm>
#include <queue>
#include <string>
#include <map>
#include <stack>
#include <list>
using namespace std;
const double eps = 1e-10;
const double PI = acos(-1.0);
int dcmp(double x)
{
if(fabs(x) < eps) return 0;
return x<0 ? -1 : 1;
}
double dis(double x1,double y1,double x2,double y2)
{
return sqrt((x2-x1)*(x2-x1)+(y2-y1)*(y2-y1));
}
int main()
{
double x1,y1,r1,x2,y2,r2;
while(scanf("%lf%lf%lf%lf%lf%lf",&x1,&y1,&r1,&x2,&y2,&r2) != EOF)
{
double r = min(r1,r2);
double d = dis(x1,y1,x2,y2);
if(dcmp(d-(r1+r2)) >= 0)
printf("0.000\n");
else if(dcmp(d-fabs(r1-r2)) <= 0)
printf("%.3lf\n",PI*r*r);
else
{
double a1 = acos((r1*r1+d*d-r2*r2)/(2*r1*d));//扇形角度的一半
double a2 = acos((r2*r2+d*d-r1*r1)/(2*r2*d));
double s1 = r1*r1*a1;//扇形面积
double s2 = r2*r2*a2;
double s = r1*d*sin(a1);//四边形面积
double area = s1+s2-s;
printf("%.3lf\n",area);
}
}
return 0;
}