http://172.20.6.3/Problem_Show.asp?id=1502
這種題用了快一天才寫出來也是真的辣雞。主要思路就是計算一下被擋住的弧度然後對弧度進行貪心。
最開始比較困擾的是求弧度值及其起始位置的部分，弧度值很好求，位置有點惡心，我的起始位置設置的是圓的十二點方向順時針到起始位置的弧度值，然後我分了四種情況討論（因為遮擋的方向有兩種不同情況，遮擋部分弧度值與180度的關系又是兩種情況），應該是有更簡單的方法的，但是我只能想出來這種了。。。

代碼

 1 #include<cstdio>
 2 #include<cstring>
 3 #include<algorithm>
 4
 
 #include<cmath>
 5 #include<iostream>
 6 using namespace std;
 7 const int maxn=1010;
 8 const double eps=0.0000001;
 9 long long n;
10 double r[maxn]={},x[maxn]={},y[maxn]={};
11 struct nod{
12     double id,v;
13 }e[maxn];
14 int tot;
15 bool mmp(nod aa,nod bb){
16     return aa.id<bb.id;
17
 
 }
18 double getit(int i){
19     tot=0;double cnt=0;
20     for(int j=n;j>i;j--){
21         double w=sqrt((x[i]-x[j])*(x[i]-x[j])+(y[i]-y[j])*(y[i]-y[j]));
22         if(r[i]+r[j]-w<=0||r[i]-r[j]-w>=0)continue;
23         if(r[j]-r[i]-w>=0) return 0;
24         
25         double z=(r[i]*r[i]-r[j]*r[j]+w*w)/(2
 
*w);
26         double h=acos(z/r[i])*2;
27         
28         double xx=x[i]+(x[j]-x[i])/w*z,yy=y[i]+(y[j]-y[i])/w*z;
29         double w1=sqrt((x[i]-xx)*(x[i]-xx)+(y[i]+z-yy)*(y[i]+z-yy));
30         double h2=asin(w1/2/z)*2;
31         if(xx<x[i]){
32             if(h2<0)h2=-h2;
33             else if(h2>0) h2=2*M_PI-h2; }
34         else if (xx>x[i]){if(h2<0) h2=2*M_PI+h2;}
35         h2=h2-h/2;
36         if(h2<0)h2+=2*M_PI;
37         h=h2+h;
38         if(h>2*M_PI){
39             e[++tot].id=h2;e[tot].v=2*M_PI;
40             e[++tot].id=0;e[tot].v=h-2*M_PI;
41         }else{e[++tot].id=h2;e[tot].v=h;}
42     }
43     sort(e+1,e+1+tot,mmp);
44     double ll=0.0,rr=0.0;
45     for(int j=1;j<=tot;j++){
46         if(e[j].id-rr<0){
47             if(e[j].v-rr>0){
48                 rr=e[j].v;
49             }
50         }
51         else if(e[j].v-rr>0){
52             cnt+=rr-ll;
53             rr=e[j].v;ll=e[j].id;
54         }
55     }
56     cnt+=rr-ll;
57     return r[i]*(M_PI*2-cnt); 
58 }
59 int main(){
60     //freopen("wtf.in","r",stdin);
61     scanf("%d",&n);
62     for(int i=1;i<=n;i++){
63         scanf("%lf%lf%lf",&r[i],&x[i],&y[i]);
64     }double ans=0;
65     for(int i=n;i>=1;i--){
66         ans+=getit(i);
67     }printf("%.3f\n",ans);
68     return 0;
69 }

JZYZOJ1502 [haoi2008]下落的圓盤 計算幾何 貪心