NOI2007.Day2.T3.追捕盗贼

这道题是看的郑暾大牛的论文《平衡思想》

里面说这道题是树的Search Number问题,有O(n)解法,恰好我这里有数据&标程,一看标程8.5k……

郑暾大牛给出了一种DP构造解

虽然不是正解,但是在大部分情况下可以保证与最优解一样

而且代码较短(我只写了120行),性价比较正解高出太多,在考试的时候不失为一种好的方法

因为是在树上,树有个性质就是每个点都是割点,所以每个点都可以把树划分成若干不连通的块,然后递归进去做

算法会枚举一个点作为根,对于子树则默认子树的初始分割点就是根(这是有时达不到正解的原因)(否则复杂度会大大增加)

因此复杂度O(n^2),实测较快,时限很宽松




//Lib
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<ctime>

#include<iostream>
#include<algorithm>
#include<vector>
#include<string>
#include<queue>
using namespace std;
//Macro
#define rep(i,a,b) for(int i=a,tt=b;i<=tt;++i)
#define drep(i,a,b) for(int i=a,tt=b;i>=tt;--i)
#define erep(i,e,x) for(int i=x;i;i=e[i].next)
#define irep(i,x) for(__typedef(x.begin()) i=x.begin();i!=x.end();i++)
#define read() (strtol(ipos,&ipos,10))
#define sqr(x) ((x)*(x))
#define pb push_back
#define PS system("pause");
typedef long long ll;
typedef pair<int,int> pii;
const int oo=~0U>>1;
const double inf=1e20;
const double eps=1e-6;
string name="",in=".in",out=".out";
//Var
struct W
{
	int x,y;char order;
}way[20008];
struct E
{
	int next,node;
}e[2008];
int cnt,n,ans=oo,as,tot;
int h[1008],f[1008];
void add(int a,int b){e[++tot].next=h[a];e[tot].node=b;h[a]=tot;}
void Land(int x){way[++cnt].order='L';way[cnt].x=x;}
void Back(int x){way[++cnt].order='B';way[cnt].x=x;}
void Move(int x,int y){way[++cnt].order='M';way[cnt].x=x;way[cnt].y=y;}
void TDP(int x,int fa)
{
	int maxx=0,msize=0,son=0,y;
	erep(i,e,h[x])
	{
		if((y=e[i].node)!=fa)
		{
			TDP(y,x);
			if(f[y]==maxx)msize++;
			if(f[y]>maxx)maxx=f[y],msize=1;
			++son;
		}
	}
	if(msize>1)f[x]=maxx+1;else f[x]=maxx;
	if(son==0)f[x]++;
}
void Solve(int x,int fa,int t,int d)
{
	int maxx=0,msize=0,son=0,mnum=-1,y;
	erep(i,e,h[x])
	{
		if((y=e[i].node)!=fa)
		{
			Solve(y,x,0,0);
			if(f[y]==maxx)msize++;
			if(f[y]>maxx)maxx=f[y],msize=1,mnum=y;
			++son;
		}
	}
	if(d==1)
	{
		erep(i,e,h[x])
			if((y=e[i].node)!=fa)
			{
				if(y!=mnum)
				{
					Land(x);
					Move(x,y);
					Solve(y,x,1,1);
				}
			}
		if(son==0){Back(x);return;}
		Move(x,mnum);Solve(mnum,x,1,1);
	}
	if(msize>1)f[x]=maxx+1;else f[x]=maxx;
	if(son==0)f[x]++;
}
void Work()
{
	scanf("%d",&n);int a,b;
	rep(i,1,n-1)scanf("%d%d",&a,&b),add(a,b),add(b,a);
	rep(i,1,n)
	{
		memset(f,0,sizeof f);
		TDP(i,0);
		if(f[i]<ans)ans=f[i],as=i;
	}
	printf("%d\n",ans);
	memset(f,0,sizeof f);
	Land(as);Solve(as,0,1,1);
}
void Outit()
{
	printf("%d\n",cnt);
	rep(i,1,cnt)
	{
		if(way[i].order=='M')printf("%c %d %d\n",way[i].order,way[i].x,way[i].y);
		else printf("%c %d\n",way[i].order,way[i].x);
	}
}
int main()
{
//	freopen((name+in).c_str(),"r",stdin);
//	freopen((name+out).c_str(),"w",stdout);
//	Init();
	Work();
	Outit();
	return 0;
}


阅读更多

没有更多推荐了,返回首页