POJ2546 Circular Area(计算几何)

                                                                            Circular Area

                                                       Time Limit: 1000MS  Memory Limit: 65536K

 

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

Source

Northeastern Europe 2000, Far-Eastern Subregion

 

 

这道题直接上模板就可以了~

 

模板：

 

struct Circle{
    double x, y, 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 ;
    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.00 / a.r / d);
    double ang2 = acos((b.r * b.r + d * d - a.r * a.r) / 2.00 / b.r / d);
    double ret = ang1 * a.r * a.r + ang2 * b.r * b.r - d * a.r * sin(ang1);
    return ret;
}
//返回值即为两圆公共部分的面积

 

 

POJ2546：

 #include <bits/stdc++.h>

 using namespace std;

 #define pi    3.1415926535897932384626

 struct Circle{
     double x, y, r;
 } 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 ;
     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.00 / a.r / d);
     double ang2 = acos((b.r * b.r + d * d - a.r * a.r) / 2.00 / 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[].x, &r[].y, &r[].r, &r[].x, &r[].y, &r[].r))
         printf("%.3f\n", solve(r[], r[]));

     return ;

 }