状态压缩dp(hdu2662)(我综合了一个人的解释和另一个人的代码）

hoj 2662


大概题意是：有一个n*m的棋盘，在这个棋盘里边放k个旗子，要求每一行每一列都不能存在一对旗子相邻，问最后总共的方案数。


这道题一看状态非常多，就一定是状压。怎么状压呢？这又是个问题。
慢慢考虑下这道题每个局面分别有哪些状态，我们很容易就能想到有以下几个状态：
第一：每一行放了多少个旗子；
第二：已经用了多少个旗子;
第三：已经放的这些旗子能不能保证合法，即上下左右均不相邻。
在这里可能还不容易想到那个状态表达式，那么我们先来考虑个简单的，假如说只有一行，要求在这一行里边填充k个旗子，要求任意两个都不相邻，这个时候的dp应该怎么表示?这就很简单了，直接就是dp[i][j][x]，代表已经到了第i列，已经使用了j个旗子，而且当前第i列的状态就是x（当然这里x只能是0和1,这里0代表这个第i列没有放旗子，1就代表这个位置放了旗子）的总方案数，递推关系是怎么写？其实也很简单，
dp[i][j][0]=dp[i-1][j][0]+dp[i-1][j][1];
dp[i][j][1]=dp[i-1][j-1][0];//这里只能是dp[i-1][j-1][0]，因为第i列已经放了，那么第i-1列就一定不能放。
当然这里你需要考虑到二维的局面，怎么考虑，把行对应于列，每一列的状态转化为每一行的状态，前i列使用了j个旗子变成前i行使用了j个旗子就这样思考。
综上考虑，我们会想到要有一个这样的dp，就是dp[i][j][x]，这里代表的是：
填充旗子已经填到第i行了，已经使用了j个旗子，而且当前第i行的状态就是x的这么一个
表示前i行的总方案数。
那么递推怎么推？
dp[i][j][x]+=dp[i-1][j-num(mark[x])][y];
解释一下，这里的x是当前第i行的状态，而这个mark[x]代表当前状态下的十进制表示，也就是说把一个状态表示成十进制之后就是mark[x]了，这里为什么是j-num(mark[x])呢？因为啊，你这样想。反过来推。如果你在前i-1行已经使用了j-num(mark[x])个旗子，而且num(mark[x])就代表第i行你使用的旗子，那么你在前i行是不是就使用了j个旗子？
就这样逆推。这个mark[x]和mark[y]分那别代表第i行和第i-1行的状态的十进制表示。
最后你只要把dp[n][k][i]其中i是从0到最多状态的那些状态，把他们加起来。

#include<iostream>  
#include<cstdio>  
#include<cstring>  
#include<string>  
#include<cmath>  
#include<algorithm>  
#include<cstdlib>  
  
using namespace std;  
  
long long n,m,k,dp[82][22][1<<9],ans,mm;//dp[i][j][x]第i行放了j个棋子当前状态为x时的方法数  
long long mark[1<<9],len;//十进制标记每一行的状态  
  
using namespace std;  
  
long long num(long long x) //记录状态x中1的个数  
{  
    long long sum=0;  
    while(x)  
    {  
        if(x&1)  
            sum++;  
        x=x>>1;  
    }  
    return sum;  
}  
  
bool judge(long long x) //判断状态x是否有相邻的棋子放在一起  
{  
  
    if(x&(x<<1)) return false;  
    return true;  
}  
  
int main()  
{  
    while(scanf("%lld %lld %lld",&n,&m,&k)!=EOF)  
    {  
           ans=0;  
           memset(dp,0,sizeof(dp));  
           memset(mark,0,sizeof(mark));  
           len = 0;  
           long long tmpp = n > m ? n : m;  
           m = n > m ? m : n; //m为小的数  
           n = tmpp;//n为大的数计为行  
             
         for(long long i=0;i<(1<<m);i++) //初始化第一行的放置方法数//剔除不合法状态  
         if(judge(i)) //若i状态没有相邻的棋子放在一起  
             {  
            dp[1][num(i)][len]=1;//则第一行状态为len（i）时1的个数为num（i）时的方法数  
            mark[len ++] = i;//标记状态             }  
               
         for(long long i=2;i<=n;i++) //第二行到第n行  
            for(long long j=0;j<=k;j++)//对于放0***n个棋子  
                for(long long x = 0 ; x < len ; x ++) //对于0***len-1个状态（第i行）//枚举  
                {  
                    for(long long y = 0 ; y < len ; y ++)//对于0***len-1个状态（第i-1行）//枚举  
                    {  
                          long long tmp = num (mark[x]);//第i行状态x中1的个数  
                        if(((mark[x] & mark[y]) == 0 ) && j >= tmp) //若上下2行没相邻的且当前的棋子数目大于此行当前状态所用的棋子  
                            dp[i][j][x] += dp[i - 1][j - tmp][y];//放法数可相加  
                            //到当前行共用了j个棋子，当前行用了tmp个棋子，状态为x，到上一行共用了j-tmp个棋子，状态为y  
                    }  
                }  
      for(long long i=0;i< len;i++) //枚举状态相加  
            ans += dp[n][k][i];  
      printf("%lld\n",ans);  
    }  
    return 0;  
}