[bzoj4987]Tree

本文介绍了一种使用树形动态规划解决寻找特定数量节点构成的联通子树,以达到最小路径总和的问题。通过定义状态转移方程,利用枚举大小法进行合并操作,最终得出最优解。

题目大意

从前有棵树。
找出K个点A1,A2,…,Ak。
使得∑dis(AiAi+1),(1<=i<=K-1)最小。

DP

结论一:一定是找一个大小为k的联通子树。
结论二:最优答案一定是所有边权和*2-直径长度。
我们把直径的两个端点叫做关键点。
然后不妨dp,设f[i,j,k]表示子树i里选择了一个包含i的联通子树,含有j个点,选择了k个关键点的最小值。
合并时用枚举size法合并复杂度是O(n^2+nm)。

#include<cstdio>
#include<algorithm>
#define fo(i,a,b) for(i=a;i<=b;i++)
using namespace std;
const int maxn=3000+10,inf=1000000000;
int f[maxn][maxn][3],g[maxn][3];
int h[maxn],go[maxn*2],nxt[maxn*2],dis[maxn*2],size[maxn];
int i,j,k,l,t,n,m,tot,top,ans;
void add(int x,int y,int z){
    go[++tot]=y;
    dis[tot]=z;
    nxt[tot]=h[x];
    h[x]=tot;
}
void merge(int x,int y,int z){
    int i,j,k,l,r,t,lx,ly;
    fo(i,0,size[x]+size[y])
        fo(j,0,2) g[i][j]=f[x][i][j];
    lx=min(size[x],m);
    ly=min(size[y],m);
    fo(i,0,lx)
        fo(j,0,2)
            fo(k,0,min(ly,m-i))
                fo(l,0,2){
                    if (j+l>=3) continue;
                    r=j+l;
                    if (l==1) t=z;else t=2*z;
                    g[i+k][r]=min(g[i+k][r],f[x][i][j]+f[y][k][l]+t);
                }
    size[x]+=size[y];
    fo(i,0,size[x])
        fo(j,0,2) f[x][i][j]=g[i][j];
}
void dfs(int x,int y){
    int i,t;
    fo(i,0,m) f[x][i][0]=f[x][i][1]=f[x][i][2]=inf;
    f[x][1][0]=f[x][1][1]=f[x][1][2]=0;
    size[x]=1;
    t=h[x];
    while (t){
        if (go[t]!=y){
            dfs(go[t],x);
            merge(x,go[t],dis[t]);
        }
        t=nxt[t];
    }
    //f[x][0][0]=0;
    ans=min(ans,f[x][m][2]);
}
int main(){
    freopen("tree.in","r",stdin);freopen("tree.out","w",stdout);
    scanf("%d%d",&n,&m);
    fo(i,1,n-1){
        scanf("%d%d%d",&j,&k,&l);
        add(j,k,l);add(k,j,l);
    }
    ans=inf;
    dfs(1,0);
    printf("%d\n",ans);
}
BZOJ 3306 是一道与树相关的编程挑战题,涉及动态维护一棵树的某些性质,并支持修改操作。题目通常要求处理以下两种操作: 1. **修改节点权值**:将某个节点的权值进行更新。 2. **查询子树中的最小值**:给定一个节点作为根,求其子树中所有节点权值的最小值。 ### 解题思路 为了解决此类问题,需要高效的树结构维护方法。以下是常用的解题策略: #### 树链剖分(Heavy-Light Decomposition) 树链剖分是一种将树结构转化为线性序列的技术,通过这种方式可以使用线段树等数据结构来高效处理区间查询和单点更新操作。具体步骤包括: - 对树进行轻重链划分,使得每个节点到根路径上的边被分成若干个连续的块。 - 使用线段树或树状数组来维护这些连续块的信息,例如最小值、最大值或总和。 #### DFS 序列 + 线段树 另一种常见的方法是利用深度优先搜索(DFS)遍历树,并记录每个节点进入和离开的时间戳。这样可以将每个子树表示为一个连续的区间。随后,使用线段树来维护这个线性序列的最小值信息: - 在预处理阶段,通过 DFS 遍历树并记录每个节点的入栈时间和出栈时间。 - 子树的范围可以通过入栈时间确定,即从该节点的入栈时间到出栈时间之间的所有节点都属于该子树。 - 使用线段树维护区间最小值,支持快速查询和更新操作[^1]。 #### 示例代码 以下是一个简单的线段树实现用于维护最小值的例子: ```cpp #include <bits/stdc++.h> using namespace std; const int MAXN = 250000 + 10; int n, m; int val[MAXN]; vector<int> G[MAXN]; int in_time[MAXN], out_time[MAXN], dfs_clock; void dfs(int u, int fa) { in_time[u] = ++dfs_clock; for (int v : G[u]) { if (v != fa) { dfs(v, u); } } out_time[u] = dfs_clock; } struct SegmentTree { int tree[MAXN << 2]; void build(int node, int l, int r) { if (l == r) { tree[node] = val[l]; return; } int mid = (l + r) / 2; build(node * 2, l, mid); build(node * 2 + 1, mid + 1, r); tree[node] = min(tree[node * 2], tree[node * 2 + 1]); } void update(int node, int l, int r, int pos, int value) { if (l == r) { tree[node] = value; return; } int mid = (l + r) / 2; if (pos <= mid) { update(node * 2, l, mid, pos, value); } else { update(node * 2 + 1, mid + 1, r, pos, value); } tree[node] = min(tree[node * 2], tree[node * 2 + 1]); } int query(int node, int l, int r, int ql, int qr) { if (ql > r || qr < l) return INT_MAX; if (ql <= l && r <= qr) return tree[node]; int mid = (l + r) / 2; return min(query(node * 2, l, mid, ql, qr), query(node * 2 + 1, mid + 1, r, ql, qr)); } } st; int main() { scanf("%d%d", &n, &m); for (int i = 1; i <= n; ++i) { scanf("%d", &val[i]); } for (int i = 1; i < n; ++i) { int u, v; scanf("%d%d", &u, &v); G[u].push_back(v); G[v].push_back(u); } dfs(1, -1); for (int i = 1; i <= n; ++i) { val[in_time[i]] = val[i]; } st.build(1, 1, n); while (m--) { int op, x; scanf("%d%d", &op, &x); if (op == 1) { int y; scanf("%d", &y); st.update(1, 1, n, in_time[x], y); } else { printf("%d\n", st.query(1, 1, n, in_time[x], out_time[x])); } } return 0; } ``` ###
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值