BZOJ3832 : [Poi2014]Rally

f[0][i]为i出发的最长路，f[1][i]为到i的最长路

新建源汇S,T，S向每个点连边，每个点向T连边

将所有点划分为两个集合S与T，一开始S中只有S，其它点都在T中

用一棵线段树维护所有连接属于两个集合的点的边，权值为f[1][u]+f[0][v]

按拓扑序依次计算去掉每个点后图中的最长路

对于当前计算的点x，先将所有连向x的边删除，此时最长路长度为线段树中的最大值

然后再将所有x出发的边加入线段树中

时间复杂度$O(m\log n)$

 

#include<cstdio>
#define N 500010
int n,m,i,j,x,y,g[2][N],nxt[2][N<<2],v[2][N<<2],ed,d[N],q[N],l,r,fin,ans,f[2][N],S,T,val[N<<2];
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 x,int y){
  d[x]++;v[0][++ed]=y;nxt[0][ed]=g[0][x];g[0][x]=ed;
  v[1][ed]=x;nxt[1][ed]=g[1][y];g[1][y]=ed;
}
inline void ins(int c,int d){
  int x=1,a=0,b=n,mid;
  while(a<=b){
    val[x]+=d;
    if(a==b)return;
    mid=(a+b)>>1;x<<=1;
    if(c<=mid)b=mid;else a=mid+1,x|=1;
  }
}
inline int ask(){
  int x=1,a=0,b=n,mid;
  while(a<b){
    mid=(a+b)>>1;x<<=1;
    if(val[x|1])a=mid+1,x|=1;else b=mid;
  }
  return a;
}
int main(){
  read(n),read(m);
  while(m--)read(x),read(y),add(x,y);
  for(i=1;i<=n;i++)if(!d[i])q[++r]=i;
  while(l<=r)for(i=g[1][q[l++]];i;i=nxt[1][i])if(!(--d[v[1][i]]))q[++r]=v[1][i];
  for(i=1;i<=n;i++)for(j=g[0][x=q[i]];j;j=nxt[0][j])if(f[0][y=v[0][j]]>=f[0][x])f[0][x]=f[0][y]+1;
  for(i=n;i;i--)for(j=g[1][x=q[i]];j;j=nxt[1][j])if(f[1][y=v[1][j]]>=f[1][x])f[1][x]=f[1][y]+1;
  for(S=n+1,T=S+1,i=1;i<=n;i++)add(S,i),add(i,T);
  for(i=1;i<=n;i++)ins(f[0][i],1);
  for(fin=i=n;i;i--){
    for(j=g[1][x=q[i]];j;j=nxt[1][j])ins(f[1][v[1][j]]+(v[1][j]<=n)+f[0][x],-1);
    if((y=ask())<fin)ans=x,fin=y;
    for(j=g[0][x];j;j=nxt[0][j])ins(f[0][v[0][j]]+(v[0][j]<=n)+f[1][x],1);
  }
  return printf("%d %d",ans,fin),0;
}