关于圆的面积 （C语言代码）

参考代码:

#include<stdio.h>
#include<math.h>

#define PI 3.1415926

int main()
{
    int t;
    double x1, y1, x2, y2, r1, r2, d, ans;
    
    scanf("%d", &t);
    while(t--){
        
        scanf("%lf%lf%lf%lf%lf%lf", &x1, &y1, &x2, &y2, &r1, &r2);
        
        d = sqrt(pow((x1 - x2), 2) + pow((y1 - y2), 2)); // 圆心距d
        
        if(d >= (r1 + r2)) ans = 0; // 相离或外切
        
        else if(d <= fabs(r1 - r2)){ // 内切或内含
            
            ans = r1 < r2 ? PI * r1 * r1 : PI * r2 * r2;
        }
        else{ // 相交
            
            double A1 = 2 * acos((r1 * r1 + d * d - r2 * r2) / (2 * r1 * d));
            double A2 = 2 * acos((r2 * r2 + d * d - r1 * r1) / (2 * r2 * d));
            double S1 = A1 / 2 * r1 * r1 + A2 / 2 * r2 * r2;
            double S2 = r1 * r1 * sin(A1) / 2 + r2 * r2 * sin(A2) / 2;
            ans = S1 - S2;
        }
        printf("%.3lf\n", ans);
    }
    return 0;
}