#include <bits/stdc++.h>
using namespace std;
const int N=110;
int n,st;
bool o[N][N];
vector <int> e[N];
bool vis[N],ver[N];
void dfs(int x){
	vis[x]=1;
	for(int i=0;i<(int)e[x].size();i++){
		int y=e[x][i];
		if(vis[y]||ver[y])
			continue;
		dfs(y);
	}
}
inline bool check(int x){
	for(int i=1;i<=n;i++)
		vis[i]=0;
	if(ver[st]==0)
		dfs(st);
	for(int i=1;i<=n;i++)
		if(ver[i]==0&&vis[i]==0&&o[x][i]==0)
			return 0;
	return 1;
}
void solve(int id){
	printf("Case #%d: ",id);
	scanf("%d%d",&n,&st);
	st++;
	for(int i=1;i<=n;i++)
		e[i].clear();
	memset(o,0,sizeof(o));
	memset(ver,0,sizeof(ver));
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++){
			char op;
			cin>>op;
			if(op=='Y'){
				e[i].push_back(j);
				o[i][j]=1;
			}
		}
	for(int i=1;i<=n;i++)
		vis[i]=0;
	dfs(st);
	for(int i=1;i<=n;i++)
		if(vis[i]==0){
			puts("IMPOSSIBLE");
			return ;
		}
	int pre=0;
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			if(ver[j])
				continue;
			if(ver[st]&&o[st][j]==0)
				continue;
			ver[j]=1;
			if(check(j)){
				printf("%d ",j-1);
				pre=o[pre][j]?pre:j;
				break;
			}
			else
				ver[j]=0;
		}
	}
	puts("");
}
int main(){
	int T;
	scanf("%d",&T);
	for(int i=1;i<=T;i++)
		solve(i);
	return 0;
}