鱼C论坛

 找回密码
 立即注册
查看: 1062|回复: 2

树链剖分代码求助

[复制链接]
发表于 2023-9-28 10:36:20 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能^_^

您需要 登录 才可以下载或查看,没有账号?立即注册

x
本帖最后由 柿子饼同学 于 2023-9-28 10:37 编辑
题目: P4315
不要用 GPT 回复....

求助求助 , 样例过了但是 0 分 , 必有鱼币重谢
多谢了...
  1. #include <bits/stdc++.h>
  2. using namespace std;

  3. const int N = 1e5 + 5;

  4. int n, u, v, w, k;
  5. vector<int> e[N], d[N], id[N];
  6. // e 是树, d 是每条边的边权, id 是每条边的编号
  7. // 因为不会链式前向星...
  8. map<int, int> dict;
  9. // dict 返回一个编号的边对应的节点编号
  10. int data[N];
  11. // 每个边的值, 为了方便都会转移到这条边上深度较大的点上
  12. int size[N], hc[N], pa[N], dep[N];
  13. // hc 是重孩子
  14. int dfn[N], rdfn[N], top[N], tot;
  15. string op;

  16. void dfs1(int u, int par){
  17.     size[u] = 1;
  18.     pa[u] = par;
  19.     dep[u] = dep[par] + 1;

  20.     for(int i = 0; i < e[u].size(); i++){
  21.         if(e[u][i] == par) continue;
  22.         data[e[u][i]] = d[u][i];
  23.         dict[id[u][i]] = e[u][i];
  24.         dfs1(e[u][i], u);
  25.         size[u] += size[e[u][i]];
  26.         if(size[e[u][i]] > size[hc[u]]) hc[u] = e[u][i];
  27.     }
  28. }

  29. void dfs2(int u, int t){   
  30.     top[u] = t;
  31.     dfn[u] = ++tot;
  32.     rdfn[tot] = u;

  33.     if(!hc[u]) return ;
  34.     dfs2(hc[u], t);

  35.     for(auto ed : e[u]){
  36.         if(ed == pa[u] || ed == hc[u]) continue;
  37.         dfs2(ed, ed);
  38.     }
  39. }

  40. struct Segment_Tree{
  41.     #define lc u<<1
  42.     #define rc u<<1|1

  43.     struct Node{
  44.         int l, r, tr, add, maxv;
  45.     } t[N << 2];

  46.     void pushup(int u){
  47.         t[u].maxv = max(t[lc].maxv, t[rc].maxv);
  48.     }

  49.     void pushdown(int u){
  50.         if(t[u].tr){
  51.             t[lc].tr = t[rc].tr = t[u].tr;
  52.             t[lc].maxv = t[rc].maxv = t[u].tr;
  53.             t[lc].add = t[rc].add = 0;
  54.             t[u].tr = 0;
  55.             
  56.         }
  57.         if(t[u].add){
  58.             t[lc].maxv += t[u].add;
  59.             t[lc].add += t[u].add;
  60.             t[rc].maxv += t[u].add;
  61.             t[rc].add += t[u].add;
  62.             t[u].add = 0;
  63.         }
  64.     }

  65.     void build(int u, int l, int r){
  66.         t[u] = {l, r};
  67.         if(l == r){
  68.             t[u].maxv = data[rdfn[l]];
  69.             t[u].add = t[u].tr = 0;
  70.             return ;
  71.         }
  72.         auto mid = (l + r) >> 1;
  73.         build(lc, l, mid);
  74.         build(rc, mid + 1, r);
  75.         pushup(u);
  76.     }

  77.     void trans(int u, int l, int r, int val){ // 全部改
  78.         if(l <= t[u].l && t[u].r <= r){
  79.             t[u].tr = val;
  80.             t[u].maxv = val;
  81.             return ;
  82.         }
  83.         auto mid = (t[u].l + t[u].r) >> 1;
  84.         pushdown(u);
  85.         if(l <= mid) trans(lc, l, r, val);
  86.         if(r > mid) trans(rc, l, r, val);
  87.         pushup(u);
  88.     }

  89.     void modify(int u, int l, int r, int val){ // 区间加
  90.         if(l <= t[u].l && t[u].r <= r){
  91.             t[u].add += val;
  92.             t[u].maxv += val;
  93.             return ;
  94.         }
  95.         auto mid = (t[u].l + t[u].r) >> 1;
  96.         pushdown(u);
  97.         if(l <= mid) modify(lc, l, r, val);
  98.         if(r > mid) modify(rc, l, r, val);
  99.         pushup(u);
  100.     }

  101.     int query(int u, int l, int r){
  102.         if(l <= t[u].l && t[u].r <= r){
  103.             return t[u].maxv;
  104.         }
  105.         auto mid = (t[u].l + t[u].r) >> 1;
  106.         int res = 0;
  107.         pushdown(u);
  108.         if(l <= mid) res = max(res, query(lc, l, r));
  109.         if(r > mid) res = max(res, query(rc, l, r));
  110.         return res;
  111.     }
  112. } sgt;

  113. int Max(int u, int v){
  114.     int res = 0;
  115.     while(top[u] != top[v]){
  116.         if(dep[top[u]] < dep[top[v]]) swap(u, v);
  117.         res = max(res, sgt.query(1, dfn[top[u]], dfn[u]));
  118.         u = pa[top[u]];
  119.     }
  120.     res = max(res, sgt.query(1, min(dfn[u], dfn[v]), max(dfn[u], dfn[v])));
  121.     return res;
  122. }

  123. void Change(int k, int w){
  124.     sgt.trans(1, dfn[k], dfn[k], w);
  125. }

  126. void Cover(int u, int v, int w){
  127.     while(top[u] != top[v]){
  128.         if(dep[top[u]] < dep[top[v]]) swap(u, v);
  129.         sgt.trans(1, dfn[top[u]], dfn[u], w);
  130.         u = pa[top[u]];
  131.     }
  132.     sgt.trans(1, min(dfn[u], dfn[v]), max(dfn[u], dfn[v]), w);
  133. }

  134. void Add(int u, int v, int w){
  135.     while(top[u] != top[v]){
  136.         if(dep[top[u]] < dep[top[v]]) swap(u, v);
  137.         sgt.modify(1, dfn[top[u]], dfn[u], w);
  138.         u = pa[top[u]];
  139.     }
  140.     sgt.modify(1, min(dfn[u], dfn[v]), max(dfn[u], dfn[v]), w);
  141. }

  142. int main(){
  143.     ios::sync_with_stdio(0);
  144.     cin.tie(0);

  145.     cin >> n;
  146.     for(int i = 1; i < n; i++){
  147.         cin >> u >> v >> w;
  148.         e[u].emplace_back(v);
  149.         e[v].emplace_back(u);
  150.         d[u].emplace_back(w);
  151.         d[v].emplace_back(w);
  152.         id[u].emplace_back(i);
  153.         id[v].emplace_back(i);
  154.     }

  155.     dfs1(1, 0);
  156.     dfs2(1, 1);

  157.     sgt.build(1, 1, n);

  158.     while(cin >> op){
  159.         if(op[0] == 'S') break;
  160.         if(op[0] == 'C'){
  161.             if(op[1] == 'h'){
  162.                 cin >> k >> w;
  163.                 Change(dict[k], w);
  164.             }
  165.             else{
  166.                 cin >> u >> v >> w;
  167.                 Cover(u, v, w);
  168.             }
  169.         }
  170.         else if(op[0] == 'A'){
  171.             cin >> u >> v >> w;
  172.             Add(u, v, w);
  173.         }
  174.         else{
  175.             cin >> u >> v;
  176.             cout << Max(u, v) << '\n';
  177.         }
  178.     }

  179.     return 0;
  180. }
复制代码




小甲鱼最新课程 -> https://ilovefishc.com
回复

使用道具 举报

发表于 2023-9-28 10:36:36 | 显示全部楼层
这是一段树链剖分的代码,用于解决树上路径问题。具体来说,它可以支持以下操作:

1. 修改某个边的权值
2. 将某个路径上的所有边的权值都改为同一个值
3. 将某个路径上的所有边的权值都加上同一个值
4. 查询某个路径上的最大权值

在这段代码中,树的节点编号从1开始,边的编号从1开始。节点的子节点存储在vector<int> e[N]中,每条边的权值存储在vector<int> d[N]中,每条边的编号存储在vector<int> id[N]中。通过dfs1函数可以预处理出每个节点的重孩子和每个节点的深度。通过dfs2函数可以预处理出每个节点所在重链的顶端节点、每个节点的dfs序和反向dfs序。通过Segment_Tree结构体可以维护每个节点所在重链的区间最大值和区间加操作。

具体操作实现如下:

1. 修改某个边的权值:通过dict映射得到这条边的两个端点,然后调用sgt.trans函数将这个端点对应的节点的权值改为w。

2. 将某个路径上的所有边的权值都改为同一个值:从u节点开始,沿着重链向上走,每次调用sgt.trans函数将当前节点到它所在重链顶端节点这条路径上的所有节点的权值都改为w。

3. 将某个路径上的所有边的权值都加上同一个值:与操作2类似,只是调用的是sgt.modify函数。

4. 查询某个路径上的最大权值:从u节点开始,沿着重链向上走,每次调用sgt.query函数查询当前节点到它所在重链顶端节点这条路径上的最大权值。最后将所有查询到的最大权值取一个最大值即可。

希望这些解释能够帮助你理解这段代码。
小甲鱼最新课程 -> https://ilovefishc.com
回复 支持 0 反对 1

使用道具 举报

发表于 2023-9-28 10:36:53 | 显示全部楼层
代码中存在一些问题,导致你的程序得不到正确的结果。以下是我发现的一些问题和相应的修改建议:

1. 在 dfs1 函数中,初始化 hc[u] 时应为 hc[u] = 0; 而不是 hc[u] = e[u][0];。因为可能存在没有孩子节点的情况。

2. 在 dfs2 函数中,递归调用 dfs2(ed, ed); 应改为 dfs2(ed, t);。因为轻边的重链顶端节点应与当前节点保持一致。

3. 在 Segment_Tree 结构体中,成员变量 addtr 的初始值应设为 0,而不是未定义的随机值。

4. 在 trans 函数中,更新左右子节点时应先推送下去再更新,即将 pushdown(u) 放在两个递归之前。

5. 在 modify 函数中,同样需要将 pushdown(u) 放在两个递归之前。

6. 在 query 函数中,同样需要将 pushdown(u) 放在两个递归之前。

7. 在 Max 函数中,获取查询范围时,应使用 max(dfn[u], dfn[v]) + 1,即将 max(dfn[u], dfn[v]) 的值加 1。

请按照上述修改建议进行代码修改,然后重新运行程序,看看是否能够得到正确的结果。

球一个最佳答案谢谢啦!这对我非常重要!
小甲鱼最新课程 -> https://ilovefishc.com
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

小黑屋|手机版|Archiver|鱼C工作室 ( 粤ICP备18085999号-1 | 粤公网安备 44051102000585号)

GMT+8, 2025-6-10 00:33

Powered by Discuz! X3.4

© 2001-2023 Discuz! Team.

快速回复 返回顶部 返回列表