Bzoj2819:Nim:树链剖分

题目链接:Nim

线段树维护异或值的沙茶题

然而dfs会爆栈所以得开手工栈(见代码)

#include
#include
#include
#include
using namespace std;
const int maxn=500010;
int tot=0,h[maxn],ind=0,s[maxn],fa[maxn];
struct edge{int to,next;}G[maxn<<1];
int dep[maxn],Belong[maxn],pos[maxn];
int v[maxn],ans=0,n,m,fz[maxn];
struct seg{
	int l,r,v;
	seg *lc,*rc;
};
seg *root=new seg();
char str[5];

void add(int x,int y){
	G[++tot].to=y;G[tot].next=h[x];h[x]=tot;
}

void dfs1(int x){
	s[x]=1; dep[x]=dep[fa[x]]+1;
	for (int i=h[x];i;i=G[i].next){
		int v=G[i].to;
		if (v==fa[x]) continue;
		fa[v]=x; dfs1(v); s[x]+=s[v]; 
	}
}

void dfs2(int x,int L){
	pos[x]=++ind; int k=0; Belong[x]=L; fz[ind]=x;
	for (int i=h[x];i;i=G[i].next)
		if (s[G[i].to]>s[k]&&dep[G[i].to]>dep[x]) k=G[i].to;
	if (!k) return; dfs2(k,L);
	for (int i=h[x];i;i=G[i].next)
		if (dep[G[i].to]>dep[x]&&k!=G[i].to) dfs2(G[i].to,G[i].to); 
}

void push_up(seg *p){
	p->v=p->lc->v^p->rc->v;
}

void build(seg *p,int l,int r){
	p->l=l; p->r=r;
	if (l+1==r){
		p->v=v[fz[l]]; p->lc=p->rc=NULL; return;
	}else if (l+1>1;
		p->lc=new seg();
		p->rc=new seg();
		if (llc,l,mid);
		if (midrc,mid,r);
		push_up(p);
	}
}

void update(seg *p,int l,int r,int val){
	if (l<=p->l&&p->r<=r){p->v=val;return;}
	int mid=(p->l+p->r)>>1;
	if (llc,l,r,val);
	if (midrc,l,r,val);
	push_up(p);
}

void ask(seg *p,int l,int r){
	if (l<=p->l&&p->r<=r){ans^=p->v;return;}
	int mid=(p->l+p->r)>>1;
	if (llc,l,r);
	if (midrc,l,r); 
}

int query(int x,int y){
	int sum=0;
	while (Belong[x]!=Belong[y]){
		if (dep[Belong[x]]


你可能感兴趣的:(OI,树链剖分)