【圆的面积并】O(n^2logn)的解法

题解完全参考的AekdyCoin牛、lishi牛和盾哥的题解，具体方法见这些人的文章，这里不罗嗦了

这里记几个值得一提的细节

1.叉积求面积记得除以2，还有扇形面积公式S = r^2 * p / 2，p是圆心角(弧度)，刚好也要除以2

2.判好两个圆完全相同的情况，我用的是一个标记数组表示i有没有曾经被包含过，这个包含的意义是要取等号的，也就是说两个完全相同的圆，出现在前面的那个会被判为被包含，后面那个则不会

3.当区间跨极角的边界时，我的做法是拆成两个区间

4.我的两圆求交点的方法，如图：

将O₁当做极点，先求出O₂的极角，再用余弦定理求出角AO₁O₂，然后就能得到A的极角，由于r₁已知，所以可以根据极角求出A的坐标，另外一个交点类似。由于实际上我们需要求的是“有向线段”，所以用这个方法就能很方便的直接确定这条边的方向——由减去角AO₁O₂的指向加上的。但是这个方法由于有反三角函数的使用，所以精度不怎么样，需要用精度高的extended才能达到6位精度的要求

另外，这个算法实际上是很好实现的，我的代码100行左右

program syj; uses math; const e=1e-10; type real=extended; point=record x,y:real; end; circle=record o:point;r:real; end; var n,i,j,t:longint; a:array[0..205]of circle; d,ans,r1,r2,p,q,z:real; l,r:array[0..405]of real; bh:array[0..205]of boolean; u,v:point; function dis(var a,b:point):real; begin dis:=sqrt(sqr(a.x-b.x)+sqr(a.y-b.y)); end; operator -(a,b:point)c:point; begin c.x:=a.x-b.x;c.y:=a.y-b.y; end; operator *(a,b:point)c:real; begin c:=a.x*b.y-a.y*b.x; end; procedure sort(s,t:longint); var i,j:longint; begin i:=s;j:=t;l[0]:=l[(s+t)>>1]; repeat while l[0]-l[i]>e do inc(i); while l[j]-l[0]>e do dec(j); if i<=j then begin z:=l[i];l[i]:=l[j];l[j]:=z; z:=r[i];r[i]:=r[j];r[j]:=z; inc(i);dec(j); end; until i>j; if s<j then sort(s,j); if i<t then sort(i,t); end; function cal(a:point;k,r:real):point; begin cal.x:=a.x+cos(k)*r;cal.y:=a.y+sin(k)*r; end; begin assign(input,'circle.in');reset(input); assign(output,'circle.out');rewrite(output); // init(10); readln(n); for i:=1 to n do begin readln(a[i].o.x,a[i].o.y,a[i].r); // putcircle(round(a[i].o.x)+5,round(a[i].o.y)+5,round(a[i].r)); end; for i:=1 to n do begin // writeln(ans:0:6); t:=0; for j:=1 to n do if (j<>i)and not bh[j] then begin d:=dis(a[i].o,a[j].o); r1:=a[i].r;r2:=a[j].r; if d-r1-r2>-e then continue; if d-abs(r1-r2)<e then begin bh[i]:=bh[i] or(r2-r1>-e);continue; end; inc(t); p:=arctan2(a[j].o.y-a[i].o.y,a[j].o.x-a[i].o.x); q:=arccos((r1*r1+d*d-r2*r2)/(2*r1*d)); l[t]:=p-q;r[t]:=p+q; if l[t]+pi<-e then begin inc(t); l[t]:=l[t-1]+pi*2;r[t]:=pi; l[t-1]:=-pi; end; if r[t]-pi>e then begin inc(t); r[t]:=r[t-1]-pi*2;l[t]:=-pi; r[t-1]:=pi; end; end; if (t=0)and not bh[i] then ans:=ans+pi*r1*r1; if (t=0)or bh[i] then continue; sort(1,t); l[0]:=0; l[t+1]:=l[1]+pi*2;r[t+1]:=r[1]+pi*2; p:=l[1]; for j:=1 to t+1 do begin if l[j]-p>e then begin u:=cal(a[i].o,p,a[i].r);v:=cal(a[i].o,l[j],a[i].r); ans:=ans+(u*v+a[i].r*a[i].r*(l[j]-p)-(u-a[i].o)*(v-a[i].o))*0.5; end; if r[j]>p then p:=r[j]; end; end; writeln(ans:0:6); close(input);close(output); end.  

Ps：调试几何题有个画图器真是爽啊！