传送门
仙人掌直径,以前好像模拟赛的时候做到过一道基环树的直径,打了个很麻烦的然而还错了……今天才发现那就是这个的弱化版啊……
如果是树的话用普通的dp即可,记\(f[u]\)表示\(u\)往下最长能伸多少。我们用一个类似tarjan的做法,在dfs的过程中记录dfn和low,如果某条边不在环内,那么这条边我们直接按树形dp转移。
否则的话,当我们做完这整个环的时候,考虑如何更新环的顶点\(f[u]\),设某个点为\(v\),那么有\(f[u]=max(f[u],f[v]+min(dep[v]-dep[u],circle\_size-dep[v]+1))\)
然后这个环上的两个顶点可能构成答案。那么维护一个单调队列,破环成链跑一遍就行了
//minamoto
#include
#define fp(i,a,b) for(register int i=a,I=b+1;iI;--i)
#define go(u) for(register int i=head[u],v=e[i].v;i;i=e[i].nx,v=e[i].v)
#define getc() (p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++)
inline int max(const int &x,const int &y){return x>y?x:y;}
inline int min(const int &x,const int &y){return x'9'||ch<'0')(ch=='-')&&(f=-1);
for(res=ch-'0';(ch=getc())>='0'&&ch<='9';res=res*10+ch-'0');
return res*f;
}
const int N=50005;
struct eg{int v,nx;}e[N<<3];int head[N],tot;
inline void add(int u,int v){e[++tot]={v,head[u]},head[u]=tot;}
int n,m,f[N],ans=1,dep[N],fa[N],dfn[N],low[N],tim;
int h,t,q1[N<<1],q[N];
void dp(int u,int y){
int top=0;
for(register int i=y;i!=u;i=fa[i])q1[++top]=i;q1[++top]=u;
reverse(q1+1,q1+1+top);fp(i,1,top)q1[i+top]=q1[i];
h=1,t=0;
fp(i,1,top+top){
while(h<=t&&i-q[h]>top/2)++h;
if(h<=t)ans=max(ans,f[q1[i]]+f[q1[q[h]]]+i-q[h]);
while(h<=t&&f[q1[i]]-i>f[q1[q[t]]]-q[t])--t;
q[++t]=i;
}
for(register int i=y;i!=u;i=fa[i])
f[u]=max(f[u],f[i]+min(dep[i]-dep[u],dep[y]-dep[i]+1));
}
void dfs(int u,int fat){
fa[u]=fat,dfn[u]=low[u]=++tim,dep[u]=dep[fat]+1;
go(u)if(v!=fat){
if(!dfn[v])dfs(v,u),low[u]=min(low[u],low[v]);
else low[u]=min(low[u],dfn[v]);
if(low[v]>dfn[u])ans=max(ans,f[u]+f[v]+1),f[u]=max(f[u],f[v]+1);
}
go(u)if(fa[v]!=u&&dfn[u]