hdu 5120 Intersection

题意：有两个内外径一样的环，给出两个环的圆心坐标，求重叠部分的面积。

思路：两圆相交求相交面积+容斥，
两圆相交求相交面积：这个可以用模版

容斥：就是大圆相交面积-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;
}