1072: 编辑距离

该博客介绍了一个计算两个字符串编辑距离的问题,编辑距离定义为通过删除、插入和修改字符将一个字符串转换为另一个字符串所需的最小操作次数。文章提供了动态规划的解决方案,并给出了状态转移方程。样例输入和输出展示了算法的应用。最后,作者提出了关于string和char[]在性能上的个人猜想。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

1072: 编辑距离 分享至QQ空间 去爱问答提问或回答

时间限制(普通/Java):1000MS/10000MS     内存限制:65536KByte

描述

假设字符串的基本操作仅为:删除一个字符、插入一个字符和将一个字符修改成另一个字符这三种操作。 
我们把进行了一次上述三种操作的任意一种操作称为进行了一步字符基本操作。 
下面我们定义两个字符串的编辑距离:对于两个字符串a和b,通过上述的基本操作,我们可以把a变成b或b变成a,那么字符串a变成字符串b需要的最少基本字符操作步数称为字符串a和字符串b的编辑距离。 
例如:a="ABC",b="CBCD",则a与b的编辑距离为2。 
你的任务就是:编一个快速的程序来计算任意两个字符串的编辑距离。

输入

输入包含多组测试数据。每组测试数据一行,为字符串A和字符串B。 
字符串的长度不大于1024,且全为字母。

输出

编辑距离。

样例输入

   
ABC CBCD

样例输出

   
2

提示

我相信动态规划算法能解决这题,因为我就是这样做的。^_^

题目来源

ZJGSU


解析:

首先定义这样一个函数——dp(i, j),它表示第一个字符串的长度为i的子串到第二个字符串的长度为j的子串的编辑距离。

显然可以有如下动态规划公式:

  • if i == 0 且 j == 0,dp(i, j) = 0
  • if i == 0 且 j > 0,dp(i, j) = j
  • if i > 0 且j == 0,dp(i, j) = i
  • if i ≥ 1  且 j ≥ 1 ,dp(i, j) == min{ dp(i-1, j) + 1, dp(i, j-1) + 1, dp(i-1, j-1) + f(i, j) },当第一个字符串的第i个字符不等于第二个字符串的第j个字符时,f(i, j) = 1;否则,f(i, j) = 0。

这道题的数据量挺大的,先给出我第一次提交的代码,结果超时。

#include <iostream>
#include <string>
using namespace std;
int dp[1025][1025];
int main()
{
    string a,b;
    int aLen,bLen;
    while(cin>>a>>b)
    {
        aLen = a.size();
        bLen = b.size();
        for(int i = 0 ; i <= aLen ; i++)
        {
            dp[i][0] = i;
        }
        for(int j = 0 ; j<=bLen ; j++)
        {
            dp[0][j] = j;
        }
        for(int i = 1 ; i <= aLen ; i++)
        {
            for(int j = 1 ; j <= bLen ; j++)
            {
                dp[i][j] = min(dp[i-1][j],dp[i][j-1])+1;
                dp[i][j] = min(dp[i][j],dp[i-1][j-1]+(a[i-1]==b[j-1]?0:1));
            }
        }
        cout<<dp[aLen][bLen]<<endl;
    }
    return 0;
}


然后改了一点点的小地方,就是定义变量的方式换了一下,将string换成char[],代码如下,结果就A了,用时906MS,可想数据量的庞大。

猜想string内部比char[]占用了一些时间,因为string要自动维持可变长度,这只是个人自己的一点猜想,没有任何考证,也没有去查阅相关资料,太懒!希望有懂得人留言告诉我一下。


#include <iostream>
#include <cstring>
using namespace std;
int dp[1025][1025];
char a[1025];
char b[1025];
int main()
{
    int aLen,bLen;
    while(cin>>a>>b)
    {
        aLen = strlen(a);
        bLen = strlen(b);
        for(int i = 0 ; i <= aLen ; i++)
        {
            dp[i][0] = i;
        }
        for(int j = 0 ; j<=bLen ; j++)
        {
            dp[0][j] = j;
        }
        for(int i = 1 ; i <= aLen ; i++)
        {
            for(int j = 1 ; j <= bLen ; j++)
            {
                dp[i][j] = min(dp[i-1][j],dp[i][j-1])+1;
                dp[i][j] = min(dp[i][j],dp[i-1][j-1]+(a[i-1]==b[j-1]?0:1));
            }
        }
        cout<<dp[aLen][bLen]<<endl;
    }
    return 0;
}



评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值