bnuoj52303 Floyd-Warshall(2016四川省赛)(生成树+倍增+SPFA)

题目链接：Floyd-Warshall

题意：

给你n个点，m条双向边，再给你q个询问，每个询问输入x,y，叫你求x到y的最短距离。

数据范围：

$1\leq n\leq 10^5,0<m-n<100,q\leq 10^5$，保证是个联通图。

题解：

1、我们先利用并查集建成一个生成树，在当前情况下，询问$x,y$的最短距离，我们利用倍增求出$x,y$的lca为$z$,他们的距离为$deep[x]+deep[y]-2*deep[z],deep[i]$表示为到1节点的距离。

2、我们再把剩下的$m-n+1$条非树边加入到树上，添加完成后，我们再对所有非树边的两个端点对剩余的所有点跑SPFA(边数较少，跑spfa更快)。时间复杂度为$O(200*m)$。

3、针对每个询问，如果答案会更优的话，我们肯定要经过非树边，所以再枚举非树边的端点，求其最短距离相加即可，最短距离已经预处理出来，所以时间复杂度为$O(q*200)$。

代码：

#include<cstdio>
#include<cstring>
#include<iostream>
#include<sstream>
#include<algorithm>
#include<vector>
#include<bitset>
#include<set>
#include<queue>
#include<stack>
#include<map>
#include<cstdlib>
#include<cmath>
#define PI 2*asin(1.0)
#define LL long long
#define pb push_back
#define pa pair<int,int>
#define clr(a,b) memset(a,b,sizeof(a))
#define lson lr<<1,l,mid
#define rson lr<<1|1,mid+1,r
#define bug(x) printf("%d++++++++++++++++++++%d\n",x,x)
#define key_value ch[ch[root][1]][0]
const int  MOD = 1000000007;
const int N = 1e5 + 1000;
const int maxn = 100+ 14;
const int letter = 130;
const int INF = 1e9;
const double pi=acos(-1.0);
const double eps=1e-8;
using namespace std;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
int n,m,q,father[N];
int u[N],v[N],nu[N],nv[N],fa[N][20],deep[N],bin[20];
int dis[220][N],ps[220],num,cnt;
vector<int>G[N];
queue<int>pq;
int Find(int x){
    if(x==father[x]) return x;
    return father[x]=Find(father[x]);
}
void init(){
    cnt=num=0;
    for(int i=1;i<=n;i++) G[i].clear(),father[i]=i;
}
void dfs(int x,int pre){
    for(int i=0;i<G[x].size();i++){
        int to=G[x][i];
        if(to==pre) continue;
        fa[to][0]=x;
        deep[to]=deep[x]+1;
        dfs(to,x);
    }
}
int lca(int x,int y){
    if(deep[x]<deep[y]) swap(x,y);
    int tmp=deep[x]-deep[y];
    for(int i=0;bin[i]<=tmp;i++){
        if(tmp&bin[i]) x=fa[x][i];
    }
    for(int i=18;i>=0;i--){
        if(fa[x][i]!=fa[y][i]) x=fa[x][i],y=fa[y][i];
    }
    if(x==y) return x;
    return fa[x][0];
}
void dij(int s){
    for(int i=1;i<=n;i++) dis[s][i]=INF;
    dis[s][ps[s]]=0;
    while(!pq.empty()) pq.pop();
    pq.push(ps[s]);
    while(!pq.empty()){
        int x=pq.front();
        pq.pop();
        for(int i=0;i<G[x].size();i++){
            int to=G[x][i];
            if(dis[s][to]>dis[s][x]+1){
                dis[s][to]=dis[s][x]+1;
                pq.push(to);
            }
        }
    }
}
int main(){
    bin[0]=1;
    for(int i=1;i<=18;i++) bin[i]=bin[i-1]<<1;
    int x,y,z;
    while(scanf("%d%d%d",&n,&m,&q)!=EOF){
        init();
        for(int i=0;i<m;i++){
            scanf("%d%d",u+i,v+i);
        }
        for(int i=0;i<m;i++){
            x=u[i],y=v[i];
            int xx=Find(x),yy=Find(y);
            if(xx!=yy){
                father[yy]=xx;
                G[x].pb(y),G[y].pb(x);
            }
            else nu[cnt]=x,nv[cnt++]=y,ps[num++]=x,ps[num++]=y;
        }
        sort(ps,ps+num);
        num=unique(ps,ps+num)-ps;
        deep[1]=0;
        dfs(1,-1);
        for(int j=1;j<=18;j++)
            for(int i=1;i<=n;i++){
                fa[i][j]=fa[fa[i][j-1]][j-1];
            }
        for(int i=0;i<cnt;i++){
            G[nu[i]].pb(nv[i]),G[nv[i]].pb(nu[i]);
        }
        for(int i=0;i<num;i++) dij(i);
        for(int i=1;i<=q;i++){
            scanf("%d%d",&x,&y);
            z=lca(x,y);
            int min1=deep[x]+deep[y]-2*deep[z];
            for(int j=0;j<num;j++) min1=min(min1,dis[j][x]+dis[j][y]);
            printf("%d\n",min1);
        }
    }
    return 0;
}