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

struct nod{
	int to,nxt,w;
}e[M*2];
int head[N],cnt=1;

void add(int u,int v,int w){
	e[++cnt].to=v;
	e[cnt].w=w;
	e[cnt].nxt=head[u];
	head[u]=cnt;
}
void qadd(int u,int v,int w){
//	cout<<"edg "<<u<<" "<<v<<" "<<w<<'\n';
	add(u,v,w);
	add(v,u,0);
}

int s,t;
int dis[N],cur[N];

bool bfs(){
	memset(dis,-1,sizeof(dis));
	queue <int> q;
	q.push(s);
	dis[s]=0;
	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].w>0 && dis[v]==-1){
				dis[v]=dis[u]+1;
				q.push(v);
			}
		}
	}
	return (dis[t]!=-1);
}

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].w && dis[v]==dis[u]+1){
			int tmp=dfs(v,min(e[i].w,delta));
			e[i].w-=tmp,e[i^1].w+=tmp;
			delta-=tmp;
			if(!delta) break;
		}
	}
	return flow-delta;
}

int dinic(){
	int ans=0;
	for(int i=2;i<=cnt;i+=2){
		e[i].w+=e[i^1].w;
		e[i^1].w=0;
	}
	while(bfs()){
		for(int i=1;i<=t;i++) cur[i]=head[i];
		int tmp;
		while(tmp=dfs(s,INF)) ans+=tmp;
	}
	return ans;
}

int T,n;
int a[N][N];
int sav[N],sav2[N],sav3[N][N];
int ans[N][N];

int main(){
	cin>>T;
	while(T--){
		cin>>n;
		for(int i=1;i<=n;i++)
			for(int j=1;j<=n;j++)
				scanf("%d",&a[i][j]);
		s=n*2+1,t=n*2+2;
		
		for(int i=1;i<=n;i++){
			qadd(s,i,0);
			sav[i]=cnt;
		}
		for(int i=1;i<=n;i++){
			qadd(i+n,t,0);
			sav2[i]=cnt;
		}
		
		for(int i=1;i<=n;i++){
			for(int j=1;j<=n;j++){
				qadd(i,a[i][j]+n,1);
				sav3[i][j]=cnt;
			}
		}
		
		for(int t=1;t<=n;t++){
			for(int i=1;i<=n;i++){
				int x=sav[i]^1;
				e[x].w=1,e[x^1].w=0;
			}
			for(int i=1;i<=n;i++){
				int x=sav2[i]^1;
				e[x].w=1,e[x^1].w=0;
			}
			int alc=dinic();
			for(int i=1;i<=n;i++){
				for(int j=1;j<=n;j++){
					int x=sav3[i][j];
					if(e[x].w){
						ans[t][i]=j;
//						cout<<j<<" ";
						e[x].w=0;
					}
				}
			}
//			cout<<'\n';
		}
		
//		for(int i=1;i<=n;i++){
//			for(int j=1;j<=n;j++){
//				cout<<ans[i][j]<<" ";
//			}
//			cout<<'\n';
//		}
		
		cout<<n*(n-1)/2<<'\n';
		for(int i=1;i<=n;i++){
			for(int j=i+1;j<=n;j++){
				cout<<i<<" "<<ans[j][i]<<" "<<j<<" "<<ans[i][j]<<'\n';
			}
		}
		
		cnt=1;
		for(int i=1;i<=t;i++) head[i]=0;
	}
	
}