机试题：查找基因序列问题

#include<iostream>
using namespace std;

const int num=101;
long score[num][num];
char a[num],//给定目标序列
b[num],//数据库序列
seq[num];//得到的最相似序列
int dp()
{
long i,j,m,n;
m=strlen(a);
n=strlen(b);
for(i=0;i<=m;++i) score[i]
=(m-i)*(-7);
for(j=0;j<=n;++j) score[m][j]=(n-j)*(-7);
for(i=m-1;i>=0;--i){
for(j=n-1;j>=0;--j){
if(a[i]!=b[j]) score[i][j]=score[i+1][j+1]-4;
else score[i][j]=score[i+1][j+1]+5;
if(score[i+1][j]-7>score[i][j]) score[i][j]=score[i+1][j]-7;
if(score[i][j+1]-7>score[i][j]) score[i][j]=score[i][j+1]-7;
}
}
return score[0][0];
}

int _tmain(int argc, _TCHAR* argv[])
{
long i,n,max,p;
scanf("%s%d",a,&n);
for(i=0;i<n;++i)
{
scanf("%s",b);
p=dp();
if(i==0||p>max||(p==max&&strcmp(b,seq)<0)){
strcpy(seq,b);
max=p;
}
}
printf("%ld\n%s\n",max,seq);
return 0;
}