【HDU6647】Bracket Sequences on Tree(树Hash 树上Dp)
大意
给出一颗树,按下列方式生成一个括号序列。
function dfs(int cur, int parent):
print('(')
for all nxt that cur is adjacent to:
dfs(nxt, cur)
print(')')
其中可以从任一点出发,且对儿子的遍历顺序是随机的。
求本质不同的括号序列个数。
思路
前置板块:树Hash
如何判断两颗有根树是否本质一样?
我们先随机生成一个\(T\)数组(随机数被卡概率小?)
令\(Siz[u]\)表示\(u\)的子树大小,\(H[u]\)表示以\(u\)为根的有根子树的Hash值。
那么\(H[u]\)的定义式就为:$$H[u]=C+T[Siz[u]]\times H[v]~~(v\in Son[u])$$其中\(C\)为随机构造的一个数。(主要是怕被卡)
若对于某对\(i,j\),有\(H[i]=H[j]\),则说明以\(i\)为根的有根树的形态与以\(j\)为根是一样的。
那么如何判断无根树呢?
我们采用换根Dp来做,考虑\(H\)数组的定义式。
那么向下转移时就有:$$H[v]=H[v]+T[Siz[u]-Siz[v]]\times(H[u]-H[v]\times T[Siz[v]])~~(v\in Son[u])$$这样,我们就可以得到以每个点为根的有根树的Hash值了。
考虑整棵无根树,我们取以每个点为根的有根树的最大Hash值就行了,即\(Max(H[i])\).
考虑固定了起点的情况。
那么对于两颗形状不同的有根子树,其所生成的括号序列集是显然没有交集的。
则对于某个点来说,我们需要记录它本质相同的儿子子树的个数。
我们设\(Cnt[u][x]\)表示\(u\)的儿子中Hash值为\(x\)的的个数。(Map好啊)
以\(u\)为根的子树内的答案\(Dp[u]\)就为
对于起点不固定的情况。
其实我们换根Dp就行了,细节有点多,注意一下。
最后我们统计答案时,每种Hash值只统计一次就行了。
# 代码
及其丑陋的代码。
可恶的出题人卡简单Hash。
```cpp
#include<map>
#include<ctime>
#include<cstdio>
#include<vector>
#include<algorithm>
using namespace std;
#define ULL unsigned long long
const int C=10007;
const int MAXT=100000;
const int MAXN=100005;
const int MOD=998244353;
const ULL LW=2305843009213693951ll;
ULL Ans;
ULL T[MAXN],H[MAXN];
int K,N,Siz[MAXN],Son[MAXN];
ULL Dp[MAXN];
ULL E[MAXN],F[MAXN];
vector<int>P[MAXN];
vector<ULL>Kind[MAXN];
map<ULL,int>Mp;
map<ULL,int>Cnt[MAXN];
ULL Mul(ULL X,ULL Y){
return (X*Y-(ULL)(X/(long double)LW*Y+1e-3)*LW+LW)%LW;
}
ULL quick_Pow(ULL x,ULL y){
if(y==0)return 1;
if(y==1)return x;
if(y%2)return (x*quick_Pow((x*x)%MOD,y/2))%MOD;
return quick_Pow((x*x)%MOD,y/2);
}
void Prepare(){
F[0]=E[0]=1;
for(int i=1;i<=MAXT;i++)
F[i]=(F[i-1]*i)%MOD;
E[MAXT]=quick_Pow(F[MAXT],MOD-2);
for(int i=MAXT-1;i>=1;i--)
E[i]=(E[i+1]*(i+1))%MOD;
}
void DFS(int u,int fa){
Siz[u]=1;H[u]=C;Dp[u]=1;
int size=P[u].size();
for(int i=0;i<size;i++){
int v=P[u][i];
if(v==fa)continue;
DFS(v,u);Son[u]++;
Dp[u]=(Dp[u]*Dp[v])%MOD;
Siz[u]+=Siz[v];
H[u]=(H[u]+Mul(H[v],T[Siz[v]]))%LW;
if(!Cnt[u][H[v]])Kind[u].push_back(H[v]);
Cnt[u][H[v]]++;
}
Dp[u]=(Dp[u]*F[Son[u]])%MOD;
size=Kind[u].size();
for(int i=0;i<size;i++){
ULL v=Kind[u][i];
Dp[u]=(Dp[u]*E[Cnt[u][v]])%MOD;
}
}
void DFS2(int u,int fa){
int size=P[u].size();
for(int i=0;i<size;i++){
int v=P[u][i];
if(v==fa)continue;
ULL Hu=(H[u]-Mul(T[Siz[v]],H[v])+LW)%LW;
ULL Dpu=Dp[u]*E[Son[u]]%MOD*F[Son[u]-1]%MOD*F[Cnt[u][H[v]]]%MOD*E[Cnt[u][H[v]]-1]%MOD*quick_Pow(Dp[v],MOD-2)%MOD;
H[v]=(H[v]+Mul(T[Siz[u]-Siz[v]],(H[u]-Mul(T[Siz[v]],H[v])+LW)))%LW;
Dp[v]=Dp[v]*E[Son[v]]%MOD*F[Son[v]+1]%MOD*F[Cnt[v][Hu]]%MOD*E[Cnt[v][Hu]+1]%MOD*Dpu%MOD;
Cnt[v][Hu]++;Son[v]++;
Siz[v]=Siz[u];
DFS2(v,u);
}
}
int main(){
srand(time(NULL));
scanf("%d",&K);
for(int i=1;i<=MAXT;i++)T[i]=(rand()*rand()*rand())%LW+1;
Prepare();
while(K--){
scanf("%d",&N);
Mp.clear();Ans=0;
for(int i=1;i<=N;i++){
H[i]=Siz[i]=0;
Dp[i]=0;Son[i]=0;
P[i].clear();
Kind[i].clear();
Cnt[i].clear();
}
for(int i=1,x,y;i<N;i++){
scanf("%d%d",&x,&y);
P[x].push_back(y);
P[y].push_back(x);
}
DFS(1,0);DFS2(1,0);
for(int i=1;i<=N;i++)
if(!Mp[H[i]]){
Ans=(Ans+Dp[i])%MOD;
Mp[H[i]]=1;
}
printf("%lld\n",Ans);
}
}
```\]
【HDU6647】Bracket Sequences on Tree(树Hash 树上Dp)的更多相关文章
- POJ 1741.Tree 树分治 树形dp 树上点对
Tree Time Limit: 1000MS Memory Limit: 30000K Total Submissions: 24258 Accepted: 8062 Description ...
- Bzoj3197/洛谷3296 [SDOI2013]刺客信条assassin(树的重心+树Hash+树形DP+KM)
题面 Bzoj 洛谷 题解 (除了代码均摘自喻队的博客,可是他退役了) 首先固定一棵树,枚举另一棵树,显然另一棵树只有与这棵树同构才有可能产生贡献 如果固定的树以重心为根,那么另一棵树最多就只有重心为 ...
- Codeforces Round #646 (Div. 2) E. Tree Shuffling(树上dp)
题目链接:https://codeforces.com/contest/1363/problem/E 题意 有一棵 $n$ 个结点,根为结点 $1$ 的树,每个结点有一个选取代价 $a_i$,当前 $ ...
- Codeforces Round #350 (Div. 2) E. Correct Bracket Sequence Editor 线段树模拟
E. Correct Bracket Sequence Editor Recently Polycarp started to develop a text editor that works o ...
- 树hash
判断树的同构,采用树hash的方式. 树hash定义在有根树上.判断无根树同构的时候,可以比较重心为根的hash值或者比较每个点为根的hash值. h[x]表示x为根的子树的hash,g[x]表示x为 ...
- 【题解】彩色树 51nod 1868 虚树 树上dp
Prelude 题目在这里:ο(=•ω<=)ρ⌒☆ Solution 蒟蒻__stdcall的第一道虚树题qaq. 首先很容易发现,这个排列是假的. 我们只需要求出每对点之间的颜色数量,然后求个 ...
- Codeforces E. Alyona and a tree(二分树上差分)
题目描述: Alyona and a tree time limit per test 2 seconds memory limit per test 256 megabytes input stan ...
- TZOJ 4292 Count the Trees(树hash)
描述 A binary tree is a tree data structure in which each node has at most two child nodes, usually di ...
- Codevs 2370 小机房的树 LCA 树上倍增
题目描述 Description 小机房有棵焕狗种的树,树上有N个节点,节点标号为0到N-1,有两只虫子名叫飘狗和大吉狗,分居在两个不同的节点上.有一天,他们想爬到一个节点上去搞基,但是作为两只虫子, ...
随机推荐
- awk 内置变量
awk 内置变量 FS 设置输入域分隔符,等价于命令行 -F选项 OFS 输出域分隔符 NF 浏览记录的域的个数NR 已读的记录数RS 控制记录分隔符 ORS 输出记录分隔符 substr 截取指定的 ...
- Sentry 企业级数据安全解决方案 - Relay 配置选项
Relay 的配置记录在文件 .relay/config.yml 中.要更改此位置,请将 --config 选项传递给任何 Relay 命令: ❯ ./relay run --config /path ...
- Python_string.Template的使用
Template是python string提供的一个字符串模板功能.主要用于文本处理 from string import Template s = Template('$who 在 $do') t ...
- 阿里云服务器ECS Ubuntu16.04 + Seafile 搭建私人网盘 (Seafile Pro)
原文链接:? 传送门 本文主要讲述 使用 Ubuntu 16.04 云服务器 通过脚本实现对 Seafile Pro 的安装,完成私人网盘的搭建 首先给出 Seafile 专业版的下载地址(Linux ...
- 如何让 Hexo 在服务器稳定运行
声明 本文地址:如何让 Hexo 在服务器稳定运行 背景 博客系统终于又搭建起来了(好一个又),但是每隔一段时间去访问自己的网站总是访问不到,去服务器查询 ps aux | grep hexo,发现 ...
- vivo推送平台架构演进
本文根据Li Qingxin老师在"2021 vivo开发者大会"现场演讲内容整理而成.公众号回复[2021VDC]获取互联网技术分会场议题相关资料. 一.vivo推送平台介绍 1 ...
- Collection的使用
常用方法 Modifier and Type Method and Description boolean add(E e) 确保此集合包含指定的元素(可选操作). boolean addAll(Co ...
- 很详细的FFT(快速傅里叶变换)概念与实现
FFT 首先要说明一个误区,很多人认为FFT只是用来处理多项式乘的,其实FFT是用来实现多项式的系数表示法和点值表示法的快速转换的,所以FFT的用处远不止多项式乘. FFT的前置知识:点值表示法,复数 ...
- 黑客是如何通过开放的Redis服务入侵服务器的
0x00 简要说明 百度百科:Redis(Remote Dictionary Server ),即远程字典服务,是一个开源的使用ANSI C语言编写.支持网络.可基于内存亦可持久化的日志型.Key-V ...
- linux信号 SIGINT SIGTERM SIGKILL
三者都是结束/终止进程运行. 1.SIGINT SIGTERM区别 前者与字符ctrl+c关联,后者没有任何控制字符关联. 前者只能结束前台进程,后者则不是. 2.SIGTERM SIGKILL的区别 ...