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;
}