序列比对——一种更简洁、更易懂的方法

'G','a','t','c'};
int sn[37][37];
void setsn()//将矩阵的第一行第一列值全初始化为0
{
for(int i=0;i<=37;i++)
{
sn[0][i]=0;
sn[i][0]=0;
}
}
int st(int i2,int j2)//查看sn[i2][j2]是怎么得到的
st1[count1++]=cha[t1];
st2[count2++]='-';
t1--;
break;
case 4:
st1[count1++]='-';
st2[count2++]=chb[t2];
t2--;
break;
case 5:
int max(int a,int b,int c,int d);
#endif
//header.cpp
#include"header.h"
int max(int a,int b,int c,int d)//求a,b,c,d中最大的数
{
return ((a>=b?a:b)>=c?(a>=b?a:b):c)>=d?((a>=b?a:b)>=c?(a>=b?a:b):c):d;
st1[count1++]=cha[t1];
st2[count2++]=chb[t2];
t1--;
t2--;
break;
case 6:
break;
}
}
if(sn[t1][t2]==0)
{
if(t1!=0)
{
if(t2!=0)
}
int w(char x,char y)
{
if(((x>='a'&&x<='z')||(x>='A'&&x<='Z'))&&((y>='a'&&y<='z')||(y>='A'&&y<='Z'))&&x==y)
return 8;
else
{
if(((x>='a'&&x<='z')||(x>='A'&&x<='Z'))&&((y>='a'&&y<='z')||(y>='A'&&y<='Z'))&&x!=y)
{
if(i2==0&&j2!=0)//sn[i2][j2]位于第一行
return 1;
else
{
if(i2!=0&&j2==0)//sn[i2][j2]位于第一列
return 2;
else
{
if(i2!=0&&j2!=0)//sn[i2][j2]既不位于第一行也不位于第一列
for(i=1;i<37;i++)//求矩阵中的值
{
for(j=1;j<37;j++)
{
st(i,j);
}
}
for(i=0;i<37;i++)//找到最大的值
{
for(j=0;j<37;j++)
{
if(sn[i][j]>=mber)
{
setsn();
int i,j,ppp,pt;
int t1=0;
int t2=0;
char st1[100000];//存放比对后的cha数组
char st2[100000];//存放比对后的chb数组
int count1=0;
int count2=0;
int maxnumber=0;
cout<<st1[pt];
cout<<endl;
for(ppp=count2-1;ppp>=0;ppp--)//输出比对后的局部chb数组st2数组
cout<<st2[ppp];
cout<<endl;
cout<<"最大的元素值："<<maxnumber<<endl;
/* float a,b,cc;
//main.cpp
#include<iostream.h>
#include<stdlib.h>
#include"header.h"
char cha[37]={'-','t','c','c','C','A','G','T','T','A','T','G','T','C','A','G','g','g','g','a','c','a','c','g','a','g','c','a','t','g','c','a','g',
'a','g','a','c'};
char chb[37]={'-','a','a','t','t','g','c','c','g','c','c','g','t','c','g','t','t','t','t','c','a','g','C','A','G','T','T','A','T','G','T','C','A',
{
st1[count1++]=cha[t1];
st2[count2++]=chb[t2];
}
else
{
st1[count1++]=cha[t1];
st2[count2++]='-';
}
}
else
{
if(t2>0)
if(maxn==numc)
return 5;//sn[i2][j2]是由左上角的值得到的
if(maxn==0)
return 6;//sn[i2][j2]的值不是由以上三个地方得到
}
}
}
}
void main()
int c,d;
c=2;
d=3;
a=d-c;
b=d+c;
cc=(float)d/c;
cout<<a<<" "<<b<<" "<<cc<<endl;*/
}
//header.h
#ifndef HEADERH
#define HEADERH
int w(char x,char y);//比较x和y是否相同，并计算相应的得分
return -5;
else
return -3;
}
}
{
int numa=sn[i2-1][j2]+w(cha[i2],'*');
int numb=sn[i2][j2-1]+w('*',chb[j2]);
int numc=sn[i2-1][j2-1]+w(cha[i2],chb[j2]);
int maxnum=max(0,numa,numb,numc);
st1[count1++]='-';
st2[count2++]=chb[t2];
t2--;
break;
case 2:
st1[count1++]=cha[t1];
st2[count2++]='-';
t1--;
break;