BZOJ[3075][Usaco2013]Necklace AC自動機+dp

傳送門ber~

要讓B不是A的子串，換句話說，就是不讓A匹配到頭(跑到AC自動機的最下面)
令fi,j$f_{i,j}$ 表示串上的位置i$i$ 匹配到樹上的結點j$j$ 最大長度(A串前i$i$ 個字符中最多選多少個字符正好匹配到節點j$j$ )
每步可走可不走，走的話就+1，不走就直接轉移
注意不能走到頭(代表可以匹配了)

代碼如下：

#include<cstring>
#include<ctype.h>
#include<cstdio>
#include<queue>
#define N 10050
using namespace std;
char s[N],s1[N];
int len,len1,cnt,ans;
int f[N][1010];
struct Node{
    Node *ch[26],*nex;
    int id;
    bool b;
    Node():b(false){
        id=++cnt;
        for(int i=0;i<26;i++)
            ch[i]=NULL;
    }
}*root=new Node,*pre[N];
queue<Node*>q;
struct Aho_Corasick_Automaton{
    inline void Insert(char *s){
        int len=strlen(s+1);
        Node *x=root;
        for(int i=1;i<=len;i++){
            if(!x->ch[s[i]-'a'])
                x->ch[s[i]-'a']=new Node;
            pre[x->id]=x;
            x=x->ch[s[i]-'a'];
        }
        pre[x->id]=x;
        x->b=true;
        return;
    }
    inline void GetFail(){
        for(int i=0;i<26;i++){
            if(root->ch[i]) root->ch[i]->nex=root,q.push(root->ch[i]);
            else root->ch[i]=root;
        }
        while(!q.empty()){
            Node *x=q.front();q.pop();
            for(int i=0;i<26;i++)
                if(!x->ch[i]) x->ch[i]=x->nex->ch[i];
                else x->ch[i]->nex=x->nex->ch[i],q.push(x->ch[i]);
        }
        return;
    }
    inline void Solve(){
        int len=strlen(s+1);
        memset(f,-1,sizeof f);
        f[1][root->id]=0;f[1][root->ch[s[1]-'a']->id]=1;
        for(int i=1;i<len;i++)
            for(int j=1;j<=cnt;j++)
                if(~f[i][j]){
                    Node *x=pre[j]->ch[s[i+1]-'a'];
                    if(!x->b) f[i+1][x->id]=max(f[i+1][x->id],f[i][j]+1);
                    if(!pre[j]->b) f[i+1][j]=max(f[i+1][j],f[i][j]);
                }
    }
}AC;
int main(){
    scanf("%s%s",s+1,s1+1);
    len=strlen(s+1),len1=strlen(s1+1);
    AC.Insert(s1);AC.GetFail();
    AC.Solve();
    for(int i=1;i<=cnt;i++)
        if(~f[len][i] && !pre[i]->b)
            ans=max(ans,f[len][i]);
    printf("%d\n",len-ans);
return 0;
}