社区讨论
为什么60?
P1140相似基因参与者 3已保存回复 2
讨论操作
快速查看讨论及其快照的属性,并进行相关操作。
- 当前回复
- 2 条
- 当前快照
- 1 份
- 快照标识符
- @mi6ol9gh
- 此快照首次捕获于
- 2025/11/20 08:16 4 个月前
- 此快照最后确认于
- 2025/11/20 08:16 4 个月前
CPP
#include<iostream>
#include<cstdio>
#include<algorithm>
using namespace std;
char a[10600];
char b[10600];
int aa[10600];
int bb[10600];
int f[10600][10600];//第一串的前i个和第二串的前j个
int t[5][5]={
{5,-1,-2,-1,-3},
{-1,5,-3,-2,-4},
{-2,-3,5,-2,-2},
{-1,-2,-2,5,-1},
{-3,-4,-2,-1,0}
};
int size_a,size_b;
int main()
{
scanf("%d",&size_a);scanf("%s",&a);scanf("%d",&size_b);scanf("%s",&b);
for(int i=1;i<=size_a;i++)
for(int j=1;j<=size_b;j++)
f[i][j]=-10000;
for(int i=1;i<=size_a;i++)//将基因转为数字
{
if(a[i-1]=='A')aa[i]=0;
if(a[i-1]=='C')aa[i]=1;
if(a[i-1]=='G')aa[i]=2;
if(a[i-1]=='T')aa[i]=3;
}
for(int j=1;j<=size_b;j++)
{
if(b[j-1]=='A')bb[j]=0;
if(b[j-1]=='C')bb[j]=1;
if(b[j-1]=='G')bb[j]=2;
if(b[j-1]=='T')bb[j]=3;
}
for(int i=1;i<=size_a;i++)f[i][0]=f[i-1][0]+t[aa[i]][4];//预处理
for(int j=1;j<=size_b;j++)f[0][j]=f[0][j-1]+t[bb[j]][4];//第1个基因前i个与第2个基因第0个匹配肯定是与空碱基匹配的值
for(int i=1;i<=size_a;i++)
for(int j=1;j<=size_b;j++)
{
f[i][j]=max(f[i][j],f[i-1][j]+t[aa[i]][4]);
f[i][j]=max(f[i][j],f[i][j-1]+t[bb[i]][4]);
f[i][j]=max(f[i][j],f[i-1][j-1]+t[aa[i]][bb[j]]);
}
cout<<f[size_a][size_b];
return 0;
}
回复
共 2 条回复,欢迎继续交流。
正在加载回复...