bzoj 1023: [SHOI2008]cactus仙人掌图

起了个仙人掌的名字但是它并不是仙人掌（的做法）。。。。

首先这个图确实是个仙人掌。。。然后如果只有树的话就可以只跑树上最长链（dfs两遍）。

然而有一些比较烦人的环。所以我们考虑把它优化掉。所以把树上最长链以dp的形式实现。

dp[i]代表这个点上的最长链。

然后环怎么办呢，可以首先造一个dfs树，然后把环的最高点作为这个环用来dp的点，

其余的点的互相连通情况在判环的过程中进行，之后这个点的dp值就是其中一个最长链的值。

而环怎么判定呢。。。用tarjan缩点的操作。。。。（其中的dp用单调队列维护）。

#include<cstdio>
#include<cmath>
#include<algorithm>
#include<cstring>
using namespace std;
int n,m;
struct node
{
	int to;
	int nxt;
}edge[300005];
int head[100005];
int dfn[100005],low[100005],fa[100005],dp[100005],que[100005],sta[100005];
int deep;
int cnt=1,ans;
void init()
{
	memset(head,-1,sizeof(head));
	deep=0;cnt=1;ans=0;
	memset(dfn,0,sizeof(dfn));
	memset(low,0,sizeof(low));
	memset(fa,0,sizeof(fa));
	memset(dp,0,sizeof(dp));
	memset(que,0,sizeof(que));
	memset(sta,0,sizeof(sta));
}
void add(int from,int to)
{
	edge[cnt].to=to;
	edge[cnt].nxt=head[from];
	head[from]=cnt++;
}
void dpit(int rt,int to)
{
	int cct=0;
	while(rt!=to)
	{
		sta[++cct]=dp[to];
		to=fa[to];
	}
	sta[++cct]=dp[rt];
	for(int i=1;i<cct;i++)
	{
		sta[i+cct]=sta[i];
	}
	int head=1,tail=1;
	que[1]=1;
	for(int i=2;i<=cct+cct/2;i++)//一圈半用来更新全每个点的状态。单调队列维护，查出最长的路径 
	{
		while(head<=tail&&i-que[head]>cct/2)head++;
		ans=max(ans,sta[i]+sta[que[head]]+i-que[head]);
		while(head<=tail&&sta[que[tail]]+i-que[tail]<=sta[i])tail--;
		que[++tail]=i;
	}
	for(int i=1;i<cct;i++)
	{
		dp[rt]=max(dp[rt],sta[i]+min(i,cct-i));
	}
}
void tarjan(int rt)
{
	dfn[rt]=low[rt]=++deep;
	for(int i=head[rt];i!=-1;i=edge[i].nxt)
	{
		int to=edge[i].to;
		if(to==fa[rt])continue;
		if(!dfn[to])
		{
			fa[to]=rt;
			tarjan(to);
			low[rt]=min(low[rt],low[to]);
			if(dfn[rt]<low[to])//这个点在rt的树边上
			{
				ans=max(ans,dp[rt]+1+dp[to]);
				dp[rt]=max(dp[rt],dp[to]+1);
			} 
		}else
		{
			low[rt]=min(low[rt],dfn[to]);
		}
	}
	for(int i=head[rt];i!=-1;i=edge[i].nxt)
	{
		int to=edge[i].to;
		if(fa[to]==rt||dfn[rt]>=dfn[to])continue;//这个点不是数边上的点也不是这个点的上面（其实就是父节点）。 
		dpit(rt,to);
	}
}
int main()
{
	while(scanf("%d%d",&n,&m)!=EOF)
	{
		init();
		for(int i=1;i<=m;i++)
		{
			int k,las=0;scanf("%d",&k);
			for(int j=1;j<=k;j++)
			{
				int x;scanf("%d",&x);
				if(!las)
				{
					las=x;continue;
				}
				add(las,x);add(x,las);
				las=x;
			}
		}
		tarjan(1);
		printf("%d\n",ans);
	}
	return 0;
}