无关青云路,无关诗书,无你处,无江湖
Time Limit: 3000/1000MS (Java/Others) Memory Limit: 65535/65535KB (Java/Others)
每个圆的属性可以使用三个参数来描述, x x-圆心横坐标, y y-圆心纵坐标, r r-半径,现给出两个圆分别的参数,求两个圆相交的面积。
Input
第一行是一个正整数 T T( T≤12 T≤12),表示有 T T组测试数据 接下来有T行,每行有六个正整数描述两个圆的属性 x1,y1,r1,x2,y2,r2 x1,y1,r1,x2,y2,r2。 0≤x1,y1,r1,x2,y2,r2≤500 0≤x1,y1,r1,x2,y2,r2≤500
Output
对于每组测试数据输出一行,表示两个圆相交的面积大小,结果保留6位小数
Sample input and output
Sample Input | Sample Output |
---|---|
3 0 0 10 15 0 10 -10 -10 5 0 -10 10 100 100 20 100 110 20 | 45.331175 35.076661 860.843690 |
思路:两圆情况分析一下,分为相切,相交,内含,直接套模板就行
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#define PI 3.1415926535898
struct 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;
}