zoukankan      html  css  js  c++  java
  • hdu5834 Magic boy Bi Luo with his excited tree 【树形dp】

    题目链接

    hdu5834

    题解

    思路很粗犷,实现很难受

    (f[i][0|1])表示向子树走回来或不回来的最大收益
    (g[i][0|1])表示向父亲走走回来或不回来的最大收益
    再设(h[i])(f[i][0])的次优收益

    对于(f[i][1]),贪心选择所有(f[v][1] - 2 * w ge 0)的子树即可
    对于(f[i][0]),贪心选择所有没有被选的子树的(f[v][0] - w le 0)的最大值 或者 被选子树(f[v][1] - 2 * w)改成(f[v][0] - w)后多产生收益的最大值
    同时维护次优(h[v])

    对于(g[i][1]),设父亲为(v),就等于(f[v][1] + g[v][1])再减去(i)(f[v][1])所作出的贡献【因为往父亲走要忽视(i)这课子树】
    对于(g[i][0])也是类似的,但是由于忽视(i)这课子树后(f[i][0])的决策可能发生改变,所以要在之前算好次优决策(h[v])

    这种树形dp简单题都做不出了

    #include<iostream>
    #include<cstdio>
    #include<cmath>
    #include<cstring>
    #include<algorithm>
    #define LL long long int
    #define Redge(u) for (int k = head[u],to; k; k = ed[k].nxt)
    #define REP(i,n) for (int i = 1; i <= (n); i++)
    #define BUG(s,n) for (int i = 1; i <= (n); i++) cout<<s[i]<<' '; puts("");
    using namespace std;
    const int maxn = 100005,maxm = 100005,INF = 1000000000;
    inline int read(){
    	int out = 0,flag = 1; char c = getchar();
    	while (c < 48 || c > 57){if (c == '-') flag = -1; c = getchar();}
    	while (c >= 48 && c <= 57){out = (out << 3) + (out << 1) + c - 48; c = getchar();}
    	return out * flag;
    }
    int head[maxn],ne = 2;
    struct EDGE{int to,nxt,w;}ed[maxn << 1];
    inline void build(int u,int v,int w){
    	ed[ne] = (EDGE){v,head[u],w}; head[u] = ne++;
    	ed[ne] = (EDGE){u,head[v],w}; head[v] = ne++;
    }
    int n,fa[maxn],d[maxn],w[maxn],f[maxn][2],g[maxn][2],h[maxn],way[maxn];
    //cal son
    void dfs1(int u){
    	f[u][0] = f[u][1] = w[u];
    	int mx = -INF,v,tmp,mx2 = -INF;
    	Redge(u) if ((to = ed[k].to) != fa[u]){
    		fa[to] = u; d[to] = ed[k].w; dfs1(to);
    		if (f[to][1] - 2 * d[to] >= 0){
    			f[u][1] += f[to][1] - 2 * d[to];
    			tmp = (f[to][0] - d[to]) - (f[to][1] - 2 * d[to]);
    			if (tmp > mx) mx2 = mx,mx = tmp,v = to;
    			else if (tmp > mx2) mx2 = tmp;
    		}
    		else if ((tmp = f[to][0] - d[to]) >= 0){
    			if (tmp > mx) mx2 = mx,mx = tmp,v = to;
    			else if (tmp > mx2) mx2 = tmp;
    		}
    	}
    	if (mx >= 0) f[u][0] = f[u][1] + mx,way[u] = v;
    	else f[u][0] = f[u][1],way[u] = 0;
    	if (mx2 >= 0) h[u] = f[u][1] + mx2;
    	else h[u] = f[u][1];
    }
    //cal father
    void dfs2(int u){
    	int v = fa[u];
    	//back
    	if (f[u][1] - 2 * d[u] >= 0)
    		g[u][1] = max(0,f[v][1] + g[v][1] - (f[u][1] - 2 * d[u]) - 2 * d[u]);
    	else g[u][1] = max(0,f[v][1] + g[v][1] - 2 * d[u]);
    	//not back
    	if (f[u][1] - 2 * d[u] >= 0){
    		g[u][0] = max(0,f[v][1] + g[v][0] - (f[u][1] - 2 * d[u]) - d[u]);
    		if (way[v] == u)
    			g[u][0] = max(g[u][0],h[v] + g[v][1] - (f[u][1] - 2 * d[u]) - d[u]);
    		else g[u][0] = max(g[u][0],f[v][0] + g[v][1] - (f[u][1] - 2 * d[u]) - d[u]);
    	}
    	else{
    		g[u][0] = max(0,f[v][1] + g[v][0] - d[u]);
    		if (way[v] == u)
    			g[u][0] = max(g[u][0],h[v] + g[v][1] - d[u]);
    		else g[u][0] = max(g[u][0],f[v][0] + g[v][1] - d[u]);
    	}
    	Redge(u) if ((to = ed[k].to) != fa[u])
    		dfs2(to);
    }
    int main(){
    	int T = read();
    	REP(C,T){
    		n = read(); ne = 2;
    		REP(i,n) w[i] = read(),head[i] = 0;
    		int a,b,w;
    		for (int i = 1; i < n; i++){
    			a = read(); b = read(); w = read();
    			build(a,b,w);
    		}
    		dfs1(1);
    		dfs2(1);
    		printf("Case #%d:
    ",C);
    		REP(i,n) printf("%d
    ",max(f[i][1] + g[i][0],f[i][0] + g[i][1]));
    	}
    	return 0;
    }
    
    
  • 相关阅读:
    实现第三方系统单点登录
    python爬虫
    webot设备motor的api
    webots学习
    python学习算术运算
    python快捷键与命令函数
    python学习构造和析构
    python学习对象相关的bif
    python学习对象:拾遗
    matlab基础知识
  • 原文地址:https://www.cnblogs.com/Mychael/p/9012958.html
Copyright © 2011-2022 走看看