1 #include<iostream>
2 #include<numeric>
3 #include<iomanip>
4 #include<algorithm>
5 #include<cmath>
6 using namespace std;
7
8 int main()
9 {
10 //freopen("acm.acm","r",stdin);
11 double x1;
12 double y1;
13 double r1;
14 double x2;
15 double y2;
16 double r2;
17 double A;
18 double B;
19 double d;
20 double s;
21 double area;
22 while(cin>>x1>>y1>>r1>>x2>>y2>>r2)
23 {
24 d = sqrt((x1 - x2)*(x1 - x2)+(y1 - y2)*(y1 - y2));
25 if(r1 < r2)
26 iter_swap(&r1,&r2);
27 if(d >= r1+r2)
28 {
29 cout<<setiosflags(ios::fixed)<<setprecision(3)<<0.0<<endl;
30 return 0;
31 }
32 if(d <= r1-r2)
33 {
34 cout<<setiosflags(ios::fixed)<<setprecision(3)<<r2*r2*2.0*acos(0.0000000)<<endl;
35 return 0;
36 }
37 A=acos((r2*r2+d*d-r1*r1)/2/r2/d);//相交
38 B=acos((r1*r1+d*d-r2*r2)/2/r1/d);
39 s=(r1+r2+d)/2;
40 s=sqrt(s*(s-r1)*(s-r2)*(s-d));
41 area=r1*r1*B+r2*r2*A-2*s;
42 cout<<setiosflags(ios::fixed)<<setprecision(3)<<area<<endl;
43 }
44 }
45
46
47
48
49
50
51 /*
52 分为两部分扇形来求：
53 设两圆圆心分别为O1，O2，交点为E，F，
54 设两圆半径分别为r1,r2,圆心距为d，则O1E=r1,O2E=r2,O1O2=d,
55 由此可求得三角形EO1O2面积S，以及∠EO1O2的度数α，∠EO2O1度数β，
56 公式为：三角形ABC中，a^2=b^2+c^2-2bc*cosA,具体的计算不用我说了吧.
57 最后阴影部分面积为α*r1^2+β*r2^2-2S.
58 */
