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-8;
const int maxn = 1e5 + 7;
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 0; //两圆相离
    else if(d<=(r1-r2)) //两圆内含
        s=pi*r2*r2;
    else { //两圆相交
        a1=acos((r1*r1+d*d-r2*r2)/(2*r1*d));//大圆中扇形圆心角的一半
        a2=acos((r2*r2+d*d-r1*r1)/(2*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 = 0;
    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) - 2 * cirinter(x1, yy, R, x2, y2, r) + cirinter(x1, yy, r, x2, y2, r);
        printf("Case #%d: %.6f\n", ++icase, ans);
    }
    return 0;
}