[SDOI2009]Elaxia的路线

P2149 [SDOI2009]Elaxia的路线

求无向图中，两对点间最短路的最长公共路径

喵啊…这题真心喵啊…orzorz

先spfa求出x1, y1, x2, y2的单源最短路

然后把x1到y1的最短路们建图

判断方式（精髓！）：

        if(dis[0][uu] + ww + dis[1][vv] == dis[0][y1]){
            addedge(uu, vv, ww, 1);
            ind[vv]++;
        }

蒟蒻欲膜又止【划掉

注意 建完是有向图哦

在建完图之后 用同样的方式判断这个新图上有哪些边在x2-y2的最短路上 并标记

注意这里有个坑点 就是这两位小彩虹不一定要并肩走的…他们满足于相互呼唤… (*￣︶￣* ) 

所以要判断一下正向和反向 只要有一个满足这条边就可以标记

那这里一段，那里一段怎么统计呢？

其实 同一条x2-y2最短路在新图上是连续的 因为如果中间有更短的地方 他们两个为什么不一起走另一条路呢？

所以最后按拓扑序bfs就好了

另 咱们能拓扑就别dfs 分分钟T爆……

int main(){
    init();
    spfa(x1, 0);
    spfa(y1, 1);
    spfa(x2, 2);
    spfa(y2, 3);
    int tmp = esize[0];
    for(int i = 1; i <= tmp; i++){
        int vv = edge[0][i].v, uu = edge[0][i].u, ww = edge[0][i].w;
        if(dis[0][uu] + ww + dis[1][vv] == dis[0][y1]){
            addedge(uu, vv, ww, 1);
            ind[vv]++;
        }
    }
    tmp = esize[1];
    for(int i = 1; i <= tmp; i++){
        int vv = edge[1][i].v, uu = edge[1][i].u, ww = edge[1][i].w;
        if(dis[2][uu] + ww + dis[3][vv] == dis[2][y2] || dis[2][vv] + ww + dis[3][uu] == dis[2][y2])
            flag[i] = 1;
    }
    bfs();
    printf("%d", f[y1]);
    return 0;
}

inline void bfs(){
    q.push(x1);
    while(!q.empty()){
        int fro = q.front(); q.pop();        
        for(int i = head[1][fro]; i != -1; i = edge[1][i].next){
            int vv = edge[1][i].v, ww = edge[1][i].w;
            f[vv] = max(f[vv], f[fro] + flag[i] * ww);
            ind[vv]--;
            if(!ind[vv]) q.push(vv);
        }
    }
}

inline void spfa(int s, int tp){
    q.push(s);
    dis[tp][s] = 0;
    for(int i = 1; i <= n; i++) vis[i] = 0;
    while(!q.empty()){
        int fro = q.front(); q.pop(); vis[fro] = 0;
        for(int i = head[0][fro]; i != -1; i = edge[0][i].next){
            int vv = edge[0][i].v, ww = edge[0][i].w;
            if(dis[tp][vv] > dis[tp][fro] + ww){
                dis[tp][vv] = dis[tp][fro] + ww;
                if(!vis[vv]){
                    q.push(vv);
                    vis[vv] = 1;
                }
            }
        }
    }
}