# BZOJ-1043 [HAOI2008]下落的圆盘

#include <cstdlib>
#include <cstdio>
#include <iostream>
#include <cstring>
#include <cmath>
#include <cctype>
#include <algorithm>
#define rep(i, l, r) for(int i=l; i<=r; i++)
#define pi acos(-1)
#define maxn 1234
typedef long long ll;
using namespace std;
{
int x=0, f=1; char ch=getchar();
while (!isdigit(ch)) {if (ch=='-') f=-1; ch=getchar();}
while (isdigit(ch)) x=x*10+ch-'0', ch=getchar();
return x*f;
}
struct line{double l, r;} q[maxn*2];
bool cmp(line a, line b){return a.l<b.l || (a.l==b.l && a.r<b.r);}
double r[maxn], x[maxn], y[maxn];
double dis(int a, int b){return sqrt((x[a]-x[b])*(x[a]-x[b])+(y[a]-y[b])*(y[a]-y[b]));}
int main()
{
rep(i, 1, n) scanf("%lf%lf%lf", &r[i], &x[i], &y[i]);
rep(a, 1, n)
{
int top=0;
bool end=0; rep(b, a+1, n) if (r[a]+dis(a, b)<=r[b]) end=1; if (end) continue;
rep(b, a+1, n) if (r[b]+dis(a, b)>r[a] && r[a]+r[b]>dis(a, b))
{
double r1=r[a], r2=r[b], d=dis(a, b), t, v;
t=acos((r1*r1-r2*r2+d*d)/(2*d)/r1);
v=atan2((x[a]-x[b]),(y[a]-y[b]));
q[++top]=(line){v-t, v+t};
}
rep(i, 1, top)
{
if (q[i].l<0) q[i].l+=2*pi;
if (q[i].l>2*pi) q[i].l-=2*pi;
if (q[i].r<0) q[i].r+=2*pi;
if (q[i].r>2*pi) q[i].r-=2*pi;
if (q[i].l>q[i].r)
q[++top]=(line){0,q[i].r}, q[i].r=2*pi;
}
sort(q+1, q+top+1, cmp);
double now=0, tmp=0;
rep(i, 1, top)
{
if (q[i].l>now) tmp+=q[i].l-now;
now=max(now, q[i].r);
}
tmp+=2*pi-now;
ans+=tmp*r[a];
}
printf("%.3lf\n", ans);
return 0;
}
posted @ 2015-05-05 18:21  NanoApe  阅读(160)  评论(0编辑  收藏  举报