【POJ3237】Tree(树链剖分+线段树(基于边权))

记录一个菜逼的成长。。

题目链接
题目大意:
有T组数据
给你n个点的树
有三种操作:
1.CHANGE i v :将第i条边的值改为v
2.NEGATE a b :将a点到b点路径上的边权取反
3.QUERY a b :查询a点到b点路径上最大的边权值

由于有取反操作,我们记录最大值的同时记录下最小值
如果要取反,则将最大最小值取反并交换即可
用线段树维护

注意懒惰标记的下放操作。
又加深了对懒惰标记的理解

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define fin freopen("D://in.txt","r",stdin)
#define fout freopen("D://out.txt","w",stdout)
#define cl(a,b) memset(a,b,sizeof(a))
#define lson t<<1,l,mid
#define rson t<<1|1,mid+1,r
const int maxn = 200000 + 10;
const int INF = 0x3f3f3f3f;
int eg[maxn][3];
int head[maxn],cnt;
struct Edge{
    int to,next;
}edge[maxn<<1];
void add(int u,int v)
{
    edge[cnt].to = v;
    edge[cnt].next = head[u];
    head[u] = cnt++;
}
int num;
int val[maxn],siz[maxn],dep[maxn],son[maxn];
int top[maxn],tid[maxn],pos[maxn],fa[maxn];
void init()
{
    cl(head,-1);cl(son,-1);
    cnt = num = 0;
}
///树链剖分部分
void dfs1(int u,int father,int d)
{
    dep[u] = d;
    fa[u] = father;
    siz[u] = 1;
    for(int i = head[u]; ~i; i = edge[i].next)
    {
        int v = edge[i].to;
        if(v != father)
        {
            dfs1(v,u,d+1);
            siz[u] += siz[v];
            if(son[u] == -1 || siz[v] > siz[son[u]])
                son[u] = v;
        }
    }
}

void dfs2(int u,int tp)
{
    top[u] = tp;
    tid[u] = ++num;
    pos[tid[u]] = u;
    if(son[u] == -1) return;
    dfs2(son[u],tp);
    for(int i = head[u]; ~i; i = edge[i].next)
    {
        int v = edge[i].to;
        if(v != son[u] && v != fa[u])
            dfs2(v,v);
    }
}
///线段树部分
struct Node{
    int l,r,mx,mn,lazy;
}node[maxn<<2];
void Swap(int &x,int &y){int t = x;x = -y;y = -t;}
void pushup(int t)
{
    node[t].mx = max(node[t<<1].mx, node[t<<1|1].mx);
    node[t].mn = min(node[t<<1].mn, node[t<<1|1].mn);
}
void pushdown(int t)
{
    if(node[t].lazy){
        int tmp = node[t].lazy;
        node[t<<1].lazy += tmp;
        node[t<<1|1].lazy += tmp;
        //注意这里,两个子区间是否取反是依靠当前区间的懒惰标记判断的
        if(tmp&1){
            Swap(node[t<<1].mx,node[t<<1].mn);
            Swap(node[t<<1|1].mx,node[t<<1|1].mn);
        }
        node[t].lazy = 0;
    }
}
void build(int t,int l,int r)
{
    node[t].l = l;
    node[t].r = r;
    node[t].lazy = 0;
    if(l == r){
        node[t].mx = node[t].mn = 0;
        return ;
    }
    int mid = (l + r) >> 1;
    build(lson);
    build(rson);
    pushup(t);
}
void update(int t,int ind,int v)
{
    if(node[t].l == node[t].r && node[t].l == ind){
        node[t].mx = v;
        node[t].mn = v;
        node[t].lazy = 0;
        return ;
    }
    pushdown(t);
    int mid = (node[t].l + node[t].r) >> 1;
    if(ind <= mid)update(t<<1,ind,v);
    else update(t<<1|1,ind,v);
    pushup(t);
}
void Negate(int t,int l,int r)
{
    if(node[t]. l >= l && node[t].r <= r){
        node[t].lazy ++;
        Swap(node[t].mx,node[t].mn);//这里要马上取反
        return ;
    }
    pushdown(t);
    int mid = (node[t].l + node[t].r) >> 1;
    if(l <= mid)Negate(t<<1,l,r);
    if(r > mid)Negate(t<<1|1,l,r);
    pushup(t);
}
int query(int t,int l,int r)
{
    if(node[t].l >= l && node[t].r <= r){
        return node[t].mx;
    }
    pushdown(t);
    int ret = -INF;
    int mid = (node[t].l + node[t].r) >> 1;
    if(l <= mid)ret = max(ret,query(t<<1,l,r));
    if(r > mid)ret = max(ret,query(t<<1|1,l,r));
    pushup(t);
    return ret;
}
int solveMax(int x,int y)
{
    int ret = -INF;
    while(top[x] != top[y]){
        if(dep[top[x]] < dep[top[y]])swap(x,y);
        ret = max(ret,query(1,tid[top[x]],tid[x]));
        x = fa[top[x]];
    }
    if(x == y)return ret;
    if(dep[x] > dep[y])swap(x,y);
    return max(ret,query(1,tid[son[x]],tid[y]));
}
void change(int x,int y)
{
    while(top[x] != top[y]){
        if(dep[top[x]] < dep[top[y]])swap(x,y);
        Negate(1,tid[top[x]],tid[x]);
        x = fa[top[x]];
    }
    if(x == y)return ;
    if(dep[x] > dep[y])swap(x,y);
    Negate(1,tid[son[x]],tid[y]);
}
int main()
{
    //fin,fout;
    char ope[10];
    int T,n;scanf("%d",&T);
    while(T--){
        scanf("%d",&n);
        init();
        for( int i = 1; i < n; i++ ){
            scanf("%d%d%d",&eg[i][0],&eg[i][1],&eg[i][2]);
            add(eg[i][0],eg[i][1]);
            add(eg[i][1],eg[i][0]);
        }
        dfs1(1,0,0);
        dfs2(1,1);
        build(1,1,n);
        for( int i = 1; i < n; i++ ){
            if(dep[eg[i][0]] < dep[eg[i][1]])swap(eg[i][0],eg[i][1]);
            update(1,tid[eg[i][0]],eg[i][2]);
        }
        int a,b;
        while(scanf("%s",ope) && ope[0] != 'D'){
            scanf("%d%d",&a,&b);
            if(ope[0] == 'Q'){
                printf("%d\n",solveMax(a,b));
            }
            if(ope[0] == 'C'){
                update(1,tid[eg[a][0]],b);
            }
            if(ope[0] == 'N'){
                change(a,b);
            }
        }
    }
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值