bzoj3694

http://www.elijahqi.win/archives/858
Description
给出一个n个点m条边的无向图，n个点的编号从1~n，定义源点为1。定义最短路树如下：从源点1经过边集T到任意一点i有且仅有一条路径，且这条路径是整个图1到i的最短路径，边集T构成最短路树。 给出最短路树，求对于除了源点1外的每个点i，求最短路，要求不经过给出的最短路树上的1到i的路径的最后一条边。

Input
第一行包含两个数n和m，表示图中有n个点和m条边。
接下来m行，每行有四个数ai，bi，li，ti，表示图中第i条边连接ai和bi权值为li，ti为1表示这条边是最短路树上的边，ti为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

对于100%的数据，n≤4000，m≤100000，1≤li≤100000

这题的做法大概真的好神奇吧

首先我们把已经存在属于树上的点存成图，不属于树上的点，单独另存

对于一条不在最短路树的有向边（无向可看成两条有向）u-v，长度L，设t=lca(u,v)

那么对于t-v的路径上所有点x，都可通过1-t-u-v-x

所以总长度等于d[u]+d[v]+L(新边的长度)-d[x]

#include<cstdio>
#include<algorithm>
#define N 110000
#define ll long long
#define inf 0x3f3f3f3f
using namespace std;
inline ll read(){
    ll x=0,f=1;char ch=getchar();
    while (ch<'0'||ch>'9') {if (ch=='-')f=-1;ch=getchar();}
    while (ch<='9'&&ch>='0'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
struct node{
    int z,y,next;
}data[N<<1];
struct node1{
    int x,y,z;
}data1[N<<1];
struct node2{
    int l,r,left,right,min,lazy;
}tree[N<<2];
int size[N],h[N],id[N],dep[N],son[N],tp[N],fa[N],d[N],num,root,n,m,tot;
void dfs1(int x){
    size[x]=1;
    for (int i=h[x];i;i=data[i].next){
        int y=data[i].y,z=data[i].z;
        if (fa[x]==y) continue;
        dep[y]=dep[x]+1;fa[y]=x ;d[y]=d[x]+z;dfs1(y);size[x]+=size[y];
        if (size[y]>size[son[x]]) son[x]=y;
    }
}
void dfs2(int x,int top){
    id[x]=++num;tp[x]=top;
    if (son[x]) dfs2(son[x],top);
    for (int i=h[x];i;i=data[i].next){
        int y=data[i].y;
        if (fa[x]==y||son[x]==y) continue;
        dfs2(y,y);
    }
}
inline void update(int x){
    int l=tree[x].left,r=tree[x].right;
    tree[x].min=min(tree[l].min,tree[r].min);
}
void build(int &x,int l,int r){
    x=++num;tree[x].l=l;tree[x].r=r;tree[x].lazy=inf;
    if (l==r){tree[x].min=inf;return;}
    int mid=l+r>>1;
    build(tree[x].left,l,mid);build(tree[x].right,mid+1,r);
    update(x);
}
inline void pushdown(int x){
    if (tree[x].lazy==inf) return;
    int l=tree[x].left,r=tree[x].right;
    tree[l].min=min(tree[l].min,tree[x].lazy);
    tree[r].min=min(tree[r].min,tree[x].lazy);
    tree[l].lazy=min(tree[l].lazy,tree[x].lazy);
    tree[r].lazy=min(tree[r].lazy,tree[x].lazy);tree[x].lazy=inf;
}
int query(int x,int l){
    if (tree[x].l==tree[x].r) return tree[x].min;
    int mid=tree[x].l+tree[x].r>>1;pushdown(x);
    if (l<=mid) return query(tree[x].left,l);else return query(tree[x].right,l);
}
void change(int x,int l,int r,int v){
    if (l<=tree[x].l&&r>=tree[x].r){tree[x].min=min(v,tree[x].min);tree[x].lazy=min(tree[x].lazy,v);return;}
    int mid=tree[x].l+tree[x].r>>1;pushdown(x);
    if (l<=mid) change(tree[x].left,l,r,v);
    if (r>mid) change(tree[x].right,l,r,v);
    update(x);
}
int main(){
    freopen("bzoj3694.in","r",stdin);
    n=read();m=read();num=1;tot=0;
    for (int i=1;i<=m;++i){
        int x=read(),y=read(),v=read(),z=read();
        if (z){
            data[++num].y=y;data[num].next=h[x];data[num].z=v;h[x]=num;
            data[++num].y=x;data[num].next=h[y];data[num].z=v;h[y]=num;
        }else{
            data1[++tot].x=x;data1[tot].y=y;data1[tot].z=v;
        }
    }dep[1]=1;num=0;
    dfs1(1);dfs2(1,1);num=0;
    build(root,1,n);
    for (int i=1;i<=tot;++i){
        int x=data1[i].x,y=data1[i].y,z=d[x]+d[y]+data1[i].z;
        while (tp[x]!=tp[y]){
            if (dep[tp[x]]<dep[tp[y]]) swap(x,y);
            change(root,id[tp[x]],id[x],z);
            x=fa[tp[x]];
        }
        if (id[x]>id[y])swap(x,y);if(id[x]==id[y]) continue;
        change(root,id[x]+1,id[y],z);
    }
    //for (int i=1;i<=n;++i) printf("%d ",id[i]);printf("\n");
    for (int i=2;i<=n;++i){
        int tmp=query(root,id[i]);
        if (tmp==inf) printf("-1 ");else printf("%d ",tmp-d[i]); 
    }
    return 0;
}