bzoj5336 [TJOI2018]party（ddp）

动态动态dp qaq
观察lcs的转移方程，我们发现f[i][j]与f[i][j+1]最多差1，于是可以差分掉状压。
复杂度$O(n2^m)$

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 33000
#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 n,m,trans[N][3],f[2][3][N],bin[20],g[20],h[20],ans[20];
//f[0][i][S],前i个串,状态为S的方案数,0-,1-N,2-NO
char s[20],op[4]={'N','O','I'};
inline void init(int S){
    g[0]=0;for(int i=1;i<=m;++i) g[i]=g[i-1]+(bool)(S&bin[i-1]);
    for(int i=0;i<3;++i){
        for(int j=1;j<=m;++j){
            h[j]=max(g[j],h[j-1]);
            if(s[j]==op[i]) h[j]=max(h[j],g[j-1]+1);
        }int ss=0;
        for(int j=1;j<=m;++j) if(h[j]-h[j-1]) ss|=bin[j-1];
        trans[S][i]=ss;
    }
}
inline void inc(int &x,int y){x+=y;x%=mod;}
int main(){
//  freopen("party.in","r",stdin);
//  freopen("party.out","w",stdout);
    n=read();m=read();scanf("%s",s+1);
    bin[0]=1;for(int i=1;i<=m;++i) bin[i]=bin[i-1]<<1;
    for(int i=0;i<=bin[m]-1;++i) init(i);int p=0;
    f[0][0][0]=1;
    for(int i=0;i<n;++i){
        for(int S=0;S<=bin[m]-1;++S){
            if(f[p][0][S]){
                inc(f[p^1][1][trans[S][0]],f[p][0][S]);
                inc(f[p^1][0][trans[S][1]],f[p][0][S]);
                inc(f[p^1][0][trans[S][2]],f[p][0][S]);
            }if(f[p][1][S]){
                inc(f[p^1][1][trans[S][0]],f[p][1][S]);
                inc(f[p^1][2][trans[S][1]],f[p][1][S]);
                inc(f[p^1][0][trans[S][2]],f[p][1][S]);
            }if(f[p][2][S]){
                inc(f[p^1][1][trans[S][0]],f[p][2][S]);
                inc(f[p^1][0][trans[S][1]],f[p][2][S]);
            }
        }memset(f[p],0,sizeof(f[p]));p^=1;
    }for(int S=0;S<=bin[m]-1;++S){
        int res=0;for(int j=0;j<m;++j) if(S&bin[j]) ++res;
        for(int j=0;j<3;++j) inc(ans[res],f[p][j][S]);
    }for(int i=0;i<=m;++i) printf("%d\n",ans[i]);
    return 0;
}