牛客 - 血压游戏


传送门:牛客 - 血压游戏

思路分析

显然只有深度一样的才会影响,所以将深度相同的节点单独拿出来建立虚树进行dp

AC代码

#pragma GCC optimize(2)
#include <bits/stdc++.h>
#define fi first
#define se second
#define ll long long
#define pb push_back
#define mp make_pair
#define fun function
#define sz(x) (x).size()
#define lowbit(x) (x)&(-x)
#define all(x) (x).begin(),(x).end()
#define mem(a,b) memset(a,b,sizeof(a))

namespace FastIO {
#define BUF_SIZE 100000
#define OUT_SIZE 100000
    bool IOerror=0;
    inline char nc() {
        static char buf[BUF_SIZE],*p1=buf+BUF_SIZE,*pend=buf+BUF_SIZE;
        if(p1==pend) {
            p1=buf;
            pend=buf+fread(buf,1,BUF_SIZE,stdin);
            if(pend==p1) {
                IOerror=1;
                return -1;
            }
        }
        return *p1++;
    }
    inline bool blank(char ch) {
        return ch==' '||ch=='\n'||ch=='\r'||ch=='\t';
    }
    template<class T> inline bool read(T &x) {
        bool sign=0;
        char ch=nc();
        x=0;
        for(; blank(ch); ch=nc());
        if(IOerror)return false;
        if(ch=='-')sign=1,ch=nc();
        for(; ch>='0'&&ch<='9'; ch=nc())x=x*10+ch-'0';
        if(sign)x=-x;
        return true;
    }
    template<class T,class... U>bool read(T& h,U&... t) {
        return read(h)&&read(t...);
    }
#undef OUT_SIZE
#undef BUF_SIZE
};
using namespace std;
using namespace FastIO;

mt19937 rng(chrono::steady_clock::now().time_since_epoch().count());

const int INF = 0x3f3f3f3f;
const int N = 2e5+10;


vector<int>e[N],tmp[N];
int dep[N],f[N],top[N],son[N],siz[N],dfn[N],tim;

void dfs1(int u,int fa) {
    dep[u]=dep[fa]+1;
    tmp[dep[u]].pb(u);
    siz[u]=1;
    son[u]=0;
    f[u]=fa;
    for(auto v:e[u]) {
        if(v==fa) continue;
        dfs1(v,u);
        siz[u]+=siz[v];
        if(siz[v]>siz[son[u]]) son[u]=v;
    }
}

void dfs2(int u,int t) {
    top[u]=t;
    dfn[u]=++tim;
    if(!son[u]) return ;
    dfs2(son[u],t);
    for(auto v:e[u]) {
        if(v==f[u] || v==son[u]) continue;
        dfs2(v,v);
    }
}

int getlca(int x,int y) {
    while(top[x]!=top[y]) {
        if(dep[top[x]]<dep[top[y]]) swap(x,y);
        x=f[top[x]];
    }
    if(dep[x]>dep[y]) swap(x,y);
    return x;
}

vector<int>new_e[N];
int st[N],Top;

void add(int u,int v){
    new_e[u].pb(v);
}

void insert(int x) {
    if(Top <= 1) {
        st[++Top] = x;
        return ;
    }
    int lca = getlca(x, st[Top]);
    if(lca == st[Top]) {
        st[++Top] = x;
        return;
    }
    while(Top > 1 && dfn[st[Top - 1]] >= dfn[lca]) {
        add(st[Top - 1], st[Top]);
        --Top;
    }
    if(lca != st[Top]) add(lca, st[Top]), st[Top] = lca;
    st[++Top] = x;
}

ll a[N]; 
int n,s;

ll dfs(int u){
    if(sz(new_e[u])==0) return a[u];
    ll ans=0;
    for(auto v:new_e[u]){
        ll sum=dfs(v);
        if(sum){
            ans+=max(1LL,sum-(dep[v]-dep[u]));
        }
    }
    return ans;
}

void del(int u){
    for(auto v:new_e[u]) del(v);
    new_e[u].clear();
}

ll solve(int pos){
    if(pos>1) tmp[pos].pb(s);
    sort(all(tmp[pos]),[](int x,int y){
        return dfn[x]<dfn[y];
    });
    Top=0;
    for(auto i:tmp[pos]) insert(i);
    while(Top>1){
        add(st[Top-1],st[Top]);
        --Top;
    }
    
    ll ans=dfs(s);
    del(s);
    return ans>1?ans-1:ans;
}


signed main() {

#ifdef xiaofan
    freopen("1.in","r",stdin);
    freopen("1.out","w",stdout);
#endif
    
    
    read(n,s);
    for(int i=1;i<=n;i++) read(a[i]);
    for(int i=1;i<n;i++){
        int u,v;
        read(u,v);
        e[u].pb(v);
        e[v].pb(u);
    }
    dfs1(s,0);
    dfs2(s,s);
    int maxdep=0;
    for(int i=1;i<=n;i++) maxdep=max(maxdep,dep[i]);
    ll ans=0;
    for(int i=1;i<=maxdep;i++) ans+=solve(i);
    printf("%lld\n",ans);



    return 0;
}


文章作者: 小凡
版权声明: 本博客所有文章除特別声明外,均采用 CC BY 4.0 许可协议。转载请注明来源 小凡 !
评论
  目录
隐藏
{% if theme.sakura.enable %}{% endif %}