POJ2236 Wireless Network

　　原题传送：http://poj.org/problem?id=2236

　　并查集。本来看到N是10^3级别，操作数P是10^6级别，那么N*P就是O(10^9)级别，一阵后怕，不停寻找Plog(N)算法，后来实在想不出，看看时间限制是10s，直接上O(10^9)碰碰运气，没想到1Y，800多ms。

　　这样的话实现就很容易了，但是1Y后又交了几次，用C++编译器有时ac，有时wa；而用G++虽然时间会慢一点，但是一直都是AC，神马情况？！

#include <stdio.h>
#include <string.h>
#define N 1002

int r[N][2], n, d, f[N];
bool vis[N];

int find_set(int x)
{
    return x == f[x] ? f[x] : f[x] = find_set(f[x]);
}

inline bool inside(int i, int j)
{
    return (r[i][0] - r[j][0]) * (r[i][0] - r[j][0]) + (r[i][1] - r[j][1]) * (r[i][1] - r[j][1]) <= d * d;
}

void work()
{
    int a, b, x, y, i, j;
    char op[3];
    while(scanf("%s", op) != EOF)
    {
        switch(op[0])
        {
            case 'O':
                scanf("%d", &j);
                vis[j] = true;
                for(i = 1; i <= n; i ++)
                {
                    if(vis[i] && inside(i, j))
                    {
                        x = find_set(i);
                        y = find_set(j);
                        x < y ? f[y] = x : f[x] = y;
                    }
                }
                break;
            case 'S':
                scanf("%d%d", &a, &b);
                if(vis[a] && vis[b])
                {
                    x = find_set(a);
                    y = find_set(b);
                    if(x == y)
                        puts("SUCCESS");
                    else 
                        puts("FAIL");
                }
                else
                    puts("FAIL");
                break;
        }
    }
    
}

void read()
{
    scanf("%d%d", &n, &d);
    for(int i = 1; i <= n; i ++)
        scanf("%d%d", &r[i][0], &r[i][1]);
}

void init()
{
    memset(vis, false, sizeof vis);
    for(int i = 0; i <= N; i ++)
        f[i] = i;
}

int main()
{
    init();
    read();
    work();
    return 0;
}