BZOJ3864 hero meet devil

题目链接：https://darkbzoj.cf/problem/3864

这个是DP套DP的典型例题。

常规LCS的求法：
设dp[i][j]表示到a[i],b[j]的LCS。
那么if (a[i]!=b[j]),then (dp[i][j]=max(dp[i][j-1],dp[i-1][j]))
if (a[i]==b[j]),then (dp[i][j]=dp[i-1][j-1]+1)

现在我们知道字符串a，但是不知道b。但是很显然的，对于b字符串，每增加一位（枚举），必然会有对相应LCS的贡献。
我们设f[i][s]表示当前处理到字符串b的第i位，且当前它与a的LCS的状态为s。g[s][0/1/2/3]表示在当前LCS状态为s的时候，在后面添加'A','T','C','G’之后的LCS情况。
那么可以列出递推式(f[i][g[s][0/1/2/3]]+=f[i-1][s])
(g)数组的处理可以写到函数里，进行预处理。

然后就是要注意。。。。空间的问题，这题空间卡得太丧心病狂了。。。。。qwqwqwq

代码如下：

#include<cstring>
#include<cstdio>
#include<algorithm>
#define mod 1000000007
using namespace std;
int n,m,t;
int g[(1<<15)+2][4],f[1001][(1<<15)+2],ans[111],pre[16],now[16],a[16];
char ss[16];
inline int calc(int x)
{
    int cur_ans=0;
    for(int i=20;i>=0;i--)
        if(x&(1<<i))
            cur_ans++;
    return cur_ans;
}
inline int trans(char x)
{
    if(x=='A') return 0;
    else if(x=='T') return 1;
    else if(x=='C') return 2;
    else return 3;
}
inline int solve(int s,int k)
{
    memset(pre,0,sizeof(pre));memset(now,0,sizeof(now));
    for(int i=0;i<n;i++) pre[i+1]=pre[i]+((s>>i)&1);
    for(int i=1;i<=n;i++)
    {
        if(k==a[i]) now[i]=max(now[i],pre[i-1]+1);
        now[i]=max(max(now[i],now[i-1]),pre[i]);
    }
    int cur_ans=0;
    for(int i=0;i<n;i++) cur_ans+=(now[i+1]-now[i])*(1<<i);
    return cur_ans;
}
int main()
{

    scanf("%d",&t);
    while(t--)
    {
        memset(f,0,sizeof(f));
        memset(g,0,sizeof(g));
        memset(ans,0,sizeof(ans));
        scanf("%s",ss+1);
        n=strlen(ss+1);
        for(int i=1;i<=n;i++) a[i]=trans(ss[i])+1;
        int maxx=1<<n;
        scanf("%d",&m);
        for(int i=0;i<=maxx;i++)
            for(int j=1;j<=4;j++)
                g[i][j]=solve(i,j);
        f[0][0]=1;
        for(int i=1;i<=m;i++)
            for(int j=0;j<=maxx;j++)
                for(int k=1;k<=4;k++)
                    f[i][g[j][k]]=(f[i][g[j][k]]+f[i-1][j])%mod;
        for(int i=0;i<=maxx;i++)
            ans[calc(i)]=(ans[calc(i)]+f[m][i])%mod;
        for(int i=0;i<=n;i++)
            printf("%d
",ans[i]);
    }
    return 0;
}