POJ 2546 Circular Area [相交园面积]【计算几何】

题目链接：http://poj.org/problem?id=2546

Circular Area
Time Limit: 1000MS Memory Limit: 65536K
Total Submissions: 5662 Accepted: 2215
Description

Your task is to write a program, which, given two circles, calculates the area of their intersection with the accuracy of three digits after decimal point.
Input

In the single line of input file there are space-separated real numbers x1 y1 r1 x2 y2 r2. They represent center coordinates and radii of two circles.
Output

The output file must contain single real number - the area.
Sample Input

20.0 30.0 15.0 40.0 30.0 30.0
Sample Output

608.366

题目大意：
就是求两个圆的相交部分的面积

解题思路：
就是简单的几何题目。
两个圆的位置一共有五种情况
1.内含 2.内切 3.相交 4.外切 5.外离
在计算几何中 内含与内切 可以当做同一种情况 外切与外离 可以当做同一种情况
Ps：如果有疑问可以自己画一画图 就会明白了

所以在计算几何中只要分成3种情况就行了
先定义 d为两圆圆心距离

1.外离

这种情况下d>=r1+r2
两圆间没有相交的部分
所以相交面积是0

2.内含

这种情况下 大圆把小圆全部覆盖
相交的面积就是校园的面积 PI * r1 * r1
此情况下 d<=r大-r小
写成代码就是d <= fabs(r1 - r2)

3.相交

相交的部分可以分成 大小两个圆上弧ab与直线ab围成的图形
即可看作
S扇形abc1-S△abc1+S扇形abc2-S△abc2<=>
S扇形abc1+S扇形abc2-S□ac1bc2;
根据中学学的扇形的面积公式可以推到如下

ang1 = acos((r1 * r1 + d * d - r2 * r2) / 2. / r1 / d)
ang2 = acos((r2 * r2 + d * d - r1 * r1) / 2. / r2 / d)
S扇形abc1=ang1 * r1 * r1
S扇形abc2=ang2 * r2 * r2
S□ac1bc2=d * r1 * sin(ang1)

S = ang1 * r1 * r1+ang2 * r2 * r2 - d * r1 * sin(ang1)

附本题代码

#include <iostream> 
#include <cstdio> 
#include <cstring> 
#include <cmath> 
#include <algorithm> 
using namespace std; 
#define repf(i,a,b) for(int i=(a);i<=(b);i++) 

typedef long long ll; 
const double PI = 3.141592653; 

struct Round {  
    double x, y; 
    double r; 
} r[2]; 

double dis(Round a, Round b) {  
    return sqrt((a.x - b.x) * (a.x - b.x) + (a.y - b.y) * (a.y - b.y)); 
}  

double solve(Round a, Round 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() {  
    while (~scanf("%lf%lf%lf%lf%lf%lf", &r[0].x, &r[0].y, &r[0].r, &r[1].x, &r[1].y, &r[1].r)) {  
        printf("%.3f\n", solve(r[0], r[1])); 
    }  
    return 0; 
}