Codeforces 786 B. Legacy

题目链接：http://codeforces.com/contest/786/problem/B

---

　　典型线段树优化连边，线段树上的每一个点表示这个区间的所有点，然后边数就被优化为了至多${nlogn}$条，注意要区分$2$,$3$操作(建两棵线段树)，然后最短路即可。

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<vector>
#include<cstdlib>
#include<cmath>
#include<cstring>
using namespace std;
#define maxn 1001000
#define llg long long
#define INC 300000
#define SIZE 2000000 
#define inf (llg)1e18
#define yyj(a) freopen(a".in","r",stdin),freopen(a".out","w",stdout);
llg n,dis[maxn],dl[SIZE+10],head,tail,pos[maxn];
bool bj[maxn];

vector<llg>a[maxn],val[maxn];

inline int getint()
{
       int w=0,q=0; char c=getchar();
       while((c<'0' || c>'9') && c!='-') c=getchar(); if(c=='-') q=1,c=getchar(); 
       while (c>='0' && c<='9') w=w*10+c-'0', c=getchar(); return q ? -w : w;
}

void link(llg x,llg y,llg z) {a[x].push_back(y),val[x].push_back(z);}

void build(llg o,llg l,llg r)
{
    if (l==r)
    {
        pos[l]=o;
        return ;
    }
    llg mid=(l+r)>>1,lc=o<<1,rc=o<<1|1;
    build(lc,l,mid); build(rc,mid+1,r);
}

void find(llg o,llg l,llg r,llg L,llg R,llg val,bool f,llg v)
{
    if (l>=L && r<=R)
    {
        if (f) link(pos[v],o,val);else link(o+INC,pos[v],val);
        return ;
    }
    llg mid=(l+r)>>1,lc=o<<1,rc=o<<1|1;
    if (mid>=L) find(lc,l,mid,L,R,val,f,v);
    if (mid<R) find(rc,mid+1,r,L,R,val,f,v);
}



void link_fa(llg o,llg l,llg r)
{
    for (llg i=l;i<=r;i++) link(o,pos[i],0),link(pos[i],o+INC,0);
    if (l==r) return ;
    llg mid=(l+r)>>1,lc=o<<1,rc=o<<1|1;
    link_fa(lc,l,mid); link_fa(rc,mid+1,r);
}

void init()
{
    llg t,v,u,w,l,r,s,q;
    n=getint(),q=getint(),s=getint();
    build(1,1,n);
    for (llg i=1;i<=q;i++)
    {
        t=getint();
        if (t==1) 
        {
            v=getint(),u=getint(),w=getint();
            link(pos[v],pos[u],w);
        }
        if (t==2)
        {
            v=getint(),l=getint(),r=getint(),w=getint();
            find(1,1,n,l,r,w,1,v);
        }
        if (t==3)
        {
            v=getint(),l=getint(),r=getint(),w=getint();
            find(1,1,n,l,r,w,0,v);
        }
    }
    link_fa(1,1,n);
    for (llg i=1;i<maxn;i++) dis[i]=inf;
    dis[pos[s]]=0;
    dl[1]=pos[s],head=0,tail=1;
}

void spfa()
{
    llg w,x,v;
    do
    {
        head%=SIZE; head++;
        x=dl[head]; w=a[x].size(); bj[x]=0;
        for (llg i=0;i<w;i++)
        {
            v=a[x][i];
            if (dis[v]>dis[x]+val[x][i])
            {
                dis[v]=dis[x]+val[x][i];
                if (!bj[v])
                {
                    bj[v]=1;
                    tail%=SIZE; tail++;
                    dl[tail]=v;
                }
            }
        }
    }while (head!=tail);
}

int main()
{
    yyj("graph");
    init();
    spfa();
    for (llg i=1;i<=n;i++) if (dis[pos[i]]==inf) printf("-1 "); else printf("%lld ",dis[pos[i]]);
    return 0;
}