#include <bits/stdc++.h>
using namespace std;
using ll = long long;
int n,m,k,p;
ll qpow(ll a,ll b){
    ll ans=1;
    while(b){
        if(b&1) ans=ans*a%p;
        a=a*a%p;
        b>>=1;
    }
    return ans;
}
ll mul(ll x,ll y){
    return x*qpow(y,p-2)%p;
}
vector<ll> operator-(vector<ll> a,vector<ll> b){
    int sz=a.size();
    for(int i=0;i<sz;i++)   a[i]=((a[i]-b[i])%p+p)%p;
    return a;
}   
vector<ll> operator+(vector<ll> a,vector<ll> b){
    int sz=a.size();
    for(int i=0;i<sz;i++)   a[i]=((a[i]+b[i]%p)%p+p)%p;
    return a;
}   
vector<ll> operator*(vector<ll> a,ll t){
    int sz=a.size();
    for(int i=0;i<sz;i++)   a[i]=t*a[i]%p;
    return a;
}  
bool iszero(vector<ll> a){
    for(auto i:a)   if(i)   return false;
    return true;
}
void solve(void) {
    cin>>n>>m>>k>>p;
    vector<vector<ll>> mp(n,vector<ll>(m,0));
    for(int i=0;i<n;i++){
        for(int j=0;j<m;j++)    cin>>mp[i][j];
    }
    bool sp=0;
    if(n<m){
        sp=1;
        swap(n,m);
        vector<vector<ll>> np(n,vector<ll>(m,0));
        for(int i=0;i<n;i++){
            for(int j=0;j<m;j++)    np[i][j]=mp[j][i];
        }
        mp=np;
    }
    vector<vector<ll>> checker=mp;
    
    int ji_num=0;
    vector<vector<ll>> bt(m);
    for(int i=0;i<n;i++){
        vector<ll> jp=mp[i];
        for(int j=0;j<m;j++){
            if(jp[j]){
                if(bt[j].empty())    {bt[j]=jp;ji_num++;break;}
                else{
                    jp=jp-bt[j]*(jp[j]*qpow(bt[j][j],p-2)%p);
                }
            }
        }
    }
    // cerr<<"???"<<ji_num<<'\n';
    // cerr<<"!!!"<<ji_num<<'\n';
    // for(int i=0;i<ji_num;i++){
    //     for(auto j:bt[i])   cerr<<j<<' ';
    //     cerr<<'\n';
    // }
    // for(int i=0;i<n;i++){
    //     for(int j=0;j<m;j++)    cerr<<mp[i][j]<<' ';
    //     cerr<<'\n';
    // }




    if(ji_num>k)    {cout<<"No solution!\n";return ;}
    
    vector<vector<ll>> a(n,vector<ll>(k));
    vector<vector<ll>> b(k,vector<ll>(m));
    if(ji_num==0){
        if(k==1){
            cout<<"No solution!\n";
            return ;
        } 
        if(m==1&&n==1){
            if(p==2){
                if(k%2==0){
                    a.assign(n,vector<ll>(k,1));
                    b.assign(k,vector<ll>(m,1));
                } 
                else{
                    cout<<"No solution!\n";
                    return ;
                }
            }
            else {
                a.assign(n,vector<ll>(k,1));
                b.assign(k,vector<ll>(m,1));
                ll dw=(p-1ll*(k-1)*(k-1)%p)%p;
                if(dw==0){
                    a[0][k-2]=2;
                    dw=p-1;
                }
                a[0][k-1]=dw;
            }
        }
        else{
            b.assign(k,vector<ll>(m,1));
            for(int i=0;i<k-1;i++){
                a[0][i]=1;
            }
            a[0][k-1]=(p-(k-1)%p)%p;
            a[1][0]=a[0][k-1];
            for(int j=1;j<k;j++)    a[1][j]=a[0][j-1];
            for(int j=2;j<n;j++)    a[j]=a[0];
        }
    }
    else{
        if(m==1&&n==1&&p==2){
            if(k%2==1){
                a.assign(n,vector<ll>(k,1));
                b.assign(k,vector<ll>(m,1));
            } 
            else{
                cout<<"No solution!\n";
                return ;
            }
        }
        else{
            int tot=0;
            for(int i=0;i<m;i++){
                if(!bt[i].empty()){
                    for(int j=0;j<m;j++){
                        b[tot][j]=bt[i][j];
                    }
                    tot++;
                }
            }

            vector<int> ko;
            for(int i=0;i<n;i++){
                bool fm=false;
                for(int j=0;j<m;j++){
                    fm=(fm|(mp[i][j]!=0));
                }
                if(!fm) ko.push_back(i);  
            }
            if(ko.size()&&ji_num!=k){
                b[tot]=b[0];
                for(auto j:ko){
                    a[j][0]=1;a[j][tot]=(p-1);
                }
                tot++;
            }
            auto get=[&](vector<ll> x)->vector<ll> {
                vector<ll> res;
                for(int i=0;i<m;i++){
                    if(bt[i].size()){
                        res.push_back(x[i]*qpow(bt[i][i],p-2)%p);
                        x=x-bt[i]*res.back();
                    }
                }
                return res;
            };
            for(int j=tot;j<k;j++){
                b[j]=b[0];
                a[0][j]=1;
                mp[0]=mp[0]-b[0];
            }
            vector<ll> happy=get(mp[0]);
            if(ko.empty()&&iszero(happy)){
                for(int j=tot;j<k;j++){
                    a[0][j]=0;
                    a[1][j]=1;
                    mp[0]=mp[0]+b[0];
                    mp[1]=mp[1]-b[0];
                }
            }



            for(int i=0;i<n;i++){
                bool fm=false;
                for(int j=0;j<m;j++){
                    fm=(fm|(mp[i][j]!=0));
                }
                if(!fm) continue;
                // cerr<<"?"<<i<<'\n';
                vector<ll> gt=get(mp[i]); 
                for(int j=0;j<gt.size();j++)    a[i][j]=gt[j];

            }
        }
    }

    
    if(sp){
        for(int i=0;i<k;i++){
            for(int j=0;j<m;j++)    cout<<b[i][j]<<' ';
            for(int j=0;j<n;j++)    cout<<a[j][i]<<' ';
            cout<<'\n';
        }
    }
    else{
        for(int i=0;i<k;i++){
            for(int j=0;j<n;j++)    cout<<a[j][i]<<' ';
            for(int j=0;j<m;j++)    cout<<b[i][j]<<' ';
            cout<<'\n';
        } 
    }


    // for(int i=0;i<n;i++){
    //     for(int j=0;j<k;j++)    cerr<<a[i][j]<<' ';
    //     cerr<<'\n';
    // }
    // for(int i=0;i<k;i++){
    //     for(int j=0;j<m;j++)    cerr<<b[i][j]<<' ';
    //     cerr<<'\n';
    // }
    // cerr<<"\n\n";
    // for(int i=0;i<n;i++){
    //     for(int j=0;j<m;j++){
    //         for(int z=0;z<k;z++){
    //             checker[i][j]=(checker[i][j]-a[i][z]*b[z][j]%p)%p;
    //         }
    //         if(checker[i][j])   assert(0);
    //     }
    // }
}

int main() {
    ios::sync_with_stdio(false), cin.tie(nullptr);

    int _ = 1;
    while (_--) solve();

    return 0;
}
/*
3 3 2 7
1 2 1
2 4 2
1 1 1
*/