【题目链接】

点击打开链接

【算法】

倍增法求最近公共祖先

【代码】

#include <algorithm>
#include <bitset>
#include <cctype>
#include <cerrno>
#include <clocale>
#include <cmath>
#include <complex>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <ctime>
#include <deque>
#include <exception>
#include <fstream>
#include <functional>
#include <limits>
#include <list>
#include <map>
#include <iomanip>
#include <ios>
#include <iosfwd>
#include <iostream>
#include <istream>
#include <ostream>
#include <queue>
#include <set>
#include <sstream>
#include <stdexcept>
#include <streambuf>
#include <string>
#include <utility>
#include <vector>
#include <cwchar>
#include <cwctype>
#include <stack>
#include <limits.h>
using namespace std;
#define MAXN 10010
#define MAXLOG 20 struct Edge
{
int to,nxt;
} e[MAXN]; int i,T,tot,n,root,x,y;
int fa[MAXN],dep[MAXN],head[MAXN],anc[MAXN][MAXLOG]; template <typename T> inline void read(T &x)
{
int f = ; x = ;
char c = getchar();
for (; !isdigit(c); c = getchar()) { if (c == '-') f = -f; }
for (; isdigit(c); c = getchar()) x = (x << ) + (x << ) + c - '';
x *= f;
}
template <typename T> inline void write(T x)
{
if (x < )
{
putchar('-');
x = -x;
}
if (x > ) write(x/);
putchar(x%+'');
}
template <typename T> inline void writeln(T x)
{
write(x);
puts("");
}
inline void add(int x,int y)
{
tot++;
e[tot] = (Edge){y,head[x]};
head[x] = tot;
}
inline void lca_init(int u)
{
int i,v;
anc[u][] = fa[u];
for (i = ; i < MAXLOG; i++)
{
if (dep[u] <= (i << )) break;
anc[u][i] = anc[anc[u][i-]][i-];
}
for (i = head[u]; i; i = e[i].nxt)
{
v = e[i].to;
dep[v] = dep[u] + ;
lca_init(v);
}
}
inline int lca(int u,int v)
{
int i,k;
if (dep[u] > dep[v]) swap(u,v);
k = dep[v] - dep[u];
for (i = ; i < MAXLOG; i++)
{
if (k & ( << i))
v = anc[v][i];
}
if (u == v) return u;
for (i = MAXLOG - ; i >= ; i--)
{
if (anc[u][i] != anc[v][i])
{
u = anc[u][i];
v = anc[v][i];
}
}
return fa[u];
} int main() { read(T);
while (T--)
{
memset(fa,,sizeof(fa));
memset(anc,,sizeof(anc));
read(n);
tot = ;
for (i = ; i <= n; i++) head[i] = ;
for (i = ; i < n; i++)
{
read(x); read(y);
add(x,y);
fa[y] = x;
}
for (i = ; i <= n; i++)
{
if (!fa[i])
root = i;
}
dep[root] = ;
lca_init(root);
read(x); read(y);
writeln(lca(x,y));
}
return ; }

最新文章

  1. WebService技术(一)
  2. jython安装
  3. 详解APM数据采样与端到端
  4. 【bzoj2179】FFT快速傅立叶 FFT模板
  5. mfc_随机数生成器
  6. 从源代码的角度聊聊java中StringBuffer、StringBuilder、String中的字符串拼接
  7. PAT-乙级-1047. 编程团体赛(20)
  8. (转载)根据数据字典表定义的表结构,生成创建表的SQL语句
  9. bzoj 1923 [Sdoi2010]外星千足虫(高斯消元+bitset)
  10. Xamarin Android 中Acitvity如何传递数据
  11. OpenCV, MatBGR2ARGB, ARGB2MatBGR
  12. Java实现杨辉三角
  13. dhtmlx Gantt知识点1
  14. activemq消息队列的使用及应用docker部署常见问题及注意事项
  15. 线程---local数据隔离
  16. 正睿 2018 提高组十连测 Day4 T3 碳
  17. 采用RedisLive监控Redis服务——安装手册
  18. JDK中ConcurrentHashMap效率测试
  19. MYSQL中group_concat有长度限制!默认1024(转载)
  20. Java设计模式(7)装饰模式(Decorator模式)

热门文章

  1. winform ComboBox/TextBox自动提示
  2. Python 字典(2)
  3. UvaLive 4917 Abstract Extract (模拟)
  4. List&lt;&gt; 集合 删除指定行
  5. 《 阿Q正传》-鲁迅 词语解释 | 经典语录
  6. mongodb shell 无法删除问题
  7. jQuery的观察者模式详解 转载
  8. PKCS填充方式
  9. HDU——2063 过山车
  10. hdu5618(cdq分治求三维偏序)