Problem 遗产
题目大意
一个平面上n个圆，任两个圆只会相离或包含，给出每个圆位置与半径。
alice&&bob轮流取圆，每取一个就可以取出这个圆以及被这个圆包含的圆。
没圆取的人输，alice先取，问谁有必胜策略。
Solution
Method #1
首先我们考虑一个暴力一点的写法：
先将圆半径从小到大排序，然后枚举两个圆，判断是否包含关系。
如果包含的话就将大圆连一条边到小圆。
很容易发现，这样执行完以后得到了一棵树。
我们接下来要做的事情就变成了：
每个人可以轮流从树上删除一条边，所有与根节点不相连的点全都被删除。
谁不能取谁就输。
预备知识：树上删边
叶子节点值为0；
除叶子节点外任意节点值为其所有儿子的值加一的异或和。
根节点值为0则后取者胜。
dfs一遍树即为答案。
接下来很高兴啊！这题A掉啦！
但是——这smg？
$O(n^2)$的算法a掉了20000？？还是700ms？？
hdu的数据怎么不上天？
自己的oj交一发，妥妥的30pts tle。
Method #2
好了我们来换思路吧。
首先还是对半径进行排序。
我们对于每一个圆，从左到右扫描这个圆的范围。
如果在这个范围内找到了点，那么就判断是否包含，若包含则连边。
大致方法与上面一样。
我们对于每一个x坐标开一个vector记录该坐标线上的点。
提交hdu，ac100ms。
自己oj，90pts，还是被卡掉了。
仔细思考扫了一下，构造数据的确是可以卡掉这个算法的。
1 #include <iostream>
2 #include <cstdio>
3 #include <algorithm>
4 #include <cstring>
5 #include <vector>
6 using namespace std;
7 vector<int> f[80010];
8 int n,T,h[80010],tot=0,t;
9 //bool cheat;
10 struct Circle{
11 int x,y,r;
12 bool inc;
13 }c[80010];
14 struct node{
15 int next,to;
16 }a[80010];
17 bool cmp(Circle a,Circle b){
18 return a.r<b.r;
19 }
20 void add(int u,int v){
21 a[++tot].to=v;
22 a[tot].next=h[u];
23 h[u]=tot;
24 }
25 int dfs_sg(int x,int fa){
26 int ret=0;
27 for(int i=h[x];~i;i=a[i].next)
28 if(a[i].to!=fa)ret^=(dfs_sg(a[i].to,x)+1);
29 return ret;
30 }
31 int main(){
32 // freopen("08b.in","r",stdin);
33 scanf("%d",&T);
34 while(T--){
35 tot=0,cheat=1;
36 memset(a,0,sizeof(a));
37 memset(c,0,sizeof(c));
38 memset(h,-1,sizeof(h));
39 for(int i=0;i<=80001;i++)f[i].clear();
40 scanf("%d",&n);
41 for(int i=1;i<=n;i++)
42 scanf("%d%d%d",&c[i].x,&c[i].y,&c[i].r),
43 c[i].x+=40000,c[i].y+=40000;
44 // for(int i=1;i<=100;i++)
45 // if(c[i].r!=1){
46 // cheat=0;
47 // break;
48 // }
49 // if(cheat){
50 // printf("Bob\n");
51 // continue;
52 // }
53 sort(c+1,c+1+n,cmp);
54 for(int i=1;i<=n;i++){
55 for(int p=c[i].x-c[i].r+1;p<=c[i].x+c[i].r-1;p++){
56 t=0;
57 while(t<f[p].size()){
58 if((c[f[p][t]].x-c[i].x)*
59 (c[f[p][t]].x-c[i].x)+
60 (c[f[p][t]].y-c[i].y)*
61 (c[f[p][t]].y-c[i].y)<=
62 (long long)((c[i].r)*(c[i].r))){
63 c[f[p][t]].inc=1;
64 add(i,f[p][t]);
65 f[p][t]=f[p][f[p].size()-1];
66 f[p].pop_back();
67 t--;
68 }
69 t++;
70 }
71 }
72 f[c[i].x].push_back(i);
73 }
74 int ans=0;
75 for(int i=1;i<=n;i++)
76 if(!c[i].inc)ans^=(dfs_sg(i,-1)+1);
77 if(!ans)printf("Bob\n");
78 else printf("Alice\n");
79 }
80 }
Method #3
我们将每个圆的左端点&&右端点出现的时间排序。
对于每一个圆的端点来说，左端点sign为+1，表示该圆开始出现，右端点sign为-1，表示该圆结束。
我们开一个set来维护这些圆与扫描线的交点。set的键值为交点的y坐标值，如下计算：
$c_i\cap(x=now)={c_i.y\pm \sqrt{(c_i.r)^2-(c_i.x-now)^2}} $
只要给set开一个operater就可以了。
对于圆$c_i$的起始交点，我们get到它的upper bound
upperbound交点有三种可能，
一种是没有该交点，这说明$c_i$一定是不被其他圆包含的圆。
一种是该交点所属的圆$c_j$包含$c_i$，这说明$c_i$的父亲是圆$c_j$，而且不会有更近的父亲。
还有是圆$c_j$不包含圆$c_i$。这说明$c_i$圆的父亲是圆$c_j$的父亲的儿子。
第三种可能我们需要特判一下，如果圆$c_j$的父亲为0，那么不就不用连边。否则会爆空间。
判断sign为正的时候insert，sign为负的时候erase即可。
建树完毕后dfs即可。
hdu 78ms，
oj终于a掉了。
强烈吐槽hdu的数据。
最垃圾的暴力都能a掉，出数据的人可能比较强吧不屑于出这种题目的数据。
1 #include <iostream>
2 #include <cstdio>
3 #include <cstring>
4 #include <algorithm>
5 #include <cmath>
6 #include <set>
7 using namespace std;
8 int n,h[80010],fa[80010],tot=0,now,T,ndtot=0;
9 struct Node{
10 int next,to;
11 }a[80010];
12 struct Circle{
13 int x,y,r;
14 }c[80010];
15 struct Time{
16 int x,id;
17 friend inline bool operator <(Time x,Time y){
18 return x.x<y.x;
19 }
20 }t[80010];
21 struct node{
22 int sign,id;
23 inline double y(){
24 return c[id].y+sign*sqrt((c[id].r)*(c[id].r)
25 -(now-c[id].x)*(now-c[id].x));
26 }
27 friend inline bool operator <(node x,node y){
28 return x.id==y.id?x.sign<y.sign:x.y()<y.y();
29 }
30 };
31 int dfs_sg(int x,int fa){
32 int ret=0;
33 for(int i=h[x];~i;i=a[i].next)
34 if(a[i].to!=fa)ret^=(dfs_sg(a[i].to,x)+1);
35 return ret;
36 }
37 void add(int u,int v){
38 a[++ndtot].to=v;
39 a[ndtot].next=h[u];
40 h[u]=ndtot;
41 }
42 set<node> S;
43 int main(){
44 // freopen("08b.in","r",stdin);
45 scanf("%d",&T);
46 while(T--){
47 memset(a,0,sizeof(a));
48 memset(c,0,sizeof(c));
49 memset(t,0,sizeof(t));
50 memset(h,-1,sizeof(h));
51 memset(fa,0,sizeof(fa));
52 S.clear();tot=0;ndtot=0;
53 scanf("%d",&n);
54 for(int i=1;i<=n;i++)scanf("%d%d%d",&c[i].x,&c[i].y,&c[i].r);
55 for(int i=1;i<=n;i++)
56 t[++tot]=(Time){c[i].x-c[i].r,i},
57 t[++tot]=(Time){c[i].x+c[i].r,-i};
58 sort(t+1,t+tot+1);
59 for(int i=1;i<=tot;i++){
60 now=t[i].x;
61 if(t[i].id>0){
62 set<node>::iterator it=S.upper_bound((node){+1,t[i].id});
63 if(it==S.end()) fa[t[i].id]=0;
64 else if(it->sign==1) add(it->id,t[i].id),fa[t[i].id]=it->id;
65 else{
66 if(fa[it->id])add(fa[it->id],t[i].id);
67 fa[t[i].id]=fa[it->id];
68 }
69 S.insert((node){-1,t[i].id}),S.insert((node){+1,t[i].id});
70 }
71 else S.erase((node){-1,-t[i].id}),S.erase((node){+1,-t[i].id});
72 }
73 int ans=0;
74 for(int i=1;i<=n;i++)
75 if(!fa[i])ans^=(dfs_sg(i,-1)+1);
76 if(!ans)printf("Bob\n");
77 else printf("Alice\n");
78 }
79 return 0;
80 }