ZOJ.1597 Circular Area【两圆相交面积】 2015/09/23

Circular Area

Time Limit: 2 Seconds
Memory Limit: 65536 KB

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.

Process to the end of input.

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

计算两圆相交面积，模板

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>

using namespace std;

double overlop_area(double x1,double y1,double r1,double x2,double y2,double r2){
double d = sqrt( (x1-x2)*(x1-x2) + (y1-y2)*(y1-y2) );
double angle1,angle2,h1,h2,s1,s2;
if( d >= r1+r2 )
return 0;
if( d <= fabs(r1-r2) )
return min( acos(-1.0)*r1*r1 , acos(-1.0)*r2*r2 );
angle1 = acos( (d*d+r1*r1-r2*r2) / (2*d*r1) );
angle2 = acos( (d*d+r2*r2-r1*r1) / (2*d*r2) );
h1 = angle1 * r1 * r1;
h2 = angle2 * r2 * r2;
s1 = r1 * r1 * cos(angle1) * sin(angle1);
s2 = r2 * r2 * cos(angle2) * sin(angle2);
return h1 + h2 - (s1+s2);
}

int main(){
double x1,x2,r1,y1,y2,r2;
while( ~scanf("%lf%lf%lf%lf%lf%lf",&x1,&y1,&r1,&x2,&y2,&r2) ){
printf("%.3lf\n",overlop_area(x1,y1,r1,x2,y2,r2));
}
return 0;
}