bzoj4557 [JLoi2016]侦察守卫——DP

本文提供了一种解决洛谷编号为4557的问题的算法实现,通过使用深度优先搜索和动态规划,有效地解决了涉及树结构的最短路径问题。代码详细展示了如何初始化数据结构、进行递归调用以及更新状态。

题目:https://www.lydsy.com/JudgeOnline/problem.php?id=4557

见这位的博客:https://www.cnblogs.com/Narh/p/9140326.html

代码如下:

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
typedef long long ll;
int const maxn=500005,inf=1e9;
ll w[maxn],f[maxn][25],g[maxn][25];
int n,d,m,head[maxn],ct;//f向下 g向上 
bool mark[maxn];
struct N{
    int to,next;
    N(int t=0,int n=0):to(t),next(n) {}
}edge[maxn<<1];
void dp(int x,int fa)
{
    if(mark[x])f[x][0]=w[x],g[x][0]=w[x];
    else f[x][0]=0,g[x][0]=0;
    for(int i=1;i<=d;i++)g[x][i]=w[x];
    g[x][d+1]=inf;
    for(int i=head[x];i;i=edge[i].next)
    {
        int v=edge[i].to;
        if(v==fa)continue;
        dp(v,x);
        for(int i=0;i<=d;i++)g[x][i]=min(g[x][i]+f[v][i],g[v][i+1]+f[x][i+1]);
        for(int i=d;i;i--)g[x][i-1]=min(g[x][i],g[x][i-1]);//倒序需要
        f[x][0]=g[x][0];
        for(int i=1;i<=d;i++)f[x][i]+=f[v][i-1];
        for(int i=1;i<=d;i++)f[x][i]=min(f[x][i],f[x][i-1]); 
    }
}
int main()
{
    scanf("%d%d",&n,&d);
    for(int i=1;i<=n;i++)scanf("%lld",&w[i]);
    scanf("%d",&m);
    for(int i=1,x;i<=m;i++)
    {
        scanf("%d",&x);
        mark[x]=1;
    }
    for(int i=1,x,y;i<n;i++)
    {
        scanf("%d%d",&x,&y);
        edge[++ct]=N(y,head[x]);head[x]=ct;
        edge[++ct]=N(x,head[y]);head[y]=ct;
    }
    dp(1,0);
    printf("%lld",f[1][0]);
    return 0;
}

 

转载于:https://www.cnblogs.com/Zinn/p/9140893.html

评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符  | 博主筛选后可见
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值