[HDU 5886] Tower Defence (树形DP)

2024-06-21 19:38
文章标签 dp hdu 树形 tower defence 5886

本文主要是介绍[HDU 5886] Tower Defence (树形DP),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

HDU - 5886

给定一棵树,每条边都有一个权值
定义一条边的战术价值为割去这条边后
剩下的图中最长链的长度
求随机割掉一条边,战术价值的期望


题目要求乘以 N1 ,所以直接把概率消掉了
所以只要求割去每条边后的最长链长度即可
这就转化为了一个树形dp
首先割去一条边 u>v 后,最长链可能在 v 为根的子树中
这个简单地做一次树形dp就能得出
另外最长链还可能是在 u 所在的联通块,这个就比较麻烦了
在 u联通块的最长链可能有几种情况

  1. 最长链在 u 之上且不经过 u
  2. 最长链在 u 之下且不经过 u,即最长链在 v的兄弟子树里
  3. 最长链经过 u

其中最长链经过 u又分为两种情况

  1. 最长链由 u 之上到 u 最长的一条链和 u 之下到 u 的最长的一条链连接而得
  2. 最长链由 v 的两个兄弟子树各出一条最长链在 u <script id="MathJax-Element-162" type="math/tex">u</script> 处连接而得
    分类讨论所有情况,然后取 max即可
#pragma comment(linker, "/STACK:102400000,102400000")
#include <cstdio>
#include <iostream>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <cctype>
#include <map>
#include <set>
#include <queue>
#include <bitset>
#include <string>
#include <complex>
using namespace std;
typedef pair<int,int> Pii;
typedef long long LL;
typedef unsigned long long ULL;
typedef double DBL;
typedef long double LDBL;
#define MST(a,b) memset(a,b,sizeof(a))
#define CLR(a) MST(a,0)
#define SQR(a) ((a)*(a))
#define PCUT puts("\n----------")const int maxn=1e5+10;;
struct Graph
{int ndn, edn, last[maxn];int u[2*maxn], v[2*maxn], w[2*maxn], nxt[2*maxn];void init(int _n){ndn=_n; edn=0; MST(last,-1);}void adde(int _u, int _v, int _w){u[edn]=_u; v[edn]=_v; w[edn]=_w;nxt[edn]=last[_u];last[_u] = edn++;}
};int N;
Graph G;
int dp[maxn][4], pos[maxn][4], dp2[maxn][3], pos2[maxn][3], res[maxn];
LL ans;void add(int,int,int);
void add2(int,int,int);
int dfs0(int,int);
int dfs1(int,int,int,int);int main()
{#ifdef LOCALfreopen("in.txt", "r", stdin);
//  freopen("out.txt", "w", stdout);#endifint T;scanf("%d", &T);for(int ck=1; ck<=T; ck++){scanf("%d", &N);G.init(N);for(int i=1, u, v, w; i<N; i++){scanf("%d%d%d", &u, &v, &w);G.adde(u,v,w); G.adde(v,u,w);}ans = 0;dfs0(1,0);dfs1(1,0,0,0);printf("%lld\n", ans);}return 0;
}void add(int u, int v, int val)
{dp[u][3] = val; pos[u][3] = v;for(int i=2; i>=0; i--) if(dp[u][i] < dp[u][i+1]){swap(dp[u][i], dp[u][i+1]);swap(pos[u][i], pos[u][i+1]);}
}void add2(int u, int v, int val)
{dp2[u][2] = val; pos2[u][2] = v;for(int i=1; i>=0; i--) if(dp2[u][i] < dp2[u][i+1]){swap(dp2[u][i], dp2[u][i+1]);swap(pos2[u][i], pos2[u][i+1]);}
}int dfs0(int u, int f)
{MST(dp[u],-1); MST(pos[u],-1); MST(dp2[u], -1); MST(pos2[u], -1);res[u] = 0;bool leaf=1;for(int e=G.last[u], v; ~e; e=G.nxt[e]) if((v=G.v[e])!=f){leaf=0;add(u, v, dfs0(v,u)+G.w[e]);add2(u, v, res[v]);res[u] = max(res[u], res[v]);}if(leaf) add(u, u, 0);res[u] = max(res[u], max(dp[u][0], dp[u][0]+dp[u][1]));return dp[u][0];
}int dfs1(int u, int f, int one, int two)
{for(int e=G.last[u], v; ~e; e=G.nxt[e]) if((v=G.v[e])!=f){int tmax=0, tp=-1, tmax2=0, none, ntwo;for(int i=0; i<3 && ~dp[u][i]; i++){if(pos[u][i] == v) continue;if(dp[u][i]>tmax) tmax = dp[u][i], tp=i;}none = max(tmax, one);ntwo = max(two, one+tmax);for(int i=0; i<3 && ~dp[u][i]; i++){if(pos[u][i] == v || i==tp) continue;tmax2 = max(tmax2, dp[u][i]);}ntwo = max(ntwo, tmax + tmax2);if(pos[u][0] == v) ntwo = max(ntwo, dp2[u][1]);else ntwo = max(ntwo, dp2[u][0]);ans += max(ntwo, res[v]);dfs1(v, u, none+G.w[e], ntwo);}}

这篇关于[HDU 5886] Tower Defence (树形DP)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1082150

相关文章

hdu4826(三维DP)

这是一个百度之星的资格赛第四题 题目链接:http://acm.hdu.edu.cn/contests/contest_showproblem.php?pid=1004&cid=500 题意:从左上角的点到右上角的点,每个点只能走一遍,走的方向有三个:向上,向下,向右,求最大值。 咋一看像搜索题,先暴搜,TLE,然后剪枝,还是TLE.然后我就改方法,用DP来做,这题和普通dp相比,多个个向上

hdu1011(背包树形DP)

没有完全理解这题, m个人,攻打一个map,map的入口是1,在攻打某个结点之前要先攻打其他一个结点 dp[i][j]表示m个人攻打以第i个结点为根节点的子树得到的最优解 状态转移dp[i][ j ] = max(dp[i][j], dp[i][k]+dp[t][j-k]),其中t是i结点的子节点 代码如下: #include<iostream>#include<algorithm

hdu4865(概率DP)

题意:已知前一天和今天的天气概率,某天的天气概率和叶子的潮湿程度的概率,n天叶子的湿度,求n天最有可能的天气情况。 思路:概率DP,dp[i][j]表示第i天天气为j的概率,状态转移如下:dp[i][j] = max(dp[i][j, dp[i-1][k]*table2[k][j]*table1[j][col] )  代码如下: #include <stdio.h>#include

usaco 1.3 Mixing Milk (结构体排序 qsort) and hdu 2020(sort)

到了这题学会了结构体排序 于是回去修改了 1.2 milking cows 的算法~ 结构体排序核心: 1.结构体定义 struct Milk{int price;int milks;}milk[5000]; 2.自定义的比较函数,若返回值为正,qsort 函数判定a>b ;为负,a<b;为0,a==b; int milkcmp(const void *va,c

usaco 1.1 Broken Necklace(DP)

直接上代码 接触的第一道dp ps.大概的思路就是 先从左往右用一个数组在每个点记下蓝或黑的个数 再从右到左算一遍 最后取出最大的即可 核心语句在于: 如果 str[i] = 'r'  ,   rl[i]=rl[i-1]+1, bl[i]=0 如果 str[i] = 'b' ,  bl[i]=bl[i-1]+1, rl[i]=0 如果 str[i] = 'w',  bl[i]=b

poj 3974 and hdu 3068 最长回文串的O(n)解法(Manacher算法)

求一段字符串中的最长回文串。 因为数据量比较大,用原来的O(n^2)会爆。 小白上的O(n^2)解法代码:TLE啦~ #include<stdio.h>#include<string.h>const int Maxn = 1000000;char s[Maxn];int main(){char e[] = {"END"};while(scanf("%s", s) != EO

hdu 2093 考试排名(sscanf)

模拟题。 直接从教程里拉解析。 因为表格里的数据格式不统一。有时候有"()",有时候又没有。而它也不会给我们提示。 这种情况下,就只能它它们统一看作字符串来处理了。现在就请出我们的主角sscanf()! sscanf 语法: #include int sscanf( const char *buffer, const char *format, ... ); 函数sscanf()和

hdu 2602 and poj 3624(01背包)

01背包的模板题。 hdu2602代码: #include<stdio.h>#include<string.h>const int MaxN = 1001;int max(int a, int b){return a > b ? a : b;}int w[MaxN];int v[MaxN];int dp[MaxN];int main(){int T;int N, V;s

hdu 1754 I Hate It(线段树,单点更新,区间最值)

题意是求一个线段中的最大数。 线段树的模板题,试用了一下交大的模板。效率有点略低。 代码: #include <stdio.h>#include <string.h>#define TREE_SIZE (1 << (20))//const int TREE_SIZE = 200000 + 10;int max(int a, int b){return a > b ? a :

hdu 1166 敌兵布阵(树状数组 or 线段树)

题意是求一个线段的和,在线段上可以进行加减的修改。 树状数组的模板题。 代码: #include <stdio.h>#include <string.h>const int maxn = 50000 + 1;int c[maxn];int n;int lowbit(int x){return x & -x;}void add(int x, int num){while