欧拉序+差分代替树链剖分维护树上信息

题目链接 [HAOI2015]树上操作

---

  如果用树链剖分来写，这当然是一道简单的维护问题，但是今天的主题是欧拉序（莫名给题加难度了）。

  先看一张图，我们列写一下本题的欧拉序的列写方法，这道题的欧拉序的跑法是进出各只记录一次。

对应欧拉序：A-B-D-D-E-E-F-F-B-C-C-A

  所以说，欧拉序是很活的，对应不同的问题，有时候需要用不一样的方法来跑。

  那么，我们也很容易的确定每个点的进欧拉序和出欧拉序。这个很有用。

  现在，就拿刚才这张图举例，如果我们用差分的思想，给进点“+1”，给出点“-1”，那么，譬如说，我现在要修改“B--E”这段枝上的点，那么实际上，可以集体的给“B-D-D-E”和“E-F-F-B”这两段区间都“+val”，“val”是要增加的值。因为不在这条链上的点都会根据“进一出一”这样的差分原则，最终抵消掉了，因为我们查询的方式，所以最终只会留下前段的值，而减去后段的值是为了查询更大子树的时候不影响外面的答案。

刚才的这段话在代码上的表述：

update(1, 1, _UP, L[x], L[y], val); //update x to y(x if father of y)
update(1, 1, _UP, R[y], R[x], val);

这样基础的操作有了之后，我们仍然未说明查询的操作。

譬如说，我们现在要查询某一条链上的权值和了。

那么，我们实际上只看他们的进欧拉序就可以了，因为出欧拉序是减去的部分，我们只留下要的部分，不要的部分不管，代码上是这样表述：

query(1, 1, _UP, L[x], L[y]); //x if father of y

这么写是正确的原因还是在于非这条链但是却在欧拉序中出现的点，他们的出欧拉点也一定在查询区间上。这也是差分的神奇之处。

对了，注释中的father都是指的祖先，我们用深度判一下就可以区别的了。

那么，问题来了，如果两个点不在由根结点下来的一条链上，而是在两条链上，中间由LCA结点分割呢？这都提到LCA了，自然求一个LCA，分两段不就好了嘛？

LCA怎么求？RMQ维护一下（码量++）。

还有个问题，对子树全体更新呢？

因为查询的链是查到它对应的“进欧拉序”，所以我们对x的子树全体都增加a的话，实际上，可以看成：

update(1, 1, _UP, L[x], R[x], a);   //update x is root

 

给出例题的代码：

#include <iostream>
#include <cstdio>
#include <cmath>
#include <string>
#include <cstring>
#include <algorithm>
#include <limits>
#include <vector>
#include <stack>
#include <queue>
#include <set>
#include <map>
#include <bitset>
//#include <unordered_map>
//#include <unordered_set>
#define lowbit(x) ( x&(-x) )
#define pi 3.141592653589793
#define e 2.718281828459045
#define INF 0x3f3f3f3f
#define eps 1e-8
#define HalF (l + r)>>1
#define lsn rt<<1
#define rsn rt<<1|1
#define Lson lsn, l, mid
#define Rson rsn, mid+1, r
#define QL Lson, ql, qr
#define QR Rson, ql, qr
#define myself rt, l, r
#define MP(a, b) make_pair(a, b)
using namespace std;
typedef unsigned long long ull;
typedef unsigned int uit;
typedef long long ll;
const int maxN = 2e5 + 7;
int N, Q, head[maxN], cnt, _UP;
struct Eddge
{
    int nex, to;
    Eddge(int a=-1, int b=0):nex(a), to(b) {}
}edge[maxN << 1];
inline void addEddge(int u, int v)
{
    edge[cnt] = Eddge(head[u], v);
    head[u] = cnt++;
}
inline void _add(int u, int v) { addEddge(u, v); addEddge(v, u); }
ll a[maxN], work[maxN], sum[maxN];
int L[maxN], R[maxN], tot;
void Euler_dfs(int u, int fa)
{
    L[u] = ++tot; work[tot] = a[u]; sum[tot] = 1;
    for(int i=head[u], v; ~i; i=edge[i].nex)
    {
        v = edge[i].to;
        if(v == fa) continue;
        Euler_dfs(v, u);
    }
    R[u] = ++tot; work[tot] = a[u]; sum[tot] = -1;
}
ll lazy[maxN << 2], tree[maxN << 2];
inline void pushup(int rt) { tree[rt] = tree[lsn] + tree[rsn]; }
inline void pushdown(int rt, int l, int r)
{
    if(lazy[rt])
    {
        int mid = HalF;
        tree[lsn] += (sum[mid] - sum[l - 1]) * lazy[rt]; tree[rsn] += (sum[r] - sum[mid]) * lazy[rt];
        lazy[lsn] += lazy[rt]; lazy[rsn] += lazy[rt];
        lazy[rt] = 0;
    }
}
void build(int rt, int l, int r)
{
    lazy[rt] = 0;
    if(l == r)
    {
        tree[rt] = (sum[l] - sum[l - 1]) * work[l];
        return;
    }
    int mid = HalF;
    build(Lson); build(Rson);
    pushup(rt);
}
void update(int rt, int l, int r, int ql, int qr, ll val)
{
    if(ql <= l && qr >= r)
    {
        tree[rt] += (sum[r] - sum[l - 1]) * val;
        lazy[rt] += val;
        return;
    }
    pushdown(myself);
    int mid = HalF;
    if(qr <= mid) update(QL, val);
    else if(ql > mid) update(QR, val);
    else { update(QL, val); update(QR, val); }
    pushup(rt);
}
ll query(int rt, int l, int r, int ql, int qr)
{
    if(ql <= l && qr >= r) return tree[rt];
    pushdown(myself);
    int mid = HalF;
    if(qr <= mid) return query(QL);
    else if(ql > mid) return query(QR);
    else return query(QL) + query(QR);
}
inline void init()
{
    cnt = tot = 0;
    for(int i=1; i<=N; i++) head[i] = -1;
}
int main()
{
    scanf("%d%d", &N, &Q);
    init();
    for(int i=1; i<=N; i++) scanf("%lld", &a[i]);
    for(int i=1, u, v; i<N; i++)
    {
        scanf("%d%d", &u, &v);
        _add(u, v);
    }
    Euler_dfs(1, 0);
    _UP = N << 1;
    for(int i=1; i<=_UP; i++) sum[i] += sum[i - 1];
    build(1, 1, _UP);
    int op, x, a;
    while(Q--)
    {
        scanf("%d%d", &op, &x);
        if(op == 1)
        {
            scanf("%d", &a);
            update(1, 1, _UP, L[x], L[x], a);
            update(1, 1, _UP, R[x], R[x], a);
        }
        else if(op == 2)
        {
            scanf("%d", &a);
            update(1, 1, _UP, L[x], R[x], a);   //update x is root
        }
        else
        {
            printf("%lld\n", query(1, 1, _UP, 1, L[x]));
        }
    }
    return 0;
}
/*
5 5
1 2 3 4 5
1 2
1 4
2 3
2 5
3 3
1 2 1
3 5
2 1 2
3 3
*/