AC自动机+DP小结 （一）

    好久没有更新博客了，最近真是懒到家了，南京赛前重点复习了下AC自动机+DP方面的题，写下来总结一下。

HDU 2457  DNA repair

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

题目大意：给你一个字符主串和很多病毒串，要求更改最少的字符使得没有一个病毒串是主串的子串。

 

思路：本人的第一道AC自动机DP，很初等的题目了，首先根据病毒串建立AC自动机，病毒串的尾节点设置一个标记，表示这个节点“危险”不能走，问题转化到求用主串在建立后的AC自动机上走len步（len为主串的长度）且不能经过“危险”节点需要更改的最小字符个数，设dp[i][j]表示长度为i的字符串到达j节点时所需改变的最小字符个数，然后根据AC自动机进行转移。建立失败指针的时候要注意的是如果一个节点的失败指针指向的节点是“危险的”，那么我们也要把该节点设为“危险”的，然后就直接状态转移即可。

 

#include <iostream>
#include <string.h>
#include <stdio.h>
#include <algorithm>
#include <queue>
#define maxn 1010
using namespace std;
struct node
{
    int next[4];
    int fail;
    int flag;
    void init()
    {
        memset(next,0,sizeof(next));
        fail=0;
        flag=0;
    }
}trie[maxn];
int n,tot,inf;
int dp[maxn][maxn];
void init()
{
    tot=0;
    trie[0].init();
    memset(dp,1,sizeof(dp));
    inf=dp[0][0];
}
int getnum(char x)
{
    switch(x)
    {
        case 'A':return 0;
        case 'C':return 1;
        case 'G':return 2;
        case 'T':return 3;
    }
    return 0;
}
void insert(char *str)
{
    int p=0,index;
    for(;*str!='\0';str++)
    {
        index=getnum(*str);
        if(trie[p].next[index]==0)
        {
            trie[++tot].init();
            trie[p].next[index]=tot;
        }
        p=trie[p].next[index];
    }
    trie[p].flag=1;
}
void build_fail()
{
    queue<int>Q;
    int p,son,cur,i;
    Q.push(0);
    while(!Q.empty())
    {
        p=Q.front();
        Q.pop();
        for(i=0;i<4;i++)
        {
            if(trie[p].next[i]!=0)
            {
                son=trie[p].next[i];
                cur=trie[p].fail;
                if(p==0)
                    trie[son].fail=0;
                else
                {
                    while(cur&&trie[cur].next[i]==0)
                        cur=trie[cur].fail;
                    trie[son].fail=trie[cur].next[i];
                    if(trie[trie[son].fail].flag)
                    trie[son].flag=1;
                }
                Q.push(son);
            }
            else
                {
                  //  if(p==0)
                   // trie[p].next[i]=0;
                  //  else
                    trie[p].next[i]=trie[trie[p].fail].next[i];
                }
        }
    }
}
char str[maxn],tmp[25];
void solve()
{
    dp[0][0]=0;
    int i,j,k,len=strlen(str+1);
    for(i=1;i<=len;i++)
    {
        for(j=0;j<=tot;j++)
        {
            if(dp[i-1][j]!=inf)
            {
                for(k=0;k<4;k++)
                {
                    int son=trie[j].next[k];
                    if(trie[son].flag==0)
                    {
                        dp[i][son]=min(dp[i][son],dp[i-1][j]+(getnum(str[i])!=k));
                    }
                }
            }
        }
    }
    int ans=inf;
    for(i=0;i<=tot;i++)
    {
        if(!trie[i].flag)
        ans=min(an