题目:https://www.luogu.org/problemnew/show/P1140

分析:

本题一看就知道是一道动归,其实和字串距离非常的像,只不过多了题目规定的匹配相似度罢了。

匹配的相似度我们之间用一个二维数组读入即可

int shuzu[6][6]={{0,0,0,0,0,0},{0,5,-1,-2,-1,-3},{0,-1,5,-3,-2,-4},{0,-2,-3,5,-2,-2},{0,-1,-2,-2,5,-1},{0,-3,-4,-2,-1,0}};

PS:换行效果更佳。

然后要把字符串转化成刚才数组中的下标,便于读写

int start(char c)
{
	if(c=='A')return 1;
	if(c=='C')return 2;
	if(c=='G')return 3;
	if(c=='T')return 4;
}

然后直接在主函数中调用即可

之后就是调用问题

for(int i=0;i<n;i++)
	{
		a[i+1]=start(s1[i]);
	}
	for(int j=0;j<m;j++)
	{
		b[j+1]=start(s2[j]);
	}

其中s1,s2均为读入的字符串,我们把它们分别逐字符转化放入a,b数组。

然后由于有负值出现,我们需要把动归的核心方程初始化为一个极小的负数

for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=m;j++)
		{
			f[i][j]=-2147483647;
		}
	}

然后就是关键的动归核心部分了。

首先,看到数据范围:100

受到启发:二位数组开的起,再加上字串距离的引导,我们很容易想到 f [ i ] [ j ] f[i][j] f[i][j]可以表示第一个字符串的前i个字符与第二个字符串的前j个字符匹配的最大值。

状态想出来,那么方程如何转移呢?

根据可以加入空碱基,我们能想到

f [ i ] [ j ] = m a x ( f [ i ] [ j ] , f [ i 1 ] [ j ] + s h u z u [ a [ i ] ] [ 5 ] , f [ i ] [ j 1 ] + s h u z u [ b [ j ] ] [ 5 ] , f [ i 1 ] [ j 1 ] + s h u z u [ a [ i ] ] [ b [ i ] ] f[i][j]=max(f[i][j],f[i-1][j]+shuzu[a[i]][5],f[i][j-1]+shuzu[b[j]][5],f[i-1][j-1]+shuzu[a[i]][b[i]] f[i][j]=max(f[i][j],f[i1][j]+shuzu[a[i]][5],f[i][j1]+shuzu[b[j]][5],f[i1][j1]+shuzu[a[i]][b[i]]
分别是s1串的最后一个字符对应一个空字符,s2串的最后一个字符对应一个空字符,s1串个s2串的最后一个字符直接对应。

显而易见的,初始化f数组就是

for(int i=1;i<=n;i++)f[i][0]=f[i-1][0]+shuzu[a[i]][5];
for(int i=1;i<=m;i++)f[0][i]=f[0][i-1]+shuzu[b[i]][5];

然后把它们拼凑起来,就完工喽!

完结撒花~