动态规划串行算法

写了两天。。感觉确实需要多些代码了。T T

void give_global_traceback(sequence &seq1,sequence &seq2,sequence &alseq1,sequence &alseq2,int db[][k])
{
    memset(alseq1.bp, 0, sizeof(char)*k);
    memset(alseq2.bp, 0, sizeof(char)*k);
    int a = 0, b = 0, c = 0;
    //从右下角开始找起，seq2是纵向序列，seq1是横向序列
    int tb_i = seq1.size;   //矩阵中横向的记录
    int tb_j = seq2.size;   //矩阵中纵向的记录
    while(tb_i >= 0 || tb_j >= 0)
    {
        if(tb_i == 0 && tb_j == 0)   //到左上角了，跳出
        {
            break;
        }
        else if(tb_i == 0 && tb_j > 0)   //到左边沿了，往上走
        {
            alseq1.bp[al_i] = 'x';
            tb_j--;
            alseq2.bp[al_j] = seq2.bp[tb_j];
        }
        else if(tb_j == 0 && tb_i > 0)   //到上边沿了，往左走
        {
            alseq1.bp[al_i] = seq1.bp[tb_i];
            tb_i--;
            alseq2.bp[al_j] = 'x';
        }
        else
        {
            a = db[tb_j][tb_i-1];
            b = db[tb_j-1][tb_i];
            c = db[tb_j-1][tb_i-1];
            if (a > b && a > c)        //对纵向序列插空
            {
                alseq1.bp[al_i] = seq1.bp[tb_i - 1];
                tb_i--;
                alseq2.bp[al_j] = 'x';
            }
            else if (b > a && b > c)    //对横向序列插空
            {
                alseq1.bp[al_i] = 'x';
                tb_j--;
                alseq2.bp[al_j] = seq2.bp[tb_j - 1];
            }
            else                        //不插空
            {
                alseq1.bp[al_i] = seq1.bp[tb_i - 1];
                alseq2.bp[al_j] = seq2.bp[tb_j - 1];
                tb_i--;
                tb_j--;
            }
        }
        al_i++;
        al_j++;
    }
}

转载于:https://www.cnblogs.com/ubiwind/p/5043375.html