#include<iostream>
#include<cstdio>
#include<cmath>
#define maxn 200005
using namespace std;
int T,n,m,u,v,dep[maxn],maxpos,s[maxn],top,ans[maxn],num=0; bool flag=0;
struct node{int to,nex;}a[maxn]; int head[maxn],cnt=0;
void add(int from,int to){a[++cnt].to=to; a[cnt].nex=head[from]; head[from]=cnt;}
int f[maxn]; int getfa(int x){return f[x]==x?x:f[x]=getfa(f[x]);}
void dfs1(int p,int fa){
	dep[p]=dep[fa]+1; maxpos=(dep[p]>dep[maxpos]?p:maxpos);
	for(int i=head[p];i;i=a[i].nex) if(a[i].to!=fa) dfs1(a[i].to,p);
}
void dfs2(int p,int fa){
	s[++top]=p; if(p==maxpos){flag=1; return;}
	for(int i=head[p];i;i=a[i].nex){if(a[i].to!=fa) dfs2(a[i].to,p); if(flag) return;}
}
int main(){
	scanf("%d",&T); while(T--){
		scanf("%d%d",&n,&m); for(int i=1;i<=n;i++) f[i]=i,head[i]=0; cnt=num=maxpos=flag=top=0;
		for(int i=1;i<=m;i++){
			scanf("%d%d",&u,&v); int r1=getfa(u),r2=getfa(v); if(r1==r2) continue;
			f[r2]=r1; add(u,v);
		} dfs1(1,0); dfs2(1,0);
		int thres=int(sqrt(n)); if(thres*thres!=n) thres++; while(top>0){
			for(int i=1;i<=thres;i++){top--; if(top==0) top++;} ans[++num]=s[top]; if(top==1) break;
		} printf("%d\n",num); for(int i=1;i<=num;i++) printf("%d ",ans[i]); printf("\n");
	} return 0;
}