bzoj 3864 Hero meet devil [dp套dp]

Description:
给一个长为$n$的串$s$($n\leq15$),求$lcs=[0,n]$长度为$m$的串的个数

---

Solution:
考虑求$lcs$的过程，$dp[i][j]=dp[i-1][j-1]+1 (a[i]==b[j])$
那么可以考虑状压当前$lcs$,设计出状态$dp[i][s]$表示当前到了$i$,$lcs$为$s$的方案数。
预处理出$trans[s][k]$表示状态$s$加上字符$k$到达的状态，那么就可以转移了$dp[i][trans[s][k]]+=dp[i-1][s]$
预处理$trans$枚举状态$s$,计算当前$lcs$数组f,然后枚举每位的下一个字符,计算$g[i]$表示添加了新的字符后的$lcs$数组。具体转移为$g[i]=max(g[i-1],f[i])),g[i]=max(g[i],f[i-1]+1)$。$f$相当于$dp[1...i][j-1]$,$g$相当于$dp[1...i][1...j]$
最后枚举 g g <script type="math/tex" id="MathJax-Element-1139">g</script>数组差分计算状态。

---

#include <bits/stdc++.h>
using namespace std;
const int P = 1e9 + 7;
int n, m;
char s[20], c[5] = "ACGT";
int dp[2][1 << 16], trans[1 << 16][4], ans[20];
void ini() {
    static int f[20], g[20];
    for(int i = 0; i < 1 << n; ++i) {
        for(int j = 0; j < n; ++j) {
            f[j + 1] = f[j] + (i >> j & 1);
        }   
        for(int k = 0; k < 4; ++k) {
            for(int j = 1; j <= n; ++j) {
                g[j] = max(g[j - 1], f[j]);
                if(s[j] == c[k]) {
                    g[j] = max(g[j], f[j - 1] + 1);
                }
            }
            trans[i][k] = 0;
            for(int j = 0; j < n; ++j) {
                if(g[j + 1] - g[j]) {
                    trans[i][k] |= 1 << j;
                }
            }
        }
    }
}
int main() {
    int T;
    scanf("%d", &T);
    while(T--) {
        scanf("%s", s + 1);
        n = strlen(s + 1);
        scanf("%d", &m);
        ini();
        memset(dp, 0, sizeof(dp));  
        memset(ans, 0, sizeof(ans));
        dp[1][0] = 1;   
        for(int k = 1, cur = 0; k <= m; ++k, cur ^= 1) {
            memset(dp[cur], 0, sizeof(dp[cur]));
            for(int i = 0; i < 1 << n; ++i) {
                for(int j = 0; j < 4; ++j) {
                    dp[cur][trans[i][j]] = (dp[cur][trans[i][j]] + dp[cur ^ 1][i]) % P;
                }   
            }   
        }
        for(int i = 0; i < 1 << n; ++i) {
            int t = __builtin_popcount(i);
            ans[t] = (ans[t] + dp[(m & 1) ^ 1][i]) % P;
        }
        for(int i = 0; i <= n; ++i) {
            printf("%d\n", ans[i]);
        }
    }
    return 0;
}