#include <bits/stdc++.h>
using namespace std;
long long mod = 1e9 + 7;
long long max(long long a, long long b) {
  if (a > b)
    return a;
  else
    return b;
}
long long min(long long a, long long b) {
  if (a < b)
    return a;
  else
    return b;
}
long long gcd(long long a, long long b, long long& x, long long& y) {
  if (b == 0) {
    x = 1;
    y = 0;
    return a;
  }
  long long x1, y1;
  long long d = gcd(b, a % b, x1, y1);
  x = y1;
  y = x1 - y1 * (a / b);
  return d;
}
long long inv(long long a) {
  long long x, y;
  long long g = gcd(a, mod, x, y);
  x = (x % mod + mod) % mod;
  return x;
}
long long binpow(long long a, long long b, long long m) {
  a %= m;
  long long res = 1;
  while (b > 0) {
    if (b & 1) res = res * a % m;
    a = a * a % m;
    b >>= 1;
  }
  return res;
}
int main() {
  int t;
  t = 1;
  while (t--) {
    string s;
    cin >> s;
    long long l = s.length();
    vector<pair<int, char>> v;
    for (int i = 0; i < l; i++) {
      int f = 0;
      if (v.size() > 0) {
        for (int j = 0; j < v.size(); j++) {
          if (v[j].second == s[i]) {
            v[j].first++;
            f = 1;
            break;
          }
        }
      }
      if (f == 0) v.push_back(make_pair(1, s[i]));
    }
    sort(v.begin(), v.end());
    vector<char> is_prime(l + 1, true);
    is_prime[0] = is_prime[1] = false;
    for (int i = 2; i <= l; i++) {
      if (is_prime[i] && (long long)i * i <= l) {
        for (int j = i * i; j <= l; j += i) is_prime[j] = false;
      }
    }
    int ct = 0;
    int fill[l + 1];
    memset(fill, 0, sizeof(fill));
    for (int i = 2; i <= l / 2; i++) {
      ct++;
    }
    for (int i = l / 2 + 1; i <= l; i++) {
      if (!is_prime[i]) ct++;
    }
    if (ct > v.back().first) {
      cout << "NO";
    } else {
      for (int i = 2; i <= l / 2; i++) {
        s[i - 1] = v.back().second;
        v.back().first--;
        fill[i] = 1;
      }
      for (int i = l / 2 + 1; i <= l; i++) {
        if (!is_prime[i]) {
          s[i - 1] = v.back().second;
          v.back().first--;
          fill[i] = 1;
        }
      }
      int k = 0;
      if (v[k].first == 0) k++;
      for (int i = 1; i <= l; i++) {
        if (!fill[i]) {
          s[i - 1] = v[k].second;
          v[k].first--;
          if (v[k].first == 0) k++;
        }
      }
      cout << "YES \n";
      cout << s;
    }
  }
  return 0;
}
