poj 2546 Circular Area （两圆相交面积）

链接：poj 2546

题意：已知两圆的圆心和半径，求两圆相交部分的面积

分析：两圆的位置关系有三种：相离，相交，内含

相离时：相交面积为0

相交时，大扇形面积+小扇形面积-四边形面积

内含时，相交面积为小圆面积

#include<stdio.h>
#include<math.h>
#define PI acos(-1.0)
typedef struct stu
{
    double x,y;
}point;
double Distance(point a,point b)
{
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
double InterArea(point a,double R,point b,double r)
{
    if(R<r){
        double temp=R;
        R=r;
        r=temp;
    }
    double dis=Distance(a,b);
    if(dis>=R+r)    //两圆相离，相交面积为0
        return 0;
    if(dis<=R-r)    //两圆内含，相交面积为小圆的面积
        return PI*r*r;
    //两圆相交时
    double angle1=acos((R*R+dis*dis-r*r)/(2.0*R*dis));
    double angle2=acos((r*r+dis*dis-R*R)/(2.0*r*dis));
    double s=R*angle1*R+r*angle2*r;
    s-=R*dis*sin(angle1);
    return s;
}
int main()
{
    point a,b;
    double r1,r2;
    while(scanf("%lf%lf%lf%lf%lf%lf",&a.x,&a.y,&r1,&b.x,&b.y,&r2)!=EOF)
        printf("%.3lf\n",InterArea(a,r1,b,r2));
    return 0;
}