#include<bits/stdc++.h>
using namespace std;
const int N =505;
int n,vst[N],ecnt;
char s[N];
vector<int> to[N];
int mp[N][N],ori[N][N];
void dfs(int x,int dep){
  if(dep>3) return;
  vst[x]=1;
  for(int t:to[x]){
    if(vst[t]) continue;
    dfs(t,dep+1);
  }
}
void solve(){
  cin>>n;
  ecnt=0;
  for(int i=1;i<=n;i++){
    cin>>s;
    for(int j=0;j<n;j++){
      if(s[j]=='1'){
        to[i].push_back(j+1);
        ecnt++;
        ori[i][j+1]=mp[i][j+1]=1;
      }else{
        mp[i][j+1]=1e8;
        ori[i][j+1]=0;
      }
    }
    mp[i][i]=0;
  }
  int flag=1;
  for(int j=1;j<=n;j++){
    for(int i=1;i<=n;i++) vst[i]=0;
    dfs(j,1);
    for(int i=1;i<=n;i++) flag&=vst[i];
  }
  if(flag){
    for(int k=1;k<=n;k++){
      for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++){
          mp[i][j]=min(mp[i][j],mp[i][k]+mp[k][j]);
        }
      }
    }
    for(int i=1;i<=n;i++){
      for(int j=1;j<=n;j++){
        if(ori[i][j]!=mp[i][j]%2){
          exit(1);
        }
      }
    }
    cout<<"YES\n";
    cout<<ecnt/2<<endl;
    for(int i=1;i<=n;i++){
      for(int j:to[i]){
        if(i<j){
          cout<<i<<' '<<j<<endl;
        }
      }
    }
  }else{
    cout<<"NO\n";
  }
  for(int i=1;i<=n;i++) to[i].clear();
}
signed main(){
  ios::sync_with_stdio(false);
  int t;
  cin>>t;
  while(t--)solve();
}