洛谷P2607 [ZJOI2008]骑士(树形dp)
洛谷P2607 [ZJOI2008]骑士(树形dp)
题目描述
Z国的骑士团是一个很有势力的组织，帮会中汇聚了来自各地的精英。他们劫富济贫，惩恶扬善，受到社会各界的赞扬。
最近发生了一件可怕的事情，邪恶的Y国发动了一场针对Z国的侵略战争。战火绵延五百里，在和平环境中安逸了数百年的Z国又怎能抵挡的住Y国的军队。于是人们把所有的希望都寄托在了骑士团的身上，就像期待有一个真龙天子的降生，带领正义打败邪恶。
骑士团是肯定具有打败邪恶势力的能力的，但是骑士们互相之间往往有一些矛盾。每个骑士都有且仅有一个自己最厌恶的骑士（当然不是他自己），他是绝对不会与自己最厌恶的人一同出征的。
战火绵延，人民生灵涂炭，组织起一个骑士军团加入战斗刻不容缓！国王交给了你一个艰巨的任务，从所有的骑士中选出一个骑士军团，使得军团内没有矛盾的两人（不存在一个骑士与他最痛恨的人一同被选入骑士军团的情况），并且，使得这支骑士军团最具有战斗力。
为了描述战斗力，我们将骑士按照1至N编号，给每名骑士一个战斗力的估计，一个军团的战斗力为所有骑士的战斗力总和。
输入输出格式
输入格式：
输入文件knight.in第一行包含一个正整数N，描述骑士团的人数。
接下来N行，每行两个正整数，按顺序描述每一名骑士的战斗力和他最痛恨的骑士。
输出格式：
输出文件knight.out应包含一行，包含一个整数，表示你所选出的骑士军团的战斗力。
输入输出样例
输入样例#1： 复制
3
10 2
20 3
30 1输出样例#1： 复制
30说明
对于30%的测试数据，满足N ≤ 10；
对于60%的测试数据，满足N ≤ 100；
对于80%的测试数据，满足N ≤ 10 000。
对于100%的测试数据，满足N ≤ 1 000 000，每名骑士的战斗力都是不大于 1 000 000的正整数。
看到标签是树形DP就点进来了
可没想到这题给了一个图？？
不过冷静下来，我们不难发现，这张图实际上只有一个环，也就是传说中的基环树
因此我们按照套路，把一条环上的边破坏掉，然后对两棵独立的树做树形DP
设$f[i][0/1]$表示该节点是否选择时的最大价值
转移的时候枚举孩子是否选择
不过在BZOJ上死活RE
// luogu-judger-enable-o2
#include<cstdio>
#include<vector>
#include<cstring>
#define LL long long
#define Pair pair<int,int>
using namespace std;
const int MAXN=3*1e6+10;
#define getchar() (p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<22,stdin),p1==p2)?EOF:*p1++)
char buf[1<<22],*p1=buf,*p2=buf;
inline int read() {
char c=getchar();int x=0,f=1;
while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}
return x*f;
}
struct Edge {
int u, v, nxt;
}edge[MAXN];
int head[MAXN], num=0;
void AddEdge(int x,int y) {
edge[num] = (Edge){x, y, head[x]};
head[x] = num++;
}
int val[MAXN], vis[MAXN], BeginEdge;
LL f[MAXN][2];
Pair Begin;
void FindCircle(int now,int fa) {
vis[now] = 1;
for(int i = head[now];i != -1;i = edge[i].nxt) {
if( (i ^ 1) == fa) continue;
if(vis[ edge[i].v ]) {BeginEdge = i;Begin.first = now;Begin.second = edge[i].v;continue ;}
FindCircle(edge[i].v, i);
}
}
LL Dp(int now,int fa) {
f[now][0] = 0;
f[now][1] = val[now];
for(int i = head[now];i != -1;i = edge[i].nxt) {
if((i ^ 1) == fa) continue;
if((i == BeginEdge) || ((i ^ 1) == BeginEdge)) continue;
Dp(edge[i].v, i);
f[now][0] += max(f[edge[i].v][1], f[edge[i].v][0]);
f[now][1] += f[edge[i].v][0];
}
return f[now][0];
}
int main() {
#ifdef WIN32
freopen("a.in","r",stdin);
#endif
int N = read();
memset(head, -1, sizeof(int) * N *4);
for(register int i=1;i<=N;i++) {
val[i] = read();
int x = read();
AddEdge(i,x),AddEdge(x,i);
}
long long ans = 0;
for(register int i=1;i<=N;i++) {
if(!vis[i]) {
FindCircle(i,-2);
ans += max(Dp(Begin.first,-1),Dp(Begin.second,-1));
}
}
printf("%lld", ans);
return 0;
}