BZOJ 4127 Abs 树链剖分

题目大意:给定一棵树,每个点有一个整数权值(可以是负数),要求支持两种操作:

1.链上加

2.链上绝对值之和

由于加的数保证非负,因此一个负数变成一个正数最多有n次

树链剖分,在线段树中维护一下区间最大负数即可

不知道为何 写了两个线段树就TLE 把两个线段树合并成一个就7s过了

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define M 100100
using namespace std;

struct Segtree{
    Segtree *ls,*rs;
    long long sum,mark;
    int cnt;//区间中负数的个数
    pair<int,int> max_pos;
    void* operator new (size_t)
    {
        static Segtree mempool[M<<1],*C=mempool;
        return C++;
    }
    void Add(int x,int y,long long val)
    {
        sum+=(y-x+1-2*cnt)*val;
        mark+=val;
        if(max_pos.first!=(signed)0xefefefef)
            max_pos.first+=val;
    }
    void Push_Up(int x,int y)
    {
        sum=ls->sum+rs->sum;
        cnt=ls->cnt+rs->cnt;
        max_pos=max(ls->max_pos,rs->max_pos);
    }
    void Push_Down(int x,int y)
    {
        int mid=x+y>>1;
        if(mark)
        {
            ls->Add(x,mid,mark);
            rs->Add(mid+1,y,mark);
            mark=0;
        }
    }
    void Build_Tree(int x,int y,int a[])
    {
        int mid=x+y>>1;
        if(x==y)
        {
            if(a[mid]<0)
                max_pos=pair<int,int>(a[mid],mid),cnt=1;
            else
                max_pos=pair<int,int>(0xefefefef,0);
            sum=abs(a[mid]);
            return ;
        }
        (ls=new Segtree)->Build_Tree(x,mid,a);
        (rs=new Segtree)->Build_Tree(mid+1,y,a);
        Push_Up(x,y);
    }
    void Reverse(int x,int y,int pos)
    {
        int mid=x+y>>1;
        if(x==y)
        {
            cnt=0;
            sum=-sum;
            max_pos=pair<int,int>(0xefefefef,0);
            return ;
        }
        Push_Down(x,y);
        if(pos<=mid)
            ls->Reverse(x,mid,pos);
        else
            rs->Reverse(mid+1,y,pos);
        Push_Up(x,y);
    }
    void Add(int x,int y,int l,int r,long long val)
    {
        int mid=x+y>>1;
        if(x==l&&y==r)
        {
            Add(x,y,val);
            return ;
        }
        Push_Down(x,y);
        if(r<=mid)
            ls->Add(x,mid,l,r,val);
        else if(l>mid)
            rs->Add(mid+1,y,l,r,val);
        else
            ls->Add(x,mid,l,mid,val) , rs->Add(mid+1,y,mid+1,r,val) ;
        Push_Up(x,y);
    }
    long long Query(int x,int y,int l,int r)
    {
        int mid=x+y>>1;
        if(x==l&&y==r)
            return sum;
        Push_Down(x,y);
        if(r<=mid)
            return ls->Query(x,mid,l,r);
        if(l>mid)
            return rs->Query(mid+1,y,l,r);
        return ls->Query(x,mid,l,mid) + rs->Query(mid+1,y,mid+1,r) ;
    }
    pair<int,int> Get_Pos(int x,int y,int l,int r)
    {
        int mid=x+y>>1;
        if(x==l&&y==r)
            return max_pos;
        Push_Down(x,y);
        if(r<=mid)
            return ls->Get_Pos(x,mid,l,r);
        if(l>mid)
            return rs->Get_Pos(mid+1,y,l,r);
        return max( ls->Get_Pos(x,mid,l,mid) , rs->Get_Pos(mid+1,y,mid+1,r) );
    }
}*tree=new Segtree;

struct abcd{
    int to,next;
}table[M<<1];
int head[M],tot;

int n,m,a[M];
int fa[M],son[M],dpt[M],size[M];
int pos[M],top[M],_a[M];

void Add(int x,int y)
{
    table[++tot].to=y;
    table[tot].next=head[x];
    head[x]=tot;
}

void DFS1(int x)
{
    int i;
    dpt[x]=dpt[fa[x]]+1;
    size[x]=1;
    for(i=head[x];i;i=table[i].next)
        if(table[i].to!=fa[x])
        {
            fa[table[i].to]=x;
            DFS1(table[i].to);
            size[x]+=size[table[i].to];
            if(size[table[i].to]>size[son[x]])
                son[x]=table[i].to;
        }
}

void DFS2(int x)
{
    static int T;
    int i;
    _a[pos[x]=++T]=a[x];
    if(son[fa[x]]==x)
        top[x]=top[fa[x]];
    else
        top[x]=x;
    if(son[x])
        DFS2(son[x]);
    for(i=head[x];i;i=table[i].next)
        if(table[i].to!=fa[x]&&table[i].to!=son[x])
            DFS2(table[i].to);
}

void Add(int x,int y,int z)
{
    int fx=top[x],fy=top[y];
    while(fx!=fy)
    {
        if(dpt[fx]<dpt[fy])
            swap(fx,fy),swap(x,y);
        while(1)
        {
            pair<int,int> temp=tree->Get_Pos(1,n,pos[fx],pos[x]);
            if(temp.first+z>=0)
                tree->Reverse(1,n,temp.second);
            else
                break;
        }
        tree->Add(1,n,pos[fx],pos[x],z);
        x=fa[fx];fx=top[x];
    }
    if(dpt[x]<dpt[y])
        swap(x,y);
    while(1)
    {
        pair<int,int> temp=tree->Get_Pos(1,n,pos[y],pos[x]);
        if(temp.first+z>=0)
            tree->Reverse(1,n,temp.second);
        else
            break;
    }
    tree->Add(1,n,pos[y],pos[x],z);
}

long long Query(int x,int y)
{
    int fx=top[x],fy=top[y];
    long long re=0;
    while(fx!=fy)
    {
        if(dpt[fx]<dpt[fy])
            swap(fx,fy),swap(x,y);
        re+=tree->Query(1,n,pos[fx],pos[x]);
        x=fa[fx];fx=top[x];
    }
    if(dpt[x]<dpt[y])
        swap(x,y);
    return re+tree->Query(1,n,pos[y],pos[x]);
}

int main()
{
    int i,p,x,y,z;
    cin>>n>>m;
    for(i=1;i<=n;i++)
        scanf("%d",&a[i]);
    for(i=1;i<n;i++)
    {
        scanf("%d%d",&x,&y);
        Add(x,y);Add(y,x);
    }
    DFS1(1);DFS2(1);
    tree->Build_Tree(1,n,_a);
    for(i=1;i<=m;i++)
    {
        scanf("%d%d%d",&p,&x,&y);
        if(p==1)
        {
            scanf("%d",&z);
            if(z) Add(x,y,z);
        }
        else
        {
            #ifdef PoPoQQQ
                printf("%I64d\n",Query(x,y));
            #else
                printf("%lld\n",Query(x,y));
            #endif
        }
    }
    return 0;
}
时间: 2024-11-08 22:39:58

BZOJ 4127 Abs 树链剖分的相关文章

BZOJ 3083 遥远的国度 树链剖分

3083: 遥远的国度 Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 797  Solved: 181[Submit][Status] Description 描述 zcwwzdjn在追杀十分sb的zhx,而zhx逃入了一个遥远的国度.当zcwwzdjn准备进入遥远的国度继续追杀时,守护神RapiD阻拦了zcwwzdjn的去路,他需要zcwwzdjn完成任务后才能进入遥远的国度继续追杀. 问题是这样的:遥远的国度有n个城市,这些城市之间由一些路连

BZOJ 2243 染色 | 树链剖分模板题进阶版

BZOJ 2243 染色 | 树链剖分模板题进阶版 这道题呢~就是个带区间修改的树链剖分~ 如何区间修改?跟树链剖分的区间询问一个道理,再加上线段树的区间修改就好了. 这道题要注意的是,无论是线段树上还是原树上,把两个区间的信息合并的时候,要注意中间相邻两个颜色是否相同. 这代码好长啊啊啊啊 幸好一次过了不然我估计永远也De不出来 #include <cstdio> #include <cstring> #include <algorithm> using namesp

【BZOJ-4127】Abs 树链剖分 + 线段树 (有趣的姿势)

4127: Abs Time Limit: 40 Sec  Memory Limit: 256 MBSubmit: 381  Solved: 132[Submit][Status][Discuss] Description 给定一棵树,设计数据结构支持以下操作 1 u   v d 表示将路径 (u,v) 加d 2 u v 表示询问路径 (u,v) 上点权绝对值的和 Input 第一行两个整数n和m,表示结点个数和操作数 接下来一行n个整数a_i,表示点i的权值 接下来n-1行,每行两个整数u,v

BZOJ 2157 旅行(树链剖分码农题)

写了5KB,1发AC... 题意:给出一颗树,支持5种操作. 1.修改某条边的权值.2.将u到v的经过的边的权值取负.3.求u到v的经过的边的权值总和.4.求u到v的经过的边的权值最大值.5.求u到v经过的边的权值最小值. 基于边权的树链剖分,放在线段树上变成了区间维护问题了,线段树维护4个量min,max,sum,tag就可以了. # include <cstdio> # include <cstring> # include <cstdlib> # include

bzoj 2243 染色(树链剖分)

题外话 首先这是个挺裸的题,由于太久没写剖分导致调了好久,前天调了一下午,一直查不到错 昨晚在看春晚的时候突然灵机一动,发现合并的时候出了问题,开电脑把它A掉了= = 感觉自己也蛮拼的 Description 给定一棵有n个节点的无根树和m个操作,操作有2类: 1:将节点a到节点b路径上所有点都染成颜色c 2:询问节点a到节点b路径上的颜色段数量(连续相同颜色被认为是同一段),如"112221"由3段组成:"11"."222"和"1&q

BZOJ 3083 遥远的国度 树链剖分+线段树

有换根的树链剖分的裸题. 在换根的时候注意讨论. 注意数据范围要开unsigned int或longlong #include<iostream> #include<cstdio> #include<cstdlib> #include<cstring> #include<cmath> #include<ctime> #include<string> #include<iomanip> #include<a

BZOJ 2819 Nim 树链剖分

题目大意:两个小人在树上玩Nim游戏,问有没有必胜策略. 思路:尼姆游戏:如果所有石子的异或值为0就是必败局面.异或有如下性质:x ^ y ^ z = x ^ (y ^ z),所以就可以进行树链剖分了.题目中还好心提醒有30%的点会使dfs爆栈..第一次写不用dfs的树链剖分,扒的别人的代码,有些丑陋. CODE: #include <queue> #include <cstdio> #include <cstring> #include <iostream>

BZOJ 2819 Nim 树链剖分/DFS序+LCA+树状数组

题意:给定一棵树,每个节点是一堆石子,给定两种操作: 1.改变x号节点的石子数量 2.用从x到y的路径上的所有堆石子玩一次Nim游戏,询问是否有必胜策略 Nim游戏有必胜策略的充要条件是所有堆的石子数异或起来不为零 这题首先一看就是树链剖分 然后题目很善良地告诉我们深搜会爆栈 于是我们可以选择广搜版的树链剖分 BFS序从左到右是深搜,从右到左是回溯,一遍BFS就够 单点修改区间查询还可以套用ZKW线段树 不过这题其实不用这么麻烦 有更简单的办法 详见 http://dzy493941464.is

bzoj 5355 kdtree 树链剖分

https://www.lydsy.com/JudgeOnline/problem.php?id=5355 想在b站搜query on a tree系列不小心看到了这题 扑鼻而来的浓浓的OI风格的题面,6个操作,放在ACM界读完题就可以喷了(误 看到前三个操作...kdtree板子题,一维dfs序,一维dep,限制区间显然 后面三个操作...考虑树链剖分,这样点到根的路径就是几条不连续的链了 就把前面kdtree的一维dfs序换成树链剖分的方法得到的dfs序就行了 对于后面三个操作相当于把常写的