ZOJ-1027-Human Gene Functions

这道题的AC率很高啊，有50%多，觉得有些奇怪，感觉应该在30%才对，不管了。

需要说明的是这一题中使用的分值矩阵map[][]是经过我变形的，这样变形后，使用起来会更加方便，尤其是对于'X'对'-'的情况可以直接减去表示的值就是它们之间的分值，当然，相应在输入时也已经把字母变成了数字方便后面直接用做map矩阵的下标查表。

这道题是很正统的DP动态规划，所以采用了标志矩阵check表示是否计算过，value矩阵存储已经算好的值。

刚开始出WA我反复检查，直到感觉消灭干净了逻辑错误，却还是WA，百思不得其解，随手把check矩阵的初始化范围改成MAXN而不是仅仅初始使用的部分，结果就AC了。。。结果速度还真是快，只是C++比C用的内存大，大就大，反正习惯了。

时间：00:00.00
内存：904K

#include  < iostream >
using   namespace  std;

#define  MAXN 101
/**/ //////////////////////-   T   G   A   C
int  map[ 5 ][ 5 ]  =   ... {...{0, -1, -2, -3, -4},//-
                 ...{-1, 5, -2, -1, -2},//T
                 ...{-2, -2, 5, -2, -3},//G
                 ...{-3, -1, -2, 5, -1},//A
                 ...{-4, -2, -3, -1, 5}} ; // C
bool  check[MAXN][MAXN]; // 标志矩阵，标志是否计算过
int  value[MAXN][MAXN]; // 存储已经计算过的值                                

void  readd( int  x[],  int   & len)
... {
    char temp[MAXN];
    cin >> len;
    cin >> temp;
    for(int i=0;i<len;i++)
    ...{
        if('T' == temp[i])      x[i] = 1;
        else if('G' == temp[i]) x[i] = 2;
        else if('A' == temp[i]) x[i] = 3;
        else                    x[i] = 4;//C    
    }
}

int  sim( int  a[],  int  b[],  int  alen,  int  blen)
... {
    if(check[alen][blen])//已经算过的
    ...{
        return value[alen][blen];
    }
    else//未算过的
    ...{            
        if(alen>0 && blen>0)
        ...{        
            int kv,temp;        
            
            kv =   sim(a,b,alen,blen-1) - b[blen-1];
            temp = sim(a,b,alen-1,blen) - a[alen-1];
            if(temp>kv) 
                kv = temp;
            temp = sim(a,b,alen-1,blen-1) + map[a[alen-1]][b[blen-1]];
            if(temp>kv) 
                kv = temp;
            
            check[alen][blen] = true;
            value[alen][blen] = kv;
            return kv;
        }
        else if((0==alen) && blen>0)
        ...{
            int temp = sim(a,b,0,blen-1)-b[blen-1];
            check[0][blen] = true;
            value[0][blen] = temp;
            return temp;            
        }
        else if((0==blen) && alen>0)
        ...{
            int temp = sim(a,b,alen-1,0)-a[alen-1];
            check[alen][0] = true;
            value[alen][0] = temp;
            return temp;
        }
    }
}

void  solve()
... {
    int len1,len2;
    int x1[MAXN],x2[MAXN];
    readd(x1,len1);
    readd(x2,len2);
    for(int i=0;i<MAXN;i++)
    ...{
        for(int j=0;j<MAXN;j++)
            check[i][j] = false;
    }
    check[0][0] = true;
    value[0][0] = 0;
    cout << sim(x1,x2,len1,len2) << endl;
}

int  main()
... {
    int t;
//    freopen("1027.txt","r",stdin);
    cin >> t;
    while(t--)
        solve();
//    fclose(stdin);
    return 0;
}