【ACMclub周赛5】Problem C——求两圆的交叉面积

题目：点击打开链接

参考了别人的推导过程，粘贴如下：

• 假设半径小的圆为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)

照着打出来，wa了两次，因为小圆减大圆的情况没有考虑，也就是说忽略了负数。。要r1<r2的

#include <iostream>#include <cstring>#include <algorithm>#include <iomanip>#include <cstdio>#include <cmath>using namespace std;const double PI=acos(-1);int main(){double rx1,ry1,rx2,ry2,r1,r2;while(cin>>rx1>>ry1>>r1>>rx2>>ry2>>r2){double drma=sqrt((rx1-rx2)*(rx1-rx2)+(ry1-ry2)*(ry1-ry2));if(r1>r2){double tmp=r1;r1=r2;r2=tmp;}if(r1+r2<=drma){cout<<"0.000"<<endl;}else if(r2-r1>=drma){double res=PI*r1*r1;cout<<setiosflags(ios::fixed)<<setprecision(3)<<res<<endl;}else {   double a1=acos((r1*r1+drma*drma-r2*r2)/(2.0*r1*drma));  double a2=acos((r2*r2+drma*drma-r1*r1)/(2.0*r2*drma)); double res=(a1*r1*r1+a2*r2*r2-r1*drma*sin(a1));   cout<<setiosflags(ios::fixed)<<setprecision(3)<<res<<endl;}}return 0;}