BZOJ 2178: 圆的面积并【辛普森积分】

最新推荐文章于 2020-01-13 20:23:43 发布

Master.Yi

最新推荐文章于 2020-01-13 20:23:43 发布

阅读量843

点赞数

分类专栏：计算几何

本文链接：https://blog.csdn.net/C20181220_xiang_m_y/article/details/89036271

版权

计算几何专栏收录该内容

21 篇文章

订阅专栏

本文介绍了一种利用辛普森积分方法求多个圆的面积并的技术方案。通过预先过滤掉被其他圆完全覆盖的圆来优化计算效率，显著减少了计算时间。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题面：

给定n个圆，求其面积并，n<=1000

题目分析：

直接辛普森积分上。。
有一个优化，预先把被其它圆包含的圆滤去。
（辛普森积分会求很多次一条线上的并，这个优化的效果很明显，从10000ms+ -> 2500ms）

Code：

#include<cstdio>
#include<cmath>
#include<algorithm>
#define maxn 1005
using namespace std;
struct Point{
	double x,y; Point(){}
	Point(double x,double y):x(x),y(y){}
	bool operator < (const Point &p)const{return x<p.x;}
}q[maxn];
struct circle{
	double x,y,r;
	bool operator < (const circle &p)const{return r>p.r;}
}a[maxn];
const double eps = 1e-8, inf = 1e20, Pi = acos(-1);
int n;
#define sqr(x) (x)*(x)
double f(double x){
	int top=0;
	for(int i=1;i<=n;i++) if(fabs(x-a[i].x)<a[i].r){
		double len=sqrt(sqr(a[i].r)-sqr(x-a[i].x));
		q[++top]=Point(a[i].y-len,a[i].y+len);
	}
	double l=-inf,r=-inf,ret=0;
	sort(q+1,q+1+top);
	for(int i=1;i<=top;i++)
		if(q[i].x>r) ret+=r-l,l=q[i].x,r=q[i].y;
		else r=max(r,q[i].y);
	return ret+=r-l;
}
double Simpson(double L,double M,double R,double fL,double fM,double fR,int dep){
	double M1=(L+M)/2,M2=(M+R)/2,fM1=f(M1),fM2=f(M2);
	double g1=(M-L)*(fL+4*fM1+fM)/6, g2=(R-M)*(fM+4*fM2+fR)/6, g=(R-L)*(fL+4*fM+fR)/6;
	if(dep>11&&fabs(g-g1-g2)<1e-8) return g;
	return Simpson(L,M1,M,fL,fM1,fM,dep+1)+Simpson(M,M2,R,fM,fM2,fR,dep+1);
}
inline double dist(circle &a,circle &b){return sqrt(sqr(a.x-b.x)+sqr(a.y-b.y));}
void prework(){
	int cnt=0;
	sort(a+1,a+1+n);
	for(int i=1;i<=n;i++){
		for(int j=1;j<=cnt;j++) if(a[i].r+dist(a[i],a[j])<=a[j].r) goto he;
		a[++cnt]=a[i]; he:;
	}
	n=cnt;
}
int main()
{
	scanf("%d",&n);
	double l=inf,r=-inf,mid;
	for(int i=1;i<=n;i++) scanf("%lf%lf%lf",&a[i].x,&a[i].y,&a[i].r),l=min(l,a[i].x-a[i].r),r=max(r,a[i].x+a[i].r);
	mid=(l+r)/2; prework();
	printf("%.3f",Simpson(l,mid,r,f(l),f(mid),f(r),0));
}