LightOJ 1331 Agent J - Basic Geometry

题意：已知三个圆（两两相切）的半径r1,r2,r3，求中间部分面积

思路：三个圆心O1,O2,O3构成三角形，则ans=三角形面积-3个扇形面积

AC Code：

//AcerGY
//#pragma comment (linker, "/STACK:1024000000,1024000000")
#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <queue>
#include <vector>
#include <cmath>
#define ll long long
#define BUG puts("BUG!")
#define ci(val) scanf("%d",&val)
#define pi(val) printf("%d\n",val)
#define clr(a,val) memset(a,val,sizeof(a))
#define lowbit(x) ((x)&(-x))
#define inf 0x3f3f3f3f
#define inf64 1LL<<60
#define PI acos(-1.0)
#define eps 1e-8

//Game Theory
#define P "First"
#define N "Second"

#define maxn 10000
#define mod 10009

using namespace std;

void readfile(char s[])
{
#ifndef ONLINE_JUDGE
	freopen(s,"r",stdin);
#endif
}

double area_sector(double r,double a,double b,double c)
{
	double angle = acos((a*a+b*b-c*c)/(2*a*b));
	return r*r*angle/2;
}

int main()
{
	//readfile("in.txt");
	int T , cas = 1;
	cin>>T;
	while(T--)
	{
		double r1,r2,r3;
		scanf("%lf%lf%lf",&r1,&r2,&r3);
		double a = r1+r2 , b = r2+r3 , c = r1+r3;
		double tri = sqrt((r1+r2+r3)*r1*r2*r3);
		double sec = area_sector(r1,a,c,b)+
			area_sector(r2,a,b,c)+area_sector(r3,c,b,a);
		double area = tri-sec;
		printf("Case %d: %lf\n",cas++,area);
	}
	return 0;
}