[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=;
 const int sizeOfType=;
 const int sizeOfMemory=;

 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=;
         ret->fail=NULL;
         memset(ret->ch, , sizeof(ret->ch));
         return ret;
     }
     inline void clear() {port=; dfa=newnode();}

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

         dfa->fail=dfa;
         for (int i=;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=;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=;i<=n;i++)
         {
             scanf("%s", str);
             insert(str);
         }
         buildDfa();
         scanf("%s", str);
         printf("Case %d: %d\n", ++cases, dp(str));
     }

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

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

     return ret==inf?-:ret;
 }

本傻装B系列