Intersection(HDU5120 + 圆交面积）

题目链接：

　　http://acm.hdu.edu.cn/showproblem.php?pid=5120

题目：

题意：

　　求两个圆环相交的面积。

思路：

　　两个大圆面积交-2×大圆与小圆面积交+两小圆面积交。

代码实现如下：

 #include <set>
 #include <map>
 #include <deque>
 #include <ctime>
 #include <stack>
 #include <cmath>
 #include <queue>
 #include <string>
 #include <cstdio>
 #include <vector>
 #include <iomanip>
 #include <cstring>
 #include <iostream>
 #include <algorithm>
 using namespace std;

 typedef long long LL;
 typedef pair<LL, LL> pll;
 typedef pair<LL, int> pli;
 typedef pair<int, int> pii;
 typedef unsigned long long uLL;

 #define lson rt<<1
 #define rson rt<<1|1
 #define name2str(name)(#name)
 #define bug printf("**********\n");
 #define IO ios::sync_with_stdio(false);
 #define debug(x) cout<<#x<<"=["<<x<<"]"<<endl;
 #define FIN freopen("/home/dillonh/CLionProjects/in.txt","r",stdin);

 const double eps = 1e-;
 const int maxn = 1e5 + ;
 const int inf = 0x3f3f3f3f;
 const double pi = acos(-1.0);
 const LL INF = 0x3f3f3f3f3f3f3f3fLL;

 int t, r, R;
 int x1, x2, yy, y2;

 double cirinter(int x1,int y1,int r1,int x2,int y2,int r2)//圆交面积公式
 {
     double d,s,t,a1,a2,s1,s2,s3;
     if(r1<r2) swap(r1, r2);
     d=sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2)); //两圆心距离
     if(d>=r1+r2)return ; //两圆相离
     else if(d<=(r1-r2)) //两圆内含
         s=pi*r2*r2;
     else { //两圆相交
         a1=acos((r1*r1+d*d-r2*r2)/(*r1*d));//大圆中扇形圆心角的一半
         a2=acos((r2*r2+d*d-r1*r1)/(*r2*d));//小圆中扇形圆心角的一半
         s1=a1*r1*r1;//大圆中的那个扇形面积
         s2=a2*r2*r2;//小圆中的那个扇形面积
         s3=r1*sin(a1)*d;//两圆心与两交点组成的四边形面积
         s=s1+s2-s3;//圆交面积
     }
     return s;
 }

 int main() {
 #ifndef ONLINE_JUDGE
     FIN;
 #endif
     scanf("%d", &t);
     int icase = ;
     while(t--) {
         scanf("%d%d", &r, &R);
         scanf("%d%d%d%d", &x1, &yy, &x2, &y2);
         double ans = cirinter(x1, yy, R, x2, y2, R) -  * cirinter(x1, yy, R, x2, y2, r) + cirinter(x1, yy, r, x2, y2, r);
         printf("Case #%d: %.6f\n", ++icase, ans);
     }
     return ;
 }