• <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>
            隨筆-65  評論-6  文章-0  trackbacks-0
             1 #include <iostream>
             2 using namespace std;
             3 #define MaxSize 102
             4 char str1[MaxSize],str2[MaxSize];
             5 int dp[MaxSize][MaxSize];
             6 int cas,len1,len2;
             7 int go[5][5]={
             8     {5,-1,-2,-1,-3},
             9     {-1,5,-3,-2,-4},
            10     {-2,-3,5,-2,-2},
            11     {-1,-2,-2,5,-1},
            12     {-3,-4,-2,-1,0}
            13 };
            14 
            15 inline int max(int a,int b,int c){
            16     a=a>c?a:c;
            17     return a>b?a:b;
            18 }
            19 
            20 inline int number(char a){
            21     switch(a) {
            22         case 'A':
            23             return 0;
            24         case 'C':
            25             return 1;
            26         case 'G':
            27             return 2;
            28         case 'T':
            29             return 3;
            30         default:
            31             return 4;
            32     }
            33 }
            34 int match(char *a,char *b,int &len_a,int &len_b){
            35     int i,j;
            36     dp[0][0]=0;
            37     for(i=1;i<=len_b;i++)
            38         dp[i][0]=dp[i-1][0]+go[number(a[i])][4];//狀態(tài)最優(yōu)解
            39     for(i=1;i<=len_a;i++)
            40         dp[0][i]=dp[0][i-1]+go[4][number(b[i])];
            41     for(i=1;i<=len_a;i++)
            42         for(j=1;j<=len_b;j++)
            43             dp[i][j]=max(dp[i-1][j-1]+go[number(a[i])][number(b[j])],
            44                          dp[i][j-1]+go[4][number(b[j])],
            45                          dp[i-1][j]+go[number(a[i])][4]);
            46     return dp[len_a][len_b];
            47 }
            48 int main(){
            49     //freopen("in.txt","r",stdin);
            50     scanf("%d",&cas);
            51     while (cas--){
            52         scanf("%d %s",&len1,str1+1);
            53         scanf("%d %s",&len2,str2+1);
            54         int ans=match(str1,str2,len1,len2);
            55         printf("%d\n",ans);
            56     }
            57     return 0;
            58 }
            posted on 2012-07-10 11:15 Leo.W 閱讀(161) 評論(0)  編輯 收藏 引用
            久久青青草原国产精品免费| 久久天天躁狠狠躁夜夜96流白浆 | 手机看片久久高清国产日韩| 国产福利电影一区二区三区,免费久久久久久久精| 久久天天躁狠狠躁夜夜avapp| 精品多毛少妇人妻AV免费久久| 精品久久久久久亚洲精品 | 亚洲国产精品无码久久SM| 久久久久高潮毛片免费全部播放| 2021国产成人精品久久| 亚洲伊人久久综合影院| 日本免费久久久久久久网站| 欧美精品丝袜久久久中文字幕 | 久久精品人人做人人爽电影蜜月| .精品久久久麻豆国产精品 | 久久久久婷婷| 狠狠干狠狠久久| 亚洲国产精品无码久久久秋霞2 | 国产欧美久久久精品影院| 国产∨亚洲V天堂无码久久久| 久久久久久一区国产精品| 国产精品久久久亚洲| 伊人久久大香线蕉综合5g| 99久久伊人精品综合观看| 亚洲精品无码久久久久去q | 久久伊人精品一区二区三区| 国产精品日韩深夜福利久久| 91精品国产9l久久久久| 亚洲女久久久噜噜噜熟女| 伊人色综合九久久天天蜜桃| 国产精品久久久久乳精品爆 | 激情伊人五月天久久综合| 国产精品99久久久精品无码| 亚洲日本va午夜中文字幕久久| 国内精品久久久久久久久电影网| 久久精品麻豆日日躁夜夜躁| 久久久久无码精品国产| 麻豆亚洲AV永久无码精品久久| 国产成人精品综合久久久| 亚洲精品高清国产一线久久| 一本一本久久A久久综合精品|