建立新圖,原圖中每條邊在新圖中是點,新圖中每個點的點權為$-e[i].c+e[i].b$,邊權為$0$。
若$e[i].d\leq e[j].a$,則連一條$i$到$j$的單向邊。
對於原圖中每個點,將所有入邊和出邊分別排序,然後建立一排虛點表示字尾,通過雙指標將邊數優化至$O(m)$。
在新圖中求出最短路,最後將答案加上$T$即可。
注意到新圖是個DAG,因此可以記搜求解。
時間複雜度$O(m\log m)$。
#include<cstdio> #include<algorithm> const int N=50010,M=100010; int n,m,cnt,P,T,i,j,k,gi[N],go[N],g[M*2],v[M*5],nxt[M*5],ed,a[M],b[M],ca,cb,vis[M*2],f[M*2],ans; struct E{int x,y,a,b,c,d;}e[M]; inline bool cmpa(int x,int y){return e[x].d<e[y].d;} inline bool cmpb(int x,int y){return e[x].a<e[y].a;} inline void read(int&a){char c;while(!(((c=getchar())>='0')&&(c<='9')));a=c-'0';while(((c=getchar())>='0')&&(c<='9'))(a*=10)+=c-'0';} inline void add(int*g,int x,int y){v[++ed]=y;nxt[ed]=g[x];g[x]=ed;} inline void up(int&x,int y){if(x>y)x=y;} int dp(int x){ if(vis[x])return f[x]; vis[x]=1; int&t=f[x]; t=1; if(x<=m)if(e[x].x==1)t=0; for(int i=g[x];i;i=nxt[i])up(t,dp(v[i])); if(t<1&&x<=m)t+=e[x].b-e[x].c; return t; } int main(){ read(n),read(m),read(P),read(T); for(i=1;i<=m;i++){ read(e[i].x),read(e[i].y),read(e[i].a),read(e[i].b),read(e[i].c),read(e[i].d); add(go,e[i].x,i); add(gi,e[i].y,i); } cnt=m; for(i=1;i<=n;i++){ for(ca=0,j=gi[i];j;j=nxt[j])a[++ca]=v[j]; if(!ca)continue; for(cb=0,j=go[i];j;j=nxt[j])b[++cb]=v[j]; if(!cb)continue; std::sort(a+1,a+ca+1,cmpa); std::sort(b+1,b+cb+1,cmpb); for(j=1;j<=cb;j++){ if(j<cb)add(g,cnt+j+1,cnt+j); add(g,b[j],cnt+j); } for(j=ca,k=cb+1;j;j--){ while(k>1&&e[a[j]].d<=e[b[k-1]].a)k--; if(k<=cb)add(g,cnt+k,a[j]); } cnt+=cb; } if(P!=1)ans=1; for(i=1;i<=m;i++)if(e[i].y==P&&e[i].d<=T){ dp(i); if(f[i]<1)up(ans,f[i]); } if(ans>0)ans=-1;else ans+=T; return printf("%d",ans),0; }