【自用复习】牛客每日一题2026.4.18 最大稳定数值


参考@此在Dasein的题解

  • l o w s m [ u ] ≤ a [ u ] ≤ u p s m [ u ] lowsm[u]\leq a[u]\leq upsm[u] lowsm[u]≤a[u]≤upsm[u],u一定是支撑点
  • a [ u ] ≤ l o w s m [ u ] , a [ u ] ≤ u p s m [ u ] a[u]\leq lowsm[u] ,a[u]\leq upsm[u] a[u]≤lowsm[u],a[u]≤upsm[u],u可能是支撑点,去掉u子树中的一个点v及其子树,可能让 l o w s m [ u ] − ( l o w s m [ v ] + a [ v ] ) ≤ a [ u ] lowsm[u]-(lowsm[v]+a[v])\leq a[u] lowsm[u]−(lowsm[v]+a[v])≤a[u],让这个u成为支撑点,但是点v及其子树中的支撑点被去掉
cpp 复制代码
const int N = 1e6+10,M=1e5;
const double PI=acos(-1);
const long long mod =998244353, inf = 1e18+10 ;
int tr[N];
inline int lowbit(int x){
    return (x&(-x));
}
void add(int x,int v){
    for(int i=x;i<N;i+=lowbit(i))tr[i]+=v;
}
int query(int x){
    int res=0;
    for(int i=x;i>0;i-=lowbit(i))res+=tr[i];
    return res;
}
vector<int>g[N];
int a[N],upsm[N],lowsm[N],loss[N];// loss:子树中支撑节点的个数
map<int,int>disc;// 离散化
vector<int>val;// 记录离散化前的值
bool fg[N],pfg[N];// 是否是支撑节点 是否可能是支撑节点perhaps_fg
int cnt,mx;// cnt是都不删的答案
int getdisc(int v){
    auto num=upper_bound(val.begin(),val.end(),v);
    if(num==val.begin())return 0;
    else return disc[*prev(num)];
}
void dfs(int now,int fa){
    upsm[now]+=a[fa]+upsm[fa];
    for(auto x:g[now]){
        dfs(x,now);
        lowsm[now]+=a[x]+lowsm[x];
        loss[now]+=loss[x];
    }
    if(upsm[now]>=a[now]){
        if(lowsm[now]<=a[now])cnt++,fg[now]=1,loss[now]++;
        else{
            pfg[now]=1;
            disc[lowsm[now]-a[now]]=0;// 记录lowsm[now]-a[now]并之后离散化(因为值到1e9开不了那么大的数组)
        }
    }
}
// 处理删边部分得到答案
void dfs2(int now,int fa){
    /*
    删去fa------now
    now子树的fg都删掉 失去loss[now]
    祖先节点中的pfg(如点u)能成为fg的条件 lowsm[u]-(lowsm[now]+a[now])<=a[u]
    即lowsm[u]-a[u]<=(lowsm[now]+a[now]) 数值大小上做前缀和 
    用树状数组维护 把祖先中pfg的lowsm[u]-a[u](离散化后)放入
    */ 
//    cout<<now<<' '<<fa<<endl;
    int gain=query(getdisc(lowsm[now]+a[now]));
    mx=max(mx,cnt-loss[now]+gain);

    if(pfg[now])add(getdisc(lowsm[now]-a[now]),1);
    for(auto x:g[now]){
        dfs2(x,now);
    }
    if(pfg[now])add(getdisc(lowsm[now]-a[now]),-1);
}
void solve(){
    int n;cin>>n;
    forr(i,1,n){
        cin>>a[i];
    }
    forr(i,1,n){
        int fa;cin>>fa;
        g[fa].push_back(i);
    }
    dfs(1,0);
    mx=cnt;
    // 得到了pfg的子树和 离散化
    int id=0;
    for(auto &[k,v]:disc){
        v=++id;
        val.push_back(k);
    }
    // forr(i,1,n)cout<<upsm[i]<<' ';cout<<endl;
    // forr(i,1,n)cout<<lowsm[i]<<' ';cout<<endl;
    // cout<<cnt<<endl;
    // forr(i,1,n)cout<<fg[i]<<' ';cout<<endl;
    // forr(i,1,n)cout<<pfg[i]<<' ';cout<<endl;
    dfs2(1,0);
    cout<<mx<<endl;
    
}
相关推荐
笨笨饿2 小时前
66_C语言与微控制器底层开发
linux·c语言·网络·数据结构·算法·机器人·个人开发
yong99902 小时前
双连杆机械臂阻抗控制仿真
算法·仿真
李可以量化3 小时前
Python之如何做出交易日历(上)
人工智能·算法·qmt·量化 qmt ptrade
会编程的土豆3 小时前
【数据结构与算法】新二叉树
数据结构·算法·二叉树
jerryinwuhan3 小时前
基于结构可控性的给水管网传感器布点选择算法
数据库·算法
南宫萧幕3 小时前
车辆控制基础:从 EKF 状态估计到非线性 MPC 轨迹跟踪的闭环实现
算法·matlab·汽车·控制·pid
点PY3 小时前
《高质量的实时直接体绘制算法》论文精读
算法
X journey3 小时前
机器学习进阶(24):主成分分析PCA
人工智能·算法·机器学习
MediaTea3 小时前
AI 术语通俗词典:精确率(分类)
人工智能·算法·机器学习·分类·数据挖掘