Problem
屠龙宝刀注册即送,只需三fan钟,你就会和我一样爱上介款游戏
题意:给定一个DAG,求删掉哪个点后图上最长链最短
Solution
这道题的脑洞还是有点大的,Orz_Poles
想到不能枚举删除每个点后跑拓扑
想到可以跑拓扑时枚举删除每个点
发现对于每个点,删掉后都等价于把连接这个点的所有边删掉,于是乎枚举删除一个点等价于枚举删除一个点所连的边
那这时如何快速求得删掉这些边后的全局最长链就成了问题的关键,这时可以考虑以边的角度贡献答案
通过一条边
(u,v)
(
u
,
v
)
的最长链一定是
“以
u
u
结束的最长链长度”+1+“以开始的最长链长度”
除此之外,还要考虑到删掉这个点所连的边后会对一些边的贡献造成影响(造成影响的情况是一条最长链同时经过这个点连的一条边和上面的边,则这条边的贡献作废),所以取贡献时只取与这个点不互达的边中贡献最大值
整个流程可以用拓扑排序解决(因为拓扑时会在DAG上形成一条分割线,分割线上的边之间不可互达),用可并堆或线段树维护
考虑这种做法的正确性:对于拓扑排序的队列,一定会把这张图分割,然而一张图中的最长链一定会经过这条分割线,所以只用统计拓扑队列所构成的分割线上的边的贡献即可
整理下思路:
先拓扑求出分别以每个点作为起点和终点的最长链长度
↓
拓扑的同时对节点 x x <script type="math/tex" id="MathJax-Element-1244">x</script>进行统计:
↓
先将所有连向这个点的边删除,删除这些边的贡献
↓
统计答案,取出当前所有贡献中的最大值
↓
激活这个点连出去的边对答案的贡献,添加这些边的贡献
↓
继续拓扑
维护贡献的情况,可以用可并堆,维护边的线段树和权值线段树,蒟蒻比较懒,代码中用的是权值线段树
Code
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
#define rg register
#define cl(x) memset(x,0,sizeof(x))
#define max(x,y) ((x)>(y)?(x):(y))
#define min(x,y) ((x)<(y)?(x):(y))
#define abs(x) ((x)>0?(x):(-(x)))
#define inf (0x7fffffff)
template <typename _Tp> inline _Tp read(_Tp&x){
rg char c11=getchar(),ob=0;x=0;
while(c11^'-'&&!isdigit(c11))c11=getchar();if(c11=='-')c11=getchar(),ob=1;
while(isdigit(c11))x=x*10+c11-'0',c11=getchar();if(ob)x=-x;return x;
}
const int N=501000,M=1001000;
struct Edge{int v,nxt;}a[M],A[M];
int in[N],out[N],head[N],Head[N],top[N],f[N],g[N],mx[N<<2],cnt[N<<2];
int n,m,_,__,top_;
inline void add(int u,int v){a[++_].v=v,a[_].nxt=head[u],head[u]=_;}
inline void Add(int u,int v){A[++__].v=v,A[__].nxt=Head[u],Head[u]=__;}
void topology();
#define mid (((l)+(r))>>1)
inline void update(int l,int r,int x,int pos,int alpha){
if(l==r){
cnt[x]+=alpha;
if(cnt[x]>0)mx[x]=l;
else mx[x]=-1,cnt[x]=0;
return ;
}
if(pos<=mid)update(l,mid,x<<1,pos,alpha);
else update(mid+1,r,x<<1|1,pos,alpha);
mx[x]=max(mx[x<<1],mx[x<<1|1]);
return ;
}
#undef mid
void work(){
for(rg int i=1;i<=n;++i)update(0,n,1,g[i],1);
rg int ans=inf,pos;
for(rg int i=1;i<=n;++i){
rg int x=top[i];
for(rg int j=Head[x];j;j=A[j].nxt)
update(0,n,1,g[x]+f[A[j].v]+1,-1);
update(0,n,1,g[x],-1);
if(mx[1]<ans)ans=mx[1],pos=x;
for(rg int j=head[x];j;j=a[j].nxt)
update(0,n,1,f[x]+g[a[j].v]+1,1);
update(0,n,1,f[x],1);
}
printf("%d %d\n",pos,ans);
return ;
}
int main(){
read(n),read(m);
for(rg int i=1,x,y;i<=m;++i)
read(x),read(y),
add(x,y),++in[y],
Add(y,x),++out[x];
topology();
work();
return 0;
}
void topology(){
queue <int> q;while(!q.empty())q.pop();
for(rg int i=1;i<=n;++i)if(!in[i])q.push(i),top[++top_]=i;
while(!q.empty()){
rg int x=q.front();q.pop();
for(rg int i=head[x];i;i=a[i].nxt){
f[a[i].v]=max(f[a[i].v],f[x]+1);
if(!(--in[a[i].v]))q.push(a[i].v),top[++top_]=a[i].v;
}
}
while(!q.empty())q.pop();
for(rg int i=1;i<=n;++i)if(!out[i])q.push(i);
while(!q.empty()){
rg int x=q.front();q.pop();
for(rg int i=Head[x];i;i=A[i].nxt){
g[A[i].v]=max(g[A[i].v],g[x]+1);
if(!(--out[A[i].v]))q.push(A[i].v);
}
}
return ;
}