CF1051F The Shortest Statement

Luogu

题意

给了你一张$n$个点$m$条边的无向图,$m$最多比$n$大$20$,有$q$次询问,每次询问给定$u,v$两个点,问两点间的最短路;

题解

我的做法非常暴力,绝对不是正解,但是在四秒时限的宽松条件下跑过了;
考虑$m$只比$n$大$20$的图有什么特点,当我们求出一颗生成树之后,没用的边最多也只有$21$条,受影响的点最多也只有$42$个,我们把这些受影响的点提出来跑一次$Floyd$得到两两之间的最短路径,考虑每次询问的两个点的最短路长啥样,必然只有这么两种情况:$1.U\to V,2.U\to 受影响的点a\to 受影响的点b\to V$;那么思路就很清晰了,每次询问$42\ast 42$的去暴力枚举中间的两个点$a,b$,总距离就是树上距离$[U][a]+$受影响点之间的最短距离$[a][b]+$树上距离$[b][V]$,然后再在这些距离与树上距离$[U][V]$中取个$min$就好了;为了跑的更快选择使用差分求树上距离以及$RMQ$求$LCA$;

#include<bits/stdc++.h>
#define mem(a,b) memset(a,b,sizeof(a))
using namespace std;
typedef long long LL;
typedef unsigned int UI;
typedef unsigned long long ULL;
template<typename T> inline void read(T& x) {
    char c = getchar();
    bool f = false;
    for (x = 0; !isdigit(c); c = getchar()) {
        if (c == '-') {
            f = true;
        }
    }
    for (; isdigit(c); c = getchar()) {
        x = x * 10 + c - '0';
    }
    if (f) {
        x = -x;
    }
}
template<typename T, typename... U> inline void read(T& x, U& ... y) {
    read(x), read(y...);
}
const int N=1e5+10,INF=0x7fffffff;
int n,m,p,Q,cnt,TIME;
int head[N],vis[N],no[N],F[N],V[N<<1],dfn[N],dep[N],ST[N<<1][20],Log[N<<1],G[50];
LL dis[50][50],D[N];
struct Edge {
    int to,last,w;
    Edge () {}
    Edge (int a,int b,int c) :to(a),last(b),w(c) {}
}edge[N<<1];
void ADD(int a,int b,int c) {
    edge[++p]=Edge(b,head[a],c); head[a]=p;
    edge[++p]=Edge(a,head[b],c); head[b]=p;
}
int Find(int a) {
    return a==F[a]?a:F[a]=Find(F[a]);
}
void DFS(int u,int fa) {
    dfn[u]=++TIME; V[TIME]=u;
    for(int i=head[u];i;i=edge[i].last) {
        int v=edge[i].to;
        if(v!=fa) {
            dep[v]=dep[u]+1; D[v]=D[u]+edge[i].w;
            DFS(v,u);
            V[++TIME]=u;
        }
    }
}
int LCA(int u,int v) {
    if(dfn[u]>dfn[v]) swap(u,v);
    int k=Log[dfn[v]-dfn[u]+1],ANS;
    if(dep[ST[dfn[u]][k]]<dep[ST[dfn[v]-(1<<k)+1][k]]) ANS=ST[dfn[u]][k];
    else ANS=ST[dfn[v]-(1<<k)+1][k];
    return ANS;
}
LL Getdis(int u,int v) {
    int lca=LCA(u,v);
    return D[u]+D[v]-(D[lca]<<1);
}
int main() {
    read(n,m);
    for(int i=1;i<=n;++i) F[i]=i;
    memset(dis,0x3f,sizeof(dis));
    for(int i=1;i<=m;++i) {
        int u,v,w; read(u,v,w);
        int x=Find(u),y=Find(v);
        if(x!=y) {
            F[y]=x;
            ADD(u,v,w);
        }
        else {
            if(!vis[u]) vis[u]=true,no[u]=++cnt,G[cnt]=u;
            if(!vis[v]) vis[v]=true,no[v]=++cnt,G[cnt]=v;
            dis[no[u]][no[v]]=dis[no[v]][no[u]]=min(dis[no[u]][no[v]],1ll*w);
        }
    }
    DFS(1,0);
    for(int i=1;i<=TIME;++i) ST[i][0]=V[i];
    for(int i=1;i<=19;++i) {
        for(int j=1;j+(1<<i)-1<=TIME;++j) {
            if(dep[ST[j][i-1]]<dep[ST[j+(1<<(i-1))][i-1]]) ST[j][i]=ST[j][i-1];
            else ST[j][i]=ST[j+(1<<(i-1))][i-1];
        }
    }
    Log[0]=-1;
    for(int i=1;i<=TIME;++i) Log[i]=Log[i>>1]+1;
    for(int i=1;i<=cnt;++i) {
        for(int j=i+1;j<=cnt;++j) {
            int lca=LCA(G[i],G[j]);
            dis[i][j]=dis[j][i]=min(dis[i][j],D[G[i]]+D[G[j]]-D[lca]*2);
        }
    }
    for(int k=1;k<=cnt;++k) {
        for(int i=1;i<=cnt;++i) {
            for(int j=1;j<=cnt;++j) {
                dis[i][j]=dis[j][i]=min(dis[i][j],dis[i][k]+dis[k][j]);
            }
        }
    }
    scanf("%d",&Q);
    while(Q--) {
        int u,v; read(u,v);
        LL ANS=Getdis(u,v);
        for(int i=1;i<=cnt;++i)
            for(int j=i+1;j<=cnt;++j) {
                ANS=min(ANS,Getdis(u,G[i])+Getdis(v,G[j])+dis[i][j]);
                ANS=min(ANS,Getdis(u,G[j])+Getdis(v,G[i])+dis[i][j]);
            }
        printf("%lld\n",ANS);
    }
    return 0;
}