poj 1470（简单LCA 倍增法）

题意：给你一个树，有若干个询问，然后让你统计每个结点在询问中做了几次LCA。按照结点顺序输出。

思路：这也是简单的LCA题目，我用的是倍增法。每次查询在相应结点标记上++，最后输出即可。这道题的输入处理比较烦，而且第一个输入的结点并不是根节点。这要注意一下

代码如下：

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>
#include <algorithm>
#include <queue>
#include <stack>
#include <vector>
#include <set>
#include <map>
#define MP(a, b) make_pair(a, b)
#define PB(a) push_back(a)

using namespace std;

typedef long long ll;
typedef pair<int ,int> pii;
typedef pair<unsigned int, unsigned int> puu;
typedef pair<int ,double> pid;
typedef pair<ll, int> pli;

const int INF = 0x3f3f3f3f;
const double eps = 1e-6;
const int LEN = 1000+10;
const int LOG_LEN = 100;
vector<int> Map[LEN];
int root, parent[LOG_LEN][LEN], depth[LEN], ind[LEN];

void dfs(int v, int p, int d){
    parent[0][v] = p;
    depth[v] = d;
    for(int i=0; i<Map[v].size(); i++){
        if(Map[v][i] != p) dfs(Map[v][i], v, d+1);
    }
}

void init(int n){
    dfs(root, -1, 0);
    for(int k=0; k+1<LOG_LEN; k++){
        for(int v=1; v<=n; v++){
            if(parent[k][v] < 0) parent[k+1][v] = -1;
            else parent[k+1][v] = parent[k][parent[k][v]];
        }
    }
}

int lca(int u, int v){
    if(depth[u] > depth[v])swap(u,v);
    for(int k=0; k < LOG_LEN; k++){
        if((depth[u] - depth[v]) >> k & 1) v = parent[k][v];
    }
    if(u==v) return u;
    for(int k=LOG_LEN-1; k>=0; k--){
        if(parent[k][v] != parent[k][u]){
            u = parent[k][u];
            v = parent[k][v];
        }
    }
    return parent[0][u];
}

int main()
{
//    freopen("in.txt", "r", stdin);

    int n, a, b, tn, q;
    while(scanf("%d", &n)!=EOF){
        for(int i=0; i<LEN; i++)Map[i].clear();
        memset(ind, 0, sizeof ind);
        for(int i=1; i<=n; i++){
            scanf("%d:(%d)", &a, &tn);
            for(int j=0; j<tn; j++){
                scanf("%d", &b);
                Map[a].PB(b);
                Map[b].PB(a);
                ind[b]++;
            }
        }
        for(int i=1; i<=n; i++)if(!ind[i]){root = i;break;}
        memset(ind, 0, sizeof ind);
        init(n);
        scanf("%d", &q);
        for(int i=0; i<q; i++){
            while(getchar() != '(');
            scanf("%d%d", &a, &b);
            while(getchar() != ')');
            ind[lca(a, b)] ++;
        }
        for(int i=1; i<=n; i++){
            if(!ind[i])continue;
            printf("%d:%d
", i, ind[i]);
        }
    }
    return 0;
}