题意
先上一道模板题:SP8073 The area of the union of circles
通过点和半径的方式给出平面上的一些圆,求这些圆覆盖的总面积。
题解
和这道题类似的还有P1222 三角形 。不过这道题由于一些特殊的性质,直接扫描线就可以过,并且效率远远高于自适应辛普森。
还是来看这道题。我们可以把平面上的一堆相交的圆看成一个整体函数,我们要求的就是每一段函数的面积和。很明显可以化成积分的形式,于是用自适应辛普森积分乱搞。
那么我们怎么求一个 x x x所对应的 f ( x ) f(x) f(x)值呢?很明显我们在 x x x确定的时候要求的就是直线 y = x y=x y=x被圆覆盖的总长度。那么我们可以枚举所有的圆,只要与直线 y = x y=x y=x有交点的圆,把两个交点之间的部分看成一条线段,然后贪心求线段覆盖长度即可。
一些小细节:自适应辛普森算法中多加几个参数,可以减少重复的计算。同时可以在开始求面积之前把被一个大圆覆盖的小圆先删除。然后对于每一堆通过交点联通的圆,求出最左端和最右端,然后进行自适应辛普森积分即可。
代码
#include <bits/stdc++.h>
#define MAX 2005
#define ll long long
#define eps 1e-7
#define INF 1e16
using namespace std;
struct seg{
double l, r;
friend bool operator <(seg a, seg b){
if(a.l == b.l) return a.r < b.r;
return a.l < b.l;
}
};
struct vec{
ll x, y;
};
double dis(vec a, vec b){
return sqrt(1.0*(a.x-b.x)*(a.x-b.x) +