#include<bits/stdc++.h>
using namespace std;
const int N=440,M=5500,INF=1e9;


struct abc{
	int to,nxt,f;
}e[M*4];
int head[N],cnt=1;
void add(int u,int v,int f){
	e[++cnt]={v,head[u],f};
	head[u]=cnt;
}
void qadd(int u,int v,int f){
	add(u,v,f);
	add(v,u,0);
}

int S,T,sav[M];

int dis[N];
queue <int> q;

bool bfs(){
	for(int i=1;i<=T;i++) dis[i]=INF;
	dis[S]=0;
	q.push(S);
	while(q.size()){
		int u=q.front(); q.pop();
		for(int i=head[u];i;i=e[i].nxt){
			int v=e[i].to;
			if(e[i].f && dis[v]==INF){
				dis[v]=dis[u]+1;
				q.push(v);
			}
		}
	}
	return (dis[T]<INF);
}

int cur[N];

int dfs(int u,int flow){
	if(u==T) return flow;
	int delta=flow;
	for(int i=cur[u];i;i=e[i].nxt){
		cur[u]=i;
		int v=e[i].to;
		if(e[i].f && dis[v]==dis[u]+1){
			int x=dfs(v,min(delta,e[i].f));
			e[i].f-=x,e[i^1].f+=x;
			delta-=x;
			if(!delta) break;
		}
	}
	return flow-delta;
}
void clr(){
	for(int i=1;i<=T;i++) head[i]=0;
	cnt=1;
}
int dinic(){
	int tot=0;
	while(bfs()){
		for(int i=1;i<=T;i++) cur[i]=head[i];
		tot+=dfs(S,INF);
	}
	return tot;
}


int T_,n,m;
int U[M],V[M];
int in[N],out[N];

vector <int> G[N];
stack <int> pat;
void dfs(int u){
	for(int i=cur[u];i<(int)G[u].size();i=cur[u]){
		cur[u]++;
		dfs(G[u][i]);
	}
	pat.push(u);
}

int main(){
	ios::sync_with_stdio(false);
	cin.tie(0);
	
	cin>>T_;
	for(int T__=1;T__<=T_;T__++){
		cout<<"Case #"<<T__<<": ";
		cin>>n>>m;
		for(int i=1;i<=n;i++) in[i]=out[i]=0,G[i].clear();
		for(int i=1;i<=m;i++){
			int u,v; cin>>u>>v;
			U[i]=u,V[i]=v;
			out[u]++,in[v]++;
		}
		
		int L=-INF;
		for(int i=1;i<=n;i++) L=max({L,in[i],out[i]});
		
		bool fl=0;
		for(int d=m;d<=m;d++){
			clr();
			S=2*n+1,T=2*n+2;
			for(int i=1;i<=n;i++) qadd(S,i,d-out[i]);
			for(int i=1;i<=n;i++) qadd(i+n,T,d-in[i]);
			for(int i=1;i<=m;i++){
				int u=U[i],v=V[i];
				qadd(u,v+n,INF);
				sav[i]=cnt;
			}
			
			int val=dinic();
			if(val==d*n-m){
//				cout<<'\n';
//				cout<<" "<<"d: "<<d<<'\n';
//				for(int i=1;i<=m;i++){
//					cout<<e[sav[i]].f+1<<" ";
//				}
				for(int i=1;i<=m;i++){
					int u=U[i],v=V[i],c=e[sav[i]].f+1;
					while(c--) G[u].push_back(v);//,cout<<"edg "<<u<<" "<<v<<'\n';
				}
				for(int i=1;i<=n;i++) cur[i]=0;
				dfs(1);
				if(pat.size()!=d*n+1){
					while(pat.size()) pat.pop();
					break;
				}
				cout<<d*n+1<<'\n';
				while(pat.size()) cout<<pat.top()<<" ",pat.pop();
				cout<<'\n';
				fl=1;
				break;
			}
			
		}
		if(!fl) cout<<"IMPOSSIBLE\n";
	}
	
}