首页 > 其他分享 >[HDU6647] Bracket Sequences on Tree 题解

[HDU6647] Bracket Sequences on Tree 题解

时间:2024-03-07 22:03:24浏览次数:27  
标签:Hs int 题解 Tree Bracket mp res deg mod

[HDU6647] Bracket Sequences on Tree 题解

一道纯靠自己推出来的换根\(dp+\)树哈希,写篇题解庆祝一下~~

题意:给定一棵无根树,你可以任意选择根节点和遍历顺序,每次遍历时进入一个节点就标记一个(,离开一个节点就标记一个),问所有存在的括号序列有多少种,对998244353取模。

先考虑根固定的情况如何统计答案。

容易发现两个序列形态相同当且仅当两棵树形态相同。这个玩意可以用树\(Hash\)维护。根据乘法原理,容易知道对于一个点\(x\),以及它的儿子 \(y\) ,其贡献 \(res_x=\prod_{}res[y]\).由于遍历顺序不固定,因此实际对答案的贡献要乘上 \((deg_x-1)!\) .但是考虑有重复子树的情况。根据全排列知识,我们计每一个 \(y\) 形态的出现次数是 \(num_i\),则 \(res_x=\prod_{}res[y]\times \frac{(deg_x-1)!}{\sum num_i!}\).

那么现在考虑换根。在第二次遍历中,根据换根 \(dp\) 的一般形式,我们令当前在 \(x\) 节点,我们要更新其儿子节点 \(y\) 的贡献。具体地,我们考虑记录每个点 \(x\) 成为根节点时它的所有子节点形态的个数 \(mp[x][y]\),按照套路先算出以 \(x\) 为根节点时除去 \(y\) 的子树部分的贡献 \(kes\),再用它更新 \(y\) 点的最终答案。

根据上面的基础式子,我们容易推出:

\[kes=\frac{res_x\times mp[x][Hs_y]}{res_y\times deg_x} \]

其中 \(Hs_y\) 是以1为根时的 \(y\) 子树 \(Hash\) 值

然后我们计算出根据题单里前几道题的套路以 \(x\) 为根节点时除去 \(y\) 的子树部分的\(Hash\) 值 \(ks\),并更新 \(Hs_y\) 为 \(y\) 为根时的 \(Hash\) 值。最后将 \(ks\) 加入 \(y\) 的桶里,并计算出新的 \(res_y\).

\[res_y=\frac{res_y\times kes\times deg_y}{mp[y][ks]} \]

详情见代码。

#include <bits/stdc++.h>
#define mod 998244353
#define N 100005
#define ull unsigned long long
#define int long long
using namespace std;
int fac[N];
int T;
int n;
struct Node {
	int to, nxt;
}e[N << 1];
int head[N], cnt;
void add(int u, int v) {
	e[++cnt].to = v;
	e[cnt].nxt = head[u];
	head[u] = cnt;
}
int qpow(int x, int y) {
	int ans = 1;
	while (y) {
		if (y & 1)
			ans = ans * x % mod;
		x = x * x % mod;
		y >>= 1;
	}
	return ans;
}
mt19937 XYCYX(time(0));
const ull mask = XYCYX();
ull shift(ull x) {
	x ^= mask;
	x ^= x << 13;
	x ^= x >> 7;
	x ^= x << 17;
	x ^= mask;
	return x;
} 
ull Hs[N];
int res[N];
int inv[N];
int INV[N];
int deg[N];
unordered_map<ull, int>mp[N];
void dfs1(int x, int fa) {
	Hs[x] = res[x] = 1;
	for (int i = head[x]; i; i = e[i].nxt) {
		int y = e[i].to;
		if (y == fa)
			continue;
		dfs1(y, x);
		Hs[x] += shift(Hs[y]);
		mp[x][Hs[y]]++;
		res[x] = res[x] * res[y] % mod;
	}
	if (x != 1)
		res[x] = res[x] * fac[deg[x] - 1] % mod;
	else
		res[x] = res[x] * fac[deg[x]] % mod;
	for (auto it : mp[x])
		res[x] = res[x] * inv[it.second] % mod;
}
void dfs2(int x, int fa) {
	for (int i = head[x]; i; i = e[i].nxt) {
		int y = e[i].to;
		if (y == fa)
			continue;
		int kes = res[x] * mp[x][Hs[y]] % mod * qpow(res[y], mod - 2) % mod * qpow(deg[x], mod - 2) % mod;
		ull ks = Hs[x] - shift(Hs[y]);
		mp[y][ks]++;
		Hs[y] = Hs[y] + shift(Hs[x] - shift(Hs[y]));
		res[y] = res[y] * kes % mod * deg[y] % mod * INV[mp[y][ks]] % mod;
		dfs2(y, x);
	}
}
unordered_map<ull, int>MP;
int u[N], v[N];
signed main() {
	ios::sync_with_stdio(0);
	cin.tie(0);
	cin >> T;
	fac[0] = inv[0] = 1;
	for (int i = 1; i < N; i++) {
		fac[i] = fac[i - 1] * i % mod;
		inv[i] = qpow(fac[i], mod - 2);
		INV[i] = qpow(i, mod - 2);
	}
	while (T--) {
		cnt = 0;
		for (int i = 0; i <= n; i++) {
			e[i].to = e[i].nxt = e[i * 2].to = e[i * 2].nxt = 0;
			head[i] = head[i * 2] = 0;
			Hs[i] = 0;
			res[i] = 0;
			deg[i] = 0;
			mp[i].clear();
		}
		MP.clear();
		cin >> n;
		for (int i = 1; i < n; i++) {
			int x, y;
			cin >> x >> y;
			u[i] = x, v[i] = y;
			add(x, y);
			add(y, x);
			deg[x]++;
			deg[y]++;
		}
		dfs1(1, 0);
		dfs2(1, 0);
		int ans = 0;
		for (int i = 1; i <= n; i++)
			MP[Hs[i]] = res[i];
		for (auto i : MP)
			ans = (ans + i.second) % mod;
		cout << ans << "\n";
	}
	return 0;

本题的难点:

  1. 灵活运用排列组合相关知识,计算根固定时的答案贡献。

  2. 熟练掌握换根 \(dp\) 的 $\ up-down\ $思想推出转移式并进行一定程度上的转化。

  3. 各种繁琐细节的处理!

\(End.\)

标签:Hs,int,题解,Tree,Bracket,mp,res,deg,mod
From: https://www.cnblogs.com/Rock-N-Roll/p/18059865

相关文章

  • CF contest 1935 Round 932 (Div. 2) A-D题解
    CodeforcesRound932(Div.2)A-D题解CodeforcesRound932(Div.2)绪言很菜,AB速度慢,卡在C,想DP,但是时间优化不下来,说服自己\(2\times10^3\)能过\(n^3\),D稍微简单,但是没看D,遂掉分。A.EntertainmentinMAC给你一个字符串\(s\)和一个偶整数\(n\)。你可以对它进行两种运......
  • 【学习笔记】 - 基础数据结构 :Link-Cut Tree(进阶篇)
    前言LCT没题写可以去写树剖和一些线段树合并的题练手LCT的概念原本的树剖是对树进行剖分,剖分为重边和轻边LCT则是对于树分为虚边和实边,特殊的,LCT可以没有虚边(例:银河英雄传说v2)单独被包含在一个实链里的点称作孤立点在树剖中,我们使用线段树/树状数组来维护重链在Link-Cut......
  • 【教程】 iOS构建版本无效问题解决方案
     引言在进行iOS应用上架时,有时会遇到构建版本无效的问题,即通过XCode上传成功后,但在AppStoreConnect的TestFlight中无法显示构建版本,或者显示一会儿后就消失了。本文将介绍可能的原因分析,并提供解决问题的方法。问题描述最近一次上传新版本至AppStore后,发现在AppStoreCon......
  • 说说Vue 3.0中Treeshaking特性?举例说明一下?
    这里给大家分享我在网上总结出来的一些知识,希望对大家有所帮助 一、是什么Treeshaking 是一种通过清除多余代码方式来优化项目打包体积的技术,专业术语叫 Deadcodeelimination简单来讲,就是在保持代码运行结果不变的前提下,去除无用的代码如果把代码打包比作制作蛋糕,传统......
  • [ARC157F] XY Ladder LCS 题解
    我们尝试给这个抽象题来一篇题解。思考过程还是很重要的。首先看了这个题,一看数据范围\(n\le50\),然后就不懂了,你告诉我这玩意可以状压??然后我们一顿乱想,发现如果\(n\)除以一个\(3\),那我们是不是就可以状压了。那怎么除以\(3\)呢。接着我们手玩一下样例,发现似乎这个答案......
  • [青少年CTF训练平台]web部分题解(已完结!)
    文章管理系统首先打开环境(>ω<。人)ZZz♪♪既然要做题,就要做全面了,图上说了,既然有假flag我就先找出来:假flag:打开vmware,使用sqlmap进行处理:sqlmap-uhttp://challenge.qsnctf.com:31645/?id=1--dbs记得中间的url换成自己的看到了六个可能:{*]ctftraining[*]information......
  • 向TreeView添加自定义信息
    可在Windows窗体TreeView控件中创建派生节点,或在ListView控件中创建派生项。通过派生可添加任何所需字段,以及添加处理这些字段的自定义方法和构造函数。此功能的用途之一是将Customer对象附加到每个树节点或列表项。虽然此处的示例是关于TreeView控件的,但该方法同样......
  • CF1353E K-periodic Garland 题解
    分析考虑DP。定义状态函数\(f_i\)表示处理完前\(i\)个字符且第\(i\)个字符为\(1\)时的最小代价。则对于\(i\),有两种情况:\(i\)不是第一个\(1\),则上一个\(1\)的位置必定为\(i-k\)。\(i\)是第一个\(1\),没有上一个\(1\)。得到转移方程:\(f_i=\min(f_{\max(......
  • P10120『STA - R4』冰红茶 题解
    分析出得很好,模板套模板,希望下次再来。难点在于维护最后连续喝的DS饮料数量。设这次喝原味饮料的区间为\([l,r]\),上一次为\([l',r']\)。则有两种情况:\([l,r]\)与\([l',r']\)不相交。如:在\([l',r']\)和\([l,r]\)两个区间中的DS连续喝的同种饮料数量都会变成\(k......
  • AT_abc343_f [ABC343F] Second Largest Query 题解
    分析考虑乱搞。对于求次大值,用线段树维护就行了。记录下每个区间的最大、次大值。则两个子区间的父区间的最大值就是这四个最大的,次大值就是这四个次大的。复杂度\(O(\logn)\)。求次大值的出现次数,乱搞就行了。因为带修,带修莫队或者分块有些麻烦。其实用线段树就行。在维护区......