zoukankan      html  css  js  c++  java
  • LCA SP913 QTREE2

    SP913 QTREE2 - Query on a tree II

    给定一棵n个点的树,边具有边权。要求作以下操作:

    DIST a b 询问点a至点b路径上的边权之和

    KTH a b k 询问点a至点b有向路径上的第k个点的编号

    有多组测试数据,每组数据以DONE结尾。

    裸的LCA。

    在处理第二个操作时,我直接向上数跳了多少个。

    顾z大佬说不能这么做,要求出跳到那个点的深度再去跳。

    真的是这样,不过懒得想了,应该是+1-1的误差。 balabala。。。

    code:

    #include <iostream>
    #include <cstdio>
    #include <cstring>
    
    using namespace std;
    
    const int wx=50017;
    
    int dep[wx],dis[wx];
    int f[wx][23];
    int head[wx];
    int num,n,t;
    char opt[7];
    
    inline int read(){
    	int sum=0,f=1; char ch=getchar();
    	while(ch<'0'||ch>'9'){if(ch=='-')f=-1; ch=getchar();}
    	while(ch>='0'&&ch<='9'){sum=(sum<<1)+(sum<<3)+ch-'0'; ch=getchar();}
    	return sum*f;
    }
    
    struct e{
    	int nxt,to,dis;
    }edge[wx*2];
    
    void add(int from,int to,int dis){
    	edge[++num].nxt=head[from];
    	edge[num].to=to;
    	edge[num].dis=dis;
    	head[from]=num;
    }
    
    void dfs(int u,int fa){
    	dep[u]=dep[fa]+1;
    	for(int i=head[u];i;i=edge[i].nxt){
    		int v=edge[i].to;
    		if(v==fa)continue;
    		f[v][0]=u;dis[v]=dis[u]+edge[i].dis;
    		dfs(v,u);
    	}
    }
    
    void pre(){
    	for(int j=1;j<=21;j++)for(int i=1;i<=n;i++)f[i][j]=f[f[i][j-1]][j-1];
    }
    
    int LCA(int x,int y){
    	if(dep[x]<dep[y])swap(x,y);
    	for(int i=21;i>=0;i--){
    		if(dep[f[x][i]]>=dep[y]){
    			x=f[x][i];
    		}
    	}
    	if(x==y)return x;
    	for(int i=21;i>=0;i--){
    		if(f[x][i]!=f[y][i]){
    			x=f[x][i]; y=f[y][i];
    		}
    	}
    	return f[x][0];
    }
    
    int find(int x,int k){
    	for(int i=21;i>=0;i--){
    		if(dep[f[x][i]]>=k)x=f[x][i];
    	}
    	return x;
    }
    
    int main(){
    	t=read();
    	while(t--){
    		n=read();
    		memset(head,0,sizeof head); num=1;
    		memset(edge,0,sizeof edge);
    		for(int i=1;i<n;i++){
    			int x,y,z;
    			x=read(); y=read(); z=read();
    			add(x,y,z); add(y,x,z);
    		}
    		dfs(1,0); pre();
    		while(1){
    			scanf("%s",opt+1);
    			if(opt[2]=='O')break;
    			if(opt[2]=='I'){
    				int x,y;
    				x=read(); y=read();
    				int lca=LCA(x,y);
    				printf("%d
    ",dis[x]+dis[y]-2*dis[lca]);
    			}
    			if(opt[1]=='K'){
    				int a,b,k;
    				a=read(); b=read(); k=read();
    				int lca=LCA(a,b);
    				if(dep[a]-dep[lca]+1>=k)printf("%d
    ",find(a,dep[a]-k+1));
    				else printf("%d
    ",find(b,k-dep[a]+2*dep[lca]-1));
    			}
    		}
    	}
    }
  • 相关阅读:
    MongoDB小结25
    MongoDB小结24
    MongoDB小结23
    MongoDB小结22
    MongoDB小结21
    MongoDB小结20
    MongoDB小结19
    MongoDB小结18
    hdu 4606 Occupy Cities
    hdu 4610 Cards
  • 原文地址:https://www.cnblogs.com/wangxiaodai/p/9858833.html
Copyright © 2011-2022 走看看