POJ 2546 Circular Area（两个圆相交面积）

POJ 2546 Circular Area（两个圆相交面积）

ACM

题目地址：POJ 2546 Circular Area

题意： 
两个圆相交面积。

分析： 
注意内含与相离的情况。 
G++下只有把PI设为3.141592653才能过TAT...以后几何题是不是都要用C++交了...

代码：

#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;
}