题目链接

这道题是一道判断无根树同构的模板题,判断同构主要的思路就是哈希。

一遇到哈希题,一百个人能有一百零一种哈希方式,这篇题解随便选用了一种——类似杨弋《Hash在信息学竞赛中的一类应用》中的这种,可能不是最简洁好写的,但是能用。

我的哈希规则:子树\(u\)的哈希值由它的每一个子树\(v_i\)的哈希值得来,首先将所有\(f(v)\)排个序(防止顺序不同造成影响),然后\(f(u) = size(u) * \sum_i f(v_i)W^{i - 1} \bmod P\),\(W\)是事先选取的一个位权,\(P\)是模数,\(size(u)\)是子树\(u\)的大小。

这样DFS一遍可求出以\(1\)号节点为根时,所有子树的哈希值\(f(u)\)。

但是这是无根树,我们想求出以任意节点为根时整棵树的哈希值。

设\(fa_u\)为以\(1\)为根时\(u\)的父亲,则上面的\(f(u)\)也是以\(fa_u\)为根时子树\(u\)的哈希值。

再求一个\(g(u)\)表示以\(u\)为根时子树\(fa_u\)的哈希值。这个\(g(u)\)怎么求呢?再DFS一遍,对于每个节点,\(g(u)\)由\(g(fa_u)\)以及\(u\)的每个兄弟\(v_i\)的\(f(v_i)\)得来。但是直接暴力枚举的话在菊花图上是\(O(n^2)\)的,那怎么办呢?

对于每个节点\(u\)维护一个数组,存储它所有儿子的哈希值\(f(v)\),如果有父亲,则\(g(u)\)也在里面,把这个数组排好序,求出每个前缀的哈希值和每个后缀的哈希值。这时,以\(u\)为根时整棵树的哈希值就是整个数组的哈希值(再乘上子树大小\(n\))。

此时求每个儿子\(v\)的\(g(v)\),就是从那个数组中间去掉\(f(v)\)后的哈希值,二分查找后把前缀哈希值和后缀哈希值拼起来就可以得到。记得乘上\(v\)为根时\(u\)的\(size\)即\(n - size(v)\)。

这样就求出以每个节点为根的哈希值了。

把A的所有哈希值存到一个set里,然后枚举B的每个度为1的点\(u\),求出以\(u\)为根它的唯一子树\(v\)的哈希值,如果set里有这个值,\(u\)就是所求的点之一。

代码比较丑,见谅 ><

#include <cstdio>
#include <cmath>
#include <cstring>
#include <algorithm>
#include <iostream>
#include <vector>
#include <set>
#define space putchar(' ')
#define enter putchar('\n')
typedef long long ll;
using namespace std;
template <class T>
void read(T &x){
char c;
bool op = 0;
while(c = getchar(), c < '0' || c > '9')
if(c == '-') op = 1;
x = c - '0';
while(c = getchar(), c >= '0' && c <= '9')
x = x * 10 + c - '0';
if(op) x = -x;
}
template <class T>
void write(T x){
if(x < 0) putchar('-'), x = -x;
if(x >= 10) write(x / 10);
putchar('0' + x % 10);
} const int N = 100005, W = 1000000021, P = 999999137;
int n, m, fa[N], f[N], g[N], pw[N], Sze[N], deg[N], ans = P;
int ecnt, adj[N], nxt[2*N], go[2*N];
vector <int> son[N], sl[N], sr[N];
set <int> vis;
bool isB; void add(int u, int v){
if(isB) deg[u]++;
go[++ecnt] = v;
nxt[ecnt] = adj[u];
adj[u] = ecnt;
}
int dfs1(int u, int pre){
Sze[u] = 1;
fa[u] = pre;
son[u].clear();
for(int e = adj[u], v; e; e = nxt[e])
if((v = go[e]) != pre){
son[u].push_back(dfs1(v, u));
Sze[u] += Sze[v];
}
if(son[u].empty()) return f[u] = 1;
sort(son[u].begin(), son[u].end());
ll ret = 0;
for(int i = 0; i < (int)son[u].size(); i++)
ret = (ret * W + son[u][i]) % P;
return f[u] = Sze[u] * ret % P;
} void dfs2(int u){
if(fa[u]){
son[u].push_back(g[u]);
sort(son[u].begin(), son[u].end());
}
int sze = son[u].size();
sl[u].resize(sze);
sl[u][0] = son[u][0];
for(int i = 1; i < sze; i++)
sl[u][i] = ((ll)sl[u][i - 1] * W + son[u][i]) % P;
sr[u].resize(sze);
sr[u][sze - 1] = son[u][sze - 1];
for(int i = sze - 2; i >= 0; i--)
sr[u][i] = (sr[u][i + 1] + (ll)son[u][i] * pw[sze - i - 1]) % P;
for(int e = adj[u], v; e; e = nxt[e])
if((v = go[e]) != fa[u]){
if(sze == 1){
g[v] = 1;
dfs2(v);
break;
}
int p = lower_bound(son[u].begin(), son[u].end(), f[v]) - son[u].begin();
g[v] = 0;
if(p + 1 < sze) g[v] = sr[u][p + 1];
if(p - 1 >= 0) g[v] = (g[v] + (ll)sl[u][p - 1] * pw[sze - 1 - p]) % P;
g[v] = (ll)g[v] * (n - Sze[v]) % P;
if(isB && deg[v] == 1 && vis.find(g[v]) != vis.end()) ans = min(ans, v);
dfs2(v);
}
if(!isB) vis.insert((ll)sl[u][sze - 1] * n % P);
} int main(){ pw[0] = 1;
for(int i = 1; i < N; i++)
pw[i] = (ll)pw[i - 1] * W % P;
read(n);
for(int i = 1, u, v; i < n; i++)
read(u), read(v), add(u, v), add(v, u);
dfs1(1, 0);
dfs2(1);
ecnt = 0, isB = 1, n++;
memset(adj, 0, sizeof(adj));
for(int i = 1, u, v; i < n; i++)
read(u), read(v), add(u, v), add(v, u);
dfs1(1, 0);
if(deg[1] == 1 && vis.find(f[go[adj[1]]]) != vis.end())
ans = 1;
dfs2(1);
write(ans), enter; return 0;
}

BZOJ 4754 [JSOI2016]独特的树叶 | 树哈希判同构的更多相关文章

  1. bzoj 4754: [Jsoi2016]独特的树叶

    不得不说这是神题. %%%   http://blog.csdn.net/samjia2000/article/details/51762811 #include <cstdio> #in ...

  2. bzoj4337: BJOI2015 树的同构 树哈希判同构

    题目链接 bzoj4337: BJOI2015 树的同构 题解 树哈希的一种方法 对于每各节点的哈希值为hash[x] = hash[sonk[x]] * p[k]; p为素数表 代码 #includ ...

  3. BZOJ4754 JSOI2016独特的树叶(哈希)

    判断两棵无根树是否同构只需要把重心提作根哈希即可.由于只添加了一个叶子,重心的位置几乎不发生偏移,所以直接把两棵树的重心提起来,逐层找哈希值不同且对应的两子树即可.被一个普及组子问题卡一年. #inc ...

  4. 【BZOJ4754】独特的树叶(哈希)

    [BZOJ4754]独特的树叶(哈希) 题面 BZOJ 给定一个\(n\)个节点的树A和一个\(n+1\)个节点的树\(B\) 求\(B\)的一个编号最小的节点,使得删去这个节点后\(A,B\)同构 ...

  5. P4323 [JSOI2016]独特的树叶(树哈希)

    传送门 树哈希?->这里 反正大概就是乱搞--的吧-- //minamoto #include<bits/stdc++.h> #define R register #define l ...

  6. BZOJ4754 & 洛谷4323 & LOJ2072:[JSOI2016]独特的树叶——题解

    https://www.lydsy.com/JudgeOnline/problem.php?id=4754 https://www.luogu.org/problemnew/show/P4323 ht ...

  7. bzoj4754[JSOI2016]独特的树叶

    这个题....别人写得怎么都....那么短啊? 我怎么....WA了好几次啊....怎么去loj扒了数据才调出来啊? 这个算法...怎么我还是不知道对不对啊 怎么回事啊怎么回事啊怎么回事啊? 请无视上 ...

  8. Luogu 4323 [JSOI2016]独特的树叶

    新技能get 树哈希,考虑到两棵树相同的条件,把每一个结点的哈希值和树的siz写进哈希值里去. 做出A树每一个结点为根时的树的哈希值丢进set中,然后暴力枚举B树中度数为1的点,求出删掉这个点之后的哈 ...

  9. [JSOI2016]独特的树叶

    https://zybuluo.com/ysner/note/1177340 题面 有一颗大小为\(n\)的树\(A\),现加上一个节点并打乱编号,形成树\(B\),询问加上的节点最后编号是多少? \ ...

随机推荐

  1. linux awk 内置函数实例

    awk内置函数,主要分4种:算数函数.字符串函数.时间函数.一般函数 一.算术函数 以下算术函数执行与 C 语言中名称相同的子例程相同的操作: 函数名 说明 atan2( y, x ) 返回 y/x ...

  2. Spring+SpringMVC+MyBatis整合(easyUI、AdminLte3)

    实战篇(付费教程) 花了几天的时间,做了一个网站小 Demo,最终效果也与此网站类似.以下是这次实战项目的 Demo 演示. 登录页: 富文本编辑页: 图片上传: 退出登录: SSM 搭建精美实用的管 ...

  3. Docker网络解决方案 - Calico部署记录

    简单来说,实现docker跨主机容器间通信,常用的第三方网络方案是Flannel,Weave,Calico:Flannel会为每个host分配一个subnet,容器从这个subnet中分配ip,这些i ...

  4. tableView优化思路

    一般优化的思路: 提前计算并缓存好高度(布局),因为heightForRowAtIndexPath:是调用最频繁的方法. 复杂界面可采用异步绘制. 在大量图片展示时,可以滑动时按需加载. 尽量少用或不 ...

  5. v-for v-if || v-else

    <el-col> <div v-for="item in resultDetail" class="physical-content" v-i ...

  6. 701 D. As Fast As Possible

    ###链接 [http://codeforces.com/group/1EzrFFyOc0/contest/701/problem/D] ###题意 n个人,走lm,有车但没人只能坐一次,车容量为k, ...

  7. Echarts学习求教

    有没有人用过百度的Echarts?刚开始接触,下面这段代码怎么理解啊,新手求指教: myChart.showLoading();$.get('data/asset/data/les-miserable ...

  8. SoftwareEngineering Individual Project - Word frequency program

    说实话前面c#实在没怎么学过.这次写起来感觉非常陌生,就连怎么引用名空间都忘记了.在经过恶补后还是慢慢地适应了. 1.项目预计用时: 构建并写出大概的数据结构,程序框架及模块: 30min 实现文件夹 ...

  9. 毕业设计 之 三 mooodle及bigbluebutton使用笔记(未完成)

    毕业设计 之 三 mooodle及bigbluebutton使用笔记 作者:20135216 平台:windows10 备注:N把辛酸泪 附:在准备过程中其他的一些零碎小问题 一.关于moodle 1 ...

  10. enumerate()函数用法

    enumerate 函数用于遍历序列中的元素以及它们的下标: