POJ1037 A decorative fence 【动态规划】

Richard just finished building his new house. Now the only thing the house misses is a cute little wooden fence. He had no idea how to make a wooden fence, so he decided to order one. Somehow he got his hands on the ACME Fence Catalogue 2002, the ultimate resource on cute little wooden fences. After reading its preface he already knew, what makes a little wooden fence cute. 
A wooden fence consists of N wooden planks, placed vertically in a row next to each other. A fence looks cute if and only if the following conditions are met: 
?The planks have different lengths, namely 1, 2, . . . , N plank length units. 
?Each plank with two neighbors is either larger than each of its neighbors or smaller than each of them. (Note that this makes the top of the fence alternately rise and fall.) 
It follows, that we may uniquely describe each cute fence with N planks as a permutation a1, . . . , aN of the numbers 1, . . . ,N such that (any i; 1 < i < N) (ai ? ai?1)*(ai ? ai+1) > 0 and vice versa, each such permutation describes a cute fence. 
It is obvious, that there are many di erent cute wooden fences made of N planks. To bring some order into their catalogue, the sales manager of ACME decided to order them in the following way: Fence A (represented by the permutation a1, . . . , aN) is in the catalogue before fence B (represented by b1, . . . , bN) if and only if there exists such i, that (any j < i) aj = bj and (ai < bi). (Also to decide, which of the two fences is earlier in the catalogue, take their corresponding permutations, find the first place on which they differ and compare the values on this place.) All the cute fences with N planks are numbered (starting from 1) in the order they appear in the catalogue. This number is called their catalogue number. 

After carefully examining all the cute little wooden fences, Richard decided to order some of them. For each of them he noted the number of its planks and its catalogue number. Later, as he met his friends, he wanted to show them the fences he ordered, but he lost the catalogue somewhere. The only thing he has got are his notes. Please help him find out, how will his fences look like.

The first line of the input file contains the number K (1 <= K <= 100) of input data sets. K lines follow, each of them describes one input data set. 
Each of the following K lines contains two integers N and C (1 <= N <= 20), separated by a space. N is the number of planks in the fence, C is the catalogue number of the fence. 
You may assume, that the total number of cute little wooden fences with 20 planks fits into a 64-bit signed integer variable (long long in C/C++, int64 in FreePascal). You may also assume that the input is correct, in particular that C is at least 1 and it doesn抰 exceed the number of cute fences with N planks.

For each input data set output one line, describing the C-th fence with N planks in the catalogue. More precisely, if the fence is described by the permutation a1, . . . , aN, then the corresponding line of the output file should contain the numbers ai (in the correct order), separated by single spaces.

2
2 1
3 3

1 2
2 3 1

#include <stdio.h>
#include <string.h>
#define maxn 22
#define UP 0
#define DOWN 1

int vis[maxn], ans[maxn];
__int64 dp[maxn][maxn][2];

void Init(int n)
{
	int i, j, k;
	dp[1][1][UP] = dp[1][1][DOWN] = 1;
	for(i = 2; i <= n; ++i){ //i is the amount of total sticks
		for(j = 1; j <= i; ++j){ //j is the first stick's position
			for(k = j; k < i; ++k) //k is the subSolution of dp array
				dp[i][j][UP] += dp[i-1][k][DOWN];
			for(k = 1; k < j; ++k)
				dp[i][j][DOWN] += dp[i-1][k][UP];
		}
	}
}

void Print(int n, __int64 c)
{
	memset(vis, 0, sizeof(vis));
	int i, j, k, rank;
	__int64 skip = 0, pre = 0;
	for(i = 1; i <= n; ++i){ //i is the position to select now
		rank = 0; //the rank's min stick
		for(j = 1; j <= n; ++j){ //j is the stick to be select		
			pre = skip;	
			if(!vis[j]){
				++rank;
				if(i == 1) skip += dp[n][rank][UP] + dp[n][rank][DOWN];
				else if(j > ans[i-1] && (i == 2 || ans[i-2] > ans[i-1]))
					skip += dp[n-i+1][rank][DOWN];
				else if(j < ans[i-1] && (i == 2 || ans[i-2] < ans[i-1]))
					skip += dp[n-i+1][rank][UP];
				if(skip >= c) break; 
			}
		}
		ans[i] = j;
		vis[j] = 1;
		skip = pre;
	}
	for(i = 1; i <= n; ++i)
		if(i != n) printf("%d ", ans[i]);
		else printf("%d\n", ans[i]);
}

int main()
{
	int t, n;
	__int64 c;
	Init(20);
	scanf("%d", &t);
	while(t--){
		scanf("%d%I64d", &n, &c);
		Print(n, c);
	}
	return 0;
}

#include <stdio.h>
#include <string.h>
#define UP 0
#define DOWN 1
#define maxn 25

long long C[maxn][maxn][2];
//C[i][k][DOWN] 是S(i)中以第k短的木棒打头的DOWN方案数,C[i][k][UP] 
//是S(i)中以第k短的木棒打头的UP方案数,第k短指i根中第k短

void Init(int n)
{
    memset(C, 0, sizeof(C));
    C[1][1][UP] = C[1][1][DOWN] = 1;
    int i, k, M, N;
    for(i = 2; i <= n; ++i){
        for(k = 1; k <= i; ++k){ //枚举第一根木棒的长度
            //枚举第二根长度
            for(M = k; M < i; ++M) C[i][k][UP] += C[i-1][M][DOWN];
            //枚举第二根长度
            for(N = 1; N < k; ++N) C[i][k][DOWN] += C[i-1][N][UP];
        }
    }
    //总方案数是 Sum{ C[n][k][DOWN] + C[n][k][UP] } k = 1.. n;
}

void Print(int n, long long cc)
{
    long long skipped = 0, oldVal; //已经跳过的方案数
    int seq[maxn]; //最终要输出的答案
    int used[maxn]; //木棒是否用过
    memset(used, 0, sizeof(used));
    for(int i = 1, k; i <= n; ++i){ //依次确定位置i的木棒序号
        int No = 0;
        for(k = 1; k <= n; ++k){ //枚举位置i的木棒
            oldVal = skipped;
            if(!used[k]){
                ++No; //k是剩下木棒里的第No短的
                if(i == 1) skipped += C[n][No][UP] + C[n][No][DOWN];
                //以下寻找合法位置
                else if(k > seq[i-1] && (i == 2 || seq[i-2] > seq[i-1]))
                    skipped += C[n-i+1][No][DOWN];
                else if(k < seq[i-1] && (i == 2 || seq[i-2] < seq[i-1]))
                    skipped += C[n-i+1][No][UP];
                if(skipped >= cc) break;
            }
        }
        used[k] = 1;
        seq[i] = k;
        skipped = oldVal;
    }
    for(int i = 1; i <= n; ++i)
        if(i < n) printf("%d ", seq[i]);
        else printf("%d\n", seq[i]);
}

int main()
{
    int T, n; long long c;
    Init(20);
    scanf("%d", &T);
    while(T--){
        scanf("%d%lld", &n, &c);
        Print(n, c);
    }
    return 0;
}