POJ 3567 Cactus Reloaded（仙人掌直径）

题意

裸的仙人掌直径。

题解

先考虑基环树的直径：先算出每颗“树”的直径，再在环上跑DP

再考虑仙人掌的直径：把每个基环树缩成一条边，边长为基环树深度。

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
using namespace std;
const int N=50010;
const int M=10010;
int cnt,head[N];
int c[N*2],dp[N],q[N*2],ans;
int dfn[N],low[N],tot,fa[N];
int n,m;
struct edge{
    int to,nxt;
}e[M*20];
void add(int u,int v){
    cnt++;
    e[cnt].nxt=head[u];
    e[cnt].to=v;
    head[u]=cnt;
}
void getdp(int x,int y){
    int l,r,i,m,p;
    for(m=0;y!=x;y=fa[y])c[++m]=dp[y];
    for(c[++m]=dp[x],i=1;i<m;i++)c[i+m]=c[i];
    l=r=q[1]=1;p=m/2;
    for(i=2;i<=m+p;i++){
        while(l<=r&&i-q[l]>p)l++;
        ans=max(ans,c[i]+c[q[l]]+i-q[l]);
        while(l<=r&&c[i]>=c[q[r]]+i-q[r])r--;
        q[++r]=i;
    }
    for(int i=1;i<m;i++)dp[x]=max(dp[x],c[i]+min(i,m-i));
}
void Tarjan(int u){
    dfn[u]=low[u]=++tot;
    for(int i=head[u];i;i=e[i].nxt){
        int v=e[i].to;
        if(fa[u]==v)continue;
        if(dfn[v]==0){
            fa[v]=u;
            Tarjan(v);
            low[u]=min(low[v],low[u]);
            if(dfn[u]<low[v]){
                ans=max(ans,dp[v]+dp[u]+1);
                dp[u]=max(dp[u],dp[v]+1);
            }
        }else low[u]=min(low[u],dfn[v]);
    }
    for(int i=head[u];i;i=e[i].nxt){
        int v=e[i].to;
        if(fa[v]!=u&&dfn[u]<dfn[v])getdp(u,v);
    }
}
int main(){
    scanf("%d%d",&n,&m);
    for(int i=1,a,x;i<=m;i++){
        scanf("%d%d",&a,&x);
        for(int j=1,y;j<a;j++){
            scanf("%d",&y);
            add(x,y);
            add(y,x);
            x=y;
        }
    }
    Tarjan(1);
    printf("%d",ans);
    return 0;
}