bzoj1487: [HNOI2009]无归岛

传送门
如果原图不是仙人掌,就是一个傻逼的树形dp
可是他就是仙人掌。
考虑破环为链。
强制断开一条边。在链上跑dp,在最后状态中求出合法的最大值就ok辣。

#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define N 100005
#define inf 2000000000
using namespace std;
int n,m,x,y,tot,tim;
int f[N],g[N],val[N],head[N],dfn[N],fa[N];
struct edge{
    int to,next;
}e[N*4];
inline int read(){
    int x=0,f=1;
    char ch=getchar();
    while (ch<'0'||ch>'9'){
        if (ch=='-') f=-1;
        ch=getchar();
    }
    while (ch>='0'&&ch<='9'){
        x=x*10+ch-48;
        ch=getchar();
    }
    return x*f;
}
inline void ins(int x,int y){
    e[++tot].to=y;
    e[tot].next=head[x];
    head[x]=tot;
}
void dp(int rt,int x){
    int u0=0,u1=0,v0,v1;
    for (int i=x;i!=rt;i=fa[i]){
        v1=u0+f[i];
        v0=u1+g[i];
        u0=v0;
        u1=max(v1,v0);
    }
    g[rt]+=u1;
    u0=-inf,u1=0;
    for (int i=x;i!=rt;i=fa[i]){
        v1=u0+f[i];
        v0=u1+g[i];
        u0=v0;
        u1=max(v1,v0);
    }
    f[rt]+=u0;
}
void dfs(int x){
    dfn[x]=++tim;
    for (int i=head[x];i;i=e[i].next)
        if (!dfn[e[i].to]){
            fa[e[i].to]=x;
            dfs(e[i].to);
        }
    f[x]=val[x];
    for (int i=head[x];i;i=e[i].next)
        if (dfn[e[i].to]>dfn[x]&&fa[e[i].to]!=x) dp(x,e[i].to);
}
int main(){
    n=read();
    m=read();
    for (int i=1;i<=m;i++){
        x=read();
        y=read();
        ins(x,y);
        ins(y,x);
    }
    for (int i=1;i<=n;i++) val[i]=read();
    dfs(1);
    printf("%d",max(f[1],g[1]));
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值