BZOJ2125 最短路

题目链接
先构造圆方树,然后圆-圆边的距离就是直接距离,圆-方边的距离就是圆点到环的起始点的距离.
然后在圆方树上倍增就好了.
注意如果\(LCA\)是方点就往下跳到圆点,然后做一个环上最短路即可.

代码如下

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<vector>
#define N (30010)
#define M (400010)
#define inf (0x7f7f7f7f)
#define rg register int
#define Label puts("NAIVE")
#define spa print(' ')
#define ent print('\n')
#define rand() (((rand())<<(15))^(rand()))
typedef long double ld;
typedef long long LL;
typedef unsigned long long ull;
using namespace std;
inline char read(){
    static const int IN_LEN=1000000;
    static char buf[IN_LEN],*s,*t;
    return (s==t?t=(s=buf)+fread(buf,1,IN_LEN,stdin),(s==t?-1:*s++):*s++);
}
template<class T>
inline void read(T &x){
    static bool iosig;
    static char c;
    for(iosig=false,c=read();!isdigit(c);c=read()){
        if(c=='-')iosig=true;
        if(c==-1)return;
    }
    for(x=0;isdigit(c);c=read())x=((x+(x<<2))<<1)+(c^'0');
    if(iosig)x=-x;
}
inline char readchar(){
    static char c;
    for(c=read();!isalpha(c);c=read())
    if(c==-1)return 0;
    return c;
}
const int OUT_LEN = 10000000;
char obuf[OUT_LEN],*ooh=obuf;
inline void print(char c) {
    if(ooh==obuf+OUT_LEN)fwrite(obuf,1,OUT_LEN,stdout),ooh=obuf;
    *ooh++=c;
}
template<class T>
inline void print(T x){
    static int buf[30],cnt;
    if(x==0)print('0');
    else{
        if(x<0)print('-'),x=-x;
        for(cnt=0;x;x/=10)buf[++cnt]=x%10+48;
        while(cnt)print((char)buf[cnt--]);
    }
}
inline void flush(){fwrite(obuf,1,ooh-obuf,stdout);}
vector<int> G[N],D[N]; bool ins[N];
int fi[N],b[M],ne[M],c[M],E,f[N][18],dis[N],fa[N],qq[N];
int n,m,Q,st[N],dfn[N],low[N],ind,top,cnt,dep[N],ab[N],cir[N];
void add(int x,int y,int z){
    ne[++E]=fi[x],fi[x]=E,b[E]=y,c[E]=z;
    ne[++E]=fi[y],fi[y]=E,b[E]=x,c[E]=z;
}
void dfs(int u,int pre){
    f[u][0]=pre;
    for(int i=fi[u];i;i=ne[i])
    if(b[i]!=pre){
        int v=b[i];
        dis[v]=dis[u]+c[i];
        dep[v]=dep[u]+1,dfs(v,u);
    }
}
void work(int u,int s,int d){
    qq[u]=0,cir[++cnt]=d;
    for(int i=u;i!=s;i=fa[i])
    qq[fa[i]]=qq[i]+ab[i],cir[cnt]+=ab[i];
    for(int i=u;i!=s;i=fa[i])
    add(cnt,i,min(qq[s]-qq[i],qq[i]+d));
    add(cnt,s,0);
}
void tarjan(int u,int pre){
    dfn[u]=low[u]=++ind;
    ins[u]=1,fa[u]=pre;
    for(int i=0;i<G[u].size();i++){
        int v=G[u][i],L=D[u][i];
        if(!dfn[v])ab[v]=L,tarjan(v,u),
        low[u]=min(low[u],low[v]);
        else if(v!=pre)low[u]=min(low[u],dfn[v]);
        if(low[v]>dfn[u])add(u,v,L);
        if(fa[v]!=u&&low[v]>=dfn[u])work(v,u,L);
    }
}
void bz(){
    for(int j=1;j<=16;j++)
    for(int i=1;i<=cnt;i++)
    f[i][j]=f[f[i][j-1]][j-1];
}
int lca(int x,int y){
    if(dep[x]<dep[y])swap(x,y);
    int k=dep[x]-dep[y];
    for(int i=0;i<=16;i++)
    if(k&(1<<i))x=f[x][i];
    if(x==y)return x;
    for(int i=16;~i;i--)
    if(f[x][i]!=f[y][i])
    x=f[x][i],y=f[y][i];
    return (x==y)?x:f[x][0];
}
int findc(int x,int rt){
    int k=dep[x]-dep[rt]-1;
    for(int i=0;i<=16;i++)
    if(k&(1<<i))x=f[x][i];
    return x;
}
int getdis(int x,int y){
    int rt=lca(x,y);
    if(rt<=n)return dis[x]+dis[y]-2*dis[rt];
    int tx=findc(x,rt),ty=findc(y,rt);
    int len=dis[x]-dis[tx]+dis[y]-dis[ty];
    int L=abs(qq[tx]-qq[ty]);
    return len+min(L,cir[rt]-L);
}
int main(){
    read(n),read(m),read(Q),cnt=n;
    for(int i=1,x,y,z;i<=m;i++)
    read(x),read(y),read(z),
    G[x].push_back(y),G[y].push_back(x),
    D[x].push_back(z),D[y].push_back(z);
    E=0,tarjan(1,0),dfs(1,0),bz();
    while(Q--){
        int x,y;
        read(x),read(y);
        printf("%d\n",getdis(x,y));
    }
}

转载于:https://www.cnblogs.com/Romeolong/p/10057397.html