博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
[BJOI 2018]求和
阅读量:4575 次
发布时间:2019-06-08

本文共 2128 字,大约阅读时间需要 7 分钟。

Description

给你一棵 \(n\) 个结点的有根树, \(m\) 次询问这棵树上一段路径上所有节点深度的 \(k\) 次方和。

\(1\leq n\leq 300000,1\leq k\leq 50\)

Solution

树剖之后维护 \(50\) 个前缀和数组。

只有刷水题才能维持得了生活这样子...

Code

#include 
using namespace std;const int N = 300000+5, yzh = 998244353;int gi() { int x = 0; char ch = getchar(); while (ch < '0' || ch > '9') ch = getchar(); while (ch >= '0' && ch <= '9') x = (x<<1)+(x<<3)+ch-48, ch = getchar(); return x;}int n, m, u, v, k, sum[51][N];struct tt {int to, next; }edge[N<<1];int path[N], TP;int size[N], top[N], fa[N], dep[N], id[N], son[N], idx;void dfs1(int u, int depth, int father) { dep[u] = depth, size[u] = 1, fa[u] = father; for (int i = path[u]; i; i = edge[i].next) if (edge[i].to != father) { dfs1(edge[i].to, depth+1, u); size[u] += size[edge[i].to]; if (size[edge[i].to] > size[son[u]]) son[u] = edge[i].to; }}void dfs2(int u, int tp) { top[u] = tp, id[u] = ++idx; if (son[u]) dfs2(son[u], tp); for (int i = path[u]; i; i = edge[i].next) if (edge[i].to != fa[u] && edge[i].to != son[u]) dfs2(edge[i].to, edge[i].to);}void add(int u, int v) {edge[++TP] = (tt){v, path[u]}; path[u] = TP; }int cal(int u, int v, int k) { int ans = 0; while (top[u] != top[v]) { if (dep[top[u]] < dep[top[v]]) swap(u, v); (ans += (sum[k][id[u]]-sum[k][id[top[u]]-1]+yzh)%yzh) %= yzh; u = fa[top[u]]; } if (dep[u] < dep[v]) swap(u, v); (ans += (sum[k][id[u]]-sum[k][id[v]-1]+yzh)%yzh) %= yzh; return ans;}void work() { n = gi(); for (int i = 1; i < n; i++) { u = gi(), v = gi(); add(u, v), add(v, u); } dfs1(1, 0, 0); dfs2(1, 1); for (int i = 1; i <= n; i++) for (int j = 1, val = dep[i]; j <= 50; j++, val = 1ll*val*dep[i]%yzh) sum[j][id[i]] = val; for (int i = 1; i <= 50; i++) for (int j = 1; j <= n; j++) (sum[i][j] += sum[i][j-1]) %= yzh; m = gi(); while (m--) { u = gi(), v = gi(), k = gi(); printf("%d\n", cal(u, v, k)); }}int main() {work(); return 0; }

转载于:https://www.cnblogs.com/NaVi-Awson/p/8962891.html

你可能感兴趣的文章
Unslider.js Tiny Sample
查看>>
FPGA的学习及注意事项
查看>>
面向对象内存分析
查看>>
Dijkstra BZOJ2763 [JLOI2011]飞行路线
查看>>
前端快捷键
查看>>
重新认识成功、失败、错误、平凡、笨拙
查看>>
【模板】Hash
查看>>
洛谷 1485 火枪打怪
查看>>
Fortran编译器
查看>>
初识go
查看>>
省市县三级联动
查看>>
多IP地址--笔记
查看>>
三级菜单小练习
查看>>
vim 插件管理
查看>>
Guid函数
查看>>
java的加减乘除问题
查看>>
Android 系统常用的权限
查看>>
类型参数约束 : Controller where T : class,new()
查看>>
N卡控制面板把physx设置为cpu
查看>>
Javsssist用InsertAt()方法对语句插桩
查看>>