/*
#include <boost/functional/hash.hpp>
#include <boost/multiprecision/cpp_int.hpp>
typedef boost::multiprecision::cpp_int mlint;
*/
#include <bits/stdc++.h>
using namespace std;

#define esc(x) cout << (x) << endl,exit(0)
#define inf(T) (numeric_limits<T>::max() / 2 - 1)
#define each(i,v) for(auto i = begin(v),end_copy = end(v); i != end_copy; ++i)
#define reach(i,v) for(auto i = rbegin(v),rend_copy = rend(v); i != rend_copy; ++i)
#define urep(i,s,t) for(int64_t i = (s),t_copy = (t); i <= t_copy; ++i)
#define drep(i,s,t) for(int64_t i = (s),t_copy = (t); i >= t_copy; --i)
#define rep(i,n) urep(i,0,(n) - 1)
#define rep1(i,n) urep(i,1,(n))
#define rrep(i,n) drep(i,(n) - 1,0)
#define all(v) begin(v),end(v)
#define rall(v) rbegin(v),rend(v)
#define vct vector
#define prque priority_queue
#define umap unordered_map
#define uset unordered_set
#define l_bnd lower_bound
#define u_bnd upper_bound
#define rsz resize
#define ers erase
#define emp emplace
#define emf emplace_front
#define emb emplace_back
#define pof pop_front
#define pob pop_back
#define mkp make_pair
#define mkt make_tuple
#define fir first
#define sec second

typedef int64_t i64;
typedef double db;
typedef vct<vct<i64>> mat;
typedef pair<int,int> pii;
typedef pair<i64,int> pli;
typedef pair<db,int> pdi;

#ifdef Local
#define debug 1
#define print(x) cout << #x << " = " << (x) << " (L" << __LINE__ << ")" << endl
// #include <sys/time.h>
// struct timeval start_time;
// clock_t start = clock();
#endif

void setup(short prec = 10) {
    cin.tie(0);
    ios::sync_with_stdio(false);
    cout << fixed << setprecision(prec);

    #ifdef Local
    // gettimeofday(&start_time,NULL);
    cout << "-- Execution At Local ---\n" << "\n<Standard Output>\n\n";
    auto print_atexit = []() {
        // struct timeval end_time;
        // gettimeofday(&end_time,NULL);
        // double runtime_millis = 1000.0 * (end_time.tv_sec - start_time.tv_sec) + 0.001 * (end_time.tv_usec - start_time.tv_usec);
        cout << "\n----------------\n";
        // cout << "\nExec time : " << runtime_millis << " ms\n\n";
        // cout << (clock() - start) * 1000.0 / CLOCKS_PER_SEC << endl;
    };
    atexit(print_atexit);
    #endif
}

template <class T, class U> ostream& operator << (ostream& s, pair<T,U> p) { return s << p.fir << " " << p.sec; }
template <class T> ostream& operator << (ostream& s, vector<T>& v) { for(auto i = v.begin(); i != v.end(); i++) { if(i != begin(v)) s << " "; s << *i; } return s; }

template <class T> void hash_combine(size_t &seed, T const &key) {
    hash<T> hasher;
    seed ^= hasher(key) + 0x9e3779b9 + (seed << 6) + (seed >> 2);
}

namespace std {
    template <class T, class U> struct hash<pair<T,U>> {
        size_t operator()(pair<T,U> const &p) const {
            size_t seed(0);
            hash_combine(seed,p.first);
            hash_combine(seed,p.second);
            return seed;
        }
    };
}

bool odd(i64 x) { return x & 1; }
bool even(i64 x) { return !odd(x); }
bool parity(i64 x, i64 y) { return odd(x) ^ even(y); }
bool bit(i64 n, uint8_t e) { return n >> e & 1; }
int ilog(i64 x, uint64_t b) { if(x) return ilog(x / b,b) + 1; return -1; }
i64 qceil(i64 x, i64 y) { return x > 0 ? (x - 1) / y + 1 :  x / y; }
uint64_t gcd(uint64_t a, uint64_t b) { if(!b) return abs(a); if(a % b) return gcd(b, a % b); return b; }
uint64_t lcm(uint64_t a, uint64_t b) { if(a & b) return a / gcd(a, b) * b; return 0; }
template <class T, class U> bool chmax(T& m, U x) { if(m < x) { m = x; return 1; } return 0; }
template <class T, class U> bool chmin(T& m, U x) { if(m > x) { m = x; return 1; } return 0; }
template <class T> T cmprs(T& v) {
    T tmp = v;
    sort(all(tmp));
    tmp.erase(unique(all(tmp)),end(tmp));
    each(i,v) *i = l_bnd(all(tmp),*i) - begin(tmp) + 1;
    return v;
}

const int dir[8][2] = { {1,0},{0,1},{-1,0},{0,-1},{1,1},{-1,1},{-1,-1},{1,-1} };
const i64 mod = 1e9 + 7;
const db eps = 1e-15;

int H,W,h,w;

int main() {
    setup();
	
	cin >> H >> W >> h >> w;
	if(H % h == 0 && W % w == 0) esc("No");
	cout << "Yes\n";
	int X = 1000000000 - 1;
	rep(i,H)rep(j,W) {
		if(i % h == 0 && j % w == 0) cout << X;
		else if(i % h == h - 1 && j % w == w - 1) cout << -1-X;
		else cout << 0;
		if(j < W - 1) cout << " ";
		else cout << endl;
	}
}


Measured memory_bytes: 2944