码迷,mamicode.com
首页 > 其他好文 > 详细

题解 poj3585 Accumulation Degree (树形dp)(二次扫描和换根法)

时间:2019-02-22 21:18:37      阅读:200      评论:0      收藏:0      [点我收藏+]

标签:扫描法   ext   含义   href   image   false   通过   hhhh   algorithm   

写一篇题解,以纪念调了一个小时的经历(就是因为边的数组没有乘2 phhhh QAQ)


题目

题目大意:找一个点使得从这个点出发作为源点,流出的流量最大,输出这个最大的流量。 

以这道题来介绍二次扫描和换根法

作为一道不定根的树形DP,如果直接对每个点进行DP,可能时间会炸掉

但是,优秀的二次换根和扫描法可以再O(n^2)内解决问题。

二次扫描的含义:(来自lyd 算法竞赛进阶指南)

第一次扫描:任选一个节点为根节点(我会选1)在树上进行树形DP,在回溯时,从儿子节点向父节点(从底向上)进行状态转移

第二次扫描:从刚才选的根出发,对树进行dfs,在每次递归前进行自顶向下的推导,计算"换根"后的解

在第一次扫描时,我们可以算出以节点u为根的子树中,从u流向其子树的最大流量

技术图片

 

在第二次扫描时,我们通过dfs,可以自上而下的求出以节点u为根,从u流向整个流域(其子树)的最大流量

技术图片

 

当我们从节点u到节点v时,已经求出F[u],而从u到v的流量为min(D[v],w(u,v)),所以从u流向v的其他部分的流量就是F[u]-min(D[v],w(u,v)),所以拿它再跟w(u,v)取min

 

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#define R register int
using namespace std;
const int N=200010;
struct edge{
    int v,nxt,w;
}e[N<<1];
int t,n,ans,cnt;
int d[N],f[N],fir[N],deg[N];
bool vis[N];

inline void add(int u,int v,int w) {e[++cnt].v=v,e[cnt].w=w,e[cnt].nxt=fir[u],fir[u]=cnt;}

inline int g()
{
    R ret=0,fix=1; register char ch;
    while(!isdigit(ch=getchar())) fix=ch==-?-1:fix;
    do ret=(ret<<3)+(ret<<1)+(ch^48); while(isdigit(ch=getchar()));
    return ret*fix;
}

void dp(int u)
{
    vis[u]=true,d[u]=0;
    for(R i=fir[u];i;i=e[i].nxt)
    {
        R v=e[i].v; if(vis[v]) continue;
        dp(v);
        if(deg[v]==1) d[u]+=e[i].w;
        else d[u]+=min(d[v],e[i].w);
    }
}

void dfs(int u)
{
    vis[u]=true;
    for(R i=fir[u];i;i=e[i].nxt)
    {
        R v=e[i].v; if(vis[v]) continue;
        if(deg[u]==1) f[v]=d[v]+e[i].w;
        else f[v]=d[v]+min(f[u]-min(d[v],e[i].w),e[i].w);
        dfs(v);
    }
}

int main()
{
    t=g();
    while(t--)
    {
        memset(vis,false,sizeof(vis));
        memset(fir,0,sizeof(fir));
        memset(deg,0,sizeof(deg)); ans=0,cnt=0;
        n=g();
        if(n==0||n==1) {putchar(0),putchar(\n);continue;}
        for(R i=1;i<n;i++) {R u=g(),v=g(),w=g(); add(u,v,w),add(v,u,w); deg[u]++,deg[v]++;}
        R rt=1;
        dp(rt); f[rt]=d[rt];
        memset(vis,false,sizeof(vis));
        dfs(rt);
        for(R i=1;i<=n;i++) ans=max(ans,f[i]);
        printf("%d\n",ans);
    }
    return 0;
}

(我太菜了。。。。QAQ)

 

题解 poj3585 Accumulation Degree (树形dp)(二次扫描和换根法)

标签:扫描法   ext   含义   href   image   false   通过   hhhh   algorithm   

原文地址:https://www.cnblogs.com/Jackpei/p/10420639.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!