CF 787D 线段树套堆优化Dij

题目链接：http://codeforces.com/contest/786/problem/B
题意：就是给了n个点，给了m个加边的关系，有向边，但是加边不一定是u->v，可能是u->[l, r]就是u到[l, r]区间里面的每一个点都加一条边，或者是[l, r] -> u，就是[l, r]区间到u加一条边。然后求单源最短路。

解法:显然单源最短路用堆优化的DIJ来求。先考虑一下，如果我们把区间拆成点来加边，那么复杂度显然会高达O(n^2logn)的。所以显然不能这样来加边，这里给了一堆线段？我们可以想怎么把线段和点连接起来呢？显然线段树可以轻易的办到，由于是有向边所以这里需要两颗线段树，线段树的节点是2*n的，2颗就是4*n，j加上原来的n个点，所以节点数是5n。然后我们加边原来最坏的O(n)，现在可以变成logn了，那么我们的复杂度可以 O(nlognlogn)了。

具体讲解可以看lwt大佬的视频：http://www.bilibili.com/video/av9365298/

//CF 787D

#include <bits/stdc++.h>
using namespace std;
const int maxn = 1e5+10;
const int maxv = maxn*5;
int n, m, s;
vector <pair<int, int>> G[maxv];
void addedge(int u, int v, int w){
    G[u].push_back(make_pair(v, w));
}
int id[2][maxn<<2], idx;
void build(int l, int r, int rt, int wh)
{
    id[wh][rt] = ++idx;
    if(l == r){
        if(wh == 0) addedge(id[wh][rt], l, 0);
        else addedge(l, id[wh][rt], 0);
        return;
    }
    int m = (l + r) / 2;
    build(l, m, rt*2, wh);
    build(m+1, r, rt*2+1, wh);
    if(wh == 0){
        addedge(id[wh][rt], id[wh][rt*2], 0);
        addedge(id[wh][rt], id[wh][rt*2+1], 0);
    }
    else{
        addedge(id[wh][rt*2], id[wh][rt], 0);
        addedge(id[wh][rt*2+1], id[wh][rt], 0);
    }
}

vector <int> vs;
void get(int L, int R, int l, int r, int rt, int wh)
{
    if(L <= l && r <= R){
        vs.push_back(id[wh][rt]);
        return ;
    }
    int m = (l + r) / 2;
    if(L <= m) get(L, R, l, m, rt*2, wh);
    if(m < R) get(L, R, m+1, r, rt*2+1, wh);
}

typedef long long LL;
const LL inf = 0x3f3f3f3f3f3f3f3fLL;
LL d[maxv]; bool vis[maxv];
void dij(int s)
{
    for(int i = 1; i <= 5*n; i++) vis[i] = 0, d[i] = inf;
    d[s] = 0;
    priority_queue <pair<LL, int>> q;
    q.push({-0, s});
    while(q.size()){
        int u = q.top().second; q.pop();
        if(vis[u]) continue;
        vis[u] = 1;
        for(auto &e : G[u]){
            int v, c;
            tie(v, c) = e;
            if(d[v] > d[u] + c){
                d[v] = d[u] + c;
                q.push({-d[v], v});
            }
        }
    }
}

int main()
{
    cin >> n >> m >> s;
    for(int i = 1; i <= 5*n; i++) G[i].clear();
    idx = n;
    build(1, n, 1, 0);
    build(1, n, 1, 1);
    while(m--){
        int t, u;
        scanf("%d%d", &t, &u);
        if(t == 1){
            int v, c;
            scanf("%d%d", &v, &c);
            addedge(u, v, c);
        }
        else if(t == 2){
            vs.clear();
            int l, r, c;
            scanf("%d%d%d", &l, &r, &c);
            get(l, r, 1, n, 1, 0);
            for(int v : vs) addedge(u, v, c);
        }
        else{
            vs.clear();
            int l, r, c;
            scanf("%d%d%d", &l, &r, &c);
            get(l, r, 1, n, 1, 1);
            for(int v : vs) addedge(v, u, c);
        }
    }
    dij(s);
    for(int i = 1; i <= n; i++){
        if(d[i] == inf) d[i] = -1;
        printf("%I64d ", d[i]);
    }
    printf("\n");
}