hdu1798(圆的位置关系)

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=1798

 

题意：给出两个圆的圆心坐标与半径，求他们相交部分的大小

 

思路：有三种情况；

1. 两圆相离，area = 0；

2. 两圆内涵，area = PI×min(r1, r2)×min(r1, r2)；

3. 一般情况相交，用扇形面积减去三角形面积；

 

代码：

 1 #include <iostream>
 2 #include <stdio.h>
 3 #include <algorithm>
 4 #include <math.h>
 5 #define PI 2*asin(1)   //**注意此处不能简单的用3.1415926，精度不够
 6 #define length(x1, y1, x2, y2) sqrt((x2-x1)*(x2-x1)+(y2-y1)*(y2-y1))
 7 #define horn(a, b, c) ((a*a+b*b-c*c)/(2*a*b))
 8 using namespace std;
 9 
10 int main(void){
11     double x1, y1, r1, x2, y2, r2;
12     while(scanf("%lf%lf%lf", &x1, &y1, &r1)!=EOF){
13         scanf("%lf%lf%lf", &x2, &y2, &r2);
14         double area;
15         double gg = length(x1, y1, x2, y2);
16         if(gg>=r1+r2){
17             area = 0;
18         }else if(fabs(r2-r1)>=gg){
19             area = PI*min(r2, r1)*min(r2, r1);
20         }else{
21             double horna = acos(horn(r1, gg, r2));
22             double hornb = acos(horn(r2, gg, r1));
23             double area1 = r1*r1*horna;
24             double area2 = r2*r2*hornb;
25             double area3 = 0.5*r1*r1*sin(horna*2);
26             double area4 = 0.5*r2*r2*sin(hornb*2);
27             area = (area1 - area3) + (area2 - area4);
28         }
29         printf("%.3lf\n", area);
30     }
31     return 0;
32 }

 

转载于:https://www.cnblogs.com/geloutingyu/p/6066940.html