Matt is a big fan of logo design. Recently he falls in love with logo made up by rings. The following figures are some famous examples you may know.
A ring is a 2-D figure bounded by two circles sharing the common center. The radius for these circles are denoted by r and R (r < R). For more details, refer to the gray part in the illustration below.
Matt just designed a new logo consisting of two rings with the same size in the 2-D plane. For his interests, Matt would like to know the area of the intersection of these two rings.
A ring is a 2-D figure bounded by two circles sharing the common center. The radius for these circles are denoted by r and R (r < R). For more details, refer to the gray part in the illustration below.
Matt just designed a new logo consisting of two rings with the same size in the 2-D plane. For his interests, Matt would like to know the area of the intersection of these two rings.
Each of the following two lines contains two integers x i, y i (0 ≤ x i, y i≤ 20) indicating the coordinates of the center of each ring.
2 2 3 0 0 0 0 2 3 0 0 5 0
Case #1: 15.707963 Case #2: 2.250778
#include <cstdio>
#include <cstring>
#include <math.h>
#include <algorithm>
#include <iostream>
using namespace std;
#define maxn 100005
const double eps = 1e-5;
const double pi = acos (-1.0);
struct Point
{
double x,y;
};
double dist(Point a,Point b)
{
return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
double Area(Point c1,double r1,Point c2,double r2)
{
double d=dist(c1,c2);
if(r1+r2<d+eps)return 0;
if(d<fabs(r1-r2)+eps)
{
double r=min(r1,r2);
return pi*r*r;
}
double x=(d*d+r1*r1-r2*r2)/(2*d);
double t1=acos(x/r1);
double t2=acos((d-x)/r2);
return r1*r1*t1+r2*r2*t2-d*r1*sin(t1);
}
double R,r;
Point O1,O2;
int main()
{
int t;
scanf("%d",&t);
for(int cas=1;cas<=t;cas++)
{
scanf("%lf%lf",&r,&R);
scanf("%lf%lf%lf%lf",&O1.x,&O1.y,&O2.x,&O2.y);
double S=Area(O1,R,O2,R)-Area(O1,R,O2,r)-Area(O2,R,O1,r)+Area(O1,r,O2,r);
printf("Case #%d: %.6lf\n",cas,S);
}
return 0;
}