hdu2457AC自动机+DP

标签：

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

2
AAA
AAG
AAAG    
2
A
TG
TGAATG
4
A
G
C
T
AGT
0

Case 1: 1
Case 2: 4
Case 3: -1

/**
hdu2457 AC自动机+DP
题目大意：给定一些模式串，对于一个字符串最少改动几个字母就能保证该字符串中不含任何一个模式串
解题思路：dp[i][j]表示长度为i以状态j为结尾改动最小值，用自动机转移一下就可以了
*/
#include <stdio.h>
#include <string.h>
#include <iostream>
#include <queue>
#include <algorithm>
using namespace std;
const int INF=0x3f3f3f3f;
struct Trie
{
    int next[1010][4],fail[1010],end[1010];
    int root,L;
    int newnode()
    {
        for(int i=0;i<4;i++)
        {
            next[L][i]=-1;
        }
        end[L++]=false;
        return L-1;
    }
    void init()
    {
        L=0;
        root=newnode();
    }
    int getch(char ch)
    {
        if(ch=='A')return 0;
        if(ch=='C')return 1;
        if(ch=='G')return 2;
        return 3;
    }
    void insert(char buf[])
    {
        int len=strlen(buf);
        int now=root;
        for(int i=0;i<len;i++)
        {
            if(next[now][getch(buf[i])]==-1)
                next[now][getch(buf[i])]=newnode();
            now=next[now][getch(buf[i])];
        }
        end[now]=true;
    }
    void build()
    {
        queue<int>Q;
        fail[root]=root;
        for(int i=0;i<4;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();
            if(end[fail[now]])end[now]=true;
            for(int i=0;i<4;i++)
            {
                if(next[now][i]==-1)
                    next[now][i]=next[fail[now]][i];
                else
                {
                    fail[next[now][i]]=next[fail[now]][i];
                    Q.push(next[now][i]);
                }
            }
        }
    }
    int dp[1010][1010];
    int solve(char *buf)
    {
        int len=strlen(buf);
        for(int i=0;i<=len;i++)
        {
            for(int j=0;j<L;j++)
            {
                dp[i][j]=INF;
            }
        }
        dp[0][root]=0;
        for(int i=0;i<len;i++)
        {
            for(int j=0;j<L;j++)
            {
                if(dp[i][j]<INF)
                {
                    for(int k=0;k<4;k++)
                    {
                        int news=next[j][k];
                        if(end[news])continue;
                        int tmp;
                        if(k==getch(buf[i]))tmp=dp[i][j];
                        else tmp=dp[i][j]+1;
                        dp[i+1][news]=min(dp[i+1][news],tmp);
                    }
                }
            }
        }
        int ans=INF;
        for(int j=0;j<L;j++)
        {
            ans=min(ans,dp[len][j]);
        }
        if(ans==INF)ans=-1;
        return ans;
    }
}ac;
char buf[1010];
int main()
{
    int n,tt=0;
    while(~scanf("%d",&n))
    {
        if(n==0)break;
        ac.init();
        while(n--)
        {
            scanf("%s",buf);
            ac.insert(buf);
        }
        ac.build();
        scanf("%s",buf);
        printf("Case %d: %d\n",++tt,ac.solve(buf));
    }
    return 0;
}

2
AAA
AAG
AAAG    
2
A
TG
TGAATG
4
A
G
C
T
AGT
0

Case 1: 1
Case 2: 4
Case 3: -1

hdu2457AC自动机+DP

标签：

原文地址：http://blog.csdn.net/lvshubao1314/article/details/45485121