[树链剖分]JZOJ 5914 盟主的忧虑

news/2024/7/7 21:12:48

Description

    江湖由 N 个门派(2≤N≤100,000,编号从 1 到 N)组成,这些门派之间有 N-1 条小道将他们连接起来,每条道路都以“尺”为单位去计量,武林盟主发现任何两个门派都能够直接或者间接通过小道连接。
    虽然整个江湖是可以互相到达的,但是他担心有心怀不轨之徒破坏这个武林的安定,破坏小道,于是武林盟主又秘密地修建了 M 条密道(1≤M≤100,000),但每条小道距离都不超过10亿尺。
    果不其然,最近一个名叫“太吾”的组织意欲破坏武林的小道,请你帮盟主想想办法,如果门派 A 到门派 B 的直连小道被破坏,从 A 走到 B 的所有路径中,经过密道的距离最少是多少? 
 
 

Input

第一行数字 N M
接下来 N-1 行,每行两个整数 A B,表示 A-B 间有一条直连小道
接下来 M 行,每行三个数字 A B V,表示 A-B 间有一条代价为 V 的密道

Output

 输出 N-1 行,对应原输入的小道,每个小道被破坏后,最少需要经过多长的密道?如果不存在替换的道路,请输出-1
 
 

Sample Input

6 3
4 1
1 3
4 5
1 2
6 5
3 6 8
2 3 7
6 4 5

Sample Output

8
7
5
7
5
 

Data Constraint

30%数据:N<=300,M<=1000
50%数据:N<=1000,M<=1000
70%数据:N<=5000,M<=5000
对于另外15%的数据点:树是一条链
100%数据:N,M<=100,000

分析

这题显然密道会密道两端到LCA中经过的边的值

这不就是树链剖分嘛

写个区间修改走人

什么?并查集做法?不会

 

#include <iostream>
#include <cstdio>
#include <cstring>
#define rep(i,a,b) for (i=a;i<=b;i++)
using namespace std;
const int N=1e5+10;
struct Node {
    int sz,f,s,seg,top,dep;
}a[N];
int scnt,rev[N];
struct Edge {
    int u,v,nx;
}g[2*N];
int cnt,list[N];
int t[4*N],lz[4*N];
int mn;
int n,m;

void Add(int u,int v) {g[++cnt].u=u;g[cnt].v=v;g[cnt].nx=list[u];list[u]=cnt;}

void Dfs1(int x,int f) {
    a[x].f=f;a[x].dep=a[f].dep+1;a[x].sz=1;
    for (int i=list[x];i;i=g[i].nx)
    if (g[i].v!=f) {
       Dfs1(g[i].v,x);
       a[x].sz+=a[g[i].v].sz;
        if (a[g[i].v].sz>a[a[x].s].sz) a[x].s=g[i].v;
    }
}

void Dfs2(int x,int f) {
    int son=a[x].s;
    if (son) {
        a[son].seg=++scnt;
        rev[scnt]=son;
        a[son].top=a[x].top;
        Dfs2(son,x);
    }
    for (int i=list[x];i;i=g[i].nx)
    if (!a[g[i].v].top) {
        a[g[i].v].seg=++scnt;
        rev[scnt]=g[i].v;
        a[g[i].v].top=g[i].v;
        Dfs2(g[i].v,x);
    }
}

void Build(int x,int l,int r) {
    if (l==r) {
        t[x]=2147483647;
        return;
    }
    int mid=l+r>>1;
    Build(x<<1,l,mid);Build((x<<1)+1,mid+1,r);
    t[x]=2147483647;
}

void Pushdown(int x) {
    if (lz[x]) {
        t[x<<1]=lz[x<<1]=min(!lz[x<<1]?2147483647:lz[x<<1],lz[x]);
        t[(x<<1)+1]=lz[(x<<1)+1]=min(!lz[(x<<1)+1]?2147483647:lz[(x<<1)+1],lz[x]);
        lz[x]=0;
    }
}

void Change_In_Segment(int x,int l,int r,int val,int ll,int rr) {
    if (ll>rr) return;
    if (l>rr||ll>r) return;
    if (ll<=l&&r<=rr) {
        t[x]=min(t[x],val);
        lz[x]=min(!lz[x]?2147483647:lz[x],val);
        return;
    }
    Pushdown(x);
    int mid=l+r>>1;
    if (mid>=ll) Change_In_Segment(x<<1,l,mid,val,ll,rr);
    if (mid<rr) Change_In_Segment((x<<1)+1,mid+1,r,val,ll,rr);
    t[x]=min(t[x<<1],t[(x<<1)+1]);
}

void Query_In_Segment(int x,int l,int r,int ll,int rr) {
    if (ll>rr) return;
    if (r<ll||l>rr) return;
    if (ll<=l&&r<=rr) {
        mn=min(mn,t[x]);
        return;
    }
    Pushdown(x);
    int mid=l+r>>1;
    if (ll<=mid) Query_In_Segment(x<<1,l,mid,ll,rr);
    if (mid+1<=rr) Query_In_Segment((x<<1)+1,mid+1,r,ll,rr);
}

void Query_In_Tree(int x,int y) {
    int fx=a[x].top,fy=a[y].top;
    while (fx!=fy) {
        if (a[fx].dep<a[fy].dep) swap(x,y),swap(fx,fy);
        Query_In_Segment(1,1,scnt,a[fx].seg,a[x].seg);
        x=a[fx].f;fx=a[x].top;
    }
    if (a[x].dep>a[y].dep) swap(x,y);
    Query_In_Segment(1,1,scnt,a[x].seg+1,a[y].seg);
}

void Change_In_Tree(int x,int y,int val) {
    int fx=a[x].top,fy=a[y].top;
    while (fx!=fy) {
        if (a[fx].dep<a[fy].dep) swap(x,y),swap(fx,fy);
        Change_In_Segment(1,1,scnt,val,a[fx].seg,a[x].seg);
        x=a[fx].f;fx=a[x].top;
    }
    if (a[x].dep>a[y].dep) swap(x,y);
    Change_In_Segment(1,1,scnt,val,a[x].seg+1,a[y].seg);
}

int main() {
    freopen("worry.in","r",stdin);
    freopen("worry.out","w",stdout);
    int i;
    char s[20];
    scanf("%d%d",&n,&m);
    rep(i,1,n-1) {
        int u,v;
        scanf("%d%d",&u,&v);
        Add(u,v);Add(v,u);
    }
    Dfs1(1,0);
    scnt=1;
    a[1].seg=a[1].top=rev[1]=1;
    Dfs2(1,0);
    Build(1,1,scnt);
    for (int i=1;i<=m;i++) {
        int x,y,w;
        scanf("%d%d%d",&x,&y,&w);
        Change_In_Tree(x,y,w);
    }
    for (int i=1;i<=n-1;i++) {
        mn=2147483647;
        Query_In_Tree(g[i*2-1].u,g[i*2-1].v);
        printf("%d\n",mn==2147483647?-1:mn);
    }
}
View Code

 

转载于:https://www.cnblogs.com/mastervan/p/9818585.html


http://www.niftyadmin.cn/n/3778345.html

相关文章

什么是Java Marker Interface(标记接口)

2019独角兽企业重金招聘Python工程师标准>>> 先看看什么是标记接口&#xff1f;标记接口有时也叫标签接口&#xff08;Tag interface&#xff09;&#xff0c;即接口不包含任何方法。在Java里很容易找到标记接口的例子&#xff0c;比如JDK里的Serializable接口就是一…

JavaWeb项目架构之Kafka分布式日志队列

架构、分布式、日志队列&#xff0c;标题自己都看着唬人&#xff0c;其实就是一个日志收集的功能&#xff0c;只不过中间加了一个Kafka做消息队列罢了。 kafka介绍 Kafka是由Apache软件基金会开发的一个开源流处理平台&#xff0c;由Scala和Java编写。Kafka是一种高吞吐量的分布…

ionic 签名、打包

ionic cordova platform add androidionic cordova build android 【debug版本&#xff0c;无需签名】 ionic cordova build android --release 【发布版&#xff0c;需要签名&#xff08;要使用jarsigner签名必须用release版本&#xff09;】 ionic cordova build android --r…

Spring框架内容整理(三):AOP

2019独角兽企业重金招聘Python工程师标准>>> AOP面向切面 切面支持类 切面就是通知和切入点的组合&#xff0c;而切面是通过配置方式定义的&#xff0c;因此这定义切面前&#xff0c;我们需要定义切面支持类&#xff0c;切面支持类提供了通知实现 package cn.javass…

python并发编程之多进程理论

一、什么是进程 进程&#xff1a;正在进行的一个过程或者说一个任务。而负责执行任务则是cpu。 举例&#xff1a;单核多道&#xff0c;实现多个进程的并发执行你在一个时间段内有很多任务要做&#xff1a;写python程序&#xff0c;打游戏&#xff0c;逛淘宝  但你同一时刻只能…

python堆排序

堆排序介绍 堆排序&#xff0c;顾名思义&#xff0c;就是基于堆。因此先来介绍一下堆的概念。 堆分为最大堆和最小堆&#xff0c;其实就是完全二叉树。最大堆要求节点的元素都要大于其孩子&#xff0c;最小堆要求节点元素都小于其左右孩子&#xff0c;两者对左右孩子的大小关系…

005-统一沟通-部署-基础-环境

准备Active Directory 和域(031-Mail01)安装完成(031-Mail01)为邮箱服务器创建DNS记录&#xff08;011-DC01&#xff09;验证DNS邮箱服务器安装证书&#xff08;031-Mail01&#xff09;前端服务器验证林准备和域&#xff08;011-DC01&#xff09;前端服务器DNS验证&#xff08;…

Shiro整合springboot,freemaker,redis

为什么80%的码农都做不了架构师&#xff1f;>>> 一、导语 今天推荐给大家一个非常简单&#xff0c;实用的开源权限框架&#xff1a;Shiro&#xff0c;这也是Java官方推荐的权限框架。 尽管网上有很多类似的文章&#xff0c;但是都不够完善&#xff0c;尤其是涉及到…