#include <cstdio>
#include <cmath>
#include <ctime>
#include <cctype>
#include <cstring>
#include <cstdlib>
#include <cassert>
#include <set>
#include <map>
#include <queue>
#include <vector>
#include <bitset>
#include <complex>
#include <iostream>
#include <algorithm>
#define fi first
#define se second
#define pb push_back
#define y1 kjfasiv
#define lowbit(x) (x&-x)
#define debug(x) cout<<#x<<"="<<x<<endl
#pragma comment(linker, "/STACK:1024000000,1024000000")
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef unsigned int uint;
typedef pair<int,int> pii;
typedef pair<ll,ll> pll;
typedef vector<int> veci;
typedef complex<double> Com;
const int mod=(int)1e9+7,INF=0x7fffffff,rx[]={-1,0,1,0},ry[]={0,1,0,-1};
const double pi=acos(-1.0),eps=1e-8;
template<class T>void rd(T &res){
    res=0;
    char c;
    while(c=getchar(),c<48);
    do res=res*10+(c^48);
        while(c=getchar(),c>47);
}
template<class T>void rec_print(T x){
    if(!x)return;
    rec_print(x/10);
    putchar(x%10^48);
}
template<class T>void print(T x){
    if(!x)putchar('0');
    else rec_print(x);
}
template<class T>inline void Max(T &a,T b){
    if(b>a)a=b;
}
template<class T>inline void Min(T &a,T b){
    if(b<a)a=b;
}
inline void mod_add(int &a,int b){
    if((a+=b)>=mod)a-=mod;
}
int fast_mod_pow(int a,int b){
    int res=1;
    for(;b;b>>=1,a=1ll*a*a%mod)
        if(b&1)res=1ll*res*a%mod;
    return res;
}
const int N=(int)1e3+5;
int num[N][N],a[N][N],b[N][N],c[N][N],d[N][N];
uint sum[N][N];
pii stk[N];
inline uint sigma(uint n){
    return n*(n+1)>>1;
}
void init(){
    for(uint i=1;i<N;i++)
        for(uint j=1;j<N;j++)
            sum[i][j]=sum[i-1][j]+sigma(i+j-1)-sigma(i-1);
}
void solve(){
    int n,m;
    rd(n);rd(m);
    for(int i=1;i<=n;++i)
        for(int j=1;j<=m;++j)
            rd(num[i][j]);
    for(int i=1;i<=n;i++){
        int top=0;
        stk[top].fi=0;
        for(int j=1;j<=m;j++){
            while(top&&stk[top].se>num[i][j])--top;
            a[i][j]=j-stk[top].fi;
            stk[++top]=pii(j,num[i][j]);
        }
        stk[top=0].fi=m+1;
        for(int j=m;j>=1;j--){
            while(top&&stk[top].se>num[i][j])--top;
            b[i][j]=stk[top].fi-j;
            stk[++top]=pii(j,num[i][j]);
        }
    }
    for(int i=1;i<=m;i++){
        int top=0;
        stk[top].fi=0;
        for(int j=1;j<=n;j++){
            while(top&&stk[top].se<num[j][i])--top;
            c[j][i]=j-stk[top].fi;
            stk[++top]=pii(j,num[j][i]);
        }
        stk[top=0].fi=n+1;
        for(int j=n;j>=1;j--){
            while(top&&stk[top].se<num[j][i])--top;
            d[j][i]=stk[top].fi-j;
            stk[++top]=pii(j,num[j][i]);
        }
    }
    uint ans=0;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            ans+=num[i][j]*sum[a[i][j]][b[i][j]]*sum[c[i][j]][d[i][j]];
    print(ans);
    putchar('\n');
}
int main(){
    int cas;
    init();
    for(rd(cas);cas--;)solve();
    return 0;
}