洛谷 P5022 旅行

day2就蓝题起步不是要我命吗

目测数据：不都是基环树和树嘛？基环树我不会但是这个树的贪心不是裸的吗？
大手一挥：这部分分我打定了！

int outp[5100],len;
void dfs(int x,int fa)
{
	outp[++len]=x;
	list<int>q;//链表……不要问我为什么要用这玩意儿（因为没用过想试一试
	for(int k=last[x];k;k=a[k].next)
		if(a[k].y!=fa)
			q.push_back(a[k].y);
	q.sort();
	for(list<int>::iterator i=q.begin();i!=q.end();i++)
		dfs(*i,x);
}
void Main()
{
	dfs(1,0);
	for(int i=1;i<=n;i++)
		printf("%d ",outp[i]);
}

光荣地水到了60分

对于剩下的分，我们来好好观察一下题面：
$n\le 5000$
常数小的$n^2\log n$都可以过啊！
容 易 证 明整棵基环树的环中一定会有一条边是不需要走过的
博主只记了结论所以以后再来填坑啦
所以我们通过删边可以把基环树转化为一棵树，然后直接跑第一种做法

首先我们要找出这个环
通过拓扑排序解决
首先，由于是双向连通图，所以每个点的入度至少为1
我们就从入度为1的点开始，进行宽搜
观察一下样例：

可以看到，3和4这两个“环的入口”度数都为3
所以，任意一个点触碰到环上的点时，都无法将环上的点入度减为1
最后，我们就可以将这些入度不为零的点计入环中。
拓扑排序：

	void topsort()
	{
		queue<int>q;
		for(int i=1;i<=n;i++)
			if(in_degree[i]==1)q.push(i);
		while(q.size())
		{
			int x=q.front();
			for(int k=last[x];k;k=a[k].next)
				if(--in_degree[a[k].y]==1)
					q.push(a[k].y);
			q.pop();
		}
	}

还真就音译topsort呗
主函数中部分：

		topsort();
		for(int i=1;i<=n;i++)
			on_cycle[i]=in_degree[i]>1;

然后，我们就按照这些点排出环上的边：

		for(int i=1;i<=n;i++)
			if(on_cycle[i])
			{
				get_edges(i,i,0);
				break;
			}

get_edges：

	bool get_edges(int x,int ed,int fa)
	{
		for(int k=last[x],y=a[k].y;k;k=a[k].next,y=a[k].y)
			if(on_cycle[y]&&y!=fa)
			{
				cycle[++cycle_len]=k;
				if(y==ed)return 1;
				if(get_edges(y,ed,x))return 1;//这个地方是为了防止st从两边同时访问这个环，造成重复计数
			}
	}

然后就一根一根地断，将cut标记置为1
不要忘记双向边

		for(int i=1;i<=cycle_len;i++)
		{
			cut[cycle[i]]=cut[cycle[i]^1]=1;//^1是对向边，为了这句话要将len提前设为1
			top=len=0;dfs(1,0);
			if(outp[1]==0||comp(tmp,outp))memcpy(outp+1,tmp+1,n<<2);
			cut[cycle[i]]=cut[cycle[i]^1]=0;//别忘了回溯
		}

dfs：

	void dfs(int x,int fa)
	{
		const int last_top=top;
		tmp[++len]=x;
		for(int k=last[x];k;k=a[k].next)
			if(a[k].y!=fa&&!cut[k])
				sta[++top]=a[k].y;//sta作为一个排序用的数组，防止爆栈又不怕STL常数大
		const int now_top=top;
		sort(sta+last_top+1,sta+now_top+1);//贪心
		for(int i=last_top+1;i<=now_top;i++)
			dfs(sta[i],x);
	}

comp：

	bool comp(int a[],int b[])
	{
		for(int i=1;i<=n;i++)
			if(a[i]!=b[i])
				if(a[i]<b[i])return 1;
				else return 0;
		return 0;
	}

（本来为了省掉comp用了basic_string<int>，结果被其他地方的RE莫名劝退了）
最后outp数组就是答案了

代码：

#include<cstdio>
#include<cstring>
#include<queue>
#include<list>
#include<algorithm>
using namespace std;
inline void read(int &x)
{
	x=0;register char c=getchar();
	while(c<48||57<c) c=getchar();
	for(;48<=c&&c<=57;c=getchar()) x=x*10+(c&15);
}
int n,m;
struct edge
{
	int y,next;
}a[11000];
int last[5100],Len;
int in_degree[5100];
inline void ins(int x,int y)
{
	a[++Len]=(edge){y,last[x]};
	last[x]=Len;
	in_degree[y]++;
}
int outp[5100];
namespace solve1
{
	int len;
	void dfs(int x,int fa)
	{
		outp[++len]=x;
		list<int>q;
		for(int k=last[x];k;k=a[k].next)
			if(a[k].y!=fa)
				q.push_back(a[k].y);
		q.sort();
		for(list<int>::iterator i=q.begin();i!=q.end();i++)
			dfs(*i,x);
	}
	void Main()
	{
		dfs(1,0);
		for(int i=1;i<=n;i++)
			printf("%d ",outp[i]);
	}
};
namespace solve2
{
	bool on_cycle[5100];
	bool cut[11000];
	int cycle[5100],cycle_len;
	int tmp[5100],len;
	int sta[5100],top;
	void topsort()
	{
		queue<int>q;
		for(int i=1;i<=n;i++)
			if(in_degree[i]==1)q.push(i);
		while(q.size())
		{
			int x=q.front();
			for(int k=last[x];k;k=a[k].next)
				if(--in_degree[a[k].y]==1)
					q.push(a[k].y);
			q.pop();
		}
	}
	bool get_edges(int x,int ed,int fa)
	{
		for(int k=last[x],y=a[k].y;k;k=a[k].next,y=a[k].y)
			if(on_cycle[y]&&y!=fa)
			{
				cycle[++cycle_len]=k;
				if(y==ed)return 1;
				if(get_edges(y,ed,x))return 1;
			}
	}
	void dfs(int x,int fa)
	{
		const int last_top=top;
		tmp[++len]=x;
		for(int k=last[x];k;k=a[k].next)
			if(a[k].y!=fa&&!cut[k])
				sta[++top]=a[k].y;
		const int now_top=top;
		sort(sta+last_top+1,sta+now_top+1);
		for(int i=last_top+1;i<=now_top;i++)
			dfs(sta[i],x);
	}
	bool comp(int a[],int b[])
	{
		for(int i=1;i<=n;i++)
			if(a[i]!=b[i])
				if(a[i]<b[i])return 1;
				else return 0;
		return 0;
	}
	void Main()
	{
		topsort();
		for(int i=1;i<=n;i++)
			on_cycle[i]=in_degree[i]>1;
		for(int i=1;i<=n;i++)
			if(on_cycle[i])
			{
				get_edges(i,i,0);
				break;
			}
		for(int i=1;i<=cycle_len;i++)
		{
			cut[cycle[i]]=cut[cycle[i]^1]=1;
			top=len=0;dfs(1,0);
			if(outp[1]==0||comp(tmp,outp))memcpy(outp+1,tmp+1,n<<2);
			cut[cycle[i]]=cut[cycle[i]^1]=0;
		}
		for(int i=1;i<=n;i++)
			printf("%d ",outp[i]);
	}
};
int main()
{
	read(n);read(m);Len=1;
	int x,y;
	for(int i=1;i<=m;i++)
	{
		scanf("%d%d",&x,&y);
		ins(x,y);
		ins(y,x);
	}
	if(n==m)
		solve2::Main();
	else solve1::Main();
	return 0;
}