[BZOJ4466][Jsoi2013]超立方体
[BZOJ4466][Jsoi2013]超立方体
试题描述
超立方体是立方体在高维空间内的拓展（其在 \(2\) 维情况下退化为正方形，\(1\) 维情况下退化成线段）。在理论计算机科学领域里，超立方体往往可以和 \(2\) 进制编码联系到一起。对理论计算机科学颇有研究的 Will 自然也会对超立方体有着自己的思考。
.png)
上图就是在 \(0 \sim 4\) 维空间内超立方体所对应的图形。显然我们可以把超立方体的每个顶点看成一个点，每一条棱看成一条边，这样就会得到一个无向图，我们称之为超立方图。
\(D\) 维空间内的超立方图有 \(2^D\) 个点，我们把这些点从 \(0\) 到 \(2^D - 1\) 依次编号。
有一个有趣而重要的充要结论是：一定存在一种编号的方式，使得图中任意两个有边相连的顶点的编号的 2进制码中，恰好有一位不同。
在 \(2\) 维和 \(3\) 维空间内这个结论可以这样形象的理解：
对于 \(2\) 维空间，我们只要把这个正方形放到第一象限内，使得 \(4\) 个顶点的坐标按逆时针顺序依次为 \((0,0)\)，\((1,0)\)，\((1,1)\)，\((0,1)\)，然后再把坐标看成 \(2\) 位 \(2\) 进制数，依次将这 \(4\) 个点编号为 \(0\)，\(1\)，\(3\)，\(2\) 即可。
对于 \(3\) 维空间，同样我们可以将立方体的一个顶点与原点重合，并使得所有棱均平行于坐标轴，然后分别确定这 \(8\) 个点的坐标，最后把 \(3\) 维空间内的坐标看成一个 \(3\) 位 \(2\) 进制数即可。对于 \(D\) 维空间，以此类推。
现在对于一个 \(N\) 个点 \(M\) 条边的无向图（每个点从 \(0\) 到 \(N-1\) 编号），Will 希望知道这个图是否同构于一个超立方图。
输入
第一行包含一个整数 \(Q\) 表示此数据中一共包含 \(Q\) 个询问。
接下来 \(Q\) 组询问，每一组询问的输入格式如下：
第一行包含两个整数 \(N\) 和 \(M\)，接下来 \(M\) 行，每行 \(2\) 个不同的整数 \(x\)，\(y\)，表示图中存在一条无向边连接编号为 \(x\) 和 \(y\) 的点（\(0 \le x,y < N\)）
输出
对于每一个询问分别输出一行，内容如下：
1、如果询问中给定的图不同构于任何一个超立方图，输出 \(-1\)；
2、如果同构于某一个超立方图，那么请给图中这 \(N\) 个点重新编号，并在这一行输出 \(N\) 个用空格隔开的整数，表示原图中每个点新的编号，使得重新编号后，满足题目中所述的结论。
注意：输出文件的每一行，要么仅包含一个整数 \(-1\)，要么则应包含一个由 \(0\) 到 \(N-1\) 这 \(N\) 个数组成的排列。如果有多组解输出任意一个均可。
输入示例
3
2 2
0 1
1 0
4 4
0 1
1 2
2 0
0 3
8 12
2 3
2 6
7 6
1 7
4 1
3 4
0 2
7 3
5 6
5 1
5 0
4 0
输出示例
-1
-1
0 6 1 5 4 2 3 7
数据规模及约定
\(Q \le 3,N \le 32768,M \le 1000000\)
题解
这题看上去有点吓人，但其实就是一个模拟，当然需要分析一些性质。
首先 \(N\) 必须是 \(2\) 整数次幂；\(M\) 也要满足一定条件，你可以一维一维地递推出 \(d\) 维下的超立方体有多少条棱，令 \(e_d\) 表示那个答案，则 \(e_0 = 0, e_i = 2e_{i-1} + 2^{i-1}\)，判一下 \(M\) 是否等于 \(e_{\log N}\) 即可。（那个递推式的意思就是，\(d\) 维是由 \(d-1\) 的超立方体平移得到的，这样原来的棱数会倍增，同时每个点的平移轨迹会产生一条新棱，所以要加上 \(d-1\) 维时的点数）
然后随便规定一个点的标号为 \(0\)，将其周围的点的编号依次设为 \(2^0, 2^1, 2^2, \cdots , 2^d\)，然后开始 bfs，当搜到一个未确定标号的点时，若图合法则一定会后两个邻居标号已经确定（因为合法的图一定会有许多的四元环），那么就可以用这两个标号计算出它的标号。最后再检查一下是否符合题意。
中间有许多细节，实现时要注意。
然后还有一个问题就是为什么可以随便选一个点开始确定，这是因为立方体高度对称，它可以任意旋转；在 \(d\) 维的超立方体中，每个点出去一定都会有刚好 \(d\) 条互相垂直的棱，所以谁都能作为“坐标原点”，它连出去的边就是坐标轴了。
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cctype>
#include <algorithm>
#include <cmath>
using namespace std;
#define rep(i, s, t) for(int i = (s), mi = (t); i <= mi; i++)
#define dwn(i, s, t) for(int i = (s), mi = (t); i >= mi; i--)
int read() {
int x = 0, f = 1; char c = getchar();
while(!isdigit(c)){ if(c == '-') f = -1; c = getchar(); }
while(isdigit(c)){ x = x * 10 + c - '0'; c = getchar(); }
return x * f;
}
#define maxn 32790
#define maxm 2000010
int n, m, head[maxn], nxt[maxm], to[maxm], D;
bool check_num(int n, int m) {
int d = 0;
while(!(n & 1)) n >>= 1, d++; D = d;
if(n > 1) return 0;
int e = 0, node = 1;
rep(i, 1, d) e = (e << 1) + node, node <<= 1;
if(m != e) return 0;
return 1;
}
void AddEdge(int a, int b) {
to[m] = b; nxt[m] = head[a]; head[a] = m++;
swap(a, b);
to[m] = b; nxt[m] = head[a]; head[a] = m++;
return ;
}
int Q[maxn], hd, tl, tag[maxn];
bool vis[maxn], used[maxn];
int confirm(int u) {
int v1 = -1, v2 = -1;
for(int e = head[u]; e != -1; e = nxt[e]) if(vis[to[e]]) {
if(v1 < 0) v1 = to[e];
else{ v2 = to[e]; break; }
}
if(v2 < 0) return -1;
int Xor = tag[v1] ^ tag[v2], cdiff = 0;
rep(i, 0, D - 1) if(Xor >> i & 1) cdiff++;
if(cdiff != 2) return -1;
tag[u] = 0;
rep(i, 0, D - 1) if(!(Xor >> i & 1)) tag[u] |= tag[v1] & (1 << i);
if(!used[tag[u]]) return used[tag[u]] = 1, 0;
// attempt 1
tag[u] ^= Xor;
if(!used[tag[u]]) return used[tag[u]] = 1, 0;
tag[u] ^= Xor;
// attempt 2
int oXor = Xor;
rep(i, 0, D - 1) if(Xor >> i & 1){ Xor ^= 1 << i; break; }
tag[u] ^= Xor;
if(!used[tag[u]]) return used[tag[u]] = 1, 0;
tag[u] ^= Xor;
// attempt 3
Xor = oXor;
dwn(i, D - 1, 0) if(Xor >> i & 1){ Xor ^= 1 << i; break; }
tag[u] ^= Xor;
if(!used[tag[u]]) return used[tag[u]] = 1, 0;
// attempt 4
return -1;
}
bool check() {
rep(u, 0, n - 1) {
for(int e = head[u]; e != -1; e = nxt[e]) {
int Xor = tag[u] ^ tag[to[e]], cdiff = 0;
rep(i, 0, D - 1) if(Xor >> i & 1) cdiff++;
if(cdiff != 1) return 0;
}
}
return 1;
}
void bfs() {
hd = tl = 0; Q[++tl] = 0;
tag[0] = 0; vis[0] = used[0] = 1;
int c = 0;
for(int e = head[0]; e != -1; e = nxt[e]) {
if(!vis[to[e]]) vis[to[e]] = used[1<<c] = 1, Q[++tl] = to[e], tag[to[e]] = 1 << c;
c++; if(c > D) return (void) puts("-1");
}
if(c != D) return (void)puts("-1");
hd++;
while(hd < tl) {
int u = Q[++hd];
c = 0;
for(int e = head[u]; e != -1; e = nxt[e]) {
if(!vis[to[e]]) {
vis[to[e]] = 1; Q[++tl] = to[e];
if(confirm(to[e]) == -1) return (void)puts("-1");
}
c++;
}
if(c != D) return (void)puts("-1");
}
rep(i, 0, n - 1) if(!vis[i]) return (void)puts("-1");
if(!check()) return (void)puts("-1");
rep(i, 0, n - 1) printf("%d%c", tag[i], i < n - 1 ? ' ' : '\n');
return ;
}
void work() {
n = read(); int M = read();
m = 0; memset(head, -1, sizeof(head));
rep(i, 1, M) {
int a = read(), b = read();
AddEdge(a, b);
}
if(!check_num(n, M)) return (void)puts("-1");
if(n == 1) return (void)puts("0");
memset(vis, 0, sizeof(vis));
memset(used, 0, sizeof(used));
bfs();
return ;
}
int main() {
int T = read();
while(T--) work();
return 0;
}
