poj 1037 A decorative fence

题目链接：http://poj.org/problem?id=1037

这道题的解析见黑书257页，后面的答案搜索花了我不少时间，后来想到了从后向前还原的方法，还是相当不错的。

#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<string>
#include<queue>
#include<algorithm>
#include<vector>
#include<stack>
#include<list>
#include<iostream>
#include<map>
using namespace std;
#define inf 0x3f3f3f3f
#define Max 110
int max(int a,int b)
{
	return a>b?a:b;
}
int min(int a,int b)
{
	return a<b?a:b;
}
int t,n;
__int64 dp[22][22][2],c;
__int64 rec[22],ans;
void init()
{
    int i,j,k;
    memset(dp,0,sizeof(dp));
    dp[1][1][0]=1;
    dp[1][1][1]=1;
    for(i=2;i<=20;i++)
        for(j=1;j<=i;j++)
        {
            for(k=j;k<=i-1;k++)
                dp[i][j][0]+=dp[i-1][k][1];
            for(k=1;k<j;k++)
                dp[i][j][1]+=dp[i-1][k][0];
          //  printf("i %d j %d dp %I64d %I64d\n",i,j,dp[i][j][0],dp[i][j][1]);
        }
}
void dfs(int pos,int last,int s)
{
    int i;
    if(pos==0)
        return;
    if(pos==n)
    {
        for(i=1;i<=n;i++)
        {
            ans+=dp[pos][i][1];
            if(ans>=c)
            {
                ans-=dp[pos][i][1];
                rec[n-pos]=i;
                dfs(pos-1,i,1);
                return;
            }
            ans+=dp[pos][i][0];
            if(ans>=c)
            {
                ans-=dp[pos][i][0];
                rec[n-pos]=i;
                dfs(pos-1,i,0);
                return;
            }
        }
    }
    if(s==0)
    {
        for(i=last;i<=pos;i++)
        {
            ans+=dp[pos][i][1];
            if(ans>=c)
            {
                ans-=dp[pos][i][1];
                rec[n-pos]=i;
                dfs(pos-1,i,1);
                return ;
            }
        }
    }
    if(s==1)
    {
        for(i=1;i<last;i++)
        {
            ans+=dp[pos][i][0];
            if(ans>=c)
            {
                ans-=dp[pos][i][0];
                rec[n-pos]=i;
                dfs(pos-1,i,0);
                return;
            }
        }
    }
}
int main()
{
    int i,j,pos,last,s;
    scanf("%d",&t);
    init();
    while(t--)
    {
        scanf("%d%I64d",&n,&c);
        ans=0;
        dfs(n,0,0);
      // s=0;
        for(i=n-2;i>=0;i--)
            for(j=i+1;j<n;j++)
            {
                if(rec[j]>=rec[i])
                    rec[j]++;
            }
        for(i=0;i<n-1;i++)
            printf("%I64d ",rec[i]);
        printf("%I64d\n",rec[n-1]);
    }
}