HDU 2457 POJ 3691 DNA repair（AC自动机+DP）

给\(n\)个模式串和一个串\(s\)，要求修改\(s\)最少的字符使得没有一个模式串是\(s\)的子串

\(dp[i][j]\)表示长度为\(i\)，到达\(AC\)自动机的节点\(j\)所需修改的最少字符数量

#include
using namespace std;
typedef long long ll;
const int maxn=2000;
struct Trie{
    int next[maxn][26], fail[maxn], end[maxn];
    int dp[maxn][maxn];
    int last[maxn];
    int root, L;
    int newnode(){
        for (int i = 0; i < 26; i++)
            next[L][i] = -1;
        fail[L]=last[L]=end[L]=0;
        return L++;
    }
    void init(){
        L = 0;
        root = newnode();
    }
    void insert(char buf[]){
        int len = strlen(buf);
        int now = root;
        for(int i=0;iQ;
        fail[root]=root;
        for(int i = 0;i < 26;i++){
            if(next[root][i]==-1)
                next[root][i]=root;
            else{
                fail[next[root][i]] = root;
                Q.push(next[root][i]);
            }
        }
        while(!Q.empty()){
            int now = Q.front();
            Q.pop();
            end[now]+=end[fail[now]];
            for(int i = 0;i < 26;i++){
                if(next[now][i] == -1) 
                    next[now][i] = next[fail[now]][i];
                else{
                    int son=next[now][i];
                    fail[son]=next[fail[now]][i];
                    last[son]=end[fail[son]]?fail[son]:last[fail[son]];
                    Q.push(next[now][i]); 
                }  
            }
        }
    }
    char to[4]={'A','G','C','T'};
    void solve(char s[]){
        memset(dp,0x3f,sizeof(dp));
        dp[0][0]=0;
        int n=strlen(s);
        for(int i=1;i<=n;i++){
            for(int j=0;j