[CTSC2008]网络管理(整体二分+树剖+树状数组)

一道经典的带修改树链第 \(k\) 大的问题。

我只想出三个 \(\log\) 的解法。。。

整体二分+树剖+树状数组。

那不是暴力随便踩的吗???

不过跑得挺快的。

\(Code\ Below:\)

// luogu-judger-enable-o2
#include <bits/stdc++.h>
#define lowbit(x) ((x)&(-(x)))
using namespace std;
const int maxn=80000+10;
const int lim=1e8;
int n,m,a[maxn],c[maxn],ans[maxn],cnt,num;
int top[maxn],dep[maxn],siz[maxn],son[maxn],fa[maxn],id[maxn],tim;
int head[maxn],to[maxn<<1],nxt[maxn<<1],tot;

struct Query{
    int op,x,y,k,id;
}q[maxn*3],q1[maxn*3],q2[maxn*3];

inline int read(){
    register int x=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
    while(isdigit(ch)){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}
    return (f==1)?x:-x;
}

inline void update(int x,int y){
    for(;x<=n;x+=lowbit(x)) c[x]+=y;
}

inline int sum(int x){
    int ans=0;
    for(;x;x-=lowbit(x)) ans+=c[x];
    return ans;
}

inline void clear(int x){
    for(;x<=n;x+=lowbit(x)){
        if(c[x]) c[x]=0;
        else break;
    }
}

inline void addedge(int x,int y){
    to[++tot]=y;
    nxt[tot]=head[x];
    head[x]=tot;
}

void dfs1(int x,int f){
    siz[x]=1;fa[x]=f;
    dep[x]=dep[f]+1;
    int maxson=-1;
    for(int i=head[x],y;i;i=nxt[i]){
        y=to[i];
        if(y==f) continue;
        dfs1(y,x);siz[x]+=siz[y];
        if(siz[y]>maxson) maxson=siz[y],son[x]=y;
    }
}

void dfs2(int x,int topf){
    top[x]=topf;id[x]=++tim;
    if(son[x]) dfs2(son[x],topf);
    for(int i=head[x],y;i;i=nxt[i]){
        y=to[i];
        if(y==fa[x]||y==son[x]) continue;
        dfs2(y,y);
    }
}

inline int query(int x,int y){
    int ans=0;
    while(top[x]!=top[y]){
        if(dep[top[x]]<dep[top[y]]) swap(x,y);
        ans+=sum(id[x])-sum(id[top[x]]-1);
        x=fa[top[x]];
    }
    if(dep[x]>dep[y]) swap(x,y);
    ans+=sum(id[y])-sum(id[x]-1);
    return ans;
}

void solve(int L,int R,int l,int r){
    if(L>R) return ;
    if(l==r){
        for(int i=L;i<=R;i++){
            if(q[i].op==1) update(q[i].x,q[i].k);
            else {
                if(q[i].k<=query(q[i].x,q[i].y)) ans[q[i].id]=l;
                else ans[q[i].id]=-1;
            }
        }
        for(int i=L;i<=R;i++)
            if(q[i].op==1) clear(q[i].x);
        return ;
    }
    int mid=(l+r)>>1,cnt1=0,cnt2=0,val;
    for(int i=L;i<=R;i++){
        if(q[i].op==1){
            if(q[i].y<=mid) q1[++cnt1]=q[i];
            else update(q[i].x,q[i].k),q2[++cnt2]=q[i];
        }
        else {
            val=query(q[i].x,q[i].y);
            if(q[i].k>val) q[i].k-=val,q1[++cnt1]=q[i];
            else q2[++cnt2]=q[i];
        }
    }
    for(int i=L;i<=R;i++)
        if(q[i].op==1&&q[i].y>mid) clear(q[i].x);
    for(int i=1;i<=cnt1;i++) q[L+i-1]=q1[i];
    for(int i=1;i<=cnt2;i++) q[L+cnt1+i-1]=q2[i];
    solve(L,L+cnt1-1,l,mid);solve(L+cnt1,R,mid+1,r);
}

int main()
{
    n=read(),m=read();
    int k,x,y;
    for(int i=1;i<=n;i++) a[i]=read();
    for(int i=1;i<n;i++){
        x=read(),y=read();
        addedge(x,y);addedge(y,x);
    }
    dfs1(1,0);dfs2(1,1);
    for(int i=1;i<=n;i++) q[++cnt]=(Query){1,id[i],a[i],1,0};
    while(m--){
        k=read(),x=read(),y=read();
        if(k==0){
            q[++cnt]=(Query){1,id[x],a[x],-1,0};
            q[++cnt]=(Query){1,id[x],y,1,0};
            a[x]=y;
        }
        else {
            q[++cnt]=(Query){2,x,y,k,++num};
        }
    }
    solve(1,cnt,0,lim);
    for(int i=1;i<=num;i++){
        if(ans[i]==-1) printf("invalid request!\n");
        else printf("%d\n",ans[i]);
    }
    return 0;
}

原文地址:https://www.cnblogs.com/owencodeisking/p/10350766.html

时间: 2024-10-11 05:44:34

[CTSC2008]网络管理(整体二分+树剖+树状数组)的相关文章

[CTSC2008]网络管理 [整体二分]

题面 bzoj luogu 所有事件按时间排序 按值划分下放 把每一个修改 改成一个删除一个插入 对于一个查询 直接查这个段区间有多少合法点 如果查询值大于等于目标值 进入左区间 如果一个查询无解 那么它要求第k大无解 k > 路径长 用lca维护即可 #include <cstdio> #include <cstdlib> #include <algorithm> #include <cmath> #include <cstring> #

BZOJ 2738 矩阵乘法(整体二分+二维树状数组)

[题目链接] http://www.lydsy.com/JudgeOnline/problem.php?id=2738 [题目大意] 给出一个方格图,询问要求求出矩阵内第k小的元素 [题解] 我们对答案的大小进行整体二分,用二维树状数组维护二维区间和, 将超过数量的分治到左区间,不满足的分治到右区间即可. [代码] #include <cstdio> #include <algorithm> #include <cstring> using namespace std;

BZOJ 2738 矩阵乘法 整体二分+二维树状数组

题目大意:给定一个矩阵,多次求某个子矩阵中的第k小 分块解法见 http://blog.csdn.net/popoqqq/article/details/41356899 <论除最小割外题目解法从来与题目名称无关系列> 整体二分 Solve(x,y,S)表示处理答案在[x,y]区间内的询问集合S 预先将所有数按照大小排序 每次将[1,mid]之间的数插入树状数组 然后对于分治内部的每一个询问 去树状数组中查询相应子矩阵的数值 如果小于等于k就划分到左集合S1 否则划分到右集合S2 然后Solv

【清澄A1333】【整体二分+二维树状数组】矩阵乘法(梁盾)

试题来源 2012中国国家集训队命题答辩 问题描述 给你一个N*N的矩阵,不用算矩阵乘法,但是每次询问一个子矩形的第K小数. 输入格式 第一行两个数N,Q,表示矩阵大小和询问组数: 接下来N行N列一共N*N个数,表示这个矩阵: 再接下来Q行每行5个数描述一个询问:x1,y1,x2,y2,k表示找到以(x1,y1)为左上角.以(x2,y2)为右下角的子矩形中的第K小数. 输出格式 对于每组询问输出第K小的数. 样例输入 2 2 2 1 3 4 1 2 1 2 1 1 1 2 2 3 样例输出 1

【BZOJ1146】【CTSC2008】网络管理 [整体二分]

网络管理 Time Limit: 50 Sec  Memory Limit: 162 MB[Submit][Status][Discuss] Description M公司是一个非常庞大的跨国公司,在许多国家都设有它的下属分支机构或部门. 为了让分布在世界各地的N个部门之间协同工作,公司搭建了一个连接整个公司的通信网络. 该网络的结构由N个路由器和N-1条高速光缆组成. 每个部门都有一个专属的路由器,部门局域网内的所有机器都联向这个路由器,然后再通过这个通信子网与其他部门进行通信联络. 该网络结

Codeforces 165D Beard Graph 边权树剖+树状数组

Beard Graph 题意:给你一颗由n个结点组成的树,支持以下操作:1 i:将第i条边染成黑色(保证此时该边是白色),2 i:将第i条边染成白色(保证此时该边是黑色),3 a b:找出a,b两点之间只由黑边组成的最短路径. 思路:树链剖分+树状数组,把每条边的权值放到它指向的点中去,初始全为黑边,黑边权值为1,白边权值为-inf,黑边变白边,将点权增加-inf,白边变黑边点权增加inf,因为不可能白边边白边,所以可以这样做,查询的时候要减去2个点的最近公共祖先的点权,最近公共祖先可通过树剖的

51nod1307(暴力树剖/二分&amp;dfs/并查集)

题目链接: http://www.51nod.com/onlineJudge/questionCode.html#!problemId=1307 题意: 中文题诶~ 思路: 解法1:暴力树剖 用一个数组 num[i] 维护编号为 i 的边当前最大能承受的重量. 在加边的过程中根据给出的父亲节点将当前边所在的链上所有边的num都减去当前加的边的重量, 注意当前边也要减自重. 那么当num首次出现负数时加的边号即位答案: 事实上这个算法的时间复杂度是O(n^2)的, 不过本题并没有出那种退化成单链的

树状数组差分求前缀和的前缀和

既然我们知道树状数组可以差分求前缀和 那么就可以对其进行 前缀和 求变量值 展开可得 $(k+1)*\sum_{i}^n c[i]-\sum_{i}^n i*c[i]$ 两个 树状数组可以搞定 顺便提一下DFS序  可以将树上问题转化为区间问题 对节点重新编号 并且有从属关系  从而 利用各种数据结构 例如 树剖  树状数组解决 原文地址:https://www.cnblogs.com/OIEREDSION/p/11380443.html

【BZOJ】1146: [CTSC2008]网络管理Network(树链剖分+线段树套平衡树+二分 / dfs序+树状数组+主席树)

第一种做法(时间太感人): 这题我真的逗了,调了一下午,疯狂造数据,始终找不到错. 后来发现自己sb了,更新那里没有打id,直接套上u了.我.... 调了一下午啊!一下午的时光啊!本来说好中午A掉去学习第二种做法,噗 好吧,现在第一种做法是hld+seg+bst+二分,常数巨大,log^4级别,目前只会这种. 树剖后仍然用线段树维护dfs序区间,然后在每个区间建一颗平衡树,我用treap,(这题找最大啊,,,囧,并且要注意,这里的rank是比他大的数量,so,我们在二分时判断要判断一个范围,即要