仙人掌直径--poj3567/bzoj1023

差不多是个模板题
给一张仙人掌图求直径
要注意环上的点到最高点的距离要去两侧min值

于是和基环树一样
环上单独求，顺便更新ans和最高点max值
可以用dp+单调队列优化复杂度
高端操作get

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#define maxn 50005
using namespace std;
int n,m,cnt=1,head[maxn],dfn[maxn],low[maxn],q[maxn];
int fa[maxn],dis[maxn],num,a[maxn],b[maxn],f[maxn],ans;

inline int rd(){
  int x=0,f=1; char c=' ';
  while(c<'0' || c>'9') {if(c=='-') f=-1;c=getchar();}
  while(c<='9' && c>='0') x=x*10+c-'0',c=getchar();
  return x*f;
}

struct EDGE{
  int to,nxt;
}edge[maxn*10];

void add(int x,int y){
  edge[++cnt].to=y;
  edge[cnt].nxt=head[x];
  head[x]=cnt;
}

void dp(int x,int y){
  for(m=0;x!=y;y=fa[y]) a[++m]=dis[y];
  a[++m]=dis[x];
  for(int i=1;i<m;i++) a[i+m]=a[i];//可以把环拆成链
  int p=m/2;
  int l=1,r=1; q[1]=1;//这个是dp+单调队列优化 不过听说暴力能过？
  for(int i=2;i<=m+p;i++){
    while(l<=r && i-q[l]>p) l++;//找到小于半环的两个点
    ans=max(ans,a[q[l]]+a[i]+i-q[l]);
    while(l<=r && a[q[r]]+i-q[r]<=a[i]) r--;//单调队列优化关键所在
    q[++r]=i;
  }
  for(int i=1;i<m;i++) dis[x]=max(dis[x],a[i]+min(i,m-i));//把最高点的dis值也更新了
}

void tarjan(int u){
  dfn[u]=low[u]=++num;
  for(int i=head[u];i;i=edge[i].nxt){
    int v=edge[i].to;
    if(fa[u]==v) continue;
    if(!dfn[v]){
      fa[v]=u;
      tarjan(v);
      low[u]=min(low[u],low[v]);
      if(dfn[u]<low[v]){//这个地方是判环的常规操作，u，v不在一个环上
        ans=max(ans,dis[u]+dis[v]+1);
        dis[u]=max(dis[u],dis[v]+1);
      }
    }
    else low[u]=min(dfn[v],low[u]);//之前这个else不小心放到了if里面调了好久···
  }
  int v;
  for(int i=head[u];i;i=edge[i].nxt)
    if(fa[v=edge[i].to]!=u && dfn[u]<dfn[v]) dp(u,v);//u，v在一个环上，常规操作
}

int main(){
  n=rd(); m=rd();
  for(int i=1;i<=m;i++){
    int k=rd(),x=rd();
    for(int j=2;j<=k;j++){
      int y=rd();
      add(x,y); add(y,x);
      x=y;
    }
  }
  tarjan(1);
  printf("%d\n",ans);
  return 0;
}