【自用复习】牛客每日一题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;
    
}
相关推荐
aramae4 小时前
MySQL复合查询(8)
java·c语言·开发语言·后端·算法
荆棘鸟智能4 小时前
城市感知设备怎么统一接入?从多协议网关到设备模型的中间件架构设计
人工智能·算法·边缘计算
INGNIGHT8 小时前
270 · 电话号码的字母组合II(Trie)
linux·算法
2401_839080549 小时前
C++常见八股
数据结构·c++·算法
青山是哪个青山9 小时前
LeetCode 188:买卖股票的最佳时机 IV
算法
mikuyyds9 小时前
geo-toolbox 插件算法解析:RUSLE 与 MUSLE 的工程化落地
算法·rust·gis
a187927218319 小时前
【算法】树(二):队列与祖先——BFS 骨架三配件、短路透传与合流
算法·leetcode·二叉树··bfs·算法讲解·树遍历
syagain_zsx10 小时前
算法基础篇 · 03 枚举(C++ 题解)
c++·算法·二进制·枚举
weixin_3077791310 小时前
C++代码实现MATLAB中的ode23tb函数功能
开发语言·c++·算法·matlab
我是章汕呐11 小时前
地级市能源消耗量及消耗强度数据【2006-2023年】平衡面板
人工智能·经验分享·算法·回归