码迷,mamicode.com
首页 > 其他好文 > 详细

poj 1080 (LCS变形)

时间:2016-03-16 23:58:29      阅读:453      评论:0      收藏:0      [点我收藏+]

标签:

Human Gene Functions

题意:

LCS:

设dp[i][j]为前i,j的最长公共序列长度;

dp[i][j] = dp[i-1][j-1]+1;(a[i] == b[j])

dp[i][j] = max(dp[i][j-1],dp[i-1][j]);

边界:dp[0][j] = 0(j<b.size) ,dp[i][0] = 0(i< a.size);

LCS变形:

设dp[i][j]为前i,j的最大价值:

value(x, y)为比较价值;

dp[i][j] = max(dp[i-1][j-1] + value(a[i],b[i]),dp[i][j-1] + value(‘-‘, b[i]), dp[i-1][j] + value(a[i],‘-‘);

边界:dp[i][0] = dp[i-1][0] + value(a[i],‘-‘),(i<=a.size);      dp[0][j] = dp[0][j-1] + value(‘-‘,b[i]), (j <= b.size);

边界是个大问题!!

 a[i]跟dp[i]的关系搞错.跪了整整一个小时啊!!!

#include <iostream>
#include <algorithm>
#include <stdlib.h>
#include <time.h>
#include <cmath>
#include <cstdio>
#include <string>
#include <cstring>
#include <vector>
#include <queue>
#include <stack>
#include <set>

#define c_false ios_base::sync_with_stdio(false); cin.tie(0)
#define INF 0x3f3f3f3f
#define INFL 0x3f3f3f3f3f3f3f3f
#define zero_(x,y) memset(x , y , sizeof(x))
#define zero(x) memset(x , 0 , sizeof(x))
#define MAX(x) memset(x , 0x3f ,sizeof(x))
#define swa(x,y) {LL s;s=x;x=y;y=s;}
using namespace std ;
#define N 105

const double PI = acos(-1.0);
typedef long long LL ;
char a[N],b[N];
int dp[N][N];

int value(char x, char y){
    if(x == y)return 5;
    else if((x == A&&y == C)||(x == C&&y == A))return -1;
    else if((x == A&&y == G)||(x == G&&y == A))return -2;
    else if((x == A&&y == T)||(x == T&&y == A))return -1;
    else if((x == C&&y == G)||(x == G&&y == C))return -3;
    else if((x == C&&y == T)||(x == T&&y == C))return -2;
    else if((x == T&&y == G)||(x == G&&y == T))return -2;
    else if((x == A&&y == -)||(x == -&&y == A))return -3;
    else if((x == C&&y == -)||(x == -&&y == C))return -4;
    else if((x == G&&y == -)||(x == -&&y == G))return -2;
    else if((x == T&&y == -)||(x == -&&y == T))return -1;
    else return 0;
}
int main(){
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    int n,aL,bL;
    cin>>n;
    while(n--){
        cin>>aL;
        scanf("%s",a);
        cin>>bL;
        scanf("%s",b);
        dp[0][0] = 0;
        for(int i = 0;i < aL; i++) {dp[i+1][0] = dp[i][0] + value(a[i],-);}
        for(int i = 0;i < bL; i++) {dp[0][i+1] = dp[0][i] + value(b[i],-);}
        for(int i = 1;i <= aL; i++){
            for(int j = 1;j <= bL; j++){
                if(a[i-1] == b[j-1]) dp[i][j] = dp[i-1][j-1] + value(a[i-1], b[j-1]);
                else dp[i][j] = max(dp[i-1][j-1]+value(a[i-1],b[j-1]), max(dp[i][j-1]+value(-,b[j-1]), dp[i-1][j]+value(a[i-1],-)));
            }
        }
        cout<<dp[aL][bL]<<endl;
    }
    return 0;
}

 

poj 1080 (LCS变形)

标签:

原文地址:http://www.cnblogs.com/yoyo-sincerely/p/5285182.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!