8
算法分析:貌似这种问题:从一个序列操作得到另一个序列都可以用动态规划DP来解决,这个题目我也是先套用dp,发现有最优子结构,也能找到递归解,所以这个题目用动态规划来解决。
假设c[i][j]表示源串m的前j位,目的串n的前i位的最少修改时间,那么:
c[i][j]=c[i-1][j-1] 如果n[i]=m[j]
c[i][j]=min{c[i-1][j-1]+c,c[i-1][j-1]+a+b,c[i][j-1]+b,c[i-1][j]+a} 如果n[i]!=m[j]
(第一个+c,相当于把源串的j位改成目的串的i位,这样的代价就是c[i-1][j-1]+c
第二个+a+b,相当于把源串的j位删除再添加一位,代价也就是c[i-1][j-1]+a+b
第三个是把源串的j位去掉
第四个是把源串后面添加一位)
上面就是递归解,有了递归解,可以选择递归方法或者数组方法,因为递归方法有重叠子问题所以我采用dp的自底向上的方法,类似于求最长公共子序列LCS(见算法导论DP一章)
具体代码:
#include<iostream> #include<vector> #include<string> using namespace std; #define num 1001 int dp[num][num]; int main() { int a,b,c; while(cin>>a>>b>>c) { for(int i=0;i<num;i++) { dp[0][i]=i*b; dp[i][0]=i*a; }//end for string source,des; cin>>source>>des; source="0"+source;//为了舍弃0号元素 des="0"+des; //开始填写数组 for(int i=1;i<source.size();i++) { for(int j=1;j<des.size();j++) { if(source[i]==des[j]) dp[i][j]=dp[i-1][j-1]; else dp[i][j]=min(min(dp[i-1][j-1]+c,dp[i-1][j-1]+a+b),min(dp[i][j-1]+b,dp[i-1][j]+a)); }//end for j }//end for i //输出结果 cout<<dp[source.size()-1][des.size()-1]<<endl; }//end while }