【BZOJ-3910】火车 倍增LCA + 并查集

3910: 火车

Time Limit: 20 Sec  Memory Limit: 512 MB
Submit: 262  Solved: 90
[Submit][Status][Discuss]

Description

A 国有n 个城市，城市之间有一些双向道路相连，并且城市两两之间有唯一路径。现在有火车在城市 a，需要经过m 个城市。火车按照以下规则行驶：每次行驶到还没有经过的城市中在 m 个城市中最靠前的。现在小 A 想知道火车经过这m 个城市后所经过的道路数量。 

Input

第一行三个整数 n、m、a，表示城市数量、需要经过的城市数量，火车开始时所在位置。 

接下来 n-1 行，每行两个整数 x和y，表示 x 和y之间有一条双向道路。 

接下来一行 m 个整数，表示需要经过的城市。 

Output

一行一个整数，表示火车经过的道路数量。 

Sample Input

5 4 2
1 2
2 3
3 4
4 5
4 3 1 5

Sample Output

9

HINT

N<=500000 ,M<=400000 

Source

Solution

水题- -最多算个并查集的有趣应用

很显然直接询问用LCA统计答案即可

至于处理走过的路径，拿并查集维护一下，很简单的把起止点到LCA的点合并一下，询问的两个点如果属于一个集合显然走过

Code

#include<iostream>
#include<cmath>
#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
int read()
{
    int x=,f=; char ch=getchar();
    while (ch<'' || ch>'') {if (ch=='-') f=-; ch=getchar();}
    while (ch>='' && ch<='') {x=x*+ch-''; ch=getchar();}
    return x*f;
}
#define maxn 501000
int n,m,a;
struct EdgeNode{int to,next;}edge[maxn<<];
int head[maxn],cnt;
void add(int u,int v) {cnt++; edge[cnt].next=head[u]; head[u]=cnt; edge[cnt].to=v;}
void insert(int u,int v) {add(u,v); add(v,u);}
int deep[maxn],father[maxn][],ffff[maxn];
long long ans;
void dfs(int now)
{
    for (int i=; i<=; i++)
        if (deep[now]>=(<<i))
            father[now][i]=father[father[now][i-]][i-];
        else
            break;
    for (int i=head[now]; i; i=edge[i].next)
        if (edge[i].to!=father[now][])
            {
                deep[edge[i].to]=deep[now]+;
                father[edge[i].to][]=now;
                dfs(edge[i].to);
            }
}
int LCA(int x,int y)
{
    if (deep[x]<deep[y]) swap(x,y);
    int dd=deep[x]-deep[y];
    for (int i=; (<<i)<=dd; i++)
        if (dd&(<<i)) x=father[x][i];
    for (int i=; i>=; i--)
        if (father[x][i]!=father[y][i])
            x=father[x][i],y=father[y][i];
    if (x==y) return x;
    return father[x][];
}
int find(int x) {if (ffff[x]==x) return x; ffff[x]=find(ffff[x]); return ffff[x];}
int ff1,ff2;
int main()
{
    n=read(),m=read(),a=read();
    for (int u,v,i=; i<=n-; i++)
        u=read(),v=read(),insert(u,v);
    dfs();
    for (int i=; i<=n; i++) ffff[i]=i;
    for (int i=; i<=m; i++)
        {
            int x=read();
            int fa=find(a),fx=find(x);
            if (fa==fx) continue;
            int lca=LCA(a,x);
            ans+=deep[a]-deep[lca]+deep[x]-deep[lca];
            int ta=a,tx=x,flca; flca=find(lca);
            while (find(ta)!=flca) {ff1=find(ta); ffff[ff1]=flca; ta=father[ff1][];}
            while (find(tx)!=flca) {ff2=find(tx); ffff[ff2]=flca; tx=father[ff2][];}
            a=x;
        }
    cout<<ans;
    return ;
}

电脑炸出奇怪的错误，所以写的比较鬼畜- -