#include<bits/stdc++.h>

using namespace std;
using ld=long double;

void fix(vector<ld> &v){
  ld s=0.0;
  for(auto &nx : v){
    s+=nx;
  }
  s=sqrt(s);
  for(auto &nx : v){
    nx/=s;
  }
}

int main(){
  ios::sync_with_stdio(false);
  cin.tie(nullptr);
  int t;
  cin >> t;
  while(t--){
    int n,m;
    cin >> n >> m;
    vector<vector<ld>> a(n,vector<ld>(m,0));
    for(auto &nx : a){
      for(auto &ny : nx){cin >> ny;}
    }
    ld alpha=1.0;
    ld beta=0.99;
    vector<ld> alice(n,((ld)1.0)/((ld)n));
    vector<ld> bob(m,((ld)1.0)/((ld)m));
    for(int tr=0;tr<2000;tr++){
      {
        // alice
        vector<ld> as(n,0.0);
        for(int i=0;i<n;i++){
          for(int j=0;j<m;j++){
            as[i]+=bob[j]*a[i][j];
          }
        }
        int tg=0;
        for(int i=0;i<n;i++){
          if(as[tg]<as[i]){tg=i;}
        }
        alice[tg]+=(1.0-alice[tg])*alpha;
        for(int i=0;i<n;i++){
          if(i!=tg){alice[i]*=(1.0-alpha);}
        }
      }
      fix(alice);
      {
        // bob
        vector<ld> bs(m,0.0);
        for(int i=0;i<n;i++){
          for(int j=0;j<m;j++){
            bs[j]+=alice[i]*a[i][j];
          }
        }
        int tg=0;
        for(int i=0;i<m;i++){
          if(bs[tg]>bs[i]){tg=i;}
        }
        bob[tg]+=(1.0-bob[tg])*alpha;
        for(int i=0;i<m;i++){
          if(i!=tg){bob[i]*=(1.0-alpha);}
        }
      }
      alpha*=beta;
      fix(bob);
    }
    ld res=0.0;
    for(int i=0;i<n;i++){
      for(int j=0;j<m;j++){
        res+=alice[i]*bob[j]*a[i][j];
      }
    }
    cout << fixed << setprecision(20) << res << "\n";
  }
  return 0;
}