#include <bits/stdc++.h>
using namespace std;
int Qread()
{
    int x=0;bool f=false;char ch=getchar();
    while(ch<'0'||ch>'9') f|=(ch=='-'),ch=getchar();
    while(ch>='0'&&ch<='9') x=x*10+(ch^48),ch=getchar();
    return f?-x:x;
}
void Qwrite(int a)
{
	if(a<10) return putchar('0'+a),void();
	return Qwrite(a/10),putchar('0'+a%10),void();
}
namespace flow{
    int p[120],edge_cnt,poi_cnt,S,T;
    struct Edge{
        int nxt,poi,val,bas;
    }l[10010];
    void reset(int n)
    {
        memset(p,0,sizeof(p));
        memset(l,0,sizeof(l));
        edge_cnt=1,poi_cnt=2*n+3;
        S=2*n+1,T=2*n+2;
    }
    int add_edge(int u,int v)
    {
        l[++edge_cnt].nxt=p[u],l[edge_cnt].poi=v,l[edge_cnt].val=1,l[edge_cnt].bas=1,p[u]=edge_cnt;
        l[++edge_cnt].nxt=p[v],l[edge_cnt].poi=u,l[edge_cnt].val=0,l[edge_cnt].bas=0,p[v]=edge_cnt;
        return edge_cnt-1;
    }
    int dep[120],cur[120];
    queue<int> Q;
    void bk_fl()
    {
        for(int i=1;i<=edge_cnt;i++)
            l[i].val=l[i].bas;
    }
    bool bfs()
    {
        memset(dep,0,poi_cnt<<2);
        dep[S]=1,Q.push(S);
        while(!Q.empty())
        {
            int rea=Q.front();Q.pop();
            for(int k=cur[rea]=p[rea];k;k=l[k].nxt)
                if(!dep[l[k].poi]&&l[k].val)
                    dep[l[k].poi]=dep[rea]+1,Q.push(l[k].poi);
        }
        return dep[T];
    }
    int dfs(int a,int flo)
    {
        if(a==T||flo==0) return flo;
        int ret=0,rem=0;
        for(int k=cur[a];k&&flo;k=l[k].nxt)
        {
            cur[a]=k;
            if(dep[l[k].poi]!=dep[a]+1||!l[k].val) continue;
            rem=dfs(l[k].poi,min(flo,l[k].val));
            flo-=rem,l[k].val-=rem;
            ret+=rem,l[k^1].val+=rem;
        }
        return ret;
    }
    void Dinic()
    {
        while(bfs()) dfs(S,20070610);
    }
}
struct ANS{int num[60],a;}rem;
vector<ANS> T;
int n,tot;
long long r[60],c[60];
int a[60][60],SS[60],TT[60];
int tk[60][60];
void solve()
{
    T.clear();
    memset(r,0,sizeof(r));
    memset(c,0,sizeof(c));
    n=Qread();
    flow::reset(n);
    tot=n*n;
    for(int i=1;i<=n;i++)
    for(int j=1;j<=n;j++)
    {
        a[i][j]=Qread();
        r[i]+=a[i][j],c[j]+=a[i][j];
    }
    for(int i=1;i<=n;i++)
    for(int j=1;j<=n;j++)
        if(a[i][j]<0)
        {
            printf("-1\n");
            return;
        }
    if(r[1]!=c[1])
    {
        printf("-1\n");
        return;
    }
    for(int i=2;i<=n;i++)
        if(r[i]!=r[i-1]||c[i]!=c[i-1])
        {
            printf("-1\n");
            return;
        }
    for(int i=1;i<=n;i++)
    for(int j=1;j<=n;j++)
        tk[i][j]=flow::add_edge(i,n+j);
    for(int i=1;i<=n;i++)
        SS[i]=flow::add_edge(flow::S,i),TT[i]=flow::add_edge(n+i,flow::T);
    for(int i=1;i<=n;i++)
    for(int j=1;j<=n;j++)
        if(!a[i][j])
            flow::l[tk[i][j]].bas=0,flow::l[tk[i][j]].val=0,tot--;
    while(tot)
    {
        flow::Dinic();
        rem.a=1e9+10;
        for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
            if(flow::l[tk[i][j]].bas&&!flow::l[tk[i][j]].val)
            {
                rem.num[i]=j;
                rem.a=min(rem.a,a[i][j]);
            }
        for(int i=1;i<=n;i++)
            if(!(a[i][rem.num[i]]-=rem.a))
            {
                flow::l[tk[i][rem.num[i]]].bas=0,tot--;
            	flow::l[tk[i][rem.num[i]]^1].val=0;
            	flow::l[SS[i]].val=1,flow::l[SS[i]^1].val=0;
            	flow::l[TT[rem.num[i]]].val=1,flow::l[TT[rem.num[i]]^1].val=0;
			}
        T.push_back(rem);
    }
    printf("%d\n",(int)T.size());
    for(ANS &g:T)
    {
        for(int i=1;i<=n;i++)
        	Qwrite(g.num[i]),putchar(' ');
		Qwrite(g.a),putchar('\n');
    }
}
int main()
{
    int TTTTT=Qread();
    while(TTTTT--) solve();
    return 0;
}