• HDU 5266 pog loves szh III 线段树,lca


    Pog and Szh are playing games. Firstly Pog draw a tree on the paper. Here we define 1 as the root of the tree.Then Szh choose some nodes from the tree. He wants Pog helps to find the least common ancestor (LCA) of these node.The question is too difficult for Pog.So he decided to simplify the problems.The nodes picked are consecutive numbers from lili to riri ([li,ri])([li,ri]). 

    Input

    Several groups of data (no more than 3 groups,n10000n≥10000 or Q10000Q≥10000). 
    The following line contains ans integers,n(2n300000)n(2≤n≤300000). 
    AT The following n1n−1 line, two integers are bibi and cici at every line, it shows an edge connecting bibi and cici. 
    The following line contains ans integers,Q(Q300000)Q(Q≤300000). 
    AT The following QQ line contains two integers li and ri(1lirin1≤li≤ri≤n).

    Output

    For each case,output QQ integers means the LCA of [li,ri][li,ri].

    Sample Input

    5
    1 2
    1 3
    3 4
    4 5
    5
    1 2
    2 3
    3 4
    3 5
    1 5
    Sample Input
    1
    1
    3
    1
    1
    一颗树以1为根节点,给出区间范围,求范围内所有点的lca公共祖先,之前学的lca算法只能求两个点的,但考虑到是区间内进行操作,故考虑线段树。这道题如果了解线段树的话可以说没什么难度了...
    做这道题这道题之前,没学过线段树...搞了好久才弄清楚,在网上抄了个线段树模板,弄清楚后没有想象中的难,实际上就是简单的lca+线段树组合罢了,可以说是模板题也不为过,收获就是学会了bfs树上倍增和搞清楚了线段树吧,
    这个数据结构挺神奇的。
    该题G++提交会TLE,C++提交如果不手动扩栈就会爆栈,当然也可以用BFS倍增就不用担心爆栈了。
      1 #include <iostream>
      2 #include <string.h>
      3 #include <stdio.h>
      4 #include <math.h>
      5 #include <cstdio>
      6 #include <queue>
      7 #pragma warning ( disable : 4996 )
      8 #pragma comment(linker, "/STACK:102400000,102400000")  //防止爆栈
      9 
     10 using namespace std;
     11 
     12 int Max( int x, int y ) { return x>y?x:y; }
     13 int Min( int x, int y ) { return x>y?y:x; }
     14 
     15 const int inf = 0x3f3f3f3f;
     16 const int vspot = 3e5 + 5;
     17 const int espot = 1e5 + 5;
     18 
     19 struct node {
     20     int from, to, next;
     21 }e[vspot<<1];
     22 
     23 int N, Q, cnt, logn;
     24 int lca[vspot<<2], linjie[vspot];
     25 int fa[vspot][22], depth[vspot];
     26 
     27 void addedge(int u, int v)
     28 {
     29     e[cnt].from = u; e[cnt].to = v;
     30     e[cnt].next = linjie[u]; linjie[u] = cnt++;
     31     e[cnt].from = v; e[cnt].to = u;
     32     e[cnt].next = linjie[v]; linjie[v] = cnt++;
     33 }
     34 
     35 void init()
     36 {
     37     cnt = 0;
     38     logn = 20;
     39     memset( linjie, -1, sizeof(linjie) );
     40     memset( lca, 0, sizeof(lca) );
     41     memset( depth, 0, sizeof(depth) );
     42     memset( fa, 0, sizeof(fa) );
     43 }
     44 
     45 void dfs( int x )
     46 {
     47     for( int i = 1; i <= logn; i++ )
     48         if ( fa[x][i-1] )
     49             fa[x][i] = fa[fa[x][i-1]][i-1];
     50         else
     51             break;
     52     for ( int i = linjie[x]; i+1; i = e[i].next )
     53         if ( fa[x][0] != e[i].to )
     54         {
     55             fa[e[i].to][0] = x;
     56             depth[e[i].to] = depth[x]+1;
     57             dfs(e[i].to);
     58         }
     59 }
     60 
     61 void bfs( int x )
     62 {
     63     queue<int> Q;
     64     Q.push(x);
     65 
     66     while (!Q.empty())
     67     {
     68         int run = Q.front();
     69         Q.pop();
     70         for( int i = linjie[run]; i+1; i = e[i].next )
     71             if ( fa[run][0] != e[i].to )
     72             {
     73                 fa[e[i].to][0] = run;
     74                 depth[e[i].to] = depth[run] + 1;
     75                 Q.push(e[i].to);
     76             }
     77     }
     78 }
     79 
     80 void init2()
     81 {
     82     for ( int i = 0; i <= logn; i++ )
     83         for ( int x = 1; x <= N; x++ )
     84             fa[x][i+1] = fa[fa[x][i]][i];
     85 }
     86 
     87 int Lca( int u, int v )
     88 {
     89     if ( depth[u] < depth[v] )
     90         swap(u,v);
     91     int d = depth[u] - depth[v];
     92     for ( int i = 0; i <= logn; i++ )
     93         if( (1<<i)&d )
     94             u = fa[u][i];
     95 
     96     if ( u == v ) return u;
     97     for ( int i = logn; i >= 0; i-- )
     98         if ( fa[u][i] != fa[v][i] )
     99         {
    100             u = fa[u][i];
    101             v = fa[v][i];
    102         }
    103     return fa[u][0];
    104 }
    105 
    106 void pushUp( int pos ) { lca[pos] = Lca(lca[pos<<1], lca[pos<<1|1]); }
    107 
    108 void build( int lhs, int rhs, int id )
    109 {
    110     if ( lhs == rhs )
    111     {
    112         lca[id] = lhs;
    113         return;
    114     }
    115     int mid = ( lhs + rhs ) >> 1;
    116     build( lhs, mid, id<<1 );
    117     build( mid+1, rhs, id<<1|1 );
    118 
    119     pushUp(id);
    120 }
    121 
    122 //lhs, rhs表示当前节点的区间,l,r表示要操作的区间
    123 int query( int lhs, int rhs, int l, int r, int id )
    124 {
    125     if ( l <= lhs && r >= rhs )
    126         return lca[id];
    127     int mid = (lhs+rhs)>>1;
    128 
    129     int llca = -1, rlca = -1;
    130     if ( l <= mid ) llca = query( lhs, mid, l, r, id<<1 );
    131     if ( r >  mid ) rlca = query( mid+1, rhs, l, r, id<<1|1 );
    132 
    133     if ( llca == -1 || rlca == -1 )
    134         return Max(llca, rlca);
    135     else
    136         return Lca(llca, rlca);
    137 }
    138 
    139 int main()
    140 {
    141     while ( ~scanf("%d", &N) )
    142     {
    143         init();
    144         int x, y;
    145         for( int i = 1; i < N; i++ )
    146             { scanf("%d %d", &x, &y ); addedge(x,y); }
    147         bfs(1);
    148         init2();
    149         // dfs(1);   用dfs并去掉上面两行也行
    150         build(1,N,1);
    151         
    152         cin >> Q;
    153         int ans;
    154         for ( int i = 1; i <= Q; i++ ) 
    155         {
    156             scanf( "%d %d", &x, &y );
    157             ans = query( 1, N, x, y, 1 );
    158             printf( "%d
    ", ans );
    159         }
    160     }
    161     return 0;
    162 }
    
    
    


  • 相关阅读:
    4、数组及字符串
    3、处理数据
    2、开始学习C++
    1、C++学习预备知识
    Django【进阶篇 】
    Django【基础篇】
    Python操作MySQL
    MySQL(二)
    python日记-使用队列写出特殊字数串
    python日记-快速排序算法
  • 原文地址:https://www.cnblogs.com/chaoswr/p/8490400.html
Copyright © 2020-2023  润新知