bzoj1043[HAOI2008]下落的圆盘 计算几何

1043: [HAOI2008]下落的圆盘

Time Limit: 10 Sec  Memory Limit: 162 MB
Submit: 1598  Solved: 676
[Submit][Status][Discuss]

Description

　　有n个圆盘从天而降，后面落下的可以盖住前面的。求最后形成的封闭区域的周长。看下面这副图, 所有的红
色线条的总长度即为所求. 

Input

　　第一行为1个整数n,N<=1000
接下来n行每行3个实数,ri,xi,yi,表示下落时第i个圆盘的半径和圆心坐标.

Output

　　最后的周长，保留三位小数

Sample Input

2
1 0 0
1 1 0

Sample Output

10.472

HINT

 

Source

可以考虑后放的盘子对先放的盘子造成的覆盖影响
枚举每个盘，这个盘贡献的答案就是它的周长减去被后面的盘所覆盖的部分
圆和圆相交部分长度可以转化成线段覆盖问题来计算：算出中轴线的角度x，再计算圆心到2交点的角度y，

圆上被覆盖的弧度区间就是[x-y,x+y]
几个圆与这个圆相交，把相交覆盖的弧转化成圆上的线段，做一个线段覆盖问题就可以解决

#include<bits/stdc++.h>
#define reg register
#define N 1005
using namespace std;
int n,tp;double x[N],y[N],r[N],ans;const double pi=acos(-1);struct node{double l,r;}q[N<<1];
double dis(int i,int j){return sqrt((x[i]-x[j])*(x[i]-x[j])+(y[i]-y[j])*(y[i]-y[j]));}
bool con(int i,int j){return r[i]>=r[j]+dis(i,j);}
bool cmp(node a,node b){return a.l<b.l;}
void push(int i,int j){
    double d,k,t,mv;
    d=dis(i,j);k=atan2(y[j]-y[i],x[j]-x[i]);
    t=(d*d+r[i]*r[i]-r[j]*r[j])/2/d;
    mv=acos(t/r[i]);
    q[++tp]=(node){k-mv,k+mv};
}
double calc(int p){
    tp=0;
    for(reg int i=p+1;i<=n;i++)
    if(con(i,p))return 0;
    for(int i=p+1;i<=n;i++){
        if(con(p,i)||dis(i,p)>r[i]+r[p])continue;
        push(p,i);
    }
    for(reg int i=1;i<=tp;i++){
        if(q[i].l<0)q[i].l+=2*pi;
        if(q[i].r<0)q[i].r+=2*pi;
        if(q[i].l>q[i].r){
            q[++tp]=(node){0,q[i].r};
            q[i].r=2*pi;
        }
    }
    sort(q+1,q+1+tp,cmp);
    double cov=0,mx=q[1].l;
    for(reg int i=1;i<=tp;i++){
        mx=max(mx,q[i].l);
        if(q[i].r<=mx)continue;
        cov+=q[i].r-mx;mx=q[i].r;
    }
    double ret=r[p]*(2*pi-cov);
    return ret;
}

int main(){
    scanf("%d",&n);
    for(reg int i=1;i<=n;i++)
    scanf("%lf%lf%lf",&r[i],&x[i],&y[i]);
    for(reg int i=1;i<=n;i++)ans+=calc(i);
    printf("%.3lf
",ans);
    return 0;
}