hdu 5120 Intersection(求相交圆的面积)

两个相同的环重合部分的面积就是A大圆和B大圆相交的面积,减去A大圆和B小圆,减去A小圆和B大圆的面积,再加上两个小圆相交的面积,直接套,两圆相交的面积模板

S = A大B大 - A大B小 - A小B大 + A小B小。

#include 
#include 
#include 
#include 
#include 

using namespace std;

const double PI = acos(-1);

// 圆的结构 
struct circle {
       double x; // 圆心横坐标 
       double y; // 圆心纵坐标 
       double r; // 半径 
};

// 计算圆心距 
double dist(circle a, circle b)
{
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}

// 计算两相交圆的面积 
double area(circle a, circle b) {

       if((dist(a, b)+min(a.r,b.r))<=max(a.r,b.r)) // 内含或重合 
        { 
            if(a.r=(a.r+b.r)) // 相离或相切 
            return 0.0;
        else // 相交 
        {
            double length=dist(a, b);
            // 利用三角形余弦定理求圆心角 
            double d1=2*acos((a.r*a.r+length*length-b.r*b.r)/(2*a.r*length)); 
            double d2=2*acos((b.r*b.r+length*length-a.r*a.r)/(2*b.r*length));
            // 利用圆心角求得扇形面积再减去三角形面积后两部分相加就是相交部分面积 
            double area1=a.r*a.r*d1/2-a.r*a.r*sin(d1)/2;
            double area2=b.r*b.r*d2/2-b.r*b.r*sin(d2)/2;
            double area=area1+area2;
            return area;
        }
}
int main()
{
    int T;
    circle a, b, A, B;
    double r, R;
    scanf("%d", &T) ;
    for(int i=1; i<=T; i++)
    {
        scanf("%lf%lf", &r, &R);   
        a.r = b.r = r;
        A.r = B.r = R;  
        scanf("%lf%lf", &a.x, &a.y);
        A.x = a.x;
        A.y = a.y;
        scanf("%lf%lf", &b.x, &b.y);
        B.x = b.x;
        B.y = b.y;     
        double sec = area(A, B) - area(A, b) - area(B, a) + area(a, b);
        printf("Case #%d: %.6lf\n", i, sec);
          
    }
    
    return 0;
}

你可能感兴趣的:(hdu 5120 Intersection(求相交圆的面积))