BZOJ1395 : [Baltic2005]Trip

Claris發表於2017-01-30

建立新圖,原圖中每條邊在新圖中是點,新圖中每個點的點權為$-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;
}