BZOJ3832 : [Poi2014]Rally

时间:2023-03-08 16:38:31

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;
}