[bzoj3694]最短路

Description

给出一个\(n\)个点\(m\)条边的无向图,\(n\)个点的编号从\(1-n\),定义源点为\(1\).

定义最短路树如下:从源点\(1\)经过边集\(T\)到任意一点\(i\)有且仅有一条路径,且这条路径是整个图\(1\)到\(i\)的最短路径,边集\(T\)构成最短路树.

给出最短路树,求对于除了源点\(1\)外的每个点\(i\),求最短路,要求不经过给出的最短路树上的\(1\)到\(i\)的路径的最后一条边.

Input

第一行包含两个数\(n\)和\(m\),表示图中有\(n\)个点和\(m\)条边.

接下来\(m\)行,每行有四个数\(a_i,b_i,l_i,t_i\),表示图中第\(i\)条边连接\(a_i\)和\(b_i\)权值为\(l_i,t_i\)为\(1\)表示这条边是最短路树上的边,\(t_i\)为\(0\)表示不是最短路树上的边.

Output

输出\(n-1\)个数,第\(i\)个数表示从\(1\)到\(i+1\)的要求的最短路.无法到达输出\(-1\).

Sample Input

5 9 
3 1 3 1 
1 4 2 1 
2 1 6 0 
2 3 4 0 
5 2 3 0 
3 2 2 1 
5 3 1 1 
3 5 2 0 
4 5 4 0

Sample Output

6 7 8 5

HINT

\(n\;\leq\;4000,m\;\leq\;100000，1\;\leq\;l_i\;\leq\;100000\)

Solution

对于一条不在最短路树上的边\(e[u][v]\),显然它只对\((u,v)\)上除了\(lca(u,v)\)之外的点有影响.

对于\((lca(u,v),v)\)上除了\(lca(u,v)\)之外的点\(k\),一条可行路径为\(1->u->v->k\).

(\(lca(u,v)\)不合法,不属于\((u,v)\)的点显然更劣.)

设\(dis[i]\)为最短路树中\((1,i)\)的长度,

则这条可行路径\(1->u->v->k\)长度为\(dis[u]+e[u][v]+dis[v]-dis[k]\).

\(dis[k]\)为定值,则要求合法的最小的\(dis[u]+e[u][v]+dis[v]\).

这个可以用树链剖分+线段树维护.

#include<cmath>
#include<ctime>
#include<queue>
#include<stack>
#include<cstdio>
#include<vector>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#define N 4005
#define M 100005
#define INF 900000000
using namespace std;
struct linetree{
    int l,r,m;
}lt[M];
struct graph{
    int nxt,to,w;
}e[M<<1];
struct line{
    int l,r,w;
}a[M];
int g[N],dis[N],n,m,tot,cnt;
int f[N],p[N],dep[N],siz[N],son[N],top[N];
inline void addedge(int x,int y,int w){
    e[++cnt].nxt=g[x];g[x]=cnt;e[cnt].to=y;e[cnt].w=w;
}
inline void dfs1(int u){
    int m=0;siz[u]=1;
    for(int i=g[u];i;i=e[i].nxt)
        if(!dep[e[i].to]){
            f[e[i].to]=u;
            dep[e[i].to]=dep[u]+1;
            dis[e[i].to]=dis[u]+e[i].w;
            dfs1(e[i].to);
            siz[u]+=siz[e[i].to];
            if(siz[e[i].to]>m){
                son[u]=e[i].to;
                m=siz[e[i].to];
            }
        }
}
inline void dfs2(int u,int tp){
    top[u]=tp;p[u]=++cnt;
    if(son[u]) dfs2(son[u],tp);
    for(int i=g[u];i;i=e[i].nxt)
        if(e[i].to!=f[u]&&e[i].to!=son[u])
            dfs2(e[i].to,e[i].to);
}

inline void build(int u,int l,int r){
    lt[u].l=l;lt[u].r=r;lt[u].m=INF;
    if(lt[u].l<lt[u].r){
        int mid=(lt[u].l+lt[u].r)>>1;
        build(u<<1,l,mid);build(u<<1|1,mid+1,r);
    }
}
inline void cover(int u,int l,int r,int k){
    if(lt[u].l>=l&&lt[u].r<=r){
        lt[u].m=min(lt[u].m,k);
    }
    else if(lt[u].l<lt[u].r){
        int lef=u<<1,rig=u<<1|1;
        int mid=(lt[u].l+lt[u].r)>>1;
        lt[lef].m=min(lt[u].m,lt[lef].m);
        lt[rig].m=min(lt[u].m,lt[rig].m);
        if(l<=mid) cover(lef,l,r,k);
        if(r>mid) cover(rig,l,r,k);
    }
}
inline int ask(int u,int x){
    if(lt[u].l<lt[u].r){
        int lef=u<<1,rig=u<<1|1;
        int mid=(lt[u].l+lt[u].r)>>1;
        lt[lef].m=min(lt[u].m,lt[lef].m);
        lt[rig].m=min(lt[u].m,lt[rig].m);
        if(x<=mid) return ask(lef,x);
        return ask(rig,x);
    }
    return lt[u].m;
}
inline void cov(int u,int v,int k){
    int t;
    while(top[u]!=top[v]){
        if(dep[top[u]]<dep[top[v]]){
            t=u;u=v;v=t;
        }
        cover(1,p[top[u]],p[u],k);
        u=f[top[u]];
    }
    if(p[u]>p[v]){
        t=u;u=v;v=t;
    }
    if(u!=v) cover(1,p[u]+1,p[v],k);
}
inline void Aireen(){
    scanf("%d%d",&n,&m);
    for(int i=1,j,k,w,t;i<=m;++i){
        scanf("%d%d%d%d",&j,&k,&w,&t);
        if(t){
            addedge(j,k,w);addedge(k,j,w);
        }
        else{
            a[++tot].l=j;a[tot].r=k;a[tot].w=w;
        }
    }
    dep[1]=1;dfs1(1);cnt=0;dfs2(1,1);
    build(1,1,n);m=tot;
    for(int i=1;i<=m;++i)
        cov(a[i].l,a[i].r,a[i].w+dis[a[i].l]+dis[a[i].r]);
    for(int i=2,k;i<=n;++i){
        k=ask(1,p[i]);
        if(k==INF) printf("-1 ");
        else printf("%d ",k-dis[i]);
    }
}
int main(){
    freopen("dis.in","r",stdin);
    freopen("dis.out","w",stdout);
    Aireen();
    fclose(stdin);
    fclose(stdout);
    return 0;
}