BZOJ3864 & HDU4899：Hero meet devil——题解

https://www.lydsy.com/JudgeOnline/problem.php?id=3864

http://acm.hdu.edu.cn/showproblem.php?pid=4899

给你字符集为{A,T,G,C}的字符串，问有多少长度为m的字符串，满足其最长公共子序列长度为0,1……|S|。

太神啦，看巨佬博客吧：https://www.cnblogs.com/RabbitHu/p/BZOJ3864.html

#include<cmath>
#include<queue>
#include<vector>
#include<cstdio>
#include<cctype>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long ll;
const int p=1e9+;
const int N=;
const int M=;
inline int read(){
    int X=,w=;char ch=;
    while(!isdigit(ch)){w|=ch=='-';ch=getchar();}
    while(isdigit(ch))X=(X<<)+(X<<)+(ch^),ch=getchar();
    return w?-X:X;
}
int s[N+];
int m,n,ans[N+],f[][<<N],last[<<N],trans[<<N][];
inline int readstring(){
    int cnt=;char ch=;
    while(ch<'A'||ch>'Z')ch=getchar();
    while(ch>='A'&&ch<='Z'){
    if(ch=='A')s[++cnt]=;if(ch=='T')s[++cnt]=;
    if(ch=='G')s[++cnt]=;if(ch=='C')s[++cnt]=;
    ch=getchar();
    }
    return cnt;
}
void init(){
    static int d[N+],g[N+];
    for(int i=;i< <<n;i++){
    if(i)last[i]=last[i^(i&-i)]+;
    for(int j=;j<n;j++)d[j+]=d[j]+(bool)(i&(<<j));
    for(int k=;k<;k++){
        for(int j=;j<=n;j++){
        g[j]=max(g[j-],d[j]);
        if(s[j]==k)g[j]=max(g[j],d[j-]+);
        }
        trans[i][k]=;
        for(int j=;j<n;j++)
        if(g[j+]-g[j])trans[i][k]|=<<j;
    }
    }
    memset(ans,,sizeof(ans));
    memset(f,,sizeof(f));
}
int main(){
    int T=read();
    while(T--){
    n=readstring();m=read();
    init();
    f[][]=;int now=;
    for(int i=;i<=m;i++){
        memset(f[now],,sizeof(f[now]));
        for(int j=;j< <<n;j++){
        for(int k=;k<;k++){
            (f[now][trans[j][k]]+=f[now^][j])%=p;
        }
        }
        now^=;
    }
    for(int i=;i< <<n;i++)
        (ans[last[i]]+=f[now^][i])%=p;
    for(int i=;i<=n;i++)printf("%d\n",ans[i]);
    }
    return ;
}