dfs序专题学习

dfs序是处理树上问题很重要的一个武器，主要能够解决对于一个点，它的子树上的一些信息的维护。
就比如那天百度之星round1A的1003题，就是dfs序+线段树维护每个点到0点的距离，然后对于每个点的更新，只需要更新它和它的子树上的点到0点的距离，查询的话就是它的子树上的最大值即可
我的dfs序开的空间就是n，因为只在入的地方时间戳++，出来的地方时间戳不变，线段树的每个节点应该是时间戳

void dfs(int u,int fa){
    p1[u]=++ti;
    dfsnum[ti]=u;
    for(int i=head[u];i!=-1;i=edge[i].next){
        int v=edge[i].v;
        if(v==fa) continue;
        dfs(v,u);
    }
    p2[u]=ti;//时间戳不变，空间为O(n)
}

废话不多说，赶紧来看例题吧

---

HDU 3887
题目传送门：http://acm.hdu.edu.cn/showproblem.php?pid=3887
问你对于每个节点，它的子树上标号比它小的点有多少个
子树的问题，dfs序可以很轻松的解决，因为点在它的子树上，所以在线段树中，必定在它的两个时间戳的区间之间，所以我们只需要从小到大考虑，它的区间里有多少个点已经放了，然后再把它放进去。很容易的解决了
代码：

#include <map>
#include <set>
#include <stack>
#include <queue>
#include <cmath>
#include <string>
#include <vector>
#include <cstdio>
#include <cctype>
#include <cstring>
#include <sstream>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#pragma comment(linker, "/STACK:102400000,102400000")

using namespace std;
#define   MAX           500005
#define   MAXN          6005
#define   maxnode       15
#define   sigma_size    30
#define   lson          l,m,rt<<1
#define   rson          m+1,r,rt<<1|1
#define   lrt           rt<<1
#define   rrt           rt<<1|1
#define   middle        int m=(r+l)>>1
#define   LL            long long
#define   ull           unsigned long long
#define   mem(x,v)      memset(x,v,sizeof(x))
#define   lowbit(x)     (x&-x)
#define   pii           pair<int,int>
#define   bits(a)       __builtin_popcount(a)
#define   mk            make_pair
#define   limit         10000

//const int    prime = 999983;
const int    INF   = 0x3f3f3f3f;
const LL     INFF  = 0x3f3f;
const double pi    = acos(-1.0);
//const double inf   = 1e18;
const double eps   = 1e-8;
const LL    mod    = 1e9+7;
const ull    mx    = 133333331;

/*****************************************************/
inline void RI(int &x) {
      char c;
      while((c=getchar())<'0' || c>'9');
      x=c-'0';
      while((c=getchar())>='0' && c<='9') x=(x<<3)+(x<<1)+c-'0';
 }
/*****************************************************/

struct Edge{
    int v,next;
}edge[MAX*2];
int head[MAX];
int tot;
int p1[MAX];
int p2[MAX];
int ti;
int sum[MAX<<2];

void init(){
    mem(head,-1);
    tot=0;ti=0;
}

void add_edge(int a,int b){
    edge[tot]=(Edge){b,head[a]};
    head[a]=tot++;
}

void dfs(int u,int fa){
    p1[u]=++ti;
    for(int i=head[u];i!=-1;i=edge[i].next){
        int v=edge[i].v;
        if(v==fa) continue;
        dfs(v,u);
    }
    p2[u]=ti;
}

void build(int l,int r,int rt){
    sum[rt]=0;
    if(l==r) return;
    middle;
    build(lson);
    build(rson);
}

void pushup(int rt){
    sum[rt]=sum[lrt]+sum[rrt];
}

void update(int l,int r,int rt,int pos,int d){
    if(l==r){
        sum[rt]+=d;
        return;
    }
    middle;
    if(pos<=m) update(lson,pos,d);
    else update(rson,pos,d);
    pushup(rt);
}

int query(int l,int r,int rt,int L,int R){
    if(L<=l&&r<=R) return sum[rt];
    middle;
    int ret=0;
    if(L<=m) ret+=query(lson,L,R);
    if(R>m) ret+=query(rson,L,R);
    return ret;
}

int main(){
    int n,p;
    while(cin>>n>>p&&n){
        init();
        for(int i=1;i<n;i++){
            int a,b;
            scanf("%d%d",&a,&b);
            add_edge(a,b);
            add_edge(b,a);
        }
        dfs(p,-1);
        build(1,n,1);
        for(int i=1;i<=n;i++){
            printf("%d",query(1,n,1,p1[i],p2[i]));
            if(i==n) printf("\n");
            else printf(" ");
            update(1,n,1,p1[i],1);
        }
    }
    return 0;
}

---

poj 3321
题目传送门：http://poj.org/problem?id=3321
这题是一开始告诉你树上每个节点都有1个苹果，然后你对一个节点操作，如果有苹果，就拿走，没苹果，就放上，然后询问你以x为根的子树上有多少个苹果。
dfs序水题，代码：

#include <map>
#include <set>
#include <stack>
#include <queue>
#include <cmath>
#include <string>
#include <vector>
#include <cstdio>
#include <cctype>
#include <cstring>
#include <sstream>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#pragma comment(linker, "/STACK:102400000,102400000")

using namespace std;
#define   MAX           500005
#define   MAXN          6005
#define   maxnode       15
#define   sigma_size    30
#define   lson          l,m,rt<<1
#define   rson          m+1,r,rt<<1|1
#define   lrt           rt<<1
#define   rrt           rt<<1|1
#define   middle        int m=(r+l)>>1
#define   LL            long long
#define   ull           unsigned long long
#define   mem(x,v)      memset(x,v,sizeof(x))
#define   lowbit(x)     (x&-x)
#define   pii           pair<int,int>
#define   bits(a)       __builtin_popcount(a)
#define   mk            make_pair
#define   limit         10000

//const int    prime = 999983;
const int    INF   = 0x3f3f3f3f;
const LL     INFF  = 0x3f3f;
const double pi    = acos(-1.0);
//const double inf   = 1e18;
const double eps   = 1e-8;
const LL    mod    = 1e9+7;
const ull    mx    = 133333331;

/*****************************************************/
inline void RI(int &x) {
      char c;
      while((c=getchar())<'0' || c>'9');
      x=c-'0';
      while((c=getchar())>='0' && c<='9') x=(x<<3)+(x<<1)+c-'0';
 }
/*****************************************************/

struct Edge{
    int v,next;
}edge[MAX*2];
int head[MAX];
int tot;
int p1[MAX];
int p2[MAX];
int ti;
int sum[MAX<<2];

void init(){
    mem(head,-1);
    tot=0;ti=0;
}

void add_edge(int a,int b){
    edge[tot]=(Edge){b,head[a]};
    head[a]=tot++;
}

void dfs(int u,int fa){
    p1[u]=++ti;
    for(int i=head[u];i!=-1;i=edge[i].next){
        int v=edge[i].v;
        if(v==fa) continue;
        dfs(v,u);
    }
    p2[u]=ti;
}

void pushup(int rt){
    sum[rt]=sum[lrt]+sum[rrt];
}

void build(int l,int r,int rt){
    if(l==r){
        sum[rt]=1;
        return;
    }
    middle;
    build(lson);
    build(rson);
    pushup(rt);
}

void update(int l,int r,int rt,int pos){
    if(l==r){
        sum[rt]=sum[rt]^1;
        return;
    }
    middle;
    if(pos<=m) update(lson,pos);
    else update(rson,pos);
    pushup(rt);
}

int query(int l,int r,int rt,int L,int R){
    if(L<=l&&r<=R) return sum[rt];
    middle;
    int ret=0;
    if(L<=m) ret+=query(lson,L,R);
    if(R>m) ret+=query(rson,L,R);
    return ret;
}

int main(){
    int n;
    cin>>n;
    init();
    for(int i=1;i<n;i++){
        int a,b;
        scanf("%d%d",&a,&b);
        add_edge(a,b);
        add_edge(b,a);
    }
    dfs(1,-1);
    build(1,n,1);
    int m;
    scanf("%d",&m);
    while(m--){
        char op;
        int a;
        getchar();
        scanf("%c%d",&op,&a);
        if(op=='Q') printf("%d\n",query(1,n,1,p1[a],p2[a]));
        else update(1,n,1,p1[a]);
    }
    return 0;
}

---

CodeForces 620E
题目传送门：http://codeforces.com/problemset/problem/620/E
给你一棵树，每个节点都有颜色，然后问你子树上有多少种不同的颜色
考虑到颜色一共只有60种，所以可以直接二进制记录，每个节点记录这个节点的颜色，然后区间直接左右子树或起来，然后对x为根的子树都变成c颜色，区间赋值，需要pushdown，pushup。
有个坑点就是需要记录每个时间戳是哪个点，然后在build线段树的时候，$sum[rt]=c[dfsnum[l]]$，这个坑点我错了好久，因为线段树的节点的下标应该是时间戳。GG，仍需努力
代码：

#include <map>
#include <set>
#include <stack>
#include <queue>
#include <cmath>
#include <string>
#include <vector>
#include <cstdio>
#include <cctype>
#include <cstring>
#include <sstream>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#pragma comment(linker, "/STACK:102400000,102400000")

using namespace std;
#define   MAX           400005
#define   MAXN          6005
#define   maxnode       15
#define   sigma_size    30
#define   lson          l,m,rt<<1
#define   rson          m+1,r,rt<<1|1
#define   lrt           rt<<1
#define   rrt           rt<<1|1
#define   middle        int m=(r+l)>>1
#define   LL            long long
#define   ull           unsigned long long
#define   mem(x,v)      memset(x,v,sizeof(x))
#define   lowbit(x)     (x&-x)
#define   pii           pair<int,int>
#define   bits(a)       __builtin_popcount(a)
#define   mk            make_pair
#define   limit         10000

//const int    prime = 999983;
const int    INF   = 0x3f3f3f3f;
const LL     INFF  = 0x3f3f;
const double pi    = acos(-1.0);
//const double inf   = 1e18;
const double eps   = 1e-8;
const LL    mod    = 1e9+7;
const ull    mx    = 133333331;

/*****************************************************/
inline void RI(int &x) {
      char c;
      while((c=getchar())<'0' || c>'9');
      x=c-'0';
      while((c=getchar())>='0' && c<='9') x=(x<<3)+(x<<1)+c-'0';
 }
/*****************************************************/

struct Edge{
    int v,next;
}edge[MAX*2];
int head[MAX];
int tot;
int c[MAX];
int p1[MAX];
int p2[MAX];
int ti;
int dfsnum[MAX];
LL sum[MAX<<2];
int col[MAX<<2];
void init(){
    mem(head,-1);
    tot=0;ti=0;
}

void add_edge(int a,int b){
    edge[tot]=(Edge){b,head[a]};
    head[a]=tot++;
}

void dfs(int u,int fa){
    p1[u]=++ti;
    dfsnum[ti]=u;
    for(int i=head[u];i!=-1;i=edge[i].next){
        int v=edge[i].v;
        if(v==fa) continue;
        dfs(v,u);
    }
    p2[u]=ti;
}

void pushup(int rt){
    sum[rt]=sum[lrt]|sum[rrt];
}

void pushdown(int rt){
    if(col[rt]){
        col[lrt]=col[rrt]=col[rt];
        sum[lrt]=sum[rrt]=(1LL<<col[rt]);
        col[rt]=0;
    }
}
void build(int l,int r,int rt){
    col[rt]=0;
    if(l==r){
        sum[rt]=(1LL<<c[dfsnum[l]]);
        return;
    }
    middle;
    build(lson);
    build(rson);
    pushup(rt);
}

void update(int l,int r,int rt,int L,int R,int d){
    if(L<=l&&r<=R){
        sum[rt]=(1LL<<d);
        col[rt]=d;
        return;
    }
    middle;
    pushdown(rt);
    if(L<=m) update(lson,L,R,d);
    if(R>m) update(rson,L,R,d);
    pushup(rt);
}

LL query(int l,int r,int rt,int L,int R){
    if(L<=l&&r<=R) return sum[rt];
    middle;
    LL ret=0;
    pushdown(rt);
    if(L<=m) ret|=query(lson,L,R);
    if(R>m) ret|=query(rson,L,R);
    return ret;
}

int main(){
    //freopen("in.txt","r",stdin);
    int n,m;
    while(cin>>n>>m){
        init();
        for(int i=1;i<=n;i++) scanf("%d",&c[i]);
        for(int i=1;i<n;i++){
            int a,b;
            scanf("%d%d",&a,&b);
            add_edge(a,b);
            add_edge(b,a);
        }
        dfs(1,-1);
        build(1,n,1);
        //cout<<query(1,n,1,p1[6],p2[6]);
        while(m--){
            int op,a;
            scanf("%d%d",&op,&a);
            if(op==1){
                int cc;
                scanf("%d",&cc);
                update(1,n,1,p1[a],p2[a],cc);
            }
            else {
                LL k=query(1,n,1,p1[a],p2[a]);
                //cout<<k<<" ";
                int num=0;
                while(k){
                    int tmp=k%2;
                    k/=2;
                    num+=tmp;
                    //if(tmp) cout<<tmp<<" ";
                }
                printf("%d\n",num);
            }
        }
    }
    return 0;
}