HDU - 5299 Circles Game 题解

题意

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=5299
给出n个不会相交不会相切的圆。Alice先手，和Bob进行博弈，每次操作可以选择一个圆，把这个圆和它所包含的圆给删掉，最后不能操作的人输。

思路

就是一个树上删边游戏，但是直接做复杂度O（n²）会TLE，所以需要圆的扫描线处理，把圆的关系转换成图，由大的圆连向它所包含的圆建边。

圆的扫描线：
把一个圆拆成两个事件，最左边进入事件，最右边离开事件。
用set来维护当前的圆的集合，set中按y轴从小到大排序。
按时间点，对所有事件按照x轴排序，然后逐一处理这些事件。当碰到进入事件时，把这个圆放入set，碰到离开事件，把这个圆从set中删掉。
然后对于每一个进入事件，从set中找到离当前圆的上方最近的圆和下方最近的圆。如果有圆包含它那就由大圆连向它，否则就是和旁边的圆同一级，和他们有同一个父亲。
（这个可能讲得不太清楚，我也是刚昨天学的）

为了方便，在最外面用一个无穷大的圆把他们都包含进去。
树上删边游戏：
叶子节点的sg值为0，其他节点的sg值等于所有儿子的（sg+1）的异或和。
根节点的sg值不为0则先手胜。

AC代码

#include<iostream>
#include<algorithm>
#include<set>
#include<vector>
#include<list>
#include<queue>
#include<cmath>
#include<cstring>
#include<sstream>
#include<cstdio>
#include<ctime>
#include<map>
#include<bitset>
#include<stack>
#include<string>
using namespace std;
#define pii pair<int,int>
#define pll pair<ll,ll>
#define ll long long
#define pb push_back
#define ull unsigned ll
#define rg register
#define bug puts("bug");
const long long INF=0x3f3f3f3f3f3f3f3fLL;
const int inf = 0x3f3f3f3f;
const int maxn =2e4+10;
const double eps=1e-6;
inline ll rd()
{
    ll x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+(ch^48);ch=getchar();}
    return x*f;
}
int n;
struct Circle{
	int x,y,r;
	void read()
	{
		x=rd(),y=rd(),r=rd();
	}
	bool friend operator<(const Circle& n1,const Circle& n2)
	{
		return n1.x<n2.x;
	}
}aa[maxn];
int now;
struct Event{
	int id,x,k;
	Event(int _id=0,int _x=0,int _k=0)
	{
		id=_id;
		x=_x;
		k=_k;
	}
	bool friend operator<(const Event& n1,const Event& n2)
	{
		double x1=aa[n1.id].y+n1.k*sqrt(aa[n1.id].r*aa[n1.id].r-(now-aa[n1.id].x)*(now-aa[n1.id].x));
		double x2=aa[n2.id].y+n2.k*sqrt(aa[n2.id].r*aa[n2.id].r-(now-aa[n2.id].x)*(now-aa[n2.id].x));
		if(x1!=x2)return x1<x2;
		return n1.k<n2.k;
	}
}event[maxn*2];
bool cmp(const Event& n1,const Event& n2)
{
	return n1.x<n2.x;
}
vector<int>ve[maxn];
set<Event>se;
int fa[maxn];
int sg(int x)
{
	int ans=0;
	for(auto y:ve[x])
	{
		ans^=sg(y)+1;
	}
	return ans;
}
int main()
{
	int t=rd();
	while(t--)
	{
		n=rd();
		se.clear();
		memset(fa,0,sizeof fa);
		for(int i=0;i<=n;i++)ve[i].clear();
		for(int i=1;i<=n;i++)aa[i].read();
		for(int i=1;i<=n;i++)
		{
			event[2*i-1]=Event(i,aa[i].x-aa[i].r,1);
			event[2*i]=Event(i,aa[i].x+aa[i].r,-1);
		}
		sort(event+1,event+1+2*n,cmp);
		for(int i=1;i<=2*n;i++)
		{
			now=event[i].x;
			if(event[i].k==1)
			{
				set<Event>::iterator p=se.upper_bound(event[i]);
				if(p!=se.end())
				{
					if(p->k==1)fa[event[i].id]=p->id;
					else fa[event[i].id]=fa[p->id];
				}
				se.insert(event[i]);
				se.insert(Event(event[i].id,event[i].x+2*aa[event[i].id].r,-1));
			}
			else
			{
				se.erase(Event(event[i].id,event[i].x-2*aa[event[i].id].r,1));
				se.erase(event[i]);
			}
		}
		for(int i=1;i<=n;i++)ve[fa[i]].pb(i);
		puts(sg(0)?"Alice":"Bob");
	}
}

在所有扫描线的板子中找了个最短的）