#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
#define fi first
#define se second
int n,m;
ll k;
ll ans=0;
ll r[505][505];
int idr[505][505];
ll c[505][505];
int idc[505][505];

int ec=0;
int f[2005];ll w[2005];
int eu[2005],ev[2005];
vector<pair<int,int> >adj[505];


ll pot[505];
void add(int u,int v,ll c){
	adj[u].push_back({ec,v});
	adj[v].push_back({ec+1,u});
	w[ec/2]=c;f[ec/2]=0;
	eu[ec/2]=u;
	ev[ec/2]=v;
	ec+=2;
}
ll dist[505];
bool vis[505];
int prv[505];
void iu(){
	for(int i=1; i<=n*m ;i++){
		dist[i]=1e18;
		vis[i]=false;
		prv[i]=0;
	}
	typedef pair<ll,int> pli;
	priority_queue<pli,vector<pli>,greater<pli> >pq;
	for(int i=1; i<=n ;i++){
		dist[(i-1)*m+1]=0;
		pq.push({0,(i-1)*m+1});
	}
	while(!pq.empty()){
		int x=pq.top().se;pq.pop();
		if(vis[x]) continue;
		for(auto c:adj[x]){
			int ec=c.fi/2;
			if(c.fi%2){
				ll cw=1e18;
				if(f[ec]==0) cw=w[ec];
				if(f[ec]==1) cw=-w[ec];
				if(dist[c.se]>dist[x]+cw){
					dist[c.se]=dist[x]+cw;
					prv[c.se]=c.fi;
					pq.push({dist[c.se]-pot[c.se],c.se});
				}
			}
			else{
				ll cw=1e18;
				if(f[ec]==0) cw=w[ec];
				if(f[ec]==-1) cw=-w[ec];
				if(dist[c.se]>dist[x]+cw){
					dist[c.se]=dist[x]+cw;
					prv[c.se]=c.fi;
					pq.push({dist[c.se]-pot[c.se],c.se});
				}
			}
		}
	}
	int god=m;
	for(int i=1; i<=n ;i++){
		if(dist[i*m]<dist[god]) god=i*m;
	}
	for(int i=1; i<=n*m ;i++){
		dist[i]=min(dist[i],dist[god]);
	}
	for(int i=1; i<=n*m ;i++){
		pot[i]=0;
	}
	k=dist[god]+k;
}
bool jieun(){
	for(int i=1; i<=n*m ;i++){
		dist[i]=k;
		vis[i]=false;
		prv[i]=0;
	}
	typedef pair<ll,int> pli;
	priority_queue<pli,vector<pli>,greater<pli> >pq;
	for(int i=1; i<=n ;i++){
		dist[(i-1)*m+1]=0;
		pq.push({0,(i-1)*m+1});
	}
	while(!pq.empty()){
		int x=pq.top().se;pq.pop();
		if(vis[x]) continue;
		for(auto c:adj[x]){
			int ec=c.fi/2;
			if(c.fi%2){
				ll cw=1e18;
				if(f[ec]==0) cw=w[ec];
				if(f[ec]==1) cw=-w[ec];
				if(dist[c.se]>dist[x]+cw){
					dist[c.se]=dist[x]+cw;
					prv[c.se]=c.fi;
					pq.push({dist[c.se]-pot[c.se],c.se});
				}
			}
			else{
				ll cw=1e18;
				if(f[ec]==0) cw=w[ec];
				if(f[ec]==-1) cw=-w[ec];
				if(dist[c.se]>dist[x]+cw){
					dist[c.se]=dist[x]+cw;
					prv[c.se]=c.fi;
					pq.push({dist[c.se]-pot[c.se],c.se});
				}
			}
		}
	}
	int god=m;
	for(int i=1; i<=n ;i++){
		if(dist[i*m]<dist[god]) god=i*m;
	}
	//cout << "outputting dist" << endl;
	for(int i=1; i<=n*m ;i++){
		dist[i]=min(dist[i],dist[god]);
		//cout << dist[i]-pot[i] << ' ' << pot[i] << endl;
	}
	ll frog=dist[god]-pot[god];
	for(int e=0; e<ec/2 ;e++){
		if(f[e]==1){	
			int add=min(frog,dist[ev[e]]-pot[ev[e]])-min(frog,dist[eu[e]]-pot[eu[e]]);
			//cout << "adding+ " << eu[e] << ' ' << ev[e] << ' ' << add << endl;
			w[e]+=add;
		}
		if(f[e]==-1){
			int add=min(frog,dist[ev[e]]-pot[ev[e]])-min(frog,dist[eu[e]]-pot[eu[e]]);
			//cout << "adding- " << eu[e] << ' ' << ev[e] << ' ' << add << endl;
			add*=-1;
			w[e]+=add;
		}
	}
	if(dist[god]==k) return false;
	ans+=k-dist[god];
	while(god%m!=1){
		int ce=prv[god];
		god^=eu[ce/2]^ev[ce/2];
		if(ce%2) f[ce/2]--;
		else f[ce/2]++;
	}
	for(int i=1; i<=n*m ;i++) pot[i]+=min(frog,dist[i]-pot[i]);
	return true;
}
void solve(){
	cin >> n >> m >> k;
	ec=0;ans=0;
	for(int i=1; i<=n*m ;i++) adj[i].clear();
	//init
	for(int i=1; i<=n ;i++){
		for(int j=1; j<m ;j++){
			cin >> r[i][j];
			idr[i][j]=ec/2;
			add((i-1)*m+j,(i-1)*m+j+1,r[i][j]);
		}
	}
	for(int i=1; i<n ;i++){
		for(int j=1; j<=m ;j++){
			cin >> c[i][j];
			idc[i][j]=ec/2;
			add((i-1)*m+j,i*m+j,c[i][j]);
		}
	}
	iu();
	while(jieun());
	cout << ans << '\n';
	for(int i=1; i<=n ;i++){
		for(int j=1; j<m ;j++){
			cout << w[idr[i][j]] << ' ';
		}
		cout << '\n';
	}
	for(int i=1; i<n ;i++){
		for(int j=1; j<=m ;j++){
			cout << w[idc[i][j]] << ' ';
		}
		cout << '\n';
	}
}
int main(){
	ios::sync_with_stdio(false);cin.tie(0);
	int t;cin >> t;while(t--) solve();
}