【算法讲21：轮廓线的状压DP入门】Campus Design | HDU4804

前置

• $dp$ 思想
  状压 $dp$ 的基本内容
  一些二进制的基本运算

题目 $\lceil CampusDesign\rfloor$

• Campus Design | HDU4804

题意

• $n\times m$ 的网格，每个格子如果是 $1$ ，表示为 $\lceil$ 空地 $\rfloor$ ；是 $0$，表示为$\lceil$ 花 $\rfloor$ 。
  你可以使用无数量限制的 $1\times 2$ 的 $\lceil$ 矩形砖块 $\rfloor$ (可以旋转成 $2\times 1$ 的砖)
  你可以使用至少 $c$ 个，至多 $d$ 个 $1\times 1$ 的 $\lceil$ 方形砖块 $\rfloor$。
  问你，用这些砖块不重叠放置，且四角只能放在格点上，铺满所有的空地而不铺到花上面的方案数
  答案取模 $1e9+7$

范围

• $1\le n\le 100$
  $1\le m\le 10$
  $0\le c\le d\le 20$

$\lceil$ 轮廓线 $dp$ $\rfloor$

• 看到列数，想到 $2^m$ 在范围之内就能很快想到用状压去做。
  但是易得我们的状态数为 $dp[n][2^m][d]$
  如果我们枚举相邻的两行，时间复杂度至少为 $O(n\times m\times 2^{2m}\times d)=2e10$，肯定会 $TLE$ 了
  所以我们就想到，是不是能只枚举一行也能保存前面所有的状态呢？答案是可以的。
• 如果我们枚举到 $(i,j)$，那么枚举的状态集合 $S$ 对应的是该位置前 $m$ 个格子的放砖块的状态。
  如果$S$ 的第 $x$ 位为 $1$，就表示第 $x$ 列的白框的该位置我们放了砖块。反之就没有放砖。
  注意，不放砖可能是我们故意空着，也可能是该位置为花，没法放砖。
  
• 接下来，我们考虑状态转移。
  • 如果 $tu[i][j]=1$，表示该位置是$\lceil$ 空地 $\rfloor$。决策：竖着放$1\ast 2$ 砖，横着放 $2\ast 1$ 砖，放 $1\ast 1$ 砖或者不放砖。
    • 如果 $S$ 的第 $j$ 位为 $0$，且 $tu[i-1][j]=1$ 表示该位置必须竖着放$1\ast 2$ 砖，否则上面的位置无法被覆盖。
    • 如果 $S$ 的第 $j$ 位为 $1$，且 $S$ 的第 $j-1$ 位为 $0$，表示该位置可以横着放 $2\ast 1$ 砖
    • 如果 $S$ 的第 $j$ 位为 $1$，表示该位置可以放 $1\ast 1$ 砖
    • 如果 $S$ 的第 $j$ 位为 $1$，表示该位置可以不放砖
  • 如果 $tu[i][j]=0$，表示该位置是$\lceil$ 花 $\rfloor$。
    • 如果 $S$ 的第 $j$ 位为 $0$，表示该位置无法达成满足要求的情况。
    • 如果 $S$ 的第 $j$ 位为 $1$，表示该位置必须不放砖。
      
• 这样，经过该位置的转移之后，我们的轮廓线转移成功！
  经过一步一步的转移，就可以得到答案了。
  
• 还有一些细节的问题：
  • 这里，我们设 $dp[i][j][S][T]$ 表示枚举到第 $i$ 行，第 $j$ 列，轮廓线状态为 $S$，用的 $1\ast 1$ 砖数量为 $T$，首先这样会 $MLE$。
    容易想到，枚举完某条轮廓线之后，上一条轮廓线就没有用了。
    我们直接设 $dp[st][S][T]$ 滚动数组即可。其中 $st=0$ 表示目前轮廓线， $st=1$ 表示下一条轮廓线。
  • 初状态，我们假设第 $0$ 行全部都铺满砖，因为我们是从第 $1$ 行开始枚举的。
    所以初状态为 $dp[st][(1<<m)-1][0]=1$，方案数为 $1$。
  • 末状态，对于最后一行，我们需要最后一条轮廓线的空位置都被填满。
    答案即为 ${\sum }_{i=c}^{d}dp[st][shu][i]$，其中 $shu$ 表示最后一行空地都为 $1$的状态。

代码

• 时间复杂度：$O(n\times m\times 2^m\times d)$
  空间复杂度：$O(nm+2^m\times d)$
  $Time(Ms):1326/8000$

/*
 _            __   __          _          _
| |           \ \ / /         | |        (_)
| |__  _   _   \ V /__ _ _ __ | |     ___ _
| '_ \| | | |   \ // _` | '_ \| |    / _ \ |
| |_) | |_| |   | | (_| | | | | |___|  __/ |
|_.__/ \__, |   \_/\__,_|_| |_\_____/\___|_|
        __/ |
       |___/
*/
const ll MOD = 1e9+7;
ll dp[2][(1<<12)][25];
char tu[105][20];
int shu;
int main()
{
    int n,m,c,d;
    while(~scanf("%d%d%d%d",&n,&m,&c,&d)){
        for(int i = 1;i <= n;++i)scanf("%s",tu[i]);
        shu = 0;
        for(int j = 0;j < m;++j){
            if(tu[n][j] == '1')shu |= (1<<j);
        }
        memset(dp,0,sizeof(dp));
        dp[0][(1<<m)-1][0] = 1;
        int st = 0;
        for(int i = 1;i <= n;++i)for(int j = 0;j < m;++j){
            for(int S = 0;S < (1<<m);++S){
                for(int T = 0;T <= d;++T){
                    /// 空过
                    if(tu[i][j] == '0'){
                        if(i >=2 && tu[i-1][j] == '1' && !(S&(1<<j)))dp[st^1][S][T] = 0;
                        else dp[st^1][S&(~(1<<j))][T] = (dp[st^1][S&(~(1<<j))][T] + dp[st][S][T]) % MOD;

                    }
                    else{
                        /// 上面要填 1*2 必须填
                        if(i >=2 && tu[i-1][j] == '1' && !(S&(1<<j)))
                            dp[st^1][S|(1<<j)][T] = (dp[st^1][S|(1<<j)][T] + dp[st][S][T]) % MOD;
                        else{
                            /// 横着填 2*1
                            if(j && tu[i][j-1] == '1' && !(S&(1<<j-1)))
                                dp[st^1][S|(1<<j-1)|(1<<j)][T] = (dp[st^1][S|(1<<j-1)|(1<<j)][T] + dp[st][S][T]) % MOD;
                            /// 填 1*1
                            dp[st^1][S|(1<<j)][T+1] = (dp[st^1][S|(1<<j)][T+1] + dp[st][S][T]) % MOD;
                            /// 不填
                            dp[st^1][S&(~(1<<j))][T] = (dp[st^1][S&(~(1<<j))][T] + dp[st][S][T]) % MOD;
                        }
                    }
                }
            }
            for(int S = 0;S < (1<<m);++S)
                for(int T = 0;T <= d;++T){
                    dp[st][S][T] = dp[st^1][S][T];
                    dp[st^1][S][T] = 0;
                }
        }
        ll res = 0;
        for(int T = c;T <= d;++T)
            res = (res + dp[st][shu][T]) % MOD;
        printf("%lld\n",res);
    }
    return 0;
}