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

转载请注明出处，谢谢http://blog.csdn.net/acm_cxlove/article/details/7854526       by---cxlove

题目：给出一些不合法的模式DNA串，给出一个原串，问最少需要修改多少个字符，使得原串中不包含非法串

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

多串匹配，先想到AC自动机，需要求出最少需要修改多少字符，DP。

结合在一起

每一次沿着Trie树往下走，不能到达叶子结点罢了。不过对于为空但是合法的孩子需要进行处理。

DP方面，dp[i][j]表示前i个字符，当前为状态j的时候，需要修改的最少字符数。

从i-1的状态，找到之后的状态，如果字符与原串相同，则不变，否则+1

#include<iostream>#include<cstdio>#include<cstring>#include<cmath>#include<algorithm>#define N 100005#define MOD 100000#define inf 1<<29#define LL long longusing namespace std;struct Trie{    Trie *next[4];    Trie *fail;    int kind,isword;};Trie *que[N],s[N];int idx;int id(char ch){    if(ch=='A') return 0;    else if(ch=='T') return 1;    else if(ch=='C') return 2;    return 3;}Trie *NewNode(){    Trie *tmp=&s[idx];    for(int i=0;i<4;i++) tmp->next[i]=NULL;    tmp->isword=0;    tmp->kind=idx++;    tmp->fail=NULL;    return tmp;}void Insert(Trie *root,char *s,int len){    Trie *p=root;    for(int i=0;i<len;i++){        if(p->next[id(s[i])]==NULL)            p->next[id(s[i])]=NewNode();        p=p->next[id(s[i])];    }    p->isword=1;}void Bulid_Fail(Trie *root){    int head=0,tail=0;    que[tail++]=root;    root->fail=NULL;    while(head<tail){        Trie *tmp=que[head++];        for(int i=0;i<4;i++){            if(tmp->next[i]){                if(tmp==root) tmp->next[i]->fail=root;                else{                    Trie *p=tmp->fail;                    while(p!=NULL){                        if(p->next[i]){                           tmp->next[i]->fail=p->next[i];                           break;                        }                        p=p->fail;                    }                    if(p==NULL) tmp->next[i]->fail=root;                }                if(tmp->next[i]->fail->isword) tmp->next[i]->isword=1;                que[tail++]=tmp->next[i];            }            else if(tmp==root) tmp->next[i]=root;            else tmp->next[i]=tmp->fail->next[i];        }    }}int dp[1005][2005];int slove(char *str,int len){    for(int i=0;i<=len;i++) for(int j=0;j<idx;j++) dp[i][j]=inf;    dp[0][0]=0;    for(int i=1;i<=len;i++){        for(int j=0;j<idx;j++){            if(s[j].isword) continue;            if(dp[i-1][j]==inf) continue;            for(int k=0;k<4;k++){                int r=s[j].next[k]->kind;                if(s[r].isword) continue;                dp[i][r]=min(dp[i][r],dp[i-1][j]+(id(str[i-1])!=k));            }        }    }    int ans=inf;    for(int i=0;i<idx;i++) ans=min(ans,dp[len][i]);    return ans==inf?-1:ans;}char str[1005];int main(){    int n,cas=0;    while(scanf("%d",&n)!=EOF&&n){        idx=0;        Trie *root=NewNode();        for(int i=0;i<n;i++){            scanf("%s",str);            Insert(root,str,strlen(str));        }        Bulid_Fail(root);        scanf("%s",str);        printf("Case %d: %d\n",++cas,slove(str,strlen(str)));    }    return 0;}