hdu 5486 Difference of Clustering 2015合肥网络赛 并查集 离散化 悲伤的题

题目

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=5486

题目来源：2015合肥网络赛第五题。出线的关键题，也是很让我悲伤的题。

简要题意：有合并，分裂，一对一三种操作，求一个结果到另一个结果每种操作的次数。

数据范围：$T\leqslant 100;\quad 0\leqslant N\leqslant 10^6;\quad |a_i|,|b_i|< 2^{31}$

题记·悲伤

这是一道很令我悲伤的题。

我们队两个多小时出完了四题，然后死磕1003。

然后一直MLE，我很费解，在最后二十分钟的时候我都要放弃治疗了。

然后跑去上厕所，上厕所的时候看到了算法分析的老师，尿到一半突然发现这题好像可以用并查集。

然后越想就越觉得就是这么搞的。

回到实验室说了一下，然后开始噼里啪啦敲。

然后因为着急然后有点混乱，最终改出来了，然后一交。艹，正好压线。

然后就特么WA了，这其实也还好，然后我想稍微调一调，发现并查集初始化的时候下标好像不对，我一般并查集是1-base的，今天正好搞成0-base了，然后改了下一交，特么就过了。

要不要这么悲伤啊，结果出线的名额就这么没了，又被苏大给压了。

更悲伤的是我之前还特意检查了一下，然后看到了这个，还特意改成了0-base，然后想想好像无所谓的又改了回去，一失足成千古恨啊。

题解

设两组盒子分别为$A[N],B[M]$

相同标号的盒子之间存在关联关系，一个组内的存在关联关系，这个关系可以用并查集维护。

合并相同序号的小球所在的两个盒子。

然后对于每个盒子对其祖先去更新含有的$A,B$盒子的数量。

最后通过祖先含有的$A,B$盒子的数量来判断操作的类型。

实现

$a_i,b_i$的范围比较大，要离散化。

其实我觉得这题不是太难，但是题意比较扭曲复杂，所以一开始也没往其他想，就模拟了。

而且这场的出题人好喜欢离散化啊。

代码

#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <cstring>
#include <stack>
#include <queue>
#include <string>
#include <vector>
#include <set>
#include <map>

#define pb push_back
#define mp make_pair
#define all(x) (x).begin(),(x).end()
#define sz(x) ((int)(x).size())
#define fi first
#define se second
using namespace std;
typedef long long LL;
typedef vector<int> VI;
typedef pair<int,int> PII;
LL powmod(LL a,LL b, LL MOD) {LL res=1;a%=MOD;for(;b;b>>=1){if(b&1)res=res*a%MOD;a=a*a%MOD;}return res;}
// head
const int N = 1E6+5;
int a[N];
int b[N];
int disa[N];
int disb[N];
int res[N][2];

// 下标1-n 路径压缩 秩优化
int r[N*2+5];
int fa[N*2+5];

int find(int key) {
    return (key == fa[key]) ? key : (fa[key] = find(fa[key]));
}

void init(int n) {
    //悲伤
    for (int i = 0; i <= n; i++) {
        r[i] = 0;
        fa[i] = i;
    }
}

void joint(int u, int v) {
    int a = find(u), b = find(v);
    if (a != b) {
        if (r[a] > r[b]) swap(a, b);
        fa[a] = b;
        if (r[a] == r[b]) r[b]++;
    }
}

bool same(int u, int v) {
    return find(u) == find(v);
}

int main() {
    int t, n, cas = 1;
    scanf("%d", &t);
    while (t--) {
        scanf("%d", &n);
        int lena, lenb, mer = 0, one = 0, spi = 0;
        for (int i = 0; i < n; i++) {
            scanf("%d%d", a+i, b+i);
            disa[i] = a[i], disb[i] = b[i];
        }
        sort(disa, disa+n);
        sort(disb, disb+n);
        lena = unique(disa, disa+n)-disa;
        lenb = unique(disb, disb+n)-disb;
        for (int i = 0; i < n; i++) {
            a[i] = lower_bound(disa, disa+lena, a[i])-disa;
            b[i] = lower_bound(disb, disb+lenb, b[i])-disb+lena;
        }

        init(lena+lenb);
        for (int i = 0; i < n; i++) {
            joint(a[i], b[i]);
        }
        for (int i = 0; i < lena; i++) {
            res[find(i)][0]++;
        }
        for (int i = 0; i < lenb; i++) {
            res[find(i+lena)][1]++;
        }
        for (int i = 0; i < lena+lenb; i++) {
            int cnta = res[i][0], cntb = res[i][1];
            if (cnta == 1 && cntb == 1) {
                one++;
            } else if (cnta == 1 && cntb > 1) {
                spi++;
            } else if (cnta > 1 && cntb == 1) {
                mer++;
            }
        }

        printf("Case #%d: %d %d %d\n", cas++, spi, mer, one);
        memset(res, 0, sizeof res);
    }
    return 0;
}