每日模板一练——树链剖分求LCA(好吧是真的写错了)

大意

给一个N个节点的树和M个询问,对于每次询问输出两点的距离。

样例:

6
1 2
1 3
2 4
2 5
3 6
2
2 6
5 6

输出:

3

4

【代码~】

#include<bits/stdc++.h>
using namespace std;
const int MAXN=1e5+10;
const int MAXM=3e5+10;

int n,m,cnt;
int head[MAXN],depth[MAXN],son[MAXN],siz[MAXN],father[MAXN];
int top[MAXN];
int nxt[MAXM],to[MAXM];

int Read()
{
	int i=0,f=1;
	char c;
	for(c=getchar();(i>'9'||c<'0')&&c!='-';c=getchar());
	if(c=='-')
	  f=-1,c=getchar();
	for(;c>='0'&&c<='9';c=getchar())
	  i=(i<<3)+(i<<1)+c-'0';
	return i*f;
}

void Add(int x,int y)
{
	cnt++;
	nxt[cnt]=head[x];
	head[x]=cnt;
	to[cnt]=y;
}

void add(int x,int y)
{
	Add(x,y);
	Add(y,x);
}

void dfs1(int u,int fa)
{
	father[u]=fa;
	siz[u]=1,son[u]=0;
	for(int i=head[u];i!=-1;i=nxt[i])
	{
		int v=to[i];
		if(v!=fa)
		{
			depth[v]=depth[u]+1;
			dfs1(v,u);
			siz[u]+=siz[v];
			if(siz[v]>siz[son[u]])
			  son[u]=v;
		}
	}
}

void dfs2(int u,int fa)
{
	if(u==son[fa])
	  top[u]=top[fa];
	else
	  top[u]=u;
	for(int i=head[u];i!=-1;i=nxt[i])
	{
		int v=to[i];
		if(v!=fa)
		  dfs2(v,u);
	}
}

int lca(int x,int y)
{
	while(top[x]!=top[y])
	{
		if(depth[top[x]]<depth[top[y]])
		  swap(x,y);
		x=father[top[x]];
	}
	return depth[x]<depth[y]?x:y;
}

int main()
{
	memset(head,-1,sizeof(head));
	n=Read();
	for(int i=1;i<n;++i)
	{
		int x=Read(),y=Read();
		add(x,y);
	}
	dfs1(1,-1);
	dfs2(1,1);
	m=Read();
	for(int i=1;i<=m;++i)
	{
		int x=Read(),y=Read();
		int res=depth[x]+depth[y]-2*depth[lca(x,y)];
		cout<<res<<'\n';
	}
	return 0;
}

 

posted @ 2018-10-11 19:11  Ishtar~  阅读(108)  评论(0编辑  收藏  举报