bzoj3073 [Pa2011]Journeys（线段树优化建图+最短路）

一个区间的点向另一个区间的点连双向边。求单源最短路。
我们考虑线段树优化建图。
因为有两个方向，因此我们需要建两颗线段树。
A线段树儿子向父亲连边，表示一条有向边的起点
B线段树父亲向儿子连边，表示一条有向边的终点
B线段树的每个节点向A线段树的对应节点连边（一条有向边的终点，同时是下一条有向边的起点，其实只连叶子就行）
这些边边权均为0.

然后对于每一天，新建一个节点表示(a,b)->(c,d)，在线段树上连一连。新建一个节点表示(c,d)->(a,b)，在线段树上连一连。

A线段树中的叶子节点就可以对应为原图中的点。跑Dij即可。
点数$O(4n+2m)$,边数$O(5n+4mlogn)$

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
#include <utility>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 500010
#define M 2210000
#define pa pair<int,int>
inline char gc(){
    static char buf[1<<16],*S,*T;
    if(S==T){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=gc();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x*f;
}
int n,m,S,tot=0,h[M],num=0,dis[M],rt1=0,rt2=0,pos[N];
struct edge{
    int to,next,val;
}data[11000000];
struct node{
    int lc,rc;
}tr[N<<2];
bool vis[M];
inline void add(int x,int y,int val){
    data[++num].to=y;data[num].next=h[x];h[x]=num;data[num].val=val;
}
inline void build1(int &p,int l,int r){
    p=++tot;if(l==r){pos[l]=p;return;}int mid=l+r>>1;
    build1(tr[p].lc,l,mid);build1(tr[p].rc,mid+1,r);
    add(tr[p].lc,p,0);add(tr[p].rc,p,0);
}
inline void build2(int &p,int l,int r){
    p=++tot;if(l==r){add(p,pos[l],0);return;}int mid=l+r>>1;
    build2(tr[p].lc,l,mid);build2(tr[p].rc,mid+1,r);
    add(p,tr[p].lc,0);add(p,tr[p].rc,0);
}
inline void cover1(int p,int l,int r,int x,int y){
    if(x<=l&&r<=y){add(p,tot,1);return;}int mid=l+r>>1;
    if(x<=mid) cover1(tr[p].lc,l,mid,x,y);
    if(y>mid) cover1(tr[p].rc,mid+1,r,x,y);
}
inline void cover2(int p,int l,int r,int x,int y){
    if(x<=l&&r<=y){add(tot,p,1);return;}int mid=l+r>>1;
    if(x<=mid) cover2(tr[p].lc,l,mid,x,y);
    if(y>mid) cover2(tr[p].rc,mid+1,r,x,y);
}
inline void Dijkstra(){
    priority_queue<pa,vector<pa>,greater<pa> >q;
    memset(dis,inf,sizeof(dis));
    q.push(make_pair(0,pos[S]));dis[pos[S]]=0;
    while(!q.empty()){
        int x=q.top().second;q.pop();
        if(vis[x]) continue;vis[x]=1;
        for(int i=h[x];i;i=data[i].next){
            int y=data[i].to;
            if(dis[x]+data[i].val<dis[y]){
                dis[y]=dis[x]+data[i].val;
                q.push(make_pair(dis[y],y));
            }
        }
    }
}
int main(){
//  freopen("a.in","r",stdin);
    n=read();m=read();S=read();
    build1(rt1,1,n);build2(rt2,1,n);
    while(m--){
        int x1=read(),y1=read(),x2=read(),y2=read();
        ++tot;cover1(rt1,1,n,x1,y1);cover2(rt2,1,n,x2,y2);
        ++tot;cover1(rt1,1,n,x2,y2);cover2(rt2,1,n,x1,y1);
    }Dijkstra();
    for(int i=1;i<=n;++i) printf("%d\n",dis[pos[i]]>>1);
    return 0;
}