题目
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=5486
题目来源:2015合肥网络赛第五题。出线的关键题,也是很让我悲伤的题。
简要题意:有合并,分裂,一对一三种操作,求一个结果到另一个结果每种操作的次数。
数据范围: T⩽100;0⩽N⩽106;|ai|,|bi|<231
题记·悲伤
这是一道很令我悲伤的题。
我们队两个多小时出完了四题,然后死磕1003。
然后一直MLE,我很费解,在最后二十分钟的时候我都要放弃治疗了。
然后跑去上厕所,上厕所的时候看到了算法分析的老师,尿到一半突然发现这题好像可以用并查集。
然后越想就越觉得就是这么搞的。
回到实验室说了一下,然后开始噼里啪啦敲。
然后因为着急然后有点混乱,最终改出来了,然后一交。艹,正好压线。
然后就特么WA了,这其实也还好,然后我想稍微调一调,发现并查集初始化的时候下标好像不对,我一般并查集是1-base的,今天正好搞成0-base了,然后改了下一交,特么就过了。
要不要这么悲伤啊,结果出线的名额就这么没了,又被苏大给压了。
更悲伤的是我之前还特意检查了一下,然后看到了这个,还特意改成了0-base,然后想想好像无所谓的又改了回去,一失足成千古恨啊。
题解
设两组盒子分别为 A[N],B[M]
相同标号的盒子之间存在关联关系,一个组内的存在关联关系,这个关系可以用并查集维护。
合并相同序号的小球所在的两个盒子。
然后对于每个盒子对其祖先去更新含有的 A,B 盒子的数量。
最后通过祖先含有的 A,B 盒子的数量来判断操作的类型。
实现
ai,bi 的范围比较大,要离散化。
其实我觉得这题不是太难,但是题意比较扭曲复杂,所以一开始也没往其他想,就模拟了。
而且这场的出题人好喜欢离散化啊。
代码
#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;
}