HDU-4899 Hero meet devil(DP套DP)

HDU-4899
第一次写dp套dp

观察普通LCS方程 f[i][j]=max(f[i][j-1],f[i-1][j],f[i-1][f-1]+1)
转移过程可以当作给串增加1个字符后，f[i-1]到 f[i] 的变化
于是我们只要处理好dp数组的状态加一个字符到另一个状态的转移,从长度为0的串一直加到m，
便可以知道长度为m时的dp数组的状态。对最后每个dp数组状态计算它的贡献即可。
怎么保存dp数组状态呢？
因为LCS与后一项比最多增加1。只要记录增量就好～ 即把每一项相对前一项的增量压缩为状态。

#include<bits/stdc++.h>using namespace std;const int MAXN=(1<<15)+7;const int MOD=1e9+7;const char w[]={"AGCT"};char s[30];int ans[30];int tr[MAXN][4],tot;int cnt[MAXN],g[2][MAXN],m;int f1[20],f2[20];void init(){    for(int i=0;i<1<<tot;i++)    {        for(int j=1;j<=tot;j++)//f[x-1]数组            f1[j]=f1[j-1]+(i>>(j-1)&1);        cnt[i]=f1[tot];        for(int j=0;j<4;j++)//枚举第x个字符        {            for(int k=1;k<=tot;k++)//更新dp数组            {            //f[i][j]=max(f[i][j-1],f[i-1][j],f[i-1][f-1]+1);                f2[k]=max(f2[k-1],f1[k]);//f[x]数组                if(s[k]==w[j])                    f2[k]=f1[k-1]+1;            }            tr[i][j]=0;            for(int k=1;k<=tot;k++)                if(f2[k]-f2[k-1])                    tr[i][j]|=1<<(k-1);//            printf("%d %d %d %d\n",i,j,tr[i][j],f1[tot]);        }    }}int main(){    int T;    scanf("%d",&T);    while(T--)    {        scanf("%s%d",s+1,&m);        tot=strlen(s+1);        init();        memset(ans,0,sizeof(ans));        for(int j=0;j<1<<tot;j++)            g[0][j]=0;        g[0][0]=1;        for(int i=1;i<=m;i++)            for(int j=0;j<1<<tot;j++)            {                if(g[~i&1][j])                    for(int k=0;k<4;k++)                        (g[i&1][tr[j][k]]+=g[~i&1][j])%=MOD;                g[~i&1][j]=0;            }        for(int i=0;i<1<<tot;i++)        {            (ans[cnt[i]]+=g[m&1][i])%=MOD;        }        for(int i=0;i<=tot;i++)            printf("%d\n",ans[i]);    }}