#include<string.h>
#include<math.h>
#include<stdlib.h>
#include<stdio.h>
#include<time.h>
#include<iostream>
#include<algorithm>
#include<string>
#include<vector>
#include<queue>
#include<stack>
#include<map>
#include<set>
using namespace std;
typedef unsigned ud;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> pii;
typedef vector<int> vec;
typedef priority_queue<int> pq;

#define pb push_back
#define ph push
#define fi first
#define se second
const int M=1005;
ud dp[M][M],sum[M<<1];
ud A[M][M],a[M][M],x[M][M],b[M][M],y[M][M],pos[M],stk[M];
inline void Max(int &a,int b){if(a<b)a=b;}
inline void Min(int &a,int b){if(a>b||a==-1)a=b;}
template<class T>void rd(T &a){
    a=0;char c;
    while(c=getchar(),!isdigit(c));
    do a=a*10+(c^48);
        while(c=getchar(),isdigit(c));
}
template<class T>void nt(T x){
    if(!x)return;
    nt(x/10);
    putchar(48+x%10);
}
template<class T>void pt(T x){
    if(!x)putchar('0');
    else nt(x);
}
inline int Mod_Pow(int x,int a,int mod){
    int res=1;
    for(int i=0;(1ll<<i)<=a;++i){
        if(a&1<<i)res=1ll*res*x%mod;
        x=1ll*x*x%mod;
    }
    return res;
}
void gao(){
    int n,m;cin>>n>>m;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)rd(A[i][j]);
//////////////////////////////////////////////////////////////////////////////////////////////////////
    pos[0]=0;
    for(int i=1;i<=n;++i){
        int top=0;
        for(int j=1;j<=m;++j){
            for(;top&&stk[top]>A[i][j];)top--;
            a[i][j]=j-pos[top];
            stk[++top]=A[i][j];
            pos[top]=j;
        }
    }
//////////////////////////////////////////////////////////////////////////////////////////////////////
    pos[0]=0;
    for(int i=1;i<=m;++i){
        int top=0;
        for(int j=1;j<=n;++j){
            for(;top&&stk[top]<A[j][i];)top--;
            b[j][i]=j-pos[top];
            stk[++top]=A[j][i];
            pos[top]=j;
        }
    }
//////////////////////////////////////////////////////////////////////////////////////////////////////
    pos[0]=m+1;
    for(int i=1;i<=n;++i){
        int top=0;
        for(int j=m;j;--j){
            for(;top&&stk[top]>A[i][j];)top--;
            x[i][j]=pos[top]-j;
            stk[++top]=A[i][j];
            pos[top]=j;
        }
    }
//////////////////////////////////////////////////////////////////////////////////////////////////////
    pos[0]=n+1;
    for(int i=1;i<=m;++i){
        int top=0;
        for(int j=n;j;--j){
            for(;top&&stk[top]<A[j][i];)top--;
            y[j][i]=pos[top]-j;
            stk[++top]=A[j][i];
            pos[top]=j;
        }
    }
//////////////////////////////////////////////////////////////////////////////////////////////////////
    ud ans=0;
    for(int i=1;i<=n;++i)
        for(int j=1;j<=m;++j){
            ud t=dp[a[i][j]][x[i][j]]*dp[b[i][j]][y[i][j]];
            ans+=t*A[i][j];
        }
    pt(ans),putchar('\n');
}
inline void pret(){
    for(int i=1;i<=2000;++i)
        sum[i]+=sum[i-1]+i;
    for(int i=1;i<=1000;++i)
        for(int j=1;j<=1000;++j)
            dp[i][j]=dp[i-1][j]+sum[i+j-1]-sum[i-1];
}
int main(){
    int _;
    for(pret(),cin>>_;_--;)gao();
}