码迷,mamicode.com
首页 > 其他好文 > 详细

hdu6031 Innumerable Ancestors

时间:2017-07-06 16:48:15      阅读:180      评论:0      收藏:0      [点我收藏+]

标签:style   size   main   get   line   无向图   its   rabl   快速   

hdu6031 Innumerable Ancestors
倍增
题意
给定一张无向图,对于每组询问,给出两个集合 A 与 B
求 lca(x,y) 最深的时候的深度 x属于A y属于B

题解
首先 我们发现答案具有单调性,于是我们就可以二分这个答案 mid
然后 把 a 集合 的 这个深度的祖先 求出来 然后看 与b 集合的这个深度 祖先 是否有交集就行了
用set 判断 是否有交集
用倍增法 快速求出一个点的祖先 求一次 只要 log n

 

  1 #include <bits/stdc++.h> 
  2 #define ll long long  
  3 #define cl(a,b) memset(a,b,sizeof(a)) 
  4 using namespace std ;  
  5 
  6 const int N = 1e5 + 11 ; 
  7 struct Edge{
  8     int to,pre ; 
  9 }e[N*2] ; 
 10 int n,Q,cnt,k1,k2,ans,l,r ; 
 11 int head[N],deep[N],anc[N][21],a[N],b[N] ;   
 12 
 13 inline int read() 
 14 {
 15     int x = 0 , f = 1 ; 
 16     char ch = getchar() ; 
 17     while(ch<0||ch>9) { if(ch==-) f = -1 ; ch = getchar() ; } 
 18     while(ch>=0&&ch<=9) { x = x * 10+ch-48 ; ch = getchar() ; } 
 19     return x * f ; 
 20 }
 21 
 22 inline void add(int x,int y) 
 23 {
 24     e[++cnt].to = y ; 
 25     e[cnt].pre = head[x] ; 
 26     head[x] = cnt ; 
 27 }
 28 
 29 inline void init()  
 30 {
 31     cl( head,-1 ) ;
 32     cl( deep,0 ) ;  
 33     cl( anc,-1 ) ; 
 34     cnt = 0 ; 
 35 }
 36 
 37 inline void dfs(int u,int fa,int d) 
 38 {
 39     deep[ u ] = d ; 
 40     for(int i=head[u]; ~i ; i = e[i].pre)
 41     {
 42         int v = e[ i ].to ; 
 43         if(v==fa) continue ; 
 44         anc[ v ][ 0 ] = u ; 
 45         dfs(v,u,d+1) ; 
 46     }
 47 }
 48 
 49 inline void lca_init() 
 50 {
 51     for(int j=1; (1<<j) <= n ; j++) 
 52         for( int i=1;i<=n;i++ ) 
 53             if(anc[i][j-1]!=-1) 
 54                 anc[i][j] = anc[anc[i][j-1]][j-1] ; 
 55 }
 56 
 57 inline int query(int u,int d) 
 58 //  返回  u的 第 d 个父亲 
 59 {
 60     if( d < 0 ) return -1 ; 
 61     if( d==0 ) return u ; 
 62     int i ; 
 63     for(i=0; (1<<i) <= d ;i++) ;   //  注意 这里有分号相当于 求 log  没有的话 死循环  
 64     i-- ; 
 65     for( ;i>=0;i--)
 66         if(d-(1<<i) >= 0) 
 67         {
 68             d-=(1<<i) ; 
 69             u = anc[u][i] ; 
 70         } 
 71     return u ; 
 72 } 
 73 
 74 inline bool check(int x) 
 75 {
 76     set<int>s ; 
 77     for(int i=1;i<=k1;i++) 
 78     {
 79         int dis = deep[a[ i ]] - x ; 
 80         int ret = query(a[ i ],dis) ; 
 81         if(ret==-1) continue ; 
 82         s.insert(ret) ; 
 83     }
 84     for(int i=1;i<=k2;i++) 
 85     {
 86         int dis = deep[b[ i ]] - x ; 
 87         int ret = query( b[ i ],dis ) ; 
 88         if(s.count(ret)) return true ;   
 89     }
 90     return false ; 
 91 }
 92 
 93 int main() 
 94 { 
 95     while(~scanf("%d%d",&n,&Q)) 
 96     {
 97         init() ; 
 98         int x,y ; 
 99         for(int i=1;i<n;i++) 
100         {
101             scanf("%d%d",&x,&y) ; 
102             add(x,y) ; add(y,x) ; 
103         }
104         dfs(1,1,1) ; 
105         lca_init() ; 
106         
107         while(Q--) 
108         {
109             scanf("%d",&k1) ; 
110             l = 1 ; r = 1 ; 
111             for(int i=1;i<=k1;i++) 
112                 scanf("%d",&a[ i ]),r = max(r,deep[ a[ i ] ]) ; 
113             scanf("%d",&k2) ; 
114             for(int i=1;i<=k2;i++) 
115                 scanf("%d",&b[ i ]) ;  
116             while(l<=r) 
117             {
118                 int mid = ( l+r ) >>1 ; 
119                 if(check(mid)) 
120                 {
121                     ans = mid ; 
122                     l = mid + 1 ; 
123                 }
124                 else 
125                     r = mid - 1 ; 
126             } 
127             printf("%d\n",ans ) ; 
128         }
129     }
130     
131     return 0 ; 
132 }

 

hdu6031 Innumerable Ancestors

标签:style   size   main   get   line   无向图   its   rabl   快速   

原文地址:http://www.cnblogs.com/third2333/p/7126846.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!