JZYZOJ1502 [haoi2008]下落的圆盘 计算几何 贪心

http://172.20.6.3/Problem_Show.asp?id=1502
这种题用了快一天才写出来也是真的辣鸡。主要思路就是计算一下被挡住的弧度然后对弧度进行贪心。
最开始比较困扰的是求弧度值及其起始位置的部分，弧度值很好求，位置有点恶心，我的起始位置设置的是圆的十二点方向顺时针到起始位置的弧度值，然后我分了四种情况讨论（因为遮挡的方向有两种不同情况，遮挡部分弧度值与180度的关系又是两种情况），应该是有更简单的方法的，但是我只能想出来这种了。。。

代码

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<iostream>
using namespace std;
const int maxn=1010;
const double eps=0.0000001;
long long n;
double r[maxn]={},x[maxn]={},y[maxn]={};
struct nod{
    double id,v;
}e[maxn];
int tot;
bool mmp(nod aa,nod bb){
    return aa.id<bb.id;
}
double getit(int i){
    tot=0;double cnt=0;
    for(int j=n;j>i;j--){
        double w=sqrt((x[i]-x[j])*(x[i]-x[j])+(y[i]-y[j])*(y[i]-y[j]));
        if(r[i]+r[j]-w<=0||r[i]-r[j]-w>=0)continue;
        if(r[j]-r[i]-w>=0) return 0;
        
        double z=(r[i]*r[i]-r[j]*r[j]+w*w)/(2*w);
        double h=acos(z/r[i])*2;
        
        double xx=x[i]+(x[j]-x[i])/w*z,yy=y[i]+(y[j]-y[i])/w*z;
        double w1=sqrt((x[i]-xx)*(x[i]-xx)+(y[i]+z-yy)*(y[i]+z-yy));
        double h2=asin(w1/2/z)*2;
        if(xx<x[i]){
            if(h2<0)h2=-h2;
            else if(h2>0) h2=2*M_PI-h2; }
        else if (xx>x[i]){if(h2<0) h2=2*M_PI+h2;}
        h2=h2-h/2;
        if(h2<0)h2+=2*M_PI;
        h=h2+h;
        if(h>2*M_PI){
            e[++tot].id=h2;e[tot].v=2*M_PI;
            e[++tot].id=0;e[tot].v=h-2*M_PI;
        }else{e[++tot].id=h2;e[tot].v=h;}
    }
    sort(e+1,e+1+tot,mmp);
    double ll=0.0,rr=0.0;
    for(int j=1;j<=tot;j++){
        if(e[j].id-rr<0){
            if(e[j].v-rr>0){
                rr=e[j].v;
            }
        }
        else if(e[j].v-rr>0){
            cnt+=rr-ll;
            rr=e[j].v;ll=e[j].id;
        }
    }
    cnt+=rr-ll;
    return r[i]*(M_PI*2-cnt); 
}
int main(){
    //freopen("wtf.in","r",stdin);
    scanf("%d",&n);
    for(int i=1;i<=n;i++){
        scanf("%lf%lf%lf",&r[i],&x[i],&y[i]);
    }double ans=0;
    for(int i=n;i>=1;i--){
        ans+=getit(i);
    }printf("%.3f
",ans);
    return 0;
}