【基础练习】【倍增LCA】codevs1036 商务旅行题解

题目描述 Description

某首都城市的商人要经常到各城镇去做生意，他们按自己的路线去做，目的是为了更好的节约时间。

假设有N个城镇，首都编号为1，商人从首都出发，其他各城镇之间都有道路连接，任意两个城镇之间如果有直连道路，在他们之间行驶需要花费单位时间。该国公路网络发达，从首都出发能到达任意一个城镇，并且公路网络不会存在环。

你的任务是帮助该商人计算一下他的最短旅行时间。

输入描述 Input Description

输入文件中的第一行有一个整数N，1<=n<=30 000，为城镇的数目。下面N-1行，每行由两个整数a 和b (1<=a, b<=n; a<>b)组成，表示城镇a和城镇b有公路连接。在第N+1行为一个整数M，下面的M行，每行有该商人需要顺次经过的各城镇编号。

输出描述 Output Description

    在输出文件中输出该商人旅行的最短时间。

样例输入 Sample Input

5

1 2

1 5

3 5

4 5

4

1

3

2

5

样例输出 Sample Output

7

这是一道很裸的LCA 乍一看似乎是图 但是边权只有1 而且符合树的结构 按顺序求出每两个点之间的最近路径相加即可

自从“水果街逛水果街2”参考了胜利一中sunshine神犇的代码 有了一些改进

先把代码放上

//codevs1036 商务旅行 LCA
//copyright by ametake
//since 19:20
#include
   
   
    
    
#include
    
    
     
     
#include
     
     
      
      
#include
      
      
       
       
using namespace std;

const int maxn=30000+10;
int depth[maxn],fa[maxn];
int head[maxn],to[maxn*2],next[maxn*2];//twice again!!!
int n,m;
int et;
int up[maxn][20];

void add(int x,int y)
{
	et++;
	to[et]=y;
	next[et]=head[x];
	head[x]=et;
}

void dfs(int now)
{
	int t=head[now];
	while (t)
	{
		if (depth[to[t]])
		{
			t=next[t];
			continue;
		}
		depth[to[t]]=depth[now]+1;//不能写depth[now]++ 这样depth[now]也会+1 
		fa[to[t]]=now;
		dfs(to[t]);
		t=next[t];
	}
	return;
}

void doit()
{
	for (int i=1;i<=n;i++) up[i][0]=fa[i];
	for (int j=1;j<=14;j++)
	{
		for (int i=1;i<=n;i++)
		{
			up[i][j]=up[up[i][j-1]][j-1];
		}
	}
}

int lca(int x,int y)
{
	if (depth[x]
       
       
         =0;i--) { if (up[x][i]!=up[y][i]) { x=up[x][i]; y=up[y][i]; } } return up[x][0]; } int main() { scanf("%d",&n); int x,y; memset(head,0,sizeof(head)); memset(depth,0,sizeof(depth)); for (int i=1;i 
         
       
      
      
     
     
    
    
   
   

具体来讲学习的是这一段核心代码

int lca(int x,int y)
{
	if (depth[x]<depth[y]) swap(x,y);
	int dt=depth[x]-depth[y];
    for (int i=0;i<=14;i++)
	{
		if (dt&(1<<i)) x=up[x][i];
	}
	if (x==y) return x;
	for (int i=14;i>=0;i--)
	{
		if (up[x][i]!=up[y][i])
		{
			x=up[x][i];
			y=up[y][i];
		}
	}
	return up[x][0];
}

这里用了一个与运算，这个与运算的作用是判断某个数x的2的幂次方分解中是否含有2的i次方 

在二进制的状况下 如果有2的i次方 那么在第i位是1 否则就是0 因此如果返回true说明有这个幂

sunshine神犇还有一个值得我学习的地方是这儿：

if(deep[x]<(1<<i))break;

这个句子是放在搜索和转移中的 有效减少了层次过多问题 虽然时间似乎没节约多少 但总是好的

sunshine神犇的lca和我的在实现上有很大的不同，因为他是一边搜索一边倍增。首先倍增浅层的点，每次处理完这一层，搜索更深一层并进行深一层的状态转移倍增。当然我认为我们时间差别不大，但是sunshine神犇的似乎更高更妙一些。

马上就要闭关了。到时候一定要三更灯火五更鸡啊。看ACM前辈们每天只睡四五个小时，真是佩服。

其实我一天也就睡四五个小时，但是时间利用效率显然还有很大提升空间啊。

十一要是不能学文化课的话，就利用晚上什么的，虽然估计我没心情。如果实在不行的话 就请假 总之必须动静分离，不能混在一起。

闭关的时候，前期每天刷题练熟算法，一天至少A十道题吧。然后中午听听歌，晚上学学英语，听听歌，看看书。

《十八岁和其他》里，杨子说，他读高中的时候，每天捧着肖洛霍夫《静静的顿河》彻夜不眠。我现在几乎也是如此，不敢彻夜不眠，然而这部书实在引人入胜。

但是现阶段还是高考和OI是最重要的。下周开始，高考也放下，OI全力以赴。

所以加油吧。为了梦想，为了快乐。

——纨绔不饿死，儒冠多误身