Codeforces 814D An overnight dance in discotheque 贪心

Codeforces 814D An overnight dance in discotheque
物语场!!!!之后的可怜和命运石之门……哈哈哈哈哈哈哈……

题意

平面直角坐标系中有n个圆,每两个圆之间最多只有一个交点.要求你把圆分成两份,
每一份的圆中加上被覆盖偶数次的圆的面积,减去被覆盖奇数次的圆的面积,求两份加起来最大的面积和.

题解

由于每两个圆之间最多只有一个交点,如果一个圆$a$被另一个圆$b$包含,我们就令$b$是$a$的父亲.
那么这$n$个圆之间必定形成了一个森林的结构,每一组嵌套的圆之间都是一棵树.
显然无关的两棵树之间是没有任何关系的.我们仅仅考虑一棵树.
这棵树肯定包括很多层,每一层都有几个节点,但是每一层的圆的面积和肯定是没有上一层大的.
故此可以贪心.
那么很显然第一层的圆的面积肯定要加上去,第二层也要.
第三层开始不管放在哪里面积都要减掉,随便选择一个.
而第四层开始不能放在与第三层相异的一边,否则会被减掉,而后面不管怎么加都没有第四层面积大,所以第四层要和第三层放在一起.
而第五层和第六层也要放在一起……等等等等.
那么我们可以用$n^2$预处理出每一个点的深度.
然后把深度为0和奇数的圆的面积全部加起来,减掉深度为正偶数的圆的面积即为答案.

#include<bits/stdc++.h> //Ithea Myse Valgulious
namespace chtholly{
typedef long long ll;
#define re0 register int
#define rec register char
#define rel register ll
#define gc getchar
#define pc putchar
#define p32 pc(' ')
#define pl puts("")
/*By Citrus*/
inline int read(){
  int x=0,f=1;char c=gc();
  for (;!isdigit(c);c=gc()) f^=c=='-';
  for (;isdigit(c);c=gc()) x=(x<<3)+(x<<1)+(c^'0');
  return f?x:-x;
  }
template <typename mitsuha>
inline bool read(mitsuha &x){
  x=0;int f=1;char c=gc();
  for (;!isdigit(c)&&~c;c=gc()) f^=c=='-';
  if (!~c) return 0;
  for (;isdigit(c);c=gc()) x=(x<<3)+(x<<1)+(c^'0');
  return x=f?x:-x,1;
  }
template <typename mitsuha>
inline int write(mitsuha x){
  if (!x) return 0&pc(48);
  if (x<0) x=-x,pc('-');
  int bit[20],i,p=0;
  for (;x;x/=10) bit[++p]=x%10;
  for (i=p;i;--i) pc(bit[i]+48);
  return 0;
  }
inline char fuhao(){
  char c=gc();
  for (;isspace(c);c=gc());
  return c;
  }
}using namespace chtholly;
using namespace std;
typedef double db;
const db pi=3.141592653589;
struct yuan{
int x,y,r;
void rd(){x=read(),y=read(),r=read();}
db area(){return pi*r*r;}
}aoi[1010];
int cnt[1010];
ll dist(int a,int b,int c,int d){return 1ll*abs(c-a)*abs(c-a)+1ll*abs(d-b)*abs(d-b);}
bool judge(yuan a,yuan b){return a.r<b.r&&dist(a.x,a.y,b.x,b.y)<1ll*(a.r+b.r)*(a.r+b.r);}

int main(){
int i,j,n=read();
for (i=1;i<=n;++i) aoi[i].rd();
for (i=1;i<=n;++i,pl){
  for (j=1;j<=n;++j){
    if (judge(aoi[i],aoi[j])) cnt[i]++;
    }
  }
db llx=0;
for (i=1;i<=n;++i) llx+=!cnt[i]||cnt[i]&1?aoi[i].area():-aoi[i].area();
printf("%.12lf",llx);
}