#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int MAXN=1005,MAXV=2e6+5;
int n,m,u,v,a[MAXN][MAXN];
int x[MAXN],y[MAXN];
bitset <MAXV> f[MAXN];
bool solve() { //x_max > y_max
	int p=abs(n-m),z=0;
	for(int j=1;j<=m;++j) {
		y[j]=a[u][v]-a[u][j];
		z-=y[j];
	}
	vector <array<int,2>> ch;
	for(int i=1;i<=n;++i) {
		int lx=-a[i][v],rx=a[i][v];
		bool fl=true,fr=true;
		for(int j=1;j<=m;++j) {
			int px=y[j]+a[i][j],qx=y[j]-a[i][j];
			fl&=(px==lx||qx==lx);
			fr&=(px==rx||qx==rx);
		}
		if(!fl&&!fr) return false;
		else if(!fr) x[i]=lx;
		else if(!fl) x[i]=rx;
		else x[i]=lx,ch.push_back({i,rx-lx});
		z+=x[i];
	}
	int q=ch.size();
	f[0].reset(),f[0].set(0);
	for(int i=0;i<q;++i) {
		f[i+1]=f[i],f[i+1]<<=ch[i][1],f[i+1]|=f[i];
	}
	int S=-1;
	if(!p) S=-z;
	else {
		for(int r=(p-z%p)%p;r<MAXV;r+=p) if(f[q][r]) {
			S=r; break;
		}
	}
	if(S<0) return false;
	for(int i=q-1;~i;--i) if(!f[i][S]) {
		S-=ch[i][1],x[ch[i][0]]+=ch[i][1];
	}
	int sx=accumulate(x+1,x+n+1,0),sy=accumulate(y+1,y+m+1,0);
	if(!p) assert(sx==sy);
	else {
		assert((sx-sy)%p==0);
		int k=(sy-sx)/(n-m);
		for(int i=1;i<=n;++i) x[i]+=k;
		for(int j=1;j<=m;++j) y[j]+=k;
	}
	memset(a,0,sizeof(a));
	for(int i=2;i<=n;++i) a[i][1]=x[i];
	for(int j=2;j<=m;++j) a[1][j]=y[j];
	a[1][1]=y[1]-accumulate(x+2,x+n+1,0);
	return true;
}
signed main() {
	ios::sync_with_stdio(false);
	cin>>n>>m,u=v=1;
	for(int i=1;i<=n;++i) for(int j=1;j<=m;++j) {
		cin>>a[i][j];
		if(a[u][v]<a[i][j]) u=i,v=j;
	}
	if(solve()) {
		for(int i=1;i<=n;++i) for(int j=1;j<=m;++j) cout<<a[i][j]<<" \n"[j==m];
		return 0;
	}
	for(int i=1;i<=n;++i) for(int j=i+1;j<=m;++j) swap(a[i][j],a[j][i]);
	swap(n,m),swap(u,v);
	solve();
	for(int i=1;i<=m;++i) for(int j=1;j<=n;++j) cout<<a[j][i]<<" \n"[j==n];
	return 0;
}