这道题与下一章的数位(dp)解题思路十分一致。
把寻找答案变成按位(并且是字典序从小到大)枚举当前这一位可以填的情况。
通过(dp)预处理的信息告诉我们可行性,就可以把答案紧逼到一个更小的(子)问题,非常有趣。
考虑 (dp) 预处理的信息:
(f[i][j][0 / 1]) 表示 (i) 块木板,最左边的长度是第 (j) 小(排名为 (j) ),最左边这块是低位 / 高位的方案数。
由于木板的数量进行了变化,在加入一块新的木板后,木板的值域从 ([1, i - 1]) 变成了 ([1, i]),所以我们可以考虑把木板的长度变为一个相对的数量,从而进行转化。
考虑第一种转移
(f[i][j][0] = sum_{k = j}^{i - 1}f[i - 1][k][1])
可以看做是把后面 (i - 1) 块木板中真实长度 (>= j) 的再抬高一格,这样再拼一个 (j) 的木板就是一个合适的状态。
第二种转移也类似:
(f[i][j][1] = sum_{k = 1}^{j - 1} f[i - 1][k][0])
#include <cstdio>
#include <iostream>
#include <cstring>
using namespace std;
const int N = 21;
typedef long long LL;
int n;
bool st[N];
LL C, f[N][N][2];
void init() {
f[1][1][0] = f[1][1][1] = 1;
for (int i = 2; i < N; i++) {
for (int j = 1; j <= i; j++) {
for (int k = j; k <= i - 1; k++) f[i][j][0] += f[i - 1][k][1];
for (int k = 1; k <= j - 1; k++) f[i][j][1] += f[i - 1][k][0];
}
}
}
int main() {
init();
int T; scanf("%d", &T);
while (T--) {
memset(st, false, sizeof st);
scanf("%d%lld", &n, &C);
int k, last;
for (int i = 1; i <= n; i++) {
if (f[n][i][1] >= C) { k = 1, st[last = i] = true; printf("%d ", i); break; }
else C -= f[n][i][1];
if (f[n][i][0] >= C) { k = 0, st[last = i] = true; printf("%d ", i); break; }
else C -= f[n][i][0];
}
for (int i = n - 1; i; i--) {
k = k ^ 1;
int d = 0;
for (int j = 1; j <= n; j++) {
if (st[j]) continue;
++d;
if ((k == 0 && j < last) || (k == 1 && j > last)) {
// 要满足 j < last
if (f[i][d][k] >= C) { st[last = j] = true; printf("%d ", j); break; }
else C -= f[i][d][k];
}
}
}
puts("");
}
return 0;
}