#include <cassert>
#include <cmath>
#include <cstdint>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <bitset>
#include <complex>
#include <deque>
#include <functional>
#include <iostream>
#include <limits>
#include <map>
#include <numeric>
#include <queue>
#include <random>
#include <set>
#include <sstream>
#include <string>
#include <unordered_map>
#include <unordered_set>
#include <utility>
#include <vector>

using namespace std;

// using Int = long long;

template <class T1, class T2> ostream &operator<<(ostream &os, const pair<T1, T2> &a) { return os << "(" << a.first << ", " << a.second << ")"; };
template <class T> ostream &operator<<(ostream &os, const vector<T> &as) { const int sz = as.size(); os << "["; for (int i = 0; i < sz; ++i) { if (i >= 256) { os << ", ..."; break; } if (i > 0) { os << ", "; } os << as[i]; } return os << "]"; }
template <class T> void pv(T a, T b) { for (T i = a; i != b; ++i) cerr << *i << " "; cerr << endl; }
template <class T> bool chmin(T &t, const T &f) { if (t > f) { t = f; return true; } return false; }
template <class T> bool chmax(T &t, const T &f) { if (t < f) { t = f; return true; } return false; }
#define COLOR(s) ("\x1b[" s "m")


#ifndef LIBRA_OTHER_INT128_H_
#define LIBRA_OTHER_INT128_H_

#include <stdio.h>
#include <iostream>

constexpr unsigned __int128 toUInt128(const char *s) {
  unsigned __int128 x = 0;
  for (; *s; ++s) x = x * 10 + (*s - '0');
  return x;
}
constexpr __int128 toInt128(const char *s) {
  if (*s == '-') return -toInt128(s + 1);
  __int128 x = 0;
  for (; *s; ++s) x = x * 10 + (*s - '0');
  return x;
}
unsigned __int128 inUInt128() {
  static char buf[41];
  scanf("%s", buf);
  return toUInt128(buf);
}
__int128 inInt128() {
  static char buf[41];
  scanf("%s", buf);
  return toInt128(buf);
}

void out(unsigned __int128 x) {
  static char buf[41];
  int len = 0;
  do { buf[len++] = '0' + static_cast<int>(x % 10); } while (x /= 10);
  for (int i = len; --i >= 0; ) putchar(buf[i]);
}
void out(__int128 x) {
  if (x < 0) {
    putchar('-');
    out(-static_cast<unsigned __int128>(x));
  } else {
    out(static_cast<unsigned __int128>(x));
  }
}
std::ostream &operator<<(std::ostream &os, unsigned __int128 x) {
  static char buf[41];
  int len = 0;
  do { buf[len++] = '0' + static_cast<int>(x % 10); } while (x /= 10);
  for (int i = len; --i >= 0; ) os << buf[i];
  return os;
}
std::ostream &operator<<(std::ostream &os, __int128 x) {
  if (x < 0) {
    os << '-' << -static_cast<unsigned __int128>(x);
  } else {
    os << static_cast<unsigned __int128>(x);
  }
  return os;
}

#endif  // LIBRA_OTHER_INT128_H_


using INT = __int128;

int N;
vector<int> P;

INT dp[2][40][2010];

INT solve(int B) {
  if (B < 0) return false;
  vector<int> xs(N, 1), ys(N, 1);
  for (int i = 0; i < N; ++i) if (~P[i]) {
    xs[i] = ys[P[i]] = 0;
  }
  for (int b = 0; b <= B; ++b) {
    dp[0][0][b] = 0;
  }
  dp[0][0][0] = 1;
  int last = 0;
  int x = 0, y = 0;
  for (int i = 0; i < N; ++i) if (xs[i] || ys[i]) {
    const int z = min(x, y);
    for (int m = 0; m <= z; ++m) for (int b = 0; b <= B; ++b) {
      dp[1][m][b] = 0;
    }
    for (int m = 0; m <= z; ++m) for (int b = 0; b <= B; ++b) {
      const int bb = b + ((x-m) + (y-m)) * (i - last);
      if (bb <= B) {
        dp[1][m][bb] += dp[0][m][b];
      }
    }
// for(int m=0;m<=z;++m)for(int b=0;b<=B;++b)if(dp[1][m][b])cerr<<i<<" "<<m<<" "<<b<<": "<<dp[1][m][b]<<endl;
    for (int m = 0; m <= z + 2; ++m) for (int b = 0; b <= B; ++b) {
      dp[0][m][b] = 0;
    }
    if (xs[i]) {
      if (ys[i]) {
        for (int m = 0; m <= z; ++m) for (int b = 0; b <= B; ++b) if (dp[1][m][b]) {
          dp[0][m][b] += dp[1][m][b];
          dp[0][m + 1][b] += dp[1][m][b] * (x-m);
          dp[0][m + 1][b] += dp[1][m][b] * (y-m);
          dp[0][m + 2][b] += dp[1][m][b] * (x-m) * (y-m);
          dp[0][m + 1][b] += dp[1][m][b];
        }
      } else {
        for (int m = 0; m <= z; ++m) for (int b = 0; b <= B; ++b) if (dp[1][m][b]) {
          dp[0][m][b] += dp[1][m][b];
          dp[0][m + 1][b] += dp[1][m][b] * (y-m);
        }
      }
    } else {
      if (ys[i]) {
        for (int m = 0; m <= z; ++m) for (int b = 0; b <= B; ++b) if (dp[1][m][b]) {
          dp[0][m][b] += dp[1][m][b];
          dp[0][m + 1][b] += dp[1][m][b] * (x-m);
        }
      } else {
        assert(false);
      }
    }
    last = i;
    x += xs[i];
    y += ys[i];
  }
  assert(x == y);
// cerr<<COLOR("33")<<"[solve] "<<P<<" "<<B<<" = "<<dp[0][x][B]<<COLOR()<<endl;
  return dp[0][x][B];
}

int main() {
  int B;
  INT K;
  for (; ~scanf("%d%d", &N, &B); ) {
    K = inInt128();
    
    P.assign(N, -1);
solve(B);
    int b = B;
    INT k = K;
    for (int i = 0; i < N; ++i) {
      const set<int> used(P.begin(), P.begin() + i);
      for (int a = 0; a < N; ++a) if (!used.count(a)) {
        P[i] = a;
        const INT res = solve(b - abs(a - i));
        if (k > res) {
          k -= res;
        } else {
          goto found;
        }
      }
      assert(false);
     found:{}
      b -= abs(P[i] - i);
    }
    
    for (int i = 0; i < N; ++i) {
      if (i) printf(" ");
      printf("%d", P[i] + 1);
    }
    puts("");
  }
  return 0;
}