POJ 1470 Closest Common Ancestors B (LCA算法)

题目链接：

http://poj.org/problem?id=1470

题目大意：

也是求最近公共祖先的个数，这道题的特别之处，在于读取数据的过程（至少我是这么认为的），其他便是LCA算法的应用，如果不懂的话，看这里 http://blog.csdn.net/qq_36287702/article/details/77368464

LCA分为在线算法和离线算法，在线算法一般是有：vector+dfs 或者邻接表+dfs 离线算法 我目前只知道 并查集+dfs 离线算法速度很快，但我还不会用QAQ。。。而在线算法我现在也只是搞懂了vector+dfs （感觉LCA好难。。）

下面是代码：

 

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <climits>
#include <cmath>
#include<vector>
using namespace std;
const int N=10010;
vector<int>vec[N];
int mnn;//不断记录出现的下标
int f[N],depth[N<<1],occur[N<<1],dp[N<<1][30],in[N],vis[N];
//f[]记录结点在搜素中最先出现的位置下标
//depth[]记录结点在树中的深度 与occur相对应
//occur[]记录结点出现的数组，重复的也要记录
//dp[][]表示从第i个位置开始的2^j个元素中最小值的下标
//in[]标记数组
void dfs(int u,int deep)
{
    occur[++mnn]=u;
    f[u]=mnn;
    depth[mnn]=deep;
    for(int i=0; i<vec[u].size(); i++)
    {
        dfs(vec[u][i],deep+1);
        occur[++mnn]=u;
        depth[mnn]=deep;
    }
}
void 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]=depth[dp[i][j-1]]<depth[dp[i+(1<<(j-1))][j-1]]?dp[i][j-1]:dp[i+(1<<(j-1))][j-1];
//认真看，不要以为长就难了 其实很简单
}
int RMQ_min(int a,int b)
{
    int l=f[a];
    int r=f[b];
    if(l>r) swap(l,r);
    int k=log(r-l+1.0)/log(2.0);
    int minx=depth[dp[l][k]]<depth[dp[r-(1<<k)+1][k]]?dp[l][k]:dp[r-(1<<k)+1][k];//最小值下标
    return occur[minx];
}
int main()
{
    int n;
    while(~scanf("%d",&n))
    {
        int q,v;
        memset(in,0,sizeof(in));
        memset(f,0,sizeof(f));
        for(int i=1; i<=n; i++)
            vec[i].clear();
        for(int i=1; i<=n; i++)
        {
            scanf("%d:(%d)",&q,&v);
            while(v--)
            {
                int u;
                scanf("%d",&u);
                vec[q].push_back(u);
                in[u]=1;
            }
        }
        int i;
        for(i=1; i<=n; i++)
            if(in[i]==0)break;//寻找根节点
        mnn=0;
        dfs(i,0);
        RMQ(2*n-1);//如果重复记录的话，就会有2*n-1条边（记住就行 慢慢就理解了）
        int s;
        memset(vis,0,sizeof(vis));
        scanf("%d",&s);
        while(s--)
        {
            while(true)
            {
                char ch;
                scanf("%c",&ch);
                if(ch=='(')//对输入的限定
                {
                    int x,y;
                    scanf("%d %d)",&x,&y);
                    int ans=RMQ_min(x,y);
                    vis[ans]++;
                    break;
                }
            }
        }
        for(int i=1; i<=n; i++)
        {
            if(vis[i])
                printf("%d:%d\n",i,vis[i]);
        }
    }
    return 0;
}

其实，自己不会的还有很多，要学会坚持，这篇是为自己更好的理解LCA而写的，走了