[HNOI2016] 网络

59 篇文章 0 订阅
7 篇文章 0 订阅

题目描述:

维护数据结构
资瓷如下操作:
0 在 u - v 路径上 插入权值 Val
1 取消 第 i 次 0 操作
2 阻断 X 的通讯 即 经过 X 的权值作废 求可行的全局最大权值 若无输出 -1

题目分析:

先剖一剖 qwq
观察到 每个点会影响经过自己的路径权值
那么我们不妨每次插入的时候 对所有不属于 路径 u - v 的点插入这个权值
在阻断通讯时 我们只需要查询阻断节点上的权值最大值即可
如果没有2操作就是个那么就结束了
但是带了2操作 我们就需要在每个线段树节点上维护很多值,资瓷删除 ,动态查询最大值
考虑堆结构
我们维护两个堆,插入堆和删除堆,当两个堆的顶元素相同时,证明这个元素删除过了,就pop,
知道不是相同元素为止
要加一个标记永久化,要不删除操作的确不好维护
树剖取反大概也是个log吧
O(m*logn^3)

题目链接:

Luogu 3250
COGS 2215
BZOJ 4538
BZOJ MLE 我也不知道为啥

Ac 代码:

#include <cstdio>
#include <iostream>
#include <queue>
#include <algorithm>
const int maxm=1e5+100;
struct treeheap{
    std::priority_queue <int> a,b;
    inline void insert(int x){a.push(x);}
    inline void del(int x){b.push(x);}
    inline int asktop()
    {
        while(!b.empty()&&a.top()==b.top()) a.pop(),b.pop();
        if(a.empty()) return -1;
        else return a.top();
    }
}st[maxm<<2];
struct asknode{
    int u,v,val;
}q[maxm<<2];
struct node{
    int u,v;
}a[maxm];
inline bool comp(node x,node y)
{
    return x.u<y.u;
}
int head[maxm],to[maxm<<1],net[maxm<<1],cnt;
int fa[maxm],deep[maxm],siz[maxm],top[maxm],son[maxm],id[maxm],tot;
int n,m;
inline void addedge(int u,int v)
{
    cnt++;
    to[cnt]=v,net[cnt]=head[u],head[u]=cnt;
}
void dfs1(int now,int fax,int dep)
{
    fa[now]=fax,deep[now]=dep,siz[now]=1;
    for(int i=head[now];i;i=net[i])
    if(to[i]!=fax)
    {
        dfs1(to[i],now,dep+1);
        siz[now]+=siz[to[i]];
        if(siz[to[i]]>siz[son[now]]) son[now]=to[i];
    }
}
void dfs2(int now,int topx)
{
    id[now]=++tot;
    top[now]=topx;
    if(son[now]) dfs2(son[now],topx);
    for(int i=head[now];i;i=net[i])
    if(!id[to[i]]) dfs2(to[i],to[i]);
}
void seginsert(int o,int l,int r,int ql,int qr,int v,int f)
{
    if(ql<=l&&r<=qr)
    {
        if(f) st[o].insert(v);
        else st[o].del(v);
        return;
    }
    int mid=(l+r)>>1;
    if(ql<=mid) seginsert((o<<1),l,mid,ql,qr,v,f);
    if(qr>mid) seginsert((o<<1)|1,mid+1,r,ql,qr,v,f);
}
int segask(int o,int l,int r,int ind)
{
    if(l>=r) return st[o].asktop();
    int mid=(l+r)>>1;
    if(ind<=mid) return std::max(segask((o<<1),l,mid,ind),st[o].asktop());
    else return std::max(segask((o<<1)|1,mid+1,r,ind),st[o].asktop());
}
void treeins(int u,int v,int val,int f)
{
    int topx=0;
    while(top[u]!=top[v])
    {
        if(deep[top[u]]<deep[top[v]]) std::swap(u,v);
        a[++topx].u=id[top[u]],a[topx].v=id[u];
        u=fa[top[u]];
    } 
    if(deep[u]<deep[v]) std::swap(u,v);
    a[++topx].u=id[v],a[topx].v=id[u];
    std::sort(a+1,a+topx+1,comp);
    if(a[1].u>1) seginsert(1,1,n,1,a[1].u-1,val,f);
    if(a[topx].v<n) seginsert(1,1,n,a[topx].v+1,n,val,f);
    for(int i=2;i<=topx;i++)
     if(a[i].u-1>=a[i-1].v+1) seginsert(1,1,n,a[i-1].v+1,a[i].u-1,val,f);
}
int main()
{
    //freopen("network_tenderRun.in","r",stdin);
    //freopen("network_tenderRun.out","w",stdout);
    scanf("%d%d",&n,&m);
    for(int i=1;i<n;i++)
    {
        int u,v;
        scanf("%d%d",&u,&v);
        addedge(u,v);
        addedge(v,u);
    }
    dfs1(1,1,1);
    dfs2(1,1);
    for(int i=1;i<=m;i++)
    {
        int opt,u,v,val;
        scanf("%d",&opt);
        if(opt==0)
        {
            scanf("%d%d%d",&u,&v,&val);
            q[i].u=u,q[i].v=v,q[i].val=val;
            treeins(u,v,val,1);
        }
        if(opt==1)
        {
            scanf("%d",&u);
            treeins(q[u].u,q[u].v,q[u].val,0);
        }
        if(opt==2)
        {
            scanf("%d",&u);
            printf("%d\n",segask(1,1,n,id[u]));
        }
    }
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值