Sol
线段树+拓扑序.
先把图的拓扑序搞出来,然后统计从起点到该点最长链,从该点到终点的最长链,然后建个起点终点,这里跟网络流很像,把它统一到一个有起点的图中,这里也要注意下细节处理.S,T的一个边割掉后最长链就是答案.
然后一开始所有点都在T的集合中,一个个将点加入S集合,用线段树维护每个节点 (从起点到该点最长链+从终点到该点的最长链)的长度,其实就是一个权值线段树,然后就是加加减减的...
Code
/************************************************************** Problem: 3832 User: BeiYu Language: C++ Result: Accepted Time:13472 ms Memory:65804 kb ****************************************************************/ #include<cstdio> #include<vector> #include<queue> #include<iostream> using namespace std; const int N = 500050; #define debug(a) cout<<#a<<"="<<a<<endl inline int in(int x=0,char ch=getchar()){ while(ch>‘9‘||ch<‘0‘) ch=getchar(); while(ch>=‘0‘&&ch<=‘9‘) x=(x<<3)+(x<<1)+ch-‘0‘,ch=getchar();return x; } int n,m,l,r,s,t;int du[N],q[N],f[N][2]; vector<int> g[N],h[N]; int ans1,ans2; #define mid ((l+r)>>1) #define lc (o<<1) #define rc (o<<1|1) int d[N<<2],ma[N<<2]; void Add(int x,int v,int o=1,int l=0,int r=n+2){ if(l==r){ d[o]+=v;return; } if(x<=mid) Add(x,v,lc,l,mid); else Add(x,v,rc,mid+1,r); d[o]=d[lc]+d[rc]; } int Query(int o=1,int l=0,int r=n+2){ if(l==r) return l; if(d[rc]>0) return Query(rc,mid+1,r); else return Query(lc,l,mid); } #undef mid #undef lc #undef rc int main(){ // freopen("in.in","r",stdin); n=in(),m=in(); for(int i=1,u,v;i<=m;i++){ u=in(),v=in(),du[v]++; g[u].push_back(v),h[v].push_back(u); } for(int i=1;i<=n;i++) if(!du[i]) q[r++]=i; int cnt=r; for(int x;l<r;){ x=q[l++]; for(int i=0,lim=g[x].size();i<lim;i++) if(!--du[g[x][i]]) q[r++]=g[x][i]; }for(int i=0;i<n;i++){ for(int j=0,u=q[i],lim=h[u].size();j<lim;j++) f[u][0]=max(f[u][0],f[h[u][j]][0]+1); for(int j=0,u=q[r-i-1],lim=g[u].size();j<lim;j++) f[u][1]=max(f[u][1],f[g[u][j]][1]+1); } // for(int i=0;i<cnt;i++) Add(f[q[i]][1],1); // debug(cnt); // debug(Query()); s=n+1,t=s+1; for(int i=1;i<=n;i++){ g[s].push_back(i),h[i].push_back(s); g[i].push_back(t),h[t].push_back(i); Add(f[i][1],1); } f[s][0]=-1,f[t][1]=-1; ans2=n+2; // for(int i=1;i<=t;i++) printf("%d %d\n",f[i][0],f[i][1]); for(int i=0;i<n;i++){ int u=q[i]; for(int j=0,lim=h[u].size();j<lim;j++) Add(f[h[u][j]][0]+f[u][1]+1,-1); int tmp=Query(); if(tmp<ans2) ans2=tmp,ans1=u; for(int j=0,lim=g[u].size();j<lim;j++) Add(f[g[u][j]][1]+f[u][0]+1,1); }return printf("%d %d\n",ans1,ans2); }
时间: 2024-10-05 15:28:33