2016 东 莞 特 长 生 T 3 : 字 串 距 离 2016 ~~东莞特长生~~T3 :字串距离 2016 东莞特长生 T3:字串距离
【问题描述】
设有字符串 X,我们称在 X 的头尾及中间插入任意多个空格后构成的新字符串为 X 的扩展串,如字符串 X 为”abcbcd”,则字符串“abcb□cd”,“□a□bcbcd□”和“abcb□cd□”都是 X 的扩展串,这里“□”代表空格字符。如果 A1 是字符串 A 的扩展串,B1 是字符串 B 的扩展串,A1 与 B1 具有相同的长度,那么我扪定义字符串 A1 与 B1 的距离为相应位置上的字符的距离总和,而两个非空格字符的距离定义为它们的 ASCII 码的差的绝对值,而空格字符与其他任意字符之间的距离为已知的定值 K,空格字符与空格字符的距离为 0。在字符串 A、B 的所有扩展串中,必定存在两个等长的扩展串 A1、B1,使得 A1与 B1 之间的距离达到最小,我们将这一距离定义为字符串 A、B 的距离。
请你写一个程序,求出字符串 A、B 的距离。
【数据输入】
从文件blast.in中读入数据,输入文件第一行为字符串A,第二行为字符串B。 A、B 均由小写字母组成且长度均不超过 2000。第三行为一个整数 K(1≤K≤100),表示空格与其他字符的距离。
【数据输出】
输出到文件 blast.out 中,仅一行包含一个整数,表示所求得字符串 A、B 的距离。
【输入输出样例】
blast.in
cmc
snmn
2
blast.out
10
思路
一个诡异的动态规划,设
f
[
i
]
[
j
]
f[i][j]
f[i][j]表示A的前i个字符和B中前j个字符的匹配后的最小距离。
明显,此时有三种策略:
- 什么都不加, f [ i ] [ j ] = f [ i − 1 ] [ j − 1 ] + ∣ A [ i ] − B [ j ] ∣ f[i][j]=f[i-1][j-1]+|~A[i]-B[j]~| f[i][j]=f[i−1][j−1]+∣ A[i]−B[j] ∣
- 在 A i A_i Ai后面添加一个空格, f [ i ] [ j ] = f [ i − 1 ] [ j ] + K f[i][j]=f[i-1][j]+K f[i][j]=f[i−1][j]+K
- 在 B j B_j Bj后面添加一个空格, f [ i ] [ j ] = f [ i ] [ j − 1 ] + K f[i][j]=f[i][j-1]+K f[i][j]=f[i][j−1]+K
在三者中取最小值即可
代码
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<queue>
#include<cstring>
#include<cmath>
#include<map>
#define r register
#define rep(i,x,y) for(r ll i=x;i<=y;++i)
#define per(i,x,y) for(r ll i=x;i>=y;--i)
using namespace std;
typedef long long ll;
const int V=100100;
string s,t;
ll k,dp[2021][2021],len1,len2;
int main()
{
freopen("blast.in","r",stdin);
freopen("blast.out","w",stdout);
cin>>s>>t;
len1=s.length(),len2=t.length();
scanf("%lld",&k);
rep(i,0,len1-1) //初始化:匹配全为空格
dp[i+1][0]=(i+1)*k;
rep(i,0,len2-1)
dp[0][i+1]=(i+1)*k;
rep(i,1,len1)
rep(j,1,len2)
dp[i][j]=min(dp[i-1][j]+k,min(dp[i][j-1]+k,dp[i-1][j-1]+abs(s[i-1]-t[j-1])));//减一是为了处理下标问题
cout<<dp[len1][len2];
return 0;
}