[poj 3691]DNA repair

好久没刷 poj 了，今天练习 AC 自动机时去水了一发喵~

在 poj 上 A 题的感觉并没有 BZOJ 上那么愉悦，准确的说是痛不欲生

真是应了那句老话，你再慢也有比你慢的，你再快也有比你快的……

跪求那些 0ms 的代码啊，还有那么多人都只跑了 32ms 啊！！

果然还是我太弱了吗？一定是我还太弱了 TAT 

一道裸裸的 AC 自动机上 dp 

令 dp[i][j] 表示母串的前 i 个字母遍历 AC 自动机，使之到达 j 节点，至少要修改多少个字母

dp[i+1][k]=min(dp[i+1][k], dp[i][j]+CHAR(j->k)!=s[i])  { k 不为匹配点}

CHAR(j->k) 是指在 AC 自动机上从点 j 转移到点 k 的边时经过的字母

答案就是 min{dp[len(s)][i]} 了喵~

我才不会说我是来晒我的 AC 自动机代码的呢~  (虽然被 poj 的各大神犇搞得一点优越感都没有……)

#include <cstdio>
#include <cstring>
const int inf=0x3F3F3F3F;
const int sizeOfText=1024;
const int sizeOfType=4;
const int sizeOfMemory=1024;

namespace trieDfa
{
    struct node
    {
        int idx;
        bool end;
        node * fail;
        node * ch[sizeOfType];
    };
    node * dfa;
    node memory[sizeOfMemory]; int port;
    node * E[sizeOfMemory];
    inline node * newnode()
    {
        node * ret=memory+port;
        E[ret->idx=port++]=ret;
        ret->end=0;
        ret->fail=NULL;
        memset(ret->ch, 0, sizeof(ret->ch));
        return ret;
    }
    inline void clear() {port=0; dfa=newnode();}

    inline int ord(char ch)
    {
        switch (ch)
        {
            case 'A':return 0;
            case 'G':return 1;
            case 'C':return 2;
            case 'T':return 3;
        }
    }
    inline void insert(char * s)
    {
        int len=strlen(s);
        node * t=dfa;
        for (int i=0;i<len;i++)
        {
            if (!t->ch[ord(s[i])]) t->ch[ord(s[i])]=newnode();
            t=t->ch[ord(s[i])];
        }
        t->end=1;
    }
    inline void buildDfa()
    {
        static node * queue[sizeOfMemory];
        int l=0, r=0;

        dfa->fail=dfa;
        for (int i=0;i<sizeOfType;i++)
            if (!dfa->ch[i]) dfa->ch[i]=dfa;
            else dfa->ch[i]->fail=dfa, queue[r++]=dfa->ch[i];

        for ( ;l<r; )
        {
            node * u=queue[l++];
            u->end|=u->fail->end;
            for (int i=0;i<sizeOfType;i++)
                if (u->ch[i])
                {
                    u->ch[i]->fail=u->fail->ch[i];
                    queue[r++]=u->ch[i];
                }
                else
                    u->ch[i]=u->fail->ch[i];
        }
    }
}
using namespace trieDfa;

int cases, n;
char str[sizeOfText];
int f[sizeOfText][sizeOfMemory];
inline int min(int x, int y) {return x<y?x:y;}
inline int dp(char * );

int main()
{
    for (scanf("%d", &n);n;scanf("%d", &n))
    {
        clear();
        for (int i=1;i<=n;i++)
        {
            scanf("%s", str);
            insert(str);
        }
        buildDfa();
        scanf("%s", str);
        printf("Case %d: %d
", ++cases, dp(str));
    }

    return 0;
}
inline int dp(char * s)
{
    int len=strlen(s);
    int ret=inf;

    memset(f, inf, sizeof(f));
    f[0][0]=0;    
    for (int i=0;i<len;i++)
        for (int j=0;j<port;j++)
            for (int k=0;k<sizeOfType;k++)
                if (!E[j]->ch[k]->end)
                    f[i+1][E[j]->ch[k]->idx]=min(f[i+1][E[j]->ch[k]->idx], f[i][j]+(ord(s[i])!=k));
    for (int i=0;i<port;i++) ret=min(ret, f[len][i]);    

    return ret==inf?-1:ret;
}