青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品

pku3691 DNA repair 自動機+DP

題意:
給出一個DNA序列,和若干病毒DNA片段,問最少修改DNA的次數使得DNA序列中不包含病毒

解法:
首先構造DFA(再次犯了和福州現場賽一樣的錯誤,沒有+結束標記狀態的轉移。。。。WA了N久),然后就是在狀態機上的DP了
令dp[i]][j]為當前在自動機第i個節點上時suffix(j)成為合法子串需要修改的最少次數。下面狀態轉移就簡單了。。(這種狀態轉移通過方程和難描述,大家看程序吧,我寫的很清楚地~)

代碼:
  1import java.io.*;
  2import java.util.Arrays;
  3public class Main {
  4
  5    /**
  6     * @param args
  7     */

  8    static class node
  9    {
 10        node nxt[]=new node[4];
 11        int id=0;
 12        boolean end=false;
 13        node pre=null;
 14    }

 15    static node head=null,buffer[]=new node[1500],q[]=new node[1500];
 16    static int s=-1,e=-1,dp[][]=new int[1005][1005];
 17    static int c=0;
 18    static String str;
 19    static void clear(node pos)
 20    {
 21        Arrays.fill(pos.nxt, null);
 22        pos.end=false;
 23        pos.pre=null;
 24    }

 25    static void insert(String str)
 26    {
 27        node p=head;
 28        for(int i=0;i<str.length();i++)
 29        {
 30            int nxt=0;
 31            switch(str.charAt(i))
 32            {
 33                case 'A':nxt=0;break;
 34                case 'C':nxt=1;break;
 35                case 'G':nxt=2;break;
 36                case 'T':nxt=3;break;
 37            }
;
 38            if(p.nxt[nxt]==null)
 39            {
 40                clear(buffer[c]);
 41                p.nxt[nxt]=buffer[c++];
 42            }

 43            p=p.nxt[nxt];
 44        }

 45        p.end=true;
 46    }

 47    static void makeper()
 48    {
 49        s=e=-1;
 50        node p=head;
 51        for(int i=0;i<4;i++)
 52            if(p.nxt[i]==null)
 53                p.nxt[i]=p;
 54            else
 55            {
 56                p.nxt[i].pre=p;
 57                q[++e]=p.nxt[i];
 58            }

 59        while(s!=e)
 60        {
 61            p=q[++s];
 62            for(int i=0;i<4;i++)
 63            {
 64                node pre=p.pre;
 65                while(pre.nxt[i]==null) pre=pre.pre;
 66                if(p.nxt[i]==null)
 67                    p.nxt[i]=pre.nxt[i];
 68                else
 69                {
 70                    p.nxt[i].pre=pre.nxt[i];
 71                    p.nxt[i].end=(p.nxt[i].end||pre.nxt[i].end);
 72                    q[++e]=p.nxt[i];
 73                }

 74            }

 75        }

 76    }

 77    static int solve(int pos,node p)
 78    {
 79        if(p.end) return 2000;
 80        else if(pos==str.length()) return 0;
 81        else if(dp[p.id][pos]!=-1return dp[p.id][pos];
 82        else
 83        {
 84            dp[p.id][pos]=0;
 85            switch(str.charAt(pos))
 86            {
 87            case 'A':
 88                dp[p.id][pos]+=Math.min(Math.min(solve(pos+1,p.nxt[0]),solve(pos+1,p.nxt[1])+1),Math.min(solve(pos+1,p.nxt[2])+1,solve(pos+1,p.nxt[3])+1));
 89                break;
 90            case 'C':
 91                dp[p.id][pos]+=Math.min(Math.min(solve(pos+1,p.nxt[0])+1,solve(pos+1,p.nxt[1])),Math.min(solve(pos+1,p.nxt[2])+1,solve(pos+1,p.nxt[3])+1));
 92                break;
 93            case 'G':
 94                dp[p.id][pos]+=Math.min(Math.min(solve(pos+1,p.nxt[0])+1,solve(pos+1,p.nxt[1])+1),Math.min(solve(pos+1,p.nxt[2]),solve(pos+1,p.nxt[3])+1));
 95                break;
 96            case 'T':
 97                dp[p.id][pos]+=Math.min(Math.min(solve(pos+1,p.nxt[0])+1,solve(pos+1,p.nxt[1])+1),Math.min(solve(pos+1,p.nxt[2])+1,solve(pos+1,p.nxt[3])));
 98                break;
 99            }
;
100            return dp[p.id][pos];
101
102        }

103    }

104    public static void main(String[] args) throws IOException{
105        BufferedReader in=new BufferedReader(new InputStreamReader(System.in));
106        int count=1;
107        for(int i=0;i<buffer.length;i++)
108        {
109            buffer[i]=new node();
110            buffer[i].id=i;
111        }

112        while(true)
113        {
114            int num=Integer.parseInt(in.readLine());
115            if(num==0break;
116            c=1;
117            clear(buffer[0]);
118            head=buffer[0];
119            while((num--)!=0)
120              insert(in.readLine());
121            for(int i=0;i<c;i++)
122                Arrays.fill(dp[i],-1);
123            makeper();
124            str=in.readLine();
125            int ans=solve(0,head);
126            if(ans>=2000) System.out.println("Case "+count+++": -1");
127            else System.out.println("Case "+count+++""+ans);
128        }

129
130    }

131
132}

133

posted on 2011-01-12 13:21 yzhw 閱讀(260) 評論(0)  編輯 收藏 引用 所屬分類: DPstring algorithm

<2011年1月>
2627282930311
2345678
9101112131415
16171819202122
23242526272829
303112345

導航

統計

公告

統計系統

留言簿(1)

隨筆分類(227)

文章分類(2)

OJ

最新隨筆

搜索

積分與排名

最新評論

閱讀排行榜

青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品
  • <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>
            欧美激情91| 欧美影院精品一区| 日韩午夜激情av| 久久人人看视频| 老司机成人在线视频| 亚洲欧洲精品天堂一级| 国产亚洲精品资源在线26u| 欧美色欧美亚洲另类七区| 亚洲激情影院| 妖精视频成人观看www| 欧美日韩四区| 久久久久久999| 久久久久久精| 在线综合亚洲| 久久久综合免费视频| 日韩亚洲成人av在线| 久久se精品一区精品二区| 日韩亚洲一区在线播放| 亚洲视频一区| 亚洲国产欧美在线| 亚洲一区二区三区三| 在线日韩精品视频| 在线观看三级视频欧美| 一本久道久久综合狠狠爱| 亚洲欧美制服另类日韩| 99热免费精品在线观看| 久久精品一二三区| 欧美一区二区三区视频在线| 蜜乳av另类精品一区二区| 欧美一区二区三区视频| 国产精品劲爆视频| 亚洲激情视频在线| 性欧美video另类hd性玩具| 在线观看av一区| 久久综合久久久| 欧美fxxxxxx另类| 国产欧美日韩一区二区三区在线| 亚洲精品三级| 免费成人美女女| 亚洲美女毛片| 欧美三级电影一区| 亚洲在线一区二区| 欧美一区二区三区在线观看| 欧美视频官网| 新狼窝色av性久久久久久| 欧美在线视频播放| 国产综合在线看| 欧美1区2区3区| 亚洲午夜黄色| 国产日韩欧美在线| 久久久久久久久伊人| 欧美电影免费网站| aⅴ色国产欧美| 亚洲欧美中文日韩v在线观看| 亚洲欧美日韩第一区| 在线观看91精品国产麻豆| 欧美黄色大片网站| 欧美国产第一页| 国产精品99久久不卡二区| 久久久久国产精品www| 最新国产成人av网站网址麻豆| 欧美性猛片xxxx免费看久爱| 久久久国产精品一区二区三区| 欧美二区不卡| 欧美在线免费一级片| 欧美黑人多人双交| 亚洲午夜激情网页| 欧美本精品男人aⅴ天堂| 国产精品99久久久久久久vr| 国语自产偷拍精品视频偷 | 亚洲国产福利在线| 国产精品第2页| 欧美 日韩 国产一区二区在线视频| 亚洲片在线资源| 免费不卡亚洲欧美| 亚洲欧美日韩精品久久奇米色影视 | 亚洲欧美日韩另类精品一区二区三区| 久久久亚洲高清| 一区二区不卡在线视频 午夜欧美不卡'| 欧美成人午夜激情视频| 久久精品人人做人人爽电影蜜月| 亚洲国产精品欧美一二99| 久久精品夜夜夜夜久久| 亚洲日韩欧美视频一区| 欧美搞黄网站| 99成人在线| 一区二区视频免费在线观看 | 午夜精品久久一牛影视| 99精品99| 亚洲主播在线观看| 亚洲女同同性videoxma| 亚洲主播在线| 亚洲国产精品一区| 一区二区三区成人 | 欧美专区日韩视频| 亚洲免费在线观看视频| 久久精品日韩一区二区三区| 麻豆精品精品国产自在97香蕉| 亚洲网站视频福利| 欧美jjzz| 久久夜色精品国产噜噜av| 99ri日韩精品视频| 欧美在线视频免费| 美女主播精品视频一二三四| 亚洲福利视频在线| 亚洲欧美日本国产专区一区| 另类欧美日韩国产在线| 国产精品视频在线观看| 亚洲人精品午夜| 免费视频一区| 欧美在线免费观看亚洲| 亚洲日本久久| 乱码第一页成人| 国产午夜精品一区二区三区视频 | 狠狠色狠色综合曰曰| 亚洲综合欧美| 亚洲高清视频一区| 国产在线一区二区三区四区| 亚洲一区视频在线| 亚洲免费av网站| 欧美成人国产一区二区| 一区二区三区视频在线播放| 欧美国产欧美亚洲国产日韩mv天天看完整 | 欧美国产日韩精品| 欧美专区日韩视频| 国产亚洲一区精品| 久久久久久久久蜜桃| 欧美在线视频一区二区三区| 国产精品午夜av在线| 国产精品一区免费观看| 亚洲女优在线| 亚洲欧美国产制服动漫| 国产日本欧美一区二区三区在线 | 欧美丝袜一区二区三区| 亚洲另类春色国产| 99天天综合性| 欧美在线视频播放| 在线看欧美日韩| 日韩一级黄色大片| 欧美视频在线播放| 麻豆国产精品va在线观看不卡| 欧美国产精品v| 亚欧成人在线| 玖玖在线精品| 亚洲欧美激情一区二区| 国产亚洲毛片在线| 91久久精品国产91久久| 国产欧美一区二区三区沐欲| 欧美精品成人在线| 久久久精品动漫| 午夜精品国产| 亚洲精品综合精品自拍| 国产一区二区三区日韩| aⅴ色国产欧美| 久久久久国产精品人| 午夜精品久久久久久久男人的天堂 | 夜夜爽99久久国产综合精品女不卡| 一本色道久久综合亚洲精品小说| 国产一区二区| 亚洲欧美中文另类| 久久都是精品| 91久久视频| 久久久久久69| 欧美成人在线免费观看| 狠狠干综合网| 久久夜色精品亚洲噜噜国产mv| 久久国产欧美| 韩国亚洲精品| 免费在线观看日韩欧美| 亚洲第一在线综合网站| 亚洲三级免费电影| 欧美日韩国产91| 国产精品综合av一区二区国产馆| 欧美涩涩视频| 亚洲女人av| 欧美成人a∨高清免费观看| 伊人夜夜躁av伊人久久| 欧美黄在线观看| 亚洲伊人观看| 欧美aaa级| 亚洲欧美成aⅴ人在线观看| 国产欧美一区二区三区久久人妖| 亚洲国产精品小视频| 欧美aⅴ99久久黑人专区| 日韩视频一区二区| 噜噜噜在线观看免费视频日韩| 亚洲欧洲精品一区二区三区| 欧美日一区二区在线观看 | 欧美三级午夜理伦三级中文幕| 亚洲一区免费看| 亚洲二区在线| 久久精品国产v日韩v亚洲| 亚洲欧洲一区二区三区在线观看| 国产精品白丝黑袜喷水久久久 | 久久精品免费电影| 亚洲激情在线播放| 国产综合精品一区| 欧美日韩国产999| 女生裸体视频一区二区三区| 亚洲欧美另类国产|