• <ins id="pjuwb"></ins>
    <blockquote id="pjuwb"><pre id="pjuwb"></pre></blockquote>
    <noscript id="pjuwb"></noscript>
          <sup id="pjuwb"><pre id="pjuwb"></pre></sup>
            <dd id="pjuwb"></dd>
            <abbr id="pjuwb"></abbr>

            Drolca

            Apologize To Drolca
            隨筆 - 28, 文章 - 1, 評論 - 6, 引用 - 0
            數據加載中……

            poj 3691 AC自動機+DP

            #include <iostream>
            #include 
            <math.h>
            using namespace std;

            const int MAXN=20*55;
            const int MAXL=1005;
            const int K=4;// kind of child
            int dp[MAXL][MAXN];

            struct Node
            {
                Node 
            *next[K], *fail;
                
            int flag, id; // use for T[j].next[k]->id
                void Init(int index)
                
            {
                    id
            =index;
                    flag
            =0;
                    fail
            =NULL;
                    
            for(int i=0; i<K; i++)next[i]=NULL;
                }

            }
            * Q[MAXN/2], *root, T[MAXN];// Q for queue, root&T for tree;

            int index=0;
            Node 
            * newNode()
            {
                T[index].Init(index);
                
            return &T[index++];
            }


            int tokind(char k)
            {
                
            if(k=='A')return 0
                
            else if(k=='C')return 1;
                
            else if(k=='G')return 2
                
            else return 3;
            }


            void insert(char *str)
            {    
                
            if(root==NULL)
                    root
            =newNode();

                Node 
            *now=root;
                
            for(int i=0; str[i]; i++)
                
            {
                    
            int kind=tokind(str[i]);
                    
            if(now->next[kind]==NULL)
                        now
            ->next[kind]=newNode();
                    now
            =now->next[kind];
                }

                now
            ->flag=1;
            }


            int search(Node *p)
            {
                
            while(p!=root)
                
            {
                    p
            =p->fail;
                    
            if(p->flag)return 1;
                }

                
            return 0;
            }


            void buildAC()
            {
                
            int head=0, tail=0;
                root
            ->fail=NULL;

                Q[tail
            ++]=root;
                
            while(head<tail)
                
            {
                    Node 
            *now=Q[head++];
                    
            if(now->flag==0)now->flag=search(now);
                    
            for(int i=0; i<K; i++)
                    
            {
                        
            if(now->next[i]!=NULL)
                        
            {
                            
            if(now==root)
                                now
            ->next[i]->fail=root;
                            
            else
                                now
            ->next[i]->fail=now->fail->next[i];
                            Q[tail
            ++]=now->next[i];
                        }

                        
            else
                        
            {
                            
            if(now==root)
                                now
            ->next[i]=root;
                            
            else
                                now
            ->next[i]=now->fail->next[i];
                        }

                    }

                }

            }


            void update(int &a, int b){if(a==-1||b<a)a=b;}

            int slove(char *str)
            {
                memset(dp, 
            -1sizeof(dp));
                dp[
            0][0]=0;
                Node 
            *p=root;

                
            for(int i=0; str[i]; i++)
                
            {
                    
            for(int j=0; j<index; j++)
                    

                        
            if(dp[i][j]!=-1&&T[j].flag==0)
                        
            {
                            
            int kind=tokind(str[i]);
                            
            for(int k=0; k<K; k++)// need to match (k, kind)
                            {                
                                Node 
            *chd=T[j].next[k];
                                
            if(chd->flag==0)
                                    update(dp[i
            +1][chd->id],dp[i][j]+(k!=kind));    
                            }

                        }

                    }

                }

                
                
            int len=strlen(str);
                
            int res=-1;
                
            for(int ii=0; ii<index; ii++)
                    
            if(dp[len][ii]!=-1)update(res, dp[len][ii]);
                
            return res;
            }


            int main()
            {
                freopen (
            "in.txt""r", stdin);
                
            int N, Tcase=0;
                
            while(scanf("%d"&N)!=EOF&&N)
                
            {
                    root
            =NULL;
                    index
            =0;
                    
            for(int i=0; i<N; i++)
                    
            {
                        
            char s[25];
                        scanf(
            "%s", s);
                        insert(s);
                    }
                    
                    buildAC();
                    
            char ss[1005];
                    scanf(
            "%s", ss);    
                    printf(
            "Case %d: %d\n"++Tcase, slove(ss));        
                }

                
            return 0;
            }

            posted on 2012-04-13 19:28 Drolca 閱讀(243) 評論(0)  編輯 收藏 引用

            欧美日韩成人精品久久久免费看| 日韩AV无码久久一区二区| www.久久热.com| 久久久久国色AV免费看图片| 手机看片久久高清国产日韩| 久久SE精品一区二区| 国产V综合V亚洲欧美久久| 精品久久久久久无码中文野结衣| 亚洲国产天堂久久综合| 国内精品久久久久影院优| 久久99精品国产99久久6| 麻豆av久久av盛宴av| 99久久精品免费看国产免费| 中文无码久久精品| 国产呻吟久久久久久久92| 模特私拍国产精品久久| 色成年激情久久综合| 色偷偷久久一区二区三区| 午夜视频久久久久一区| 久久精品国产亚洲欧美| 一本一道久久综合狠狠老| 日韩美女18网站久久精品| 中文字幕亚洲综合久久| 1000部精品久久久久久久久| 午夜精品久久久久久| 国产精品va久久久久久久| 久久精品无码一区二区三区| 狠狠色噜噜色狠狠狠综合久久| 久久久网中文字幕| 国产精品午夜久久| 青青国产成人久久91网| 国产精品久久成人影院| 久久无码人妻一区二区三区午夜| 亚洲AV乱码久久精品蜜桃| 久久这里只有精品首页| 久久精品中文无码资源站| 精产国品久久一二三产区区别| 性做久久久久久久久浪潮| 久久天天躁狠狠躁夜夜躁2014| 久久久久综合国产欧美一区二区| 青青草原综合久久大伊人导航|