#include <bits/stdc++.h>
using namespace std;
const int mod = 1e9 + 7;
const int maxn = 1e6 + 100;
int a[maxn];
vector<int> num[5];
int t[5];
int main() {
  int n;
  scanf("%d", &n);
  for (int i = 0; i < n; ++i) {
    scanf("%d", &a[i]);
    if (a[i] & 1) {
      num[1].push_back(a[i]);
    } else
      num[a[i] % 4].push_back(a[i]);
  }
  if (num[2].size() > 0 && num[0].size() < num[1].size() ||
      num[1].size() > num[0].size() + 1) {
    printf("No\n");
  } else {
    printf("YES\n");
    if (t[1] < num[1].size()) {
      printf("%d %d ", num[1][t[1]++], num[0][t[0]++]);
    }
    while (t[2] < num[2].size()) {
      printf("%d ", num[2][t[2]++]);
    }
    while (t[0] < num[0].size() && t[1] < num[1].size()) {
      printf("%d %d ", num[0][t[0]++], num[1][t[1]++]);
    }
    while (t[0] < num[0].size() && t[3] < num[3].size()) {
      printf("%d %d ", num[3][t[3]++], num[0][t[0]++]);
    }
    if (t[1] < num[1].size()) {
      printf("%d ", num[1][t[1]++]);
    }
    while (t[0] < num[0].size()) {
      printf("%d ", num[0][t[0]++]);
    }
    printf("\n");
  }
  return 0;
}
