zoj 1597 Circular Area（求两圆相交的面积）

题目地址

题目大意：给出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;
}