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


参考@此在Dasein的题解

  • l o w s m u ≤ a u ≤ u p s m u lowsmu\leq au\leq upsmu lowsmu≤au≤upsmu,u一定是支撑点
  • a u ≤ l o w s m u , a u ≤ u p s m u au\leq lowsmu ,au\leq upsmu au≤lowsmu,au≤upsmu,u可能是支撑点,去掉u子树中的一个点v及其子树,可能让 l o w s m u − ( l o w s m v + a v ) ≤ a u lowsmu-(lowsmv+av)\leq au lowsmu−(lowsmv+av)≤au,让这个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;
    
}
相关推荐
KaMeidebaby13 小时前
卡梅德生物技术快报|PD1 单克隆抗体定制配套 N 糖全谱质控开发
前端·人工智能·算法·数据挖掘·数据分析
8Qi814 小时前
LeetCode 235. 二叉搜索树的最近公共祖先(LCA)
算法·leetcode·二叉树·递归·二叉搜索树·lca·迭代
bIo7lyA8v14 小时前
算法稳定性分析中的随机扰动建模的技术8
算法
科研online14 小时前
基于多源数据和XGBoost-SHAP分析中国大陆绿地碳汇空间变异影响因素的非线性相关性与尺度差异
算法·学习方法
Cthy_hy15 小时前
拓扑排序超详解:原理 + Kahn 贪心算法
python·算法·贪心算法
三品吉他手会点灯15 小时前
C语言学习笔记 - 43.运算符与表达式 - 运算符1 - 运算符的分类和简单介绍
c语言·笔记·学习·算法
VkN2X2X4b15 小时前
算法复杂度的实验验证与误差分析的技术8
算法
其利天下技术15 小时前
风扇灯无刷电机自适应算法实战指南
算法·cocos2d·无刷电机自适应算法·bldc驱动自适应算法·其利无刷电机驱动算法
8Qi816 小时前
LeetCode 494:目标和(Target Sum)—— 题解 ✅
算法·leetcode·职场和发展·动态规划·01背包
hujinyuan2016016 小时前
2026年3月 中国电子学会青少年软件编程(Python)三级考试试卷 真题及答案
java·python·算法