hdu 1798 Tell me the area（相交圓的面積）

相交圓的面積模板

純數學知識，S1，S2是兩個扇形的面積，S3是矩形的面積，

a1，a2是兩個角，用來求扇形面積

q是周長的一半，用來求矩形面積

#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
using namespace std;

struct circle
{
    double x,y,r;
};
int main()
{
    circle a,b;
    double S1,S2,S3,a1,a2,d,p,ans,mi;
    while(cin>>a.x>>a.y>>a.r>>b.x>>b.y>>b.r)
    {
        d=(double)sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
        mi=min(a.r,b.r);
        if(d>a.r+b.r||d==a.r+b.r)
            ans=0.0;
        else if(d<abs(a.r-b.r)||d==abs(a.r-b.r))
            ans=acos(-1.0)*mi*mi;
        else
        {
            p=(a.r+b.r+d)/2.0;
            a1=acos((a.r*a.r+d*d-b.r*b.r)/(2.0*a.r*d));
            a2=acos((b.r*b.r+d*d-a.r*a.r)/(2.0*b.r*d));
            S1=a1*a.r*a.r;
            S2=a2*b.r*b.r;
            S3=2*sqrt(p*(p-a.r)*(p-b.r)*(p-d));
            ans=S1+S2-S3;
        }
        printf("%.3lf\n",ans);

    }
    return 0;
}

WA：輸入double類型

相離相交內含的判斷

pi 用acos（-1）表示