The Shortest Statement //最短路+lca

题目链接:https://codeforces.com/problemset/problem/1051/F

题意:一个图由若干顶点和若干边组成,且说边的个数减去顶点的个数不会多于二十,有很多组点对,问这些点对的最短距离为多少。
解析:询问次数较多,不能每次都跑dij,由于题目给出边的个数减去顶点的个数不会多于二十,说明如果以原图建dfs树,最多要删21条边,所以可以把这些多余的边所连顶点单独拿出来跑dij,然后找任意一个顶点为根建dfs树,那么对于任意一组询问,这个点对(u,v)都在dfs树上,需在树上找一个顶点x,使得(u,x)和(x,v)的距离之和为(u,v)的最短距离,很容易想到这个点就是(u,v)的lca,不妨设dis(i)为点i到根的距离,则(u,v)的最小距离即为dis(u)+dis(v)-2*dis(lca(u,v))。对于非dfs树上的边所涉及顶点可能存在使结果更小的情况,那样的话就暴力去查看是否存在这样的情况就行了(反正这样的顶点最多才42个),存在就更新答案。
#include <bits/stdc++.h>

using namespace std;
typedef long long ll;
typedef pair<ll,int> P;
#define ft first
#define sd second
const int maxn=1e5+10;
const ll inf=1e15;
int n,m;
ll d[50][maxn];
vector<P> g[maxn];
bool vis[maxn];
int qq[maxn],top=0;
ll dist[maxn];
int dep[maxn],f[maxn][22];
void dfs(int u,int fa){//简单树形dp求所有顶点深度和距离根的距离
    f[u][0]=fa;vis[u]=1;
    for(int i=0;i<(int)g[u].size();i++){
        int v=g[u][i].sd;ll w=g[u][i].ft;
        if(v==fa)continue;
        if(vis[v]){
            qq[top++]=u,qq[top++]=v;
        }
        else {
            dep[v]=dep[u]+1;
            dist[v]=dist[u]+w;
            dfs(v,u);
        }
    }
}
int lca(int u,int v){//求u,v两点的最近公共祖先
    if(dep[u]<dep[v])swap(u,v);
    for(int i=19;i>=0;i--)if(dep[f[u][i]]>=dep[v])u=f[u][i];
    if(u==v)return u;
    for(int j=19;j>=0;j--)
        if(f[u][j]!=f[v][j])
            u=f[u][j],v=f[v][j];
    return f[u][0];
}
void dijkstra(int s){//求qq[s]到所有顶点的最短距离
    for(int i=1;i<=n;i++)d[s][i]=inf;
    d[s][qq[s]]=0;
    priority_queue<P,vector<P>,greater<P> > pq;
    pq.push(P(0,qq[s]));
    while(!pq.empty()){
        P p=pq.top();pq.pop();
        int u=p.sd;
        if(d[s][u]<p.ft)continue;
        for(int i=0;i<(int)g[u].size();i++){
            P now=g[u][i];
            int v=now.sd;
            if(d[s][v]>d[s][u]+now.ft){
                d[s][v]=d[s][u]+now.ft;
                pq.push(P(d[s][v],v));
            }
        }
    }
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=0;i<m;i++){
        int u,v;ll w;
        scanf("%d%d%lld",&u,&v,&w);
        g[u].push_back(P(w,v));
        g[v].push_back(P(w,u));
    }
    dep[1]=1;
    dfs(1,0);
    for(int j=1;j<=19;j++)//倍增思想求顶点i的2^j祖先
         for(int i=1;i<=n;i++)
             f[i][j]=f[f[i][j-1]][j-1];
    sort(qq,qq+top);top=unique(qq,qq+top)-qq;//非dfs树上的边所连顶点个数
    for(int i=0;i<top;i++)dijkstra(i);//求这些点到其他点的最短距离
    int q;
    cin>>q;
    while(q--){
        int u,v;
        scanf("%d%d",&u,&v);
        ll ans=dist[u]+dist[v]-dist[lca(u,v)]*2;//这是在dfs树上u,v的最短距离
        for(int i=0;i<top;i++){
            ans=min(ans,d[i][u]+d[i][v]);//看是否存在不是dfs树上的边使得u,v距离更短
        }
        cout<<ans<<'\n';
    }
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值