• <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 閱讀(249) 評論(0)  編輯 收藏 引用

            久久久久久午夜精品| 亚洲欧美日韩久久精品 | 99久久精品国产综合一区| 亚洲一区中文字幕久久| 伊色综合久久之综合久久| 中文字幕乱码人妻无码久久| 91精品国产高清久久久久久io| 久久精品无码专区免费| 麻豆亚洲AV永久无码精品久久| 热久久这里只有精品| 亚州日韩精品专区久久久| 国产精品久久久久久一区二区三区| 国产AⅤ精品一区二区三区久久| 久久亚洲精品国产亚洲老地址| 欧美精品一区二区精品久久| 亚洲国产精品无码久久青草 | 久久久久久无码国产精品中文字幕| 超级97碰碰碰碰久久久久最新| 色综合久久中文综合网| 伊人久久大香线焦AV综合影院 | 日韩人妻无码精品久久久不卡| 久久国产影院| 99久久国产综合精品五月天喷水| 天天躁日日躁狠狠久久| 国产69精品久久久久观看软件| 丁香久久婷婷国产午夜视频| 国产精品久久久久久久| 久久男人Av资源网站无码软件| 久久九九久精品国产免费直播| 久久精品天天中文字幕人妻| 亚洲精品无码久久久久| 久久中文字幕人妻丝袜| 亚洲欧洲久久久精品| 久久亚洲天堂| 偷窥少妇久久久久久久久| 伊人伊成久久人综合网777| 人妻少妇精品久久| 久久精品国产亚洲AV久| 亚洲综合伊人久久大杳蕉| 亚洲色欲久久久综合网东京热| 久久精品aⅴ无码中文字字幕不卡|