hdu 5299 Circles Game

http://acm.hdu.edu.cn/showproblem.php?pid=5299

题意：在坐标系中给出一些圆，这些圆可能相互包含和相离，现在两个人玩一个游戏，每个人每次选择去掉一个圆，如果这个圆中有小圆，同时去掉，没有圆可取的人输。

分析：等级制的威佐夫博奕
因为圆与圆有相互包含的关系，那么我们可以建一棵树，虚设一个无穷大的圆作为根节点，那么剩下的圆与大圆建立联系。
如果知道树的删边游戏的话，这道题应该是很好做的。

树的删边游戏规则如下：

1.给出一颗n个结点的树，其中有一个点规定为根。

2.两名游戏者轮流操作，每次从树中删去一条边，删去该边后，不与根节点连通的部分将被移走。

3.无法操作者输。

定理：

1.叶子结点的sg值为0；

2.分支结点的sg值为其所有子节点的sg值加1后的异或和。

有了这些知识以后就可以做了，设0结点为根，求其sg值即可。

建边可以用扫描线，也可以暴力建边

#include <algorithm>
#include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;

const int N = 11111;
int head[N];
int e;

struct Edge 
{
    int v, next;
} edge[N];

void addEdge( int u, int v )
{
    edge[e].v = v;
    edge[e].next = head[u];
    head[u] = e++;
}

struct Circle 
{
    int x, y, r;
    bool operator < ( const Circle & o ) const 
    {
        return r < o.r;
    }
} circle[N];

int dfs( int u )
{
    int res = 0;
    for ( int i = head[u]; i != -1; i = edge[i].next )
    {
        res ^= ( dfs( edge[i].v ) + 1 );
    }
    return res;
}

bool involve( int i, int j )
{
    int x1 = circle[i].x, x2 = circle[j].x;
    int y1 = circle[i].y, y2 = circle[j].y;
    int r1 = circle[i].r, r2 = circle[j].r;
    return ( x1 - x2 ) * ( x1 - x2 ) + ( y1 - y2 ) * ( y1 - y2 ) < ( r1 - r2 ) * ( r1 - r2 );
}

int main ()
{
    int t;
    scanf("%d", &t);
    while ( t-- )
    {
        int n;
        scanf("%d", &n);
        for ( int i = 1; i <= n; i++ )
        {
            scanf("%d%d%d", &circle[i].x, &circle[i].y, &circle[i].r);
        }
        e = 0;
        memset( head, -1, sizeof(head) );
        sort( circle + 1, circle + 1 + n );
        for ( int i = 1; i <= n; i++ )
        {
            bool flag = false;
            for ( int j = i + 1; j <= n; j++ )
            {
                if ( involve( i, j ) )
                {
                    addEdge( j, i );
                    flag = true;
                    break;
                }
            }
            if ( !flag )
            {
                addEdge( 0, i );
            }
        }
        int ans = dfs(0);
        if ( ans ) puts("Alice");
        else puts("Bob");
    }
    return 0;
}

因为仅要找到每个圆外侧的最内侧的那个，因此用set二分查找左右边界优化，以加边171ms，对比上面是889ms

#include <bits/stdc++.h>
using namespace std;
#define prt(k) cerr<<#k" = "<<k<<endl
typedef long long ll;
const int N = 20003;
double sqr(double x) { return x*x ; }
struct P
{
    int x, y, r; int id;
    P() {}
    P(int _x, int _y, int _r =0) { x=_x, y=_y, r=_r; }
    void out() { printf("%d : (%d, %d) ---%d\n", id, x,y,r); }
    bool operator < (P b) const
    {
            if (x - b.x) return x < b.x;
            if (y - b.y) return y < b.y;
            return r < b.r;
    }
    double dis(P b) {
                return sqrt(sqr(x-b.x) + sqr(y-b.y) );
        }
}p[N];
double dis(P a, P b) { return a.dis(b); }
bool cmp(P a, P b)
{
        return a.r > b.r;
}
bool vis[N];
struct Edge
{
    int v, next;
} e[N<<1] ;
int head[N], mm;
void add(int u, int v)
{
    e[mm].v = v;
    e[mm].next = head[u];
    head[u] = mm++;
}
int n;
const int Root = N - 1;
int fa[N];
int dfs(int u)
{
    int ret=  0;
    for (int i=head[u];~i;i=e[i].next) {
        ret ^= dfs(e[i].v)+1;
    }
    return ret;
}
set<P>::iterator it;
const int inf = 0x3f3f3f3f;
int main()
{
        int re, ca=1;
        cin>>re;
        while (re--)
        {
                memset(head, -1, sizeof head); mm=  0;
                cin>>n;
                set<P> se;
                for (int i=0;i<n;i++) {
                        scanf("%d%d%d", &p[i].x, &p[i].y, &p[i].r);
                }
                se.insert(P(inf, inf));
                sort(p, p+n, cmp);
                for (int i=0;i<n;i++) p[i].id= i ;
                memset(vis, 0, sizeof vis);
                for (int i=n-1;i>=0;i--) {
                        P l = P(p[i].x - p[i].r, p[i].y - p[i].r);
                        P r = P(p[i].x + p[i].r, p[i].y + p[i].r);
                        set<P>::iterator e,f ;
                        e = se.lower_bound(l);
                        f = se.lower_bound(r);
                        for (set<P>::iterator it = e; it!=se.end() && it!=f; it++)
                        {
                                P tmp = *it;
                                double d = dis(tmp, p[i]);
                                if (!vis[it->id] && d < abs(p[i].r - tmp.r))
                                {
                                        add(i , tmp.id);
                                        vis[tmp.id] = 1;
                                }
                        }
                        se.insert(p[i]);
                }
                for (int i=0;i<n;i++) if (!vis[i]) add(Root, i);
                int ans = dfs(Root);
                puts(ans ? "Alice" : "Bob");
        }
}