• <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, 評(píng)論 - 6, 引用 - 0
            數(shù)據(jù)加載中……

            poj 3691 AC自動(dòng)機(jī)+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 閱讀(252) 評(píng)論(0)  編輯 收藏 引用


            只有注冊(cè)用戶登錄后才能發(fā)表評(píng)論。
            網(wǎng)站導(dǎo)航: 博客園   IT新聞   BlogJava   博問   Chat2DB   管理


            久久久久久国产精品美女| 亚洲AV成人无码久久精品老人| 一级做a爰片久久毛片16| 久久99精品久久久久久动态图| 日本福利片国产午夜久久| 国产欧美久久久精品影院| 久久久噜噜噜久久中文福利| 99久久伊人精品综合观看| 久久久久久久久久免免费精品| 久久男人Av资源网站无码软件| 嫩草伊人久久精品少妇AV| 成人久久综合网| 亚洲精品乱码久久久久久久久久久久 | 91精品国产高清91久久久久久| 99精品久久精品| 亚洲精品美女久久久久99小说| 国产韩国精品一区二区三区久久 | 国产精品免费久久久久影院| 怡红院日本一道日本久久| 亚洲精品乱码久久久久久久久久久久| 精品久久777| 久久男人Av资源网站无码软件| 婷婷久久综合| 久久人人爽人人爽AV片| 久久最新精品国产| 国产一区二区三区久久精品| 久久精品一区二区三区AV| 青青青青久久精品国产h久久精品五福影院1421 | 久久91精品国产91久久户| 色综合久久久久综合体桃花网 | 久久精品无码一区二区app| 久久免费美女视频| 青青热久久综合网伊人| 国产成人久久精品激情| 久久精品国产清高在天天线| 无码AV中文字幕久久专区| 亚洲狠狠婷婷综合久久蜜芽| 五月丁香综合激情六月久久| 亚洲中文字幕无码久久精品1| 亚洲伊人久久精品影院| 日日噜噜夜夜狠狠久久丁香五月 |