[补档]vijos1883 月光的魔法

vijos1883 月光的魔法

题目

传送门：https://www.vijos.org/p/1883

背景

影几欺哄了众生了

天以外——

月儿何曾圆缺

 

描述

有些东西就如同月光的魔法一般.

Luke是爱着vijos的.

他想为自己心爱的东西画些什么.

就画N个圆吧.

把它们的圆心都固定在x轴上.

圆与圆.

为了爱,两两不能相交.

为了爱,它们可以互相贴在一起.

内切或外切,都是允许的.

vijos的美丽,在于人心.

vijos的孩子们,一定能告诉大家:Luke画的圆究竟把平面分割成了多少块?

月光恬美地洒在大地上.

Luke知道,如果什么都不画,平面就只有一块.多美呢!

Luke知道,只画一个圆,平面就分成了两块.也很美呢!

但Luke还是要多画一些的,因为他真的深爱着vijos呢.

INPUT

输入数据第一行:输出一个整数N,1<=N<=300,000.表示圆的个数.

之后N行,每一行有2个整数,x[i]和r[i]表示圆心固定在x[i]的位置,半径为r[i].

-1,000,000,000<=x[i]<=1,000,000,000

1<=r[i]<=1,000,000,000

所有圆都是唯一的,不会出现重叠.

OUTPUT

输出只有一行,要求输出平面被分割成了多少块.

SAMPLE

INPUT

4

7 5

-9 11

11 9

0 20

OUTPUT

6

解题报告

本次考试觉得唯一能A的一道题 （唯一能拿分的题吧啊喂），然而还是打挂了- -

正解：

我们考虑一下每个圆对答案的贡献，当它只是单独的一个圆时，它只把整个平面分割成圆内与圆外两部分，故贡献为1，同理，各种内切与外切也没有什么影响，但我们考虑，如果一个圆被沿直径一个点不差的被分开，它自己就又被分成两个部分，贡献就为2。

如何判断呢，显然10^9的坐标是无法正常处理的，所以我们需要离散，用线段树维护区间被覆盖的点，把圆压成线段，每加入一个线段，就判断该区间是否已经被覆盖，假如被覆盖，说明该圆已经被沿直径分开了，那么它的贡献为2，否则为1.

你以为这样就结束了？NONONO（莫名中二= =）。我们考虑这样一种情况，我们已经有了两个圆，一个占领了1~2，一个占领了3~4，我们现在要加入一个1~4的圆，首先，我们自然会查询1~4的总权值，得到4，是不是就会认为1~4全部覆盖然后开心的加了2？显然是错误的，2~3这一段区间并未被覆盖，但我们认为它被覆盖了，那么我们如何处理呢？

我们可以把每个点劈成两半，一个称为该点的左半点，一个称为该点的右半点，那么当两个圆切于该点时，自然可以认为，一个圆覆盖了左半点，一个圆覆盖了右半点，此时再按上面那种方法判断，就啥事也没有啦

 #include<algorithm>
 #include<iostream>
 #include<cstring>
 #include<cstdio>
 #include<map>
 using namespace std;
 inline int read(){
     int sum(),f();
     char ch(getchar());
     for(;ch<''||ch>'';ch=getchar())
         if(ch=='-')
             f=-;
     for(;ch>=''&&ch<='';sum=sum*+ch-'',ch=getchar());
     return sum*f;
 }
 map<int,int>ma;
 struct node{
     int l,r;
     friend bool operator<(const node &a,const node &b){
         return (a.r-a.l)<(b.r-b.l);
     }
 }a[];
 int n;
 int tmp[];
 int cnt;
 int sum[<<],add[<<];
 inline void pushup(int i){
     sum[i]=sum[i<<]+sum[i<<|];
 }
 inline void pushdown(int i,int len){
     if(add[i]){
         add[i<<]=add[i];
         add[i<<|]=add[i];
         sum[i<<]=add[i]*(len-(len>>));
         sum[i<<|]=add[i]*(len>>);
         add[i]=;
     }
 }
 inline void update(int ll,int rr,int c,int l,int r,int i){
     if(ll>rr)
         return;
     if(ll<=l&&r<=rr){
         add[i]=c;
         sum[i]=c*(r-l+);
         return;
     }
     pushdown(i,r-l+);
     int mid((l+r)>>);
     if(ll<=mid)
         update(ll,rr,c,l,mid,i<<);
     if(mid<rr)
         update(ll,rr,c,mid+,r,i<<|);
     pushup(i);
 }
 inline int query(int ll,int rr,int l,int r,int i){
     if(ll>rr)
         return ;
     if(ll<=l&&r<=rr)
         return sum[i];
     pushdown(i,r-l+);
     int mid((l+r)>>);
     int ret();
     if(ll<=mid)
         ret+=query(ll,rr,l,mid,i<<);
     if(mid<rr)
         ret+=query(ll,rr,mid+,r,i<<|);
     return ret;
 }
 int ans();
 int main(){
     n=read();
     for(int i=;i<=n;i++){
         int x(read()),r(read());
         a[i].l=x-r,a[i].r=x+r;
         tmp[++cnt]=a[i].l,tmp[++cnt]=a[i].r;
     }
     sort(tmp+,tmp+cnt+);
     cnt=;
     for(int i=;i<=(n<<);i++)
         if(!ma.count(tmp[i]))
             ma[tmp[i]]=++cnt;
     sort(a+,a+n+);
     for(int i=;i<=n;i++){
         a[i].l=ma[a[i].l]<<;
         a[i].r=(ma[a[i].r]<<)-;//cout<<i<<' '<<l<<' '<<r<<endl;
         if(query(a[i].l,a[i].r,,n<<,)==a[i].r-a[i].l+)
             ans+=;
         else
             ans++;
         update(a[i].l,a[i].r,,,n<<,);
     }
     printf("%d",ans);
 }

话说我当时用线段建了个图，然后dfs，现在想想也是蠢= =