【自用复习】牛客每日一题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;
    
}
相关推荐
罗西的思考11 小时前
机器人 / 强化学习】HIL-SERL:人类在环驱动的具身智能进化框架
人工智能·算法·机器学习
美团技术团队14 小时前
LongCat 开源 VitaBench 2.0:长期动态智能体基准新标杆
人工智能·算法
To_OC1 天前
LC 207 课程表:刚学图论那会儿,我连这是拓扑排序都没看出来
javascript·算法·leetcode
To_OC1 天前
LC 208 实现 Trie 前缀树:曾被名字劝退,写完发现是送分题
javascript·算法·leetcode
BadBadBad__AK1 天前
线段树维护区间 k 次方和
c++·数学·算法·stl
_清歌2 天前
DSpark 深度解读:DeepSeek-V4 如何用「半自回归」把推理速度提升 85%
算法
统计实现局2 天前
SVD 的三步走:双对角化、Givens 收敛、排序
算法
躬行见万象2 天前
《VLA 系列》UniLab 强化训练 | G1 机器人 |复现
算法
统计实现局2 天前
对称不定分解(Bunch-Kaufman):为什么 Cholesky 不够用
算法