欢迎访问~原文出处——博客园-zhouzhendong

去博客园看该题解


题意概括

有n个人,每一个人有一个最恨的人。

并且,每一个人有一个权值。

一个人不可以和他最恨的人同时被选中。

现在请你求出在这n个人中选出一些人,使得其权值和最大。

(题解在“心塞史”后面)


心塞史

注:蒟蒻第一次遇见这种基环树题QAQ。

先看样例。

3

10 2

20 3

30 1

瞬间开心。

随便选了一道,居然是水题?

这种类似的题目我貌似做过,好像是找几个环就可以了?

当然,于是我找完环,打完线性动归,一交wa掉,才发现错了。约2000B。

然而可怕的事情没这么容易结束。

我灵机一动~哦!不就是环再连一条线出去嘛……

分3种情况dp一下就可以了。

于是又wa掉了。约2500B。

不过事情总是有转机的。

终于发现这是一道树形dp题,只不过有基环树林(环套树,而且是森林)。

理了2个小时思路,证明了一个连通分量中有且仅有一个环(应该没有自环,一个人总不至于恨自己吧)。

然后开开心心的编。一开始我还以为这一题卡栈空间,为了不用手工扩栈这种卑(ji)鄙(zhi)的手段,我打了bfs序写树形dp。结果有情况没考虑,一发wa掉。3000B

第二天重看这一题。

之前的代码居然一点都看不懂,于是果断again。最终过了。


题解

这一题是一道基环树林的树形动归题。

我们按照最恨的人建边。

首先,我们可以分析,这个图虽然看似是有向的,实际上是无向的。

那么,两个不同的连通分量是互不影响的。所以我们下面只考虑一个连通分量。对于每一个连通分量,把结果加起来就可以了。

对于一个连通分量,总共有n条边。如果只有n-1条边,那么显然是一棵树。

但是多了一条边,那么这棵树中一定会出现一个环。

有,且仅有一个环!

如果不考虑环,只考虑树,那么我们明显可以写一份树形dp。

用dp[i][0]表示节点i不取的最大权值和,dp[i][1]表示取节点i的最大权值和。

那么设S为节点i的儿子集合,那么,

有dp[i][0] = Σmax(dp[x][0] , dp[x][1])  x∈S

dp[i][1] = Σdp[x][0]  x∈S

那么有环的情况怎么办?

其实就是在树形dp的基础上,再限制两个节点不能同时为1。那么实际上我们可以特殊标记。

我们可以先找出环,然后断环,然后分别以环的两个端点为根跑一次树形dp。

但是事情没这么简单。

我们认为不为根的那个端点一定不取,那么我们只要在求完该子树的时候,把它的dp[][0]赋值为-INF即可。

关键是,我们不能随随便便的不取某一个点。

比如说这样一个图:

我假设以红色边断开,上面那个节点为根。

我们不可以不访问下面那个。我们要dp。

但是,我们也不能通过两个五角星之间的连边互相访问。

于是我们要判断边。

这里又有一个问题了。

聪(lan)明(duo)的我一开始选择了记录两个端点。

然后莫名wa掉了。

呵呵,对着hzwer大佬的标算对拍了一会儿。排出一组数据。

是这样的:

两个人互相恨对方。任何我保存点的方法就wa掉了。因为他们之间有两条无向边,一条是找环的时候剪掉的,一条其实还是留着的,然而都被我当作没有了。

所以我只能写了记录边的,好像改变不是很大。

我想找个环的事情总不用说了吧?既然你来做这题了,至少应该会找环吧?。

该程序挺繁的,100w个数字,建议读入优化。

代码

#pragma comment(linker, "/STACK:1024000000,1024000000")
#include <cstring>
#include <algorithm>
#include <cstdio>
#include <cstdlib>
#include <cmath>
using namespace std;
typedef long long LL;
bool isd(char ch){return ''<=ch&&ch<='';}
void read(int &x) {x=;char ch=getchar();while (!isd(ch))ch=getchar();while (isd(ch))x=x*+ch-,ch=getchar();}
void readLL(LL &x){x=;char ch=getchar();while (!isd(ch))ch=getchar();while (isd(ch))x=x*+ch-,ch=getchar();}
LL max(LL a,LL b){
return a>b?a:b;
}
const LL N=+,M=N*,INF=N*N;
struct Gragh{
int cnt,x[M],y[M],nxt[M],fst[N];
void set(){
cnt=;
memset(fst,,sizeof fst);
}
void add(int a,int b){
y[++cnt]=b,x[cnt]=a;
nxt[cnt]=fst[a],fst[a]=cnt;
}
}g;
LL v[N],dp[N][];
int n,p[N],head,tail,fa[N],q[N],vis[N],vis_cnt,dis_use_point;
int pointx,pointy,del_point,del_e;
void dfs(int prev,int rt){
dp[rt][]=,dp[rt][]=v[rt];
for (int i=g.fst[rt];i;i=g.nxt[i]){
int y=g.y[i];
if (y==prev||i==del_e||(i^)==del_e)
continue;
dfs(rt,y);
dp[rt][]+=max(dp[y][],dp[y][]);
dp[rt][]+=dp[y][];
}
if (rt==del_point)
dp[rt][]=-INF;
}
LL solve(int root){
head=tail=;
q[++tail]=root,vis[root]=vis_cnt;
fa[root]=;
pointx=pointy=del_point=;
while (head<tail)
for (int x=q[++head],i=g.fst[x];i;i=g.nxt[i]){
int y=g.y[i];
if (y==fa[x])
continue;
if (vis[y]==vis_cnt){
pointx=x,pointy=y;
del_e=i;
continue;
}
vis[y]=vis_cnt;
fa[y]=x;
q[++tail]=y;
}
for (int i=;i<=tail;i++)
dp[q[i]][]=dp[q[i]][]=;
del_point=pointx;
dfs(,pointy);
LL ans1=max(dp[pointy][],dp[pointy][]);
for (int i=;i<=tail;i++)
dp[q[i]][]=dp[q[i]][]=;
del_point=pointy;
dfs(,pointx);
LL ans2=max(dp[pointx][],dp[pointx][]);
return max(ans1,ans2);
}
int main(){
read(n);
for (int i=;i<=n;i++)
readLL(v[i]),read(p[i]);
g.set();
for (int i=;i<=n;i++)
g.add(i,p[i]),g.add(p[i],i);
memset(vis,,sizeof vis);
vis_cnt=;
LL ans=;
for (int i=;i<=n;i++){
if (vis[i])
continue;
vis_cnt++;
ans+=solve(i);
}
printf("%lld",ans);
return ;
}

BZOJ1040 [ZJOI2008]骑士 基环树林(环套树) 树形动态规划的更多相关文章

  1. BZOJ 1040: [ZJOI2008]骑士 基环加外向树

    1040: [ZJOI2008]骑士 Time Limit: 10 Sec  Memory Limit: 162 MBSubmit: 1190  Solved: 465[Submit][Status] ...

  2. [BZOJ1040][ZJOI2008]骑士 基环树DP

    题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=1040 题目给出了$n$个点和$n$条无向边,即一棵基环树或者基环树森林. 如果题目给的关系 ...

  3. luogu2607/bzoj1040 [ZJOI2008]骑士 (基环树形dp)

    N个点,每个点发出一条边,那么这个图的形状一定是一个基环树森林(如果有重边就会出现森林) 那我做f[0][x]和f[1][x]分别表示对于x子树,x这个点选还是不选所带来的最大价值 然后就变成了这好几 ...

  4. bzoj1040(ZJOI2008)骑士——基环树

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=1040 基环树的模板. 套路就是把环断开,先把一端作为根节点,强制不选:再把另一端作为根节点, ...

  5. [bzoj1040][ZJOI2008]骑士_树形dp_基环树_并查集

    骑士 bzoj-1040 ZJOI-2008 题目大意:n个骑士,每个骑士有权值val和一个讨厌的骑士.如果一个骑士讨厌另一个骑士那么他们将不会一起出战.问出战的骑士最大atk是多少. 注释:$1\l ...

  6. BZOJ1040 骑士 【环套树 树形dp】

    1040: [ZJOI2008]骑士 Time Limit: 10 Sec  Memory Limit: 162 MB Submit: 5611  Solved: 2166 [Submit][Stat ...

  7. BZOJ 1040 [ZJOI2008]骑士 (基环树+树形DP)

    <题目链接> 题目大意: Z国的骑士团是一个很有势力的组织,帮会中汇聚了来自各地的精英.他们劫富济贫,惩恶扬善,受到社会各界的赞扬.最近发生了一件可怕的事情,邪恶的Y国发动了一场针对Z国的 ...

  8. [BZOJ1040][ZJOI2008]骑士(环套树dp)

    1040: [ZJOI2008]骑士 Time Limit: 10 Sec  Memory Limit: 162 MBSubmit: 5816  Solved: 2263[Submit][Status ...

  9. BZOJ 1040 ZJOI 2008 骑士 基环树林+树形DP

    题目大意:有一些骑士.他们每个人都有一个权值.可是因为一些问题,每个骑士都特别讨厌还有一个骑士.所以不能把他们安排在一起.求这些骑士所组成的编队的最大权值和是多少. 思路:首先貌似是有向图的样子,可是 ...

随机推荐

  1. 并行动画组QParallelAnimationGroup

    QParallelAnimationGroup会同时执行添加到该组的所有动画 import sys from PyQt5.QtGui import QPixmap from PyQt5.QtCore ...

  2. ado.net 访问excel

    1 类:OleDbConnection ... 2 connection string : Provider=Microsoft.ACE.OLEDB.12.0;Data Source={excelPa ...

  3. luogu P4568 [JLOI2011]飞行路线

    传送门 看到免费次数\(k\)最多只有10,可以考虑构建\(k+1\)层的分层图,即每一层正常连边,上下两层对应点连边权为0的单向边,最后对所有层里面的\(di_t\)取\(\max\)救星了 #in ...

  4. 2017-2018-2 165X 『Java程序设计』课程 助教总结

    2017-2018-2 165X 『Java程序设计』课程 助教总结 本学期完成的助教工作主要包括: 编写300道左右测试题,用于蓝墨云课下测试: 发布博客三篇:<2017-2018-2 165 ...

  5. Linux定时任务调度

    ⒈概述 任务调度:是指系统在某个时间执行的特定的命令或程序 分类:1)系统任务:有些重要的工作必须周而复始的执行,例如病毒扫描等 2)用户任务:个别用户可能希望定时执行某些程序,例如mysql定时备份 ...

  6. jQuery-contextMenu使用教程

    jQuery-contextMenu使用教程 效果如下图所示.在[右击菜单]处右击,会出现下面的效果. 添加引用 <script src="jQuery-contextMenu-mas ...

  7. Linq基于两个属性的分组

    1.需求 我们看下面的定义 #region 学生类 /// <summary> /// 学生类 /// </summary> class Student { /// <s ...

  8. Cpp读文件、CString转String、String转CString

    场景 C++读取文件 技术点 读取文件 fstream提供了三个类,用来实现c++对文件的操作.(文件的创建.读.写). ifstream -- 从已有的文件读入 ofstream -- 向文件写内容 ...

  9. matplotlib 画图

    matplotlib 画图 1. 画曲线图       Tompson = np.array([0, 0, 0, 0, 0.011, 0.051, 0.15, 0.251, 0.35, 0.44, 0 ...

  10. CONFIG_DEBUG_USER【转】

    转自:https://blog.csdn.net/adaptiver/article/details/12778621 关于CONFIG_DEBUG_USER 把menuconfig中查到的 CONF ...