poj-1330-Nearest Common Ancestors-LCA（RMQ）

题目传送门：http://poj.org/problem?id=1330

 

 

什么是LCA？

公共祖先简称LCA（Lowest Common Ancestor），所谓LCA，是当给定一个有根树T时，对于任意两个结点u、v，找到一个离根最远的结点x，使得x同时是u和v的祖先，x 便是u、v的最近公共祖先。

 

额￣へ￣  听说有好多种解法，但是我就先发一种，其他回来补上

 如果RMQ不会，请看http://blog.csdn.net/newproblems/article/details/77183394

开始

 

 

 

 

 

1 ．在遍历数组中找到7和 16 第一次出现的位置。

 

2 ．在节点深度中，找到涂色序列的最小值----1，则4就是7和16 的共同祖先。

 

 

 

上图解释的最快O(∩_∩)O

 

邻接表存图

 

不定长数组：vecter

定义方式：vecter<int> a

 

基本的操作：

 

a.size()     读取它的大小

a.resize()   改变大小

a.push_back()  向尾部添加元素

a.pop_back()   删除最后一个元素

a.empty()   判断数组是否为空

a.clear()   清空数组

 

另外，常用vecter<int> a[max]（相当于二维数组，第一维大小固定，第二维大小不固定）

 

邻接表是对图中的每个顶点建立起一个单链表。大白话，就是与下标，相连的点放到对应下标的数组中。用vector就可以实现。例如：上图的a[1]->14->13

 

vector数组内部的存储方式类似于链表，直接用就行了。

 

图的度是指和该顶点相关联的边数。

在有向图中，度又分为入度和出度。

顶点的入度：以某顶点为弧头，终止于该顶点的弧的数目。（指向它的有几条）

根节点的入度为0。

 

 

 

完整代码

 

代码解释：

 

遍历数组------------>  vis
节点深度------------>  depth

邻接表图------------>  G     

入度表--------------->  in

 

id[i]是记录i第一次出现的位置

 

本题之所以要记录入度，是因为题目未给出根节点的位置需要通过遍历in找出

还要注意，这是下标在找最小值

 

#include <algorithm>  
#include <iostream>  
#include <cstring>  
#include <cstdlib>  
#include <cstdio>  
#include <cmath>  
#include <queue>  
#include <vector>  
#include <set>  
#include <map>  
#define N 300010  
using namespace std;  
  
vector<int> G[N];  
int vis[N], id[N], depth[N], in[N], k;  
int dp[N][30];  
  
int  Min(int a, int b)  
{  
    if(depth[a] <= depth[b])  
        return a;  
    return b;  
}  
  
void init_rmq(int n)  
{  
    for(int i = 1; i <= n; i++)  
        dp[i][0] = i;  
    for(int j = 1; (1 << j) <= n; j++)  
        for(int i = 1; i + (1 << j) - 1 <= n; i++)  
            dp[i][j] = Min(dp[i][j - 1], dp[i + (1 << (j - 1))][j - 1]);  
}  
  
void dfs(int v, int d)  
{  
    id[v] = k;  
    vis[k] = v;  
    depth[k] = d;  
    k++;  
    for(int i = 0; i < G[v].size(); i++)  
    {  
        dfs(G[v][i], d + 1);  
        vis[k] = v;  
        depth[k] = d;  
        k++;  
    }  
}  
  
void init(int root, int n)  
{  
    k = 1;  
    dfs(root, 0);//递归生成上面的那个表格  
    init_rmq(2 * n - 1);  
}  
  
int query(int l, int r)  
{  
    int len = r - l + 1;  
    int k = (int)(log10(len * 1.0) / log10(2.0));  
    int ans = Min(dp[l][k], dp[r - (1 << k) + 1][k]);  
    return vis[ans];  
}  
  
int main()  
{  
    int T;  
    scanf("%d", &T);  
    while(T--)  
    {  
        int n;  
        scanf("%d", &n);  
  
        for(int i = 1; i < N; i++)//多组数据需清空  
            G[i].clear();  
  
        memset(in, 0, sizeof(in));  
  
        for(int i = 1; i < n; i++)// 存图  
        {  
            int a, b;  
            scanf("%d%d", &a, &b);  
            G[a].push_back(b);  
            in[b] = 1;  
        }  
  
        int i;  
        for(i = 1; i <= n; i++)  
            if(in[i] == 0)  
                break;//找根节点  
  
        init(i, n);  
  
        int x, y;  
        scanf("%d%d", &x, &y);  
  
        int ans = query(min(id[x], id[y]), max(id[x], id[y]));  
        printf("%d\n", ans);  
    }  
    return 0;  
}  

 

转载，大佬写的，感觉写的很好。