#pragma GCC optimize("Ofast","inline","-ffast-math")
#pragma GCC target("avx,sse2,sse3,sse4,mmx")
#include <bits/stdc++.h>
#include<ext/pb_ds/assoc_container.hpp>
#include<ext/pb_ds/hash_policy.hpp>
//#define int long long
using namespace __gnu_pbds;
using namespace std;

struct custom_hash {
    static uint64_t splitmix64(uint64_t x) {
        x += 0x9e3779b97f4a7c15;
        x = (x ^ (x >> 30)) * 0xbf58476d1ce4e5b9;
        x = (x ^ (x >> 27)) * 0x94d049bb133111eb;
        return x ^ (x >> 31);
    }

    size_t operator()(uint64_t x) const {
        static const uint64_t FIXED_RANDOM = chrono::steady_clock::now().time_since_epoch().count();
        return splitmix64(x + FIXED_RANDOM);
    }
};

int n,X,vis[110],num[110],nxt[110],nt[110],ok[110][110];
vector<int> v[110],g[110],res;

inline void dfs(int x)
{
	vis[x]=1;
	for(int i=0;i<(int)v[x].size();i++) if(!vis[v[x][i]]){
		dfs(v[x][i]);
	}
	res.push_back(x);
}

inline void rdfs(int x,int id)
{
	vis[x]=1;num[x]=id;
	for(int i=0;i<(int)g[x].size();i++) if(!vis[g[x][i]]){
		rdfs(g[x][i],id);
	}
}

signed main()
{
	ios::sync_with_stdio(false);cin.tie(0);
	int t;cin>>t;
	int tot=0;
	while(t--){
		cin>>n>>X;X++;tot++;
		for(int i=1;i<=n;i++) v[i].clear(),g[i].clear();
		for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) ok[i][j]=0;
		for(int i=1;i<=n;i++){
			for(int j=1;j<=n;j++){
				char ch;cin>>ch;
				if(ch=='Y'){
					ok[i][j]=1;
					v[i].push_back(j);
					g[j].push_back(i);
				}
			}
		}
		for(int i=1;i<=n;i++) vis[i]=0;
		res.clear();
		for(int i=1;i<=n;i++) if(!vis[i]) dfs(i);
		for(int i=1;i<=n;i++) vis[i]=0,num[i]=0;
		int cnt=0;
		for(int i=(int)res.size()-1;i>=0;i--) if(!vis[res[i]]){
			cnt++;rdfs(res[i],cnt);
		}
		vector<int> vec,nw;
		for(int i=1;i<=n;i++){
			if(num[i]==1) vec.push_back(i);
			else nw.push_back(i);
		}
		bool flag=false;
		for(int i=0;i<(int)vec.size();i++) if(vec[i]==X) flag=true;
		if(!flag){
			cout<<"IMPOSSIBLE\n";continue;
		}
		if((int)vec.size()==1){
			cout<<"Case #"<<tot<<": ";
			for(int i=0;i<(int)vec.size();i++) cout<<vec[i]-1<<' ';
			for(int i=0;i<(int)nw.size();i++) cout<<nw[i]-1<<' ';
			cout<<'\n';
			continue;
		}
		int l=vec[0],r=vec[0];
		for(int i=1;i<(int)vec.size();i++){
			int u=vec[i];
			if(ok[u][l]) nxt[u]=l,l=u;
			else if(ok[r][u]) nxt[r]=u,r=u;
			else{
				for(int j=l;j!=r;j=nxt[j]){
					if(ok[j][u]&&ok[u][nxt[j]]){
						nxt[u]=nxt[j];nxt[j]=u;
						break;
					}
				}
			}
		}
		for(int i=l;i!=r;i=nxt[i]) nt[i]=nxt[i];
		int now=-1;
		for(int i=l;;i=nxt[i]){
			if(now==-1){
				if(ok[i][l]) nt[i]=l,now=i;
			}
			else{
				int s=nt[now];
				if(ok[i][s]){
					for(int j=now;j!=i;j=nxt[j]) nt[j]=nxt[j];
					nt[i]=s;now=i;
				}
				else{
					for(int j=nt[s],pre=s;j!=s;j=nt[j],pre=nt[pre]){
						if(ok[i][j]){
							nt[i]=j;nt[pre]=nxt[now];now=i;
							break;
						}
					}
				}
			}
			if(i==r) break;
		}
		vector<int> cyc;
		for(int i=now;;i=nt[i]){
			cyc.push_back(i);
			if(nt[i]==now) break;
		}
		int pos=-1;
		for(int i=0;i<(int)cyc.size();i++) if(cyc[i]==X){
			pos=i;break;
		}
		cout<<"Case #"<<tot<<": ";
		for(int i=pos-1;i>=0;i--) cout<<cyc[i]-1<<' ';
		for(int i=(int)cyc.size()-1;i>=pos;i--) cout<<cyc[i]-1<<' ';
		for(int i=0;i<(int)nw.size();i++) cout<<nw[i]-1<<' ';
		cout<<'\n';
	}
	return 0;
}