UESTC1713 无关青云路，无关诗书，无你处，无江湖(计算几何)

无关青云路，无关诗书，无你处，无江湖

Time Limit: 3000/1000MS (Java/Others)     Memory Limit: 65535/65535KB (Java/Others)

Submit Status

每个圆的属性可以使用三个参数来描述，xx-圆心横坐标，yy-圆心纵坐标，rr-半径，现给出两个圆分别的参数，求两个圆相交的面积。

Input

第一行是一个正整数TT(T≤12T≤12)，表示有TT组测试数据 接下来有T行，每行有六个正整数描述两个圆的属性x1,y1,r1,x2,y2,r2x1,y1,r1,x2,y2,r2。 0≤x1,y1,r1,x2,y2,r2≤5000≤x1,y1,r1,x2,y2,r2≤500

Output

对于每组测试数据输出一行，表示两个圆相交的面积大小，结果保留6位小数

Sample input and output

Sample InputSample Output

30 0 10 15 0 10-10 -10 5 0 -10 10100 100 20 100 110 20

45.33117535.076661860.843690

题意：计算两圆相交面积

思路：两圆情况分析一下，分为相切，相交，内含，直接套模板就行

#include <stdio.h>#include <stdlib.h>#include <math.h>#define PI 3.1415926535898struct circle{    double x,y;    double r;};double dis(circle a,circle b){    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));}double solve(circle a, circle b) {    double d = dis(a, b);    if (d >= a.r + b.r)        return 0;    if (d <= fabs(a.r - b.r)) {        double r = a.r < b.r ? a.r : b.r;        return PI * r * r;    }    double ang1 = acos((a.r * a.r + d * d - b.r * b.r) / 2. / a.r / d);    double ang2 = acos((b.r * b.r + d * d - a.r * a.r) / 2. / b.r / d);    double ret = ang1 * a.r * a.r + ang2 * b.r * b.r - d * a.r * sin(ang1);    return ret;}int main(){    int t;    while(scanf("%d",&t)!=EOF)    {        for(int i=0;i<t;i++)        {            circle a,b;        scanf("%lf%lf%lf%lf%lf%lf",&a.x,&a.y,&a.r,&b.x,&b.y,&b.r);        printf("%lf\n",solve(a,b));        }    }    return 0;}