您的位置:首页 > 其它

AGTC(hdu 3356)

2015-08-15 16:39 393 查看
本题链接(单击这里)

题意:

这题是一个编辑距离题,编辑距离是指两个字串之间,由一个转成另一个所需的最少编辑操作次数(可以替换、插入、删除)。

思路:

用表格来记录状态的变化。

第一个表格是dp初始化后的状态。

dp[i][0]=dp[0][i]=i;

为什么这样初始化:假如有一个字符串什么都没有,而另一个有5个字符是不是要操作5次,11个是不是要操作11次;

也可以这样想,想要5个字符变为没有字符是不是要删除5次。

AGTAAGTAGGC
01234567891011
A1
G2
T3
C4
T5
G6
A7
C8
G9
C10
拆分一下,先假如只有一个字母的串变成11个字母的串。

再拆一下,行两个字母,列一个字母。

AG
012
A101
AGTAAGTAGGC
01234567891011
A1012345678910
然后给出DP方程:

dp[i][j]=min(dp[i-1][j]+1,dp[i][j-1]+1);

if(a[i]==b[j])

dp[i][j]=min(dp[i][j],dp[i-1][j-1]);

else

dp[i][j]=min(dp[i][j],dp[i-1][j-1]+1);

AGTAAGTAGGC
01234567891011
A1012345678910
G210123456789
T321123456678
C432112345677
T543222334567
G654333234456
A765433333456
C876544444455
G987655455445
C1098766556554
#include <iostream>
#include <cstdio>
#include <algorithm>
using namespace std;
char a[1010];
char b[1010];
int dp[1010][1010];
int main()
{
int ai,bi;
while(~scanf("%d%s",&ai,a+1))
{
scanf("%d%s",&bi,b+1);
int Max=max(ai,bi);
for(int i=0;i<=Max;i++)
dp[i][0]=dp[0][i]=i;
for(int i=1;i<=ai;i++)
for(int j=1;j<=bi;j++)
{
dp[i][j]=min(dp[i-1][j]+1,dp[i][j-1]+1);
if(a[i]==b[j])
dp[i][j]=min(dp[i][j],dp[i-1][j-1]);
else
dp[i][j]=min(dp[i][j],dp[i-1][j-1]+1);
}
printf("%d\n",dp[ai][bi]);
}
return 0;
}
内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
标签: