POJ 2546

 #include<iostream>
 #include<numeric>
 #include<iomanip>
 #include<algorithm>
 #include<cmath>
 using namespace std;

 int main()
 {
     //freopen("acm.acm","r",stdin);
     double x1;
     double y1;
     double r1;
     double x2;
     double y2;
     double r2;
     double A;
     double B;
     double d;
     double s;
     double area;
     while(cin>>x1>>y1>>r1>>x2>>y2>>r2)
     {
     d = sqrt((x1 - x2)*(x1 - x2)+(y1 - y2)*(y1 - y2));
     if(r1 < r2)
         iter_swap(&r1,&r2);
     if(d >= r1+r2)
     {
         cout<<setiosflags(ios::fixed)<<setprecision()<<0.0<<endl;
         return ;
     }
     if(d <= r1-r2)
     {
         cout<<setiosflags(ios::fixed)<<setprecision()<<r2*r2*2.0*acos(0.0000000)<<endl;
         return ;
     }
     A=acos((r2*r2+d*d-r1*r1)//r2/d);//相交
     B=acos((r1*r1+d*d-r2*r2)//r1/d);
     s=(r1+r2+d)/;
     s=sqrt(s*(s-r1)*(s-r2)*(s-d));
     area=r1*r1*B+r2*r2*A-*s;
     cout<<setiosflags(ios::fixed)<<setprecision()<<area<<endl;
     }
 }

 /*
 分为两部分扇形来求：
 设两圆圆心分别为O1，O2，交点为E，F，
 设两圆半径分别为r1,r2,圆心距为d，则O1E=r1,O2E=r2,O1O2=d,
 由此可求得三角形EO1O2面积S，以及∠EO1O2的度数α，∠EO2O1度数β，
 公式为：三角形ABC中，a^2=b^2+c^2-2bc*cosA,具体的计算不用我说了吧.
 最后阴影部分面积为α*r1^2+β*r2^2-2S.
 */