首页 > 代码库 > BZOJ1395 : [Baltic2005]Trip
BZOJ1395 : [Baltic2005]Trip
建立新图,原图中每条边在新图中是点,新图中每个点的点权为$-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;}
BZOJ1395 : [Baltic2005]Trip
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。