#include <bits/stdc++.h>
#pragma GCC optimize("Ofast")
using namespace std;
template <typename T>
bool uax(T &a, const T b) {
  if (a < b) {
    a = b;
    return true;
  }
  return false;
}
template <typename T>
bool uin(T &a, const T b) {
  if (a > b) {
    a = b;
    return true;
  }
  return false;
}
template <typename T1, typename T2>
istream &operator>>(istream &in, pair<T1, T2> &p) {
  in >> p.first >> p.second;
  return in;
}
template <typename T1, typename T2>
ostream &operator<<(ostream &out, pair<T1, T2> &p) {
  out << p.first << ' ' << p.second;
  return out;
}
int a, b, c, d;
int aa, bb, cc, dd;
int cnt = 0;
string st, fn;
void viv() {
  cout << st;
  for (int i = 0; i < aa; i++) {
    cout << "0 1 ";
  }
  for (int i = 0; i < cnt; i++) {
    cout << "2 1 ";
  }
  for (int i = 0; i < dd; i++) {
    cout << "2 3 ";
  }
  cout << fn;
}
int main() {
  ios_base::sync_with_stdio(false);
  cin.tie(0);
  st = fn = "";
  cin >> a >> b >> c >> d;
  if (a == 0 && b == 0) {
    if (d == c + 1) {
      cout << "YES"
           << "\n";
      cout << "3 ";
      for (int i = 0; i < c; i++) {
        cout << "2 3 ";
      }
      exit(0);
    }
  }
  if (c == 0 && d == 0) {
    if (a == b + 1) {
      cout << "YES"
           << "\n";
      cout << "0 ";
      for (int i = 0; i < b; i++) {
        cout << "1 0 ";
      }
      exit(0);
    }
  }
  aa = a;
  bb = b;
  cc = c;
  dd = d;
  b -= a;
  c -= d;
  if (b < 0 || c < 0) {
    cout << "NO";
    exit(0);
  }
  cnt = min(c, b);
  c -= cnt;
  b -= cnt;
  if (b <= 1 && c <= 1) {
    if (b == 1) st = "1 ";
    if (c == 1) fn = "2";
    cout << "YES"
         << "\n";
    viv();
    exit(0);
  }
  cout << "NO";
}
