bzoj 3768: spoj 4660 Binary palindrome二进制回文串

Description

给定k个长度不超过L的01串，求有多少长度为n的01串S满足：

1.该串是回文串

2.该串不存在两个不重叠的子串，在给定的k个串中。

即不存在a<=b<c<=d，S[a,b]是k个串中的一个，S[c,d]是k个串中的一个

(It does not contain two non-overlapped substrings in the given list of K binary strings.)

举个例子，若给定2(k=2)个01串：101和1001

1010001和1010101均不满足条件。前者不满足条件1，后者不满足条件2

Input

第一行两个整数n,k

以下k行，每行一个01串

Output

输出一个整数表示答案，答案mod 1000000007（10^9+7）

对原串和反串分别建ac自动机，从两边向中间处理回文串，同时在自动机上走，若匹配到一个串就记录并回到起点

dp状态表示为 当前已确定前缀、后缀在ac自动机上的位置，以及是否已匹配到一个串 的方案数

最后检查是否在跨过中点的位置匹配到一个串

#include<cstdio>
#include<cstring>
#include<algorithm>
const int P=1e9+;
int n,k,ls[];
char s[][];
inline void inc(int&a,int b){
    b+=a-P;
    a=b+(b>>&P);
}
struct acam{
    int ch[][],fa[];
    int S[][];
    bool e[];
    int ptr;
    acam(){
        memset(this,,sizeof(acam));
        ptr=;
        ch[][]=ch[][]=;
    }
    void ins(char*s,int a){
        int w=;
        for(int i=;s[i];++i){
            int c=s[i]-'',&u=ch[w][c];
            if(!u)u=++ptr;
            w=u;
            S[w][a]|=<<i;
        }
        e[w]=;
    }
    void ins(char*s,int a,int len){
        int w=;
        for(int i=len-;i>=;--i){
            int c=s[i]-'',&u=ch[w][c];
            if(!u)u=++ptr;
            w=u;
            S[w][a]|=<<i;
        }
        e[w]=;
    }
    void build(){
        int q[],ql=,qr=;
        q[++qr]=;
        while(ql!=qr){
            int w=q[++ql];
            for(int i=;i<;++i){
                int&u=ch[w][i];
                (u?fa[q[++qr]=u]:u)=ch[fa[w]][i];
            }
        }
        for(int i=;i<=qr;++i){
            int w=q[i],f=fa[w];
            e[w]|=e[f];
            for(int j=;j<k;++j)S[w][j]|=S[f][j];
        }
    }
}t1,t2;
int f[][][][];
int main(){
    scanf("%d%d",&n,&k);
    for(int i=;i<k;++i){
        scanf("%s",s[i]);
        ls[i]=strlen(s[i]);
        t1.ins(s[i],i);
        t2.ins(s[i],i,ls[i]);
    }
    t1.build();
    t2.build();
    int now=;
    f[][][][]=;
    for(int t=;t<n/;++t,now^=){
        memset(f[now^],,sizeof(f[][])*(t1.ptr+));
        for(int c=;c<;++c){
            for(int a=;a<=t1.ptr;++a)if(!t1.e[a]){
                int a1=t1.ch[a][c],is=;
                if(t1.e[a1])a1=,is=;
                int(*f1a)[]=f[now^][a1];
                int(*f0a)[]=f[now][a];
                for(int b=;b<=t2.ptr;++b)if(!t2.e[b]){
                    int b1=t2.ch[b][c];
                    int v=is;
                    if(t2.e[b1])b1=,++v;
                    for(int x=;x+v<;++x){
                        inc(f1a[b1][x+v],f0a[b][x]);
                    }
                }
            }
        }
    }
    if(n&){
        memset(f[now^],,sizeof(f[][])*(t1.ptr+));
        for(int c=;c<;++c){
            for(int a=;a<=t1.ptr;++a)if(!t1.e[a]){
                int a1=t1.ch[a][c],is=;
                if(t1.e[a1])a1=,is=;
                int(*f1a)[]=f[now^][a1];
                int(*f0a)[]=f[now][a];
                for(int b=;b<=t2.ptr;++b)if(!t2.e[b]){
                    int v=is;
                    for(int x=;x+v<;++x){
                        inc(f1a[b][x+v],f0a[b][x]);
                    }
                }
            }
        }
        now^=;
    }
    int ans=;
    for(int a=;a<=t1.ptr;++a){
        for(int b=;b<=t2.ptr;++b){
            inc(ans,f[now][a][b][]);
            for(int c=;c<k;++c)if(t1.S[a][c]<<&t2.S[b][c])goto o;
            inc(ans,f[now][a][b][]);
            o:;
        }
    }
    printf("%d\n",ans);
    return ;
}