#ifndef DEBUG
#  define NDEBUG
#endif
#include <iostream>
#include <algorithm>
#include <string>
#include <iterator>
#include <utility>
#include <type_traits>
#include <limits>
#include <cassert>
#include <numeric>
#include <vector>
#include <array>

using namespace std;

// template {{{
// misc {{{
using Z = long long;
[[maybe_unused]]
constexpr Z INF = numeric_limits<Z>::max() / 3;
// }}}

// debug {{{
template <class T>
using enable_if_sequence_t = enable_if_t<
  is_same<
    typename iterator_traits<
      typename remove_cv_t<remove_reference_t<T>>::iterator
    >::iterator_category,
    random_access_iterator_tag>{}>;
#ifdef NDEBUG
# define dump(x)
# define debug if (0)
template <class T>
inline T&& ix(T&& value) { return value; }
template <class T, class Index, class... Indices, class = enable_if_sequence_t<T>>
inline auto&& ix(T&& container, Index index, Indices... indices) {
  return ix(container[index], indices...);
}
#else
# define dump(x) do { \
    cerr << "\e[33m[dump: " << __LINE__ << "] " << (#x) << " = " << x << "\e[0m" << endl; \
  } while (0);
# define debug
template <class T>
inline auto&& ix_impl([[maybe_unused]] int line, [[maybe_unused]] const char* code, T&& value) {
  return value;
}
template <class T, class Index, class... Indices, class = enable_if_sequence_t<T>>
inline auto&& ix_impl(int line, const char* code, T&& container, Index index, Indices... indices) {
  if (index < 0) {
    cerr << "\e[31mix: negative index at line " << line << ", ix(" << code << ")\n"
      << "index = " << index << "\e[0m" << endl;
    exit(1);
  }
  int size = distance(begin(container), end(container));
  if (index >= size) {
    cerr << "\e[31mix: index out of range at line " << line << ", ix(" << code << ")\n"
      << "index = " << index << " whereas size = " << size << "\e[0m" << endl;
    exit(1);
  }
  return ix_impl(line, code, container[index], indices...);
}
#define ix(...) ix_impl(__LINE__, #__VA_ARGS__, __VA_ARGS__)
#endif
// }}}

// IO {{{
template <class T, class U>
ostream& operator<< (ostream& os, pair<T, U> p) {
  os << '(' << p.first << ", " << p.second << ')';
  return os;
}
template <class T, class U>
istream& operator>> (istream& is, pair<T, U>& p) {
  is >> p.first >> p.second;
  return is;
}
template <class T>
false_type is_container_impl(...);
template <class T, class = typename T::iterator>
true_type is_container_impl(T);
template <class T>
struct is_container : decltype(is_container_impl<T>(declval<T>())) {};
template <>
struct is_container <string> : std::false_type {};
template <class Container, class = enable_if_t<is_container<Container>{}>>
istream& operator>> (istream& is, Container& cont) {
  for (auto& v : cont) is >> v;
  return is;
}
template <class Container, class = enable_if_t<is_container<Container>{}>>
ostream& operator<< (ostream& os, Container const& cont) {
  bool first = true;
  constexpr char sep =
    is_container<typename Container::value_type>{} ? '\n' : ' ';
  for (auto const& v : cont) {
    if (!first) {
      os << sep;
    }
    first = false;
    os << v;
  }
  return os;
}
template <class T>
void conclude(T const& message) {
  cout << message << endl;
  exit(0);
}
// }}}
// }}}


template <class T, class F>
T boundary(T ok, T ng, F is_ok) {
  while (abs(ok - ng) > 1) {
    T m = (ok + ng) / 2;
    (is_ok(m) ? ok : ng) = m;
  }
  return ok;
}

int main() {
  ios_base::sync_with_stdio(false);
  cin.tie(nullptr);
  cout.precision(10);
  cout << fixed;
  Z H, W, h, w;
  cin >> H >> W >> h >> w;
  if (H % h == 0 && W % w == 0) {
    conclude("No");
  }
  cout << "Yes" << '\n';
  if (h == 1) {
    Z pos = boundary<Z>(1, INF / 1000, [&] (Z k) {
      return k * (w - 1) < 1000000000;
    });
    vector<Z> row(W);
    for (auto i = 0; i < W; i++) {
      ix(row, i) = (i + 1) % w == 0 ? -1000000000 : pos;
    }
    for (auto i = 0; i < H; i++) {
      cout << row << '\n';
    }
  }
  else if (w == 1) {
    Z pos = boundary<Z>(1, INF / 1000, [&] (Z k) {
      return k * (h - 1) < 1000000000;
    });
    for (auto i = 0; i < H; i++) {
      if ((i + 1) % h == 0) {
        cout << vector<Z>(W, -1000000000) << '\n';
      }
      else {
        cout << vector<Z>(W, pos) << '\n';
      }
    }
  }
  else {
    Z pos = boundary<Z>(1, INF / 100000000ll, [&] (Z k) {
      Z negs = 10000000ll * (h + w - 1);
      return k * (h - 1) * (w - 1) < negs;
    });
    vector<Z> usual(W), line(W, -10000000);
    for (auto i = 0; i < W; i++) {
      ix(usual, i) = (i + 1) % w == 0 ? -10000000ll : pos;
    }
    for (auto i = 0; i < H; i++) {
      cout << ((i + 1) % h == 0 ? line : usual) << '\n';
    }
  }
}



Measured memory_bytes: 2944