[SDOI2006]最短距离题解

                                      [SDOI2006]最短距离题解

                                                                                           ——HM

题目描述

一种EDIT字母编辑器，它的功能是可以通过不同的变换操作可以把一个源串X [l..m]变换为新的目标串y[1..n]。EDIT提供的变换操作有：

源串中的单个字符可被删除(delete)；

被替换 (replace)；

被复制到目标串中去(copy)；

字符也可被插入(insert)；

源串中的两个相邻字符可进行交换并复制到目标串中去(twiddle)；

在完成其它所有操作之后，源串中余下的全部后缀就可用删至行末的操作删除(kill)。

例如，将源"algorithm"转换成目标串"altruistic"的一种方法是采取下面的操作序列:

要达到这个结果还可能有其它一些操作序列。

操作delete,replace，copy，insert，twiddle和kill中每一个都有一个相联系的代价cost。例如

cost(delete)=3;
cost(replace)=6;
cost(copy)=5;
cost(insert)=4;
cost(twiddle)=4;
cost(kill)=被删除的串长*cost(delete)-1;

一个给定的操作序列的代价为序列中各操作代价之和。 例如上述操作序列的代价为

3*cost(copy)+2*cost(replace)+cost(delete)+3*cost(insert) + cost(twiddle) +cost(kill)

=3*5+2*6+3+3*4+4+1*3-1=48

编程任务：

给定两个序列x[1..m],y[1..n]和一些操作代价集合，X到Y的最短距离为将X转化为Y的最小的转换序列的代价。请给出一个算法来找出x[1..m]至y[1..n]的最短距离。

输入格式

第一行：源序列x[1..m]。（m<200）

第二行：目标序列y[1..n]。(n<200)

第三行：5个正整数（<100）：分别是：delete 、replace 、copy、 insert、 twiddle的代价。

输出格式

X到Y的最短距离（最小代价和）。

输入输出样例

输入 #1

algorithm
altruistic
3 6 5 4 4

输出 #1

48

题解

    一道挺考验码力的题。

    看完题目应该会发现这是一道dp题，但dp方程比较复杂，因为状态比较多。

    设我们有dp[i][j]表示初始串操作到第i位，目标串完成到第j位，则各个操作的状态转移方程为：

delete:dp[i][j]=min(dp[i][j],dp[i-1][j]+cost[1])

replace:dp[i][j]=min(dp[i][j],dp[i-1][j-1]+cost[2])

copy：dp[i][j]=min(dp[i][j],dp[i-1][j-1]+cost[3])

insert:dp[i][j]=min(dp[i][j],dp[i][j-1]+cost[4])

twiddle:dp[i][j]=min(dp[i][j],dp[i-2][j-2]+cost[5])

kill:dp[Len1][Len2]=min(dp[Len1][Len2],dp[i][Len2]+cost[1]*(Len1-i)-1)

    解释一下：delete操作中将初始串第i位的前一位即删除一个字符时加上删除代价便是状态，将其与当前状态比较即可。前5个操作都是如此，可以理解一下，应该比较简单吧。

    最后一个操作也很好理解，因为kill操作优于delete操作(最后-1)，当目标串已经完成时进行枚举，按照题目要求进行操作即可。

    大约就这样了，注意每个情况的条件与特殊情况即可。

附代码：

#include <bits/stdc++.h> 
using namespace std;
const int SIZE=205;
const int INF=0x3f3f3f3f;
#define ll long long

char s1[SIZE],s2[SIZE];
int dp[SIZE][SIZE],cost[10];
int Len1,Len2;

int main()
{
//	cin>>s1>>s2;
	scanf("%s%s",s1+1,s2+1);
	for (int i=1;i<=5;i++)
		scanf("%d",&cost[i]);
	
	Len1=strlen(s1+1);
	Len2=strlen(s2+1);
	
	memset(dp,INF,sizeof(dp));
	
	if (Len1!=0 && Len2==0){
		printf("%d",Len1*cost[1]);
		return 0;
	}
	
	if (Len1==0 && Len2!=0){
		printf("%d",Len2*cost[4]);
		return 0;
	}
	dp[0][0]=0;
	for (int i=1;i<=Len1;i++) dp[i][0]=i*cost[1];
	for (int i=1;i<=Len2;i++) dp[0][i]=i*cost[4];
	
	for (int i=1;i<=Len1;i++)
		for (int j=1;j<=Len2;j++){
			if (s1[i]==s2[j])
				dp[i][j]=min(dp[i][j],dp[i-1][j-1]+cost[3]);
			dp[i][j]=min(dp[i][j],min(dp[i-1][j-1]+cost[2],min(dp[i-1][j]+cost[1],dp[i][j-1]+cost[4])));
			
			if (i==1 || j==1) continue;
			if (s1[i-1]==s2[j] && s1[i]==s2[j-1])
				dp[i][j]=min(dp[i][j],dp[i-2][j-2]+cost[5]);
		}
		
	for (int i=1;i<Len1;i++)
		dp[Len1][Len2]=min(dp[Len1][Len2],(Len1-i)*cost[1]+dp[i][Len2]-1);
		
	printf("%d",dp[Len1][Len2]);
		
	return 0;
}

吐槽：貌似这题与状态压缩没有太大关系。

谢谢观看！