两个相同的环重合部分的面积就是A大圆和B大圆相交的面积,减去A大圆和B小圆,减去A小圆和B大圆的面积,再加上两个小圆相交的面积,直接套,两圆相交的面积模板
S = A大B大 - A大B小 - A小B大 + A小B小。
#include <cstdlib>
#include <iostream>
#include <cmath>
#include <cstring>
#include <algorithm>
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<b.r)
return PI*a.r*a.r;
else
return PI*b.r*b.r;
}
else if(dist(a, b)>=(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;
}
#include <cstdlib>
#include <iostream>
#include <cmath>
#include <cstring>
#include <algorithm>
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<b.r)
return PI*a.r*a.r;
else
return PI*b.r*b.r;
}
else if(dist(a, b)>=(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;
}