bzoj3864 Hero meet devil（ddp）

这大概就是动态动态规划了吧qaq

题意：给一个长度<=15的串S（仅包含ACGT），对于所有的i（0<=i<=|S|），求有多少长度为M的串T满足串T和串S的最长公共子序列长度为i。

我们考虑普通的O(n^2)dp求lcs,f[i][j]表示LCS(S[1..i],T[1..j])
我们有转移:
$f[i][j]=max(f[i][j-1],f[i-1][j],[S[i]==T[j]]*(f[i-1][j-1]+1))$
我们发现f[i][j]与f[i+1][j]的差值只可能是0或者1，因此我们可以差分掉，用一个状态来表示f[1][j]…f[n][j]，状态数是$O(2^n)$的。我们又发现f[1…n][j]只与f[1…n][j-1]有关，因此我们可以提前处理出串T再加上一个字符后状态S转移到的新状态S’.
然后我们再套上一个外层的dp：
dp[i][S]表示匹配到串T的第i位，f[1…n][i]状态为S的输入数。
需要滚动数组优化空间qaq

复杂度$O(4n2^n+4m2^n)$

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 10010
#define mod 1000000007
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
int bin[16],n,m,trans[33000][4],dp[2][33000],cnt[33000],ans[20];
char ch[4]={'A','C','G','T'},a[20];
inline void predp(){
    int f[20],g[20];
    for(int S=0;S<bin[n];++S){
        memset(f,0,sizeof(f));memset(g,0,sizeof(g));
        for(int i=0;i<n;++i) f[i+1]=f[i]+(bool)(S&bin[i]);
        cnt[S]=f[n];
        for(int k=0;k<4;++k){
            for(int i=1;i<=n;++i){
                g[i]=max(g[i-1],f[i]);
                if(ch[k]==a[i]) g[i]=max(g[i],f[i-1]+1);
            }trans[S][k]=0;
            for(int i=0;i<n;++i) if(g[i+1]-g[i]) trans[S][k]|=bin[i];
        }
    }
}
int main(){
//  freopen("a.in","r",stdin);
    int tst=read();bin[0]=1;
    for(int i=1;i<=15;++i) bin[i]=bin[i-1]<<1;
    while(tst--){
        scanf("%s",a+1);n=strlen(a+1);m=read();predp();
        memset(dp[0],0,sizeof(dp[0]));int p=0;dp[p][0]=1;
        for(int i=1;i<=m;++i){
            memset(dp[p^1],0,sizeof(dp[0]));
            for(int S=0;S<bin[n];++S){
                if(!dp[p][S]) continue;
                for(int k=0;k<4;++k)
                    (dp[p^1][trans[S][k]]+=dp[p][S])%=mod;
            }p^=1;
        }memset(ans,0,sizeof(ans));
        for(int i=0;i<bin[n];++i) (ans[cnt[i]]+=dp[p][i])%=mod;
        for(int i=0;i<=n;++i) printf("%d\n",ans[i]);
    }return 0;
}