求两圆相交的面积

转自博客：http://blog.sina.com.cn/s/blog_850498e20100w6fq.html

假设半径小的圆为c1，半径大的圆为c2。
c1的半径r1，圆心坐标（x1,y1)。c2的半径r2,圆心坐标（x2,y2)。
d为两圆圆心连线的长度。
相交面积为S
d=sqrt((x1-x2)^2+(y1-y2)^2)
（1）如果r1+r2<=d
那么两圆相离，相交面积S=0
（2）如果r2-r1>=d
那么半径小的圆内含半径大的圆，那么相交面积为小圆的面积S=pi*r1*r1
（3）既非（1）也非（2）
在图上画两个相交圆，结合图像看。
那么两圆相交，连接小圆的圆心与两个圆的交点，连接大圆的圆心和两个圆的交点。
可以发现形成的图形被两个圆心的连线平分成2个全等三角形。
由小圆圆心和交点所连两条线(长度为半径）以及在大圆之内的弧所形成的扇形为S1
由大圆圆心和交点所连两条线(长度为半径）以及在小圆之内的弧所形成的扇形为S2
由小圆圆心和交点所连两条线以及由大圆圆心和交点所连两条线所形成的四边形的面积为S3
可见相交面积S=S1+S2-S3
要求出扇形的面积，要知道扇形的圆心角。
小圆包含的扇形的圆心角为2*a1(考虑一个三角形）
a1=acos((r1^2+d^2-r2^2)/(2.0*r1*d)) 余弦定理
a2=acos((r2^2+d^2-r1^2)/(2.0*r2*d))
S1=pi*r1*r1*2*a1/(2*pi)=a1*r1*r1
同理
S2=a2*r2*r2
S3为一个三角形面积的2倍
S3=2*r1*d*sin(a1)/2=r1*d*sin(a1)
则S=a1*r1*r1+a2*r2*r2-r1*d*sin(a1)

代码：

//求两圆相交的面积
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<algorithm>
#define esp 1e-8

using namespace std;

struct Circle{
    double x,y;
    double r;
};

double calArea(Circle c1, Circle c2)
{
    double d;
    double s,s1,s2,s3,angle1,angle2,temp;

    d=sqrt((c1.x-c2.x)*(c1.x-c2.x)+(c1.y-c2.y)*(c1.y-c2.y));
    if(d>=(c1.r+c2.r))//两圆相离
        return 0;
    if((c1.r-c2.r)>=d)//两圆内含,c1大
        return acos(-1.0)*c2.r*c2.r;
    if((c2.r-c1.r)>=d)//两圆内含,c2大
        return acos(-1.0)*c1.r*c1.r;

    angle1=acos((c1.r*c1.r+d*d-c2.r*c2.r)/(2*c1.r*d));
    angle2=acos((c2.r*c2.r+d*d-c1.r*c1.r)/(2*c2.r*d));

    s1=angle1*c1.r*c1.r;s2=angle2*c2.r*c2.r;
    s3=c1.r*d*sin(angle1);
    s=s1+s2-s3;

    return s;
}