HDU 2874 Connections between cities 离线LCA

题目链接

http://acm.hdu.edu.cn/showproblem.php?pid=2874

题意

给定一个深林，问每两个点之间的最短距离。

思路

与HDU2586一样的套路,只是这要判断是不是在同一棵树内。

#include<cstdio>
#include<queue>
#include<iostream>
#include<vector>
#include<map>
#include<cstring>
#include<string>
#include<set>
#include<stack>
#include<algorithm>
#define cle(a) memset(a,0,sizeof(a))
#define inf(a) memset(a,0x3f,sizeof(a))
#define ll long long
#define Rep(i,a,n) for(int i=a;i<=n;i++)
using namespace std;
#define INF2 9223372036854775807ll
const int INF = ( 2e9 ) + 2;
const ll maxn = 1e4+10;
int n,m,c;
struct edge
{
    int v,w,next;
}e[maxn*2];
struct Q
{
    int v,id,next;
}que[1000010*2];
int ans[1000010];
int head[maxn],head2[maxn],vis[maxn],scc[maxn],f[maxn],d[maxn];
int tot,tot2,time,num;
stack<int> s;
int Find(int x)
{
    return f[x]==x?x:f[x]=Find(f[x]);
}
void addedge(int u,int v,int w)
{
    e[tot].v=v;
    e[tot].w=w;
    e[tot].next=head[u];
    head[u]=tot++;
    e[tot].v=u;
    e[tot].w=w;
    e[tot].next=head[v];
    head[v]=tot++;
}
void init(int n,int c)
{
    tot=tot2=num=0;
    memset(head,-1,sizeof(head));
    memset(head2,-1,sizeof(head2));
    memset(vis,0,sizeof(vis));
    memset(scc,0,sizeof(scc));
    memset(d,0,sizeof(d));
    memset(ans,-1,sizeof(ans));
    for(int i=1;i<=n;i++)f[i]=i;
}
void dfs(int u,int num)
{
    f[u]=u;
    vis[u]=1;
    scc[u]=num;
    for(int i=head2[u];i!=-1;i=que[i].next)
    {
        int v=que[i].v;
        int id=que[i].id;
        if(vis[v]&&scc[v]==scc[u])
        ans[id]=d[u]+d[v]-d[Find(v)]*2;
    }
    for(int i=head[u];i!=-1;i=e[i].next)
    {
        int v=e[i].v;
        if(!vis[v])
        {
            d[v]=d[u]+e[i].w;
            dfs(v,num);
            f[v]=u;
        }
    }
}
int main()
{
    while(~scanf("%d%d%d",&n,&m,&c))
    {
        init(n,c);
        for(int i=0;i<m;i++)
        {
            int u,v,w;
            scanf("%d%d%d",&u,&v,&w);
            addedge(u,v,w);
        }
        for(int i=0;i<c;i++)
        {
            int u,v;
            scanf("%d%d",&u,&v);
            que[tot2].v=v;
            que[tot2].id=i;
            que[tot2].next=head2[u];
            head2[u]=tot2++;
            que[tot2].v=u;
            que[tot2].id=i;
            que[tot2].next=head2[v];
            head2[v]=tot2++;
        }

        for(int i=1;i<=n;i++)
        {
            if(!vis[i])
            {
                num++;
                dfs(i,num);
            }
        }

        for(int i=0;i<c;i++)
        if(ans[i]!=-1)
        printf("%d\n",ans[i]);
        else
        printf("Not connected\n");
    }
}