社区讨论

为什么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 条回复,欢迎继续交流。

正在加载回复...