#include<bits/stdc++.h>
using namespace std;
namespace gza{
	#define pb push_back
	#define MT int TTT=R;while(TTT--)
	#define pc putchar
	#define R read()
	#define fo(i,a,b) for(signed i=a;i<=b;i++)
	#define rep(i,a,b) for(signed i=a;i>=b;i--)
	#define m1(a,b) memset(a,b,sizeof a)
	namespace IO
	{
		inline int read()
		{
		    int x=0;
		    char ch=getchar();
		    bool f=0;
		    while(!isdigit(ch)){if(ch=='-') f=1;ch=getchar();}
		    while(isdigit(ch)) x=(x<<1)+(x<<3)+(ch^48),ch=getchar();
		    if(f) x=-x;
		    return x;
		}
		template<typename T> inline void write(T x)
		{
		    if(x<0) pc('-'),x=-x;
		    if(x>9) write(x/10);
		    pc(x%10+'0');
		}
	};
	using namespace IO;

	const int N=5e4+10,M=16;
	int n,s,t,ans;
	vector<int> g[N];
	int deg[N];
	int fa[N],dep[N],siz[N];
	int son[N],top[N];
	int dfn[N],cnt;
	int p[N],ansp[N];
	bool vis[N];
	mt19937 myrand(time(0));
	void dfs1(int u,int f)
	{
		fa[u]=f,dep[u]=dep[f]+1,siz[u]=1,son[u]=0;
		for(auto j:g[u]) if(j!=f)
		{
			dfs1(j,u);
			siz[u]+=siz[j];
			if(siz[j]>siz[son[u]]) son[u]=j;
		}
	}
	void dfs2(int u,int f)
	{
		top[u]=f,dfn[u]=++cnt;
		if(son[u]) dfs2(son[u],f);
		for(auto j:g[u]) if(j!=son[u]&&j!=fa[u]) dfs2(j,j);
	}
	int lca(int u,int v)
	{
		while(top[u]!=top[v])
		{
			if(dep[top[u]]<dep[top[v]]) swap(u,v);
			u=fa[top[u]];
		}
		if(dep[u]>dep[v]) swap(u,v);
		return u;
	}
	int dis(int u,int v)
	{
		return dep[u]+dep[v]-2*dep[lca(u,v)];
	}
	int calc()
	{
		int res=0;
		fo(i,1,n-1) res^=dis(p[i],p[i+1]);
		return res;
	}
	int get(int x)
	{
		return dis(p[x],p[x-1])^dis(p[x],p[x+1]);
	}
	void dfs3(int x,int sum)
	{
		if(x==n)
		{
			int tt=(sum^dis(p[n-1],p[n]));
			// cout<<tt<<endl;
			if(tt<ans)
			{
				ans=tt;
				fo(i,1,n) ansp[i]=p[i];
			}
			return;
		}
		fo(i,1,n) if(!vis[i])
		{
			p[x]=i;
			vis[i]=1;
			dfs3(x+1,(sum^dis(p[x-1],p[x])));
			vis[i]=0;
		}
	}
	void solve(){
		n=R,s=R,t=R;
		fo(i,1,n) g[i].clear(),deg[i]=fa[i]=dep[i]=siz[i]=son[i]=top[i]=dfn[i]=p[i]=ansp[i]=vis[i]=0;
		cnt=0;
		fo(i,2,n)
		{
			int u=R,v=R;
			g[u].pb(v),g[v].pb(u);
			deg[u]++,deg[v]++;
		}
		dfs1(1,0),dfs2(1,1);
		int maxd=(*max_element(deg+1,deg+n+1));
		p[1]=s,p[n]=t;
		fo(i,1,n) if(i!=s&&i!=t) p[i-(i>s)-(i>t)+1]=i;
		if(maxd==n-1)
		{
			fo(i,1,n) write(p[i]),pc(' ');
			puts("");
		}
		else if(n<=9)
		{
			ans=calc();
			fo(i,1,n) ansp[i]=p[i];
			vis[s]=vis[t]=1;
			dfs3(2,0);
			fo(i,1,n) write(ansp[i]),pc(' ');
			puts("");
		}
		else
		{
			ans=calc();
			while(ans>1)
			{
				int x=myrand()%(n-2)+2,y=myrand()%(n-2)+2;
				while(x==y) y=myrand()%(n-2)+2;
				ans^=get(x)^get(y);
				swap(p[x],p[y]);
				ans^=get(x)^get(y);
			}
			fo(i,1,n) write(p[i]),pc(' ');
			puts("");
		}
	}
	void main(){
		srand(time(0));
		MT solve();
	}
}
signed main(){
	
	gza::main();
}