Codeforces787D - Legacy

Posted

V

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了Codeforces787D - Legacy相关的知识,希望对你有一定的参考价值。

Description

\\(n(n\\leq10^5)\\)个点构成的有向图,有\\(m(m\\leq10^5)\\)条连通信息,信息有三种:

  • 1 u v w,表示存在一条边权为\\(w\\)的有向边\\((u,v)\\)
  • 2 u L R w,表示\\(\\forall v\\in[L,R]\\),存在一条边权为\\(w\\)的有向边\\((u,v)\\)
  • 3 u L R w,表示\\(\\forall v\\in[L,R]\\),存在一条边权为\\(w\\)的有向边\\((v,u)\\)

其中\\(w\\leq10^9\\)。求点\\(s\\)到每个点的最短路,不存在输出\\(-1\\)

Solution

线段树优化建图。
建立两棵线段树,其上点的点权分别表示“到达这个区间内所有点的最小花费”和“到达这个区间内任意一个点的最小花费”。

第一棵线段树上,由于花费\\(v_{[L,R]}\\)能够到达\\([L,R]\\)中所有点,当然也包含\\([L,mid]\\)\\([mid+1,R]\\),所以父节点向子节点连0边;第二棵线段树上,由于花费\\(v_{[L,R]}\\)能够到达\\([L,R]\\)中的一个点,这个点当然也包含在其父节点中,所以子节点向父节点连0边。
如果不做感性理解的话,两棵线段树上的点分别用于连和被连,连向第一棵树上的\\([L,R]\\)就等价于连向\\([L,R]\\)中的每一个点,被第二棵树上的\\([L,R]\\)连就等价于被\\([L,R]\\)中的每一个点连。
由于每一条信息最多建立\\(O(logn)\\)条边,所以总边数是\\(O(mlogn+4n)\\)
建完图后直接跑一遍单源最短路就好啦。

Code

//Legacy
#include <cstdio>
#include <cstring>
#include <queue>
typedef long long lint;
inline char gc()
{
    static char now[1<<16],*s,*t;
    if(s==t) {t=(s=now)+fread(now,1,1<<16,stdin); if(s==t) return EOF;}
    return *s++;
}
inline int read()
{
    int x=0; char ch=gc();
    while(ch<\'0\'||\'9\'<ch) ch=gc();
    while(\'0\'<=ch&&ch<=\'9\') x=x*10+ch-\'0\',ch=gc();
    return x;
}
inline int min(int x,int y) {return x<y?x:y;}
const int N=1e5+10;
int n,m,s;
const int N1=3e5+110;
int cnt,rt1,rt2,ch[N1][2];
int h[N1],edCnt;
struct edge{int v,w,nxt;} ed[N*20];
inline void edAdd(int u,int v,int w)
{
    edCnt++; ed[edCnt].v=v,ed[edCnt].w=w;
    ed[edCnt].nxt=h[u],h[u]=edCnt;
}
void bldTr1(int &p,int L0,int R0)
{
    if(L0==R0) {p=L0; return;}
    p=++cnt;
    int mid=L0+R0>>1;
    bldTr1(ch[p][0],L0,mid);
    bldTr1(ch[p][1],mid+1,R0);
    edAdd(p,ch[p][0],0),edAdd(p,ch[p][1],0);
}
void bldTr2(int &p,int L0,int R0)
{
    if(L0==R0) {p=L0; return;}
    p=++cnt;
    int mid=L0+R0>>1;
    bldTr2(ch[p][0],L0,mid);
    bldTr2(ch[p][1],mid+1,R0);
    edAdd(ch[p][0],p,0),edAdd(ch[p][1],p,0);
}
int optL,optR;
void add(int p,int L0,int R0,int u,int w,int type)
{
    if(optL<=L0&&R0<=optR)
    {
        if(type==2) edAdd(u,p,w); else edAdd(p,u,w);
        return;
    }
    int mid=L0+R0>>1;
    if(optL<=mid) add(ch[p][0],L0,mid,u,w,type);
    if(mid<optR) add(ch[p][1],mid+1,R0,u,w,type);
}
const lint INF=0x3F3F3F3F3F3F3F3F;
lint dst[N1];
std::queue<int> Q;
void SPFA(int s)
{
    memset(dst,0x3F,sizeof dst);
    dst[s]=0; Q.push(s);
    while(!Q.empty())
    {
        int u=Q.front(); Q.pop();
        for(int i=h[u];i;i=ed[i].nxt)
        {
            int v=ed[i].v,w=ed[i].w;
            if(dst[u]+w<dst[v]) dst[v]=dst[u]+w,Q.push(v);
        }
    }
}
int main()
{
    n=read(),m=read(),s=read();
    cnt=n;
    bldTr1(rt1,1,n); bldTr2(rt2,1,n);
    while(m--)
    {
        int opt=read(),u,v,w;
        if(opt==1)
        {
            u=read(),v=read(),w=read();
            edAdd(u,v,w); continue;
        }
        u=read(); optL=read(),optR=read(); w=read();
        add(opt==2?rt1:rt2,1,n,u,w,opt);
    }
    SPFA(s);
    for(int i=1;i<=n;i++) printf("%lld ",dst[i]<INF?dst[i]:-1);
    puts("");
    return 0;
}

以上是关于Codeforces787D - Legacy的主要内容,如果未能解决你的问题,请参考以下文章

[CF787D] legacy

CF787D Legacy

[CF787D]遗产(Legacy)-线段树-优化Dijkstra(内含数据生成器)

Legacy CodeForces

Codeforces 787D. Legacy 线段树优化建图+最短路

codeforces 786B Legacy