#pragma GCC optimize(2)
#pragma GCC optimize("Ofast")
#pragma GCC optimize("inline","fast-math","unroll-loops","no-stack-protector")
#pragma GCC diagnostic error "-fwhole-program"
#pragma GCC diagnostic error "-fcse-skip-blocks"
#pragma GCC diagnostic error "-funsafe-loop-optimizations"
#include "perm.h"
// MagicDark
#include <bits/stdc++.h>
#define debug cerr << "\033[32m[" << __LINE__ << "]\033[0m "
#define SZ(x) ((int) x.size() - 1)
#define all(x) x.begin(), x.end()
#define ms(x, y) memset(x, y, sizeof x)
#define F(i, x, y) for (int i = (x); i <= (y); i++)
#define DF(i, x, y) for (int i = (x); i >= (y); i--)
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
template <typename T> T& chkmax(T& x, T y) {return x = max(x, y);}
template <typename T> T& chkmin(T& x, T y) {return x = min(x, y);}
template <typename T> T& read(T &x) {
	x = 0; int f = 1; char c = getchar();
	for (; !isdigit(c); c = getchar()) if (c == '-') f = -f;
	for (; isdigit(c); c = getchar()) x = (x << 1) + (x << 3) + (c ^ 48);
	return x *= f;
}
mt19937 mrand(chrono::steady_clock::now().time_since_epoch().count());
inline int rnd(int l, int r) {return mrand() % (r - l + 1) + l;}
vector <int> solve(ll n) {
	vector <int> ans;
	DF(i, 60, 1) {
		ll w = (1ll << i) - 1;
		while (n >= w) {
			for (int& j: ans) j += i;
			F(j, 1, i) ans.push_back(j - 1);
			n -= w;
		}
	}
	return ans;
}
int dp[1010], t[1010], len;
void modify(int x, int y) {
	x++;
	for (; x <= len; x += x & -x) t[x] += y;
}
int query(int x) {
	int s = 0;
	for (; x; x ^= x & -x) s += t[x];
	return s;
}
int calc(vector <int> x) {
	len = x.size();
	F(i, 1, len) t[i] = 0;
	int s = 0;
	F(i, 0, SZ(x)) {
		dp[i] = 1 + query(x[i]);
		// F(j, 0, i - 1) {
		// 	if (x[j] < x[i]) {
		// 		dp[i] += dp[j];
		// 	}
		// }
		s += dp[i];
		modify(x[i], dp[i]);
	}
	return s;
}
std::vector<int> construct_permutation(long long n) {
	n--;
	vector <int> ans;
	// if (n <= 90) {
	// 	F(i, 1, n) ans.push_back(n - i);
	// 	return ans;
	// }
	ans = solve(n);
	int lim = 1;
	F(i, 1, 30) {
		vector <int> t(i);
		F(j, 1, i) t[j - 1] = j - 1;
		lim = min(100, lim * i);
		// do {
		F(j, 1, lim) {
			shuffle(all(t), mrand);
			int v = calc(t);
			// debug << v << endl;
			if (n >= v) {
				vector <int> w = solve(n - v);
				if (w.size() + i < ans.size()) {
					ans.clear();
					for (int j: t) ans.push_back(j + w.size());
					for (int j: w) ans.push_back(j);
				}
			}
		}
		// } while (next_permutation(all(t)));
	}
	return ans;
}