BZOJ 1023 SHOI2008 cactus仙人掌图 仙人掌DP

题目大意：给定一棵仙人掌，求这棵仙人掌的直径

首先Tarjan缩点双，开vector或者链表记录每个点属于哪些点双，以及每个点双中有哪些点

有些点双可能不是环，我们可以补上一条边看成环，无伤大雅

每次DP时，首先枚举环的根节点以外的点，对这些点所在的其它点双DP一遍

然后令f[x]为以x为根的子仙人掌的所有点和x之间的最大距离

然后我们将环倍增 用单调队列来更新答案 保证决策点和被更新点的距离不超过环长度的一半

最后再用环上所有点更新环的根节点的f值即可

注意整棵仙人掌的根节点可能在多个点双中

#include <vector>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define M 50500
#define INF 0x3f3f3f3f
using namespace std;
int n,m,k,cnt,ans=-INF;
namespace Cactus_Graph{
	int size[M],f[M];
	//f[x]表示以x为根的子仙人掌中所有节点到x的最长距离
	vector<int> belong[M];
	vector<int> contain[M];
	int Cactus_DP(int scc)
	{
		vector<int>::iterator it;
		vector<int>::reverse_iterator _it;
		for(_it=contain[scc].rbegin(),_it++;_it!=contain[scc].rend();_it++)
		{
			int x=*_it;
			for(it=belong[x].begin();it!=belong[x].end();it++)
				if( *it!=scc )
					f[x]=max(f[x],Cactus_DP(*it));
		}
		static pair<int,int> q[M<<1];
		int re=-INF,r=0,h=0,limit=size[scc]>>1,pos=0;
		for(_it=contain[scc].rbegin();_it!=contain[scc].rend();_it++)
		{
			int x=*_it;++pos;

			while( r-h>=1 && f[q[r].first]+(pos-q[r].second)<f[x] )
				r--;
			q[++r]=make_pair(x,pos);
			while( r-h>=1 && pos-q[h+1].second>=limit )
				++h;
		}
		for(_it=contain[scc].rbegin();_it!=contain[scc].rend();_it++)
		{
			int x=*_it;++pos;


			if(r-h>=1)
			ans=max(ans,f[q[h+1].first]+(pos-q[h+1].second)+f[x] );

			while( r-h>=1 && f[q[r].first]+(pos-q[r].second)<f[x] )
				r--;
			q[++r]=make_pair(x,pos);
			while( r-h>=1 && pos-q[h+1].second>=limit )
				++h;
		}
		for(_it=contain[scc].rbegin(),pos=0;_it!=contain[scc].rend();_it++,pos++)
			re=max(re,f[*_it]+min(pos,size[scc]-pos) );
		return re;
	}
}
namespace Origin_Graph{
	struct abcd{
		int to,next;
	}table[M<<2];
	int head[M],tot=1;
	void Add(int x,int y)
	{
		table[++tot].to=y;
		table[tot].next=head[x];
		head[x]=tot;
	}
	void Tarjan(int x)
	{
		using namespace Cactus_Graph;

		static int dpt[M],low[M],T,stack[M],top;
		int i;
		
		dpt[x]=low[x]=++T;
		stack[++top]=x;

		for(i=head[x];i;i=table[i].next)
		{
			if(dpt[table[i].to])
				low[x]=min(low[x],dpt[table[i].to]);
			else
			{
				Tarjan(table[i].to),low[x]=min(low[x],low[table[i].to]);
				if(dpt[x]==low[table[i].to])
				{
					int t;++cnt;
					do{
						t=stack[top--];
						belong[t].push_back(cnt);
						contain[cnt].push_back(t);
						size[cnt]++;
					}while(t!=table[i].to);
					belong[x].push_back(cnt);
					contain[cnt].push_back(x);
					size[cnt]++;
				}
			}
		}
	}
}
int main()
{
	using namespace Origin_Graph;
	using namespace Cactus_Graph;
	int i,j,x,y;
	cin>>n>>m;
	for(i=1;i<=m;i++)
	{
		scanf("%d%d",&k,&x);
		for(j=2;j<=k;j++)
		{
			scanf("%d",&y);
			Add(x,y);Add(y,x);
			x=y;
		}
	}
	Tarjan(1);
	vector<int>::iterator it;
	for(it=belong[1].begin();it!=belong[1].end();it++)
		f[1]=max(f[1],Cactus_DP(*it));
	cout<<ans<<endl;
	return 0;
}