【SHOI2008】【BZOJ1023】cactus仙人掌图

Description

如果某个无向连通图的任意一条边至多只出现在一条简单回路（simple cycle）里，我们就称这张图为仙人图（cactus）。所谓简单回路就是指在图上不重复经过任何一个顶点的回路。

举例来说，上面的第一个例子是一张仙人图，而第二个不是——注意到它有三条简单回路：（4，3，2，1，6，5，4）、（7，8，9，10，2，3，7）以及（4，3，7，8，9，10，2，1，6，5，4），而（2，3）同时出现在前两个的简单回路里。另外，第三张图也不是仙人图，因为它并不是连通图。显然，仙人图上的每条边，或者是这张仙人图的桥（bridge），或者在且仅在一个简单回路里，两者必居其一。定义在图上两点之间的距离为这两点之间最短路径的距离。定义一个图的直径为这张图相距最远的两个点的距离。现在我们假定仙人图的每条边的权值都是1，你的任务是求出给定的仙人图的直径。

Input

输入的第一行包括两个整数n和m（1≤n≤50000以及0≤m≤10000）。其中n代表顶点个数，我们约定图中的顶点将从1到n编号。接下来一共有m行。代表m条路径。每行的开始有一个整数k（2≤k≤1000），代表在这条路径上的顶点个数。接下来是k个1到n之间的整数，分别对应了一个顶点，相邻的顶点表示存在一条连接这两个顶点的边。一条路径上可能通过一个顶点好几次，比如对于第一个样例，第一条路径从3经过8，又从8返回到了3，但是我们保证所有的边都会出现在某条路径上，而且不会重复出现在两条路径上，或者在一条路径上出现两次。

Output

只需输出一个数，这个数表示仙人图的直径长度。

Sample Input

15 3

9 1 2 3 4 5 6 7 8 3

7 2 9 10 11 12 13 10

5 2 14 9 15 10 8

10 1

10 1 2 3 4 5 6 7 8 9 10
Sample Output

9
HINT

对第一个样例的说明：如图，6号点和12号点的最短路径长度为8，所以这张图的直径为8。

【注意】使用Pascal语言的选手请注意：你的程序在处理大数据的时候可能会出现栈溢出。如果需要调整栈空间的大小，可以在程序的开头填加一句：{$M 5000000}，其中5000000即指代栈空间的大小，请根据自己的程序选择适当的数值。

做法我在那个仙人掌DP的学习笔记里写了
记得要拆环成两倍长度序列所以数组要开大一倍

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define MAXN 50010
using namespace std;
int n,m,top;
int ans;
int k,u,v;
int Top;
int deep[MAXN<<1],dfn[MAXN<<1],fa[MAXN<<1],minn[MAXN<<1];
int f[MAXN<<1],a[MAXN<<1];
int que[MAXN<<1],head,tail;//维护一个单调队列 
struct edge
{
    int to;
    edge *next;
}e[MAXN<<2],*prev[MAXN];
inline void insert(int u,int v)
{
    e[++top].to=v;e[top].next=prev[u];prev[u]=&e[top];
}
inline void in(int &x)
{
    char ch=getchar();x=0;int flag=1;
    while (!(ch>='0'&&ch<='9')) flag=ch=='-'?-1:1,ch=getchar();
    while (ch>='0'&&ch<='9')    x=x*10+ch-'0',ch=getchar();x*=flag;
}
inline void dp(int root,int x)
{
    int tot=deep[x]-deep[root]+1,temp=tot;
    for (int i=x;i!=root;i=fa[i])   a[temp--]=f[i];
    a[temp]=f[root];
    for (int i=1;i<=tot;i++)    a[i+tot]=a[i];//把环拆开 
    que[1]=1;head=tail=1;//head,tail分别对应区间端点 
    for (int i=2;i<=2*tot;i++)
    {
        while (head<=tail&&i-que[head]>tot/2)   head++;
        ans=max(ans,a[i]+i+a[que[head]]-que[head]);
        while (head<=tail&&a[que[tail]]-que[tail]<=a[i]-i)  tail--;
        que[++tail]=i;
    }
    for (int i=2;i<=tot;i++)
        f[root]=max(f[root],a[i]+min(i-1,tot-i+1));
}
inline void dfs(int x)
{
    dfn[x]=minn[x]=++Top;
    for (edge *i=prev[x];i;i=i->next)
    {
        if (i->to==fa[x])   continue;
        if (!dfn[i->to])//没有访问过 
        {
            fa[i->to]=x;deep[i->to]=deep[x]+1;
            dfs(i->to);
            minn[x]=min(minn[x],minn[i->to]);
        }   
        else    minn[x]=min(minn[x],dfn[i->to]);
        if (dfn[x]<minn[i->to])
        {
            ans=max(ans,f[x]+f[i->to]+1);
            f[x]=max(f[x],f[i->to]+1);
        }
    }
    for (edge *i=prev[x];i;i=i->next)
        if (fa[i->to]!=x&&dfn[x]<dfn[i->to])    dp(x,i->to);
}
int main()
{
    in(n);in(m);
    for (int i=1;i<=m;i++)
    {
        in(k);in(u);
        for (int j=2;j<=k;j++)  in(v),insert(u,v),insert(v,u),u=v;
    }
    dfs(1);
    printf("%d\n",ans);
}