第八届蓝桥杯 最大公共子串(动态规划)

2024-09-05 21:38

本文主要是介绍第八届蓝桥杯 最大公共子串(动态规划),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

标题:最大公共子串

最大公共子串长度问题就是:

求两个串的所有子串中能够匹配上的最大长度是多少。

比如:"abcdkkk" 和 "baabcdadabc",

可以找到的最长的公共子串是"abcd",所以最大公共子串长度为4。

下面的程序是采用矩阵法进行求解的,这对串的规模不大的情况还是比较有效的解法。

请分析该解法的思路,并补全划线部分缺失的代码。

#include <stdio.h>
#include <string.h>


#define N 256
int f(const char* s1, const char* s2)
{
    int a[N][N];
    int len1 = strlen(s1);
    int len2 = strlen(s2);
    int i,j;


    memset(a,0,sizeof(int)*N*N);
    int max = 0;
    for(i=1; i<=len1; i++){
        for(j=1; j<=len2; j++){
            if(s1[i-1]==s2[j-1]) {
                a[i][j] = __________________________;  //填空
                if(a[i][j] > max) max = a[i][j];
            }
        }
    }


    return max;
}


int main()
{
    printf("%d\n", f("abcdkkk", "baabcdadabc"));
    return 0;
}


注意:只提交缺少的代码,不要提交已有的代码和符号。也不要提交说明性文字。

观察可以看出该题的代码是动态规划的类型;


子问题 :第一个字符串中以i结尾的字符串和第二个字符串中以j结尾的字符串的最长公共子串的长度;

状态:dp[i][j]表示第一个字符串中以i结尾的字符串和第二个字符串中以j结尾的字符串的最长公共子串的长度;

状态转移方程:考虑状态方程要结合子问题来考虑。如果str1[i]==str2[j],这时dp[i][j]就可以从dp[i-1][j-1]转移过来。

如果str1[i],str2[j]不相等的话,状态就不转移,或者dp[i][j]为0,二者是等价的。因为我们要考虑到后来的状态可能会从当前状态转移过来,例如当str1[i]!=str2[j],str[i+1]=str[j+1]时,dp[i][j]=0,dp[i+1][j+1]=dp[i][j]+1。符合连续子串的定义。、

正是因为这个性质,所以最长公共子串的最大值不一定是在dp[len1][len2]中,所以我们要用一个变量来存这个最大值。

答案为a[i-1][j-1]+1;


——————————————————————————————————————————————————

接下来我们考虑另一个问题,如何输出最长公共子串。

因为最长公共子串的不唯一,我们可以用一个数组CoordinateMax来存最长公共子串的末尾的字符在第一个字符串和第二个字符串的下标。在动态规划过程中,如果max发生了更新,我们就更新CoordinateMax最后一个存有实际内容元素的值。

当max == dp[i][j],说明在当前子问题中存在多个长度相同的最长公共子串,我们就在CoordinateMax中加入该点信息

代码如下

#include<stdio.h>
#include<string.h>const int INF = 1005;int CoordinateMax[100][2], index = 0, max = -1;int LongestCommonSubstring(char* str1, char* str2, int len1, int len2);int dp[INF][INF] = {0};
/*
abcdkkkkedydfghfs
baabcdadabckkkksgsdghdghfsjgh
*/
int main()
{char str1[INF],  str2[INF];scanf("%s",str1);scanf("%s",str2);int len1 = strlen(str1), len2 = strlen(str2);LongestCommonSubstring(str1,str2,len1,len2) ;printf("%d\n",max);return 0;
} 
int LongestCommonSubstring(char* str1, char* str2, int len1, int len2)
{for(int i = 1; i <= len1; i++){for(int j = 1; j <= len2; j++){if(str1[i-1] == str2[j-1]){dp[i][j] = dp[i-1][j-1] + 1;if(max < dp[i][j]){max = dp[i][j];CoordinateMax[index][0] = i, CoordinateMax[index][1] = j;			//当max更新时,CoordiateMax也需要更新(注意不是添加一个元素) }else if(max == dp[i][j]){index++;CoordinateMax[index][0] = i, CoordinateMax[index][1] = j;			//当max与dp[i][j]相等时,说明在当前子问题中存在多个最长公共子串,//我们需要在CoordinateMax中添加该点的信息 }}}} for(int i = index; i >= 0; i--){if(dp[CoordinateMax[i][0]][CoordinateMax[i][1]] == max){int end1 = CoordinateMax[i][0]; for(int i = end1-max; i < end1; i++)printf("%c",str1[i]);printf("\n");}elsebreak;}return max;}




这篇关于第八届蓝桥杯 最大公共子串(动态规划)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1140142

相关文章

第10章 中断和动态时钟显示

第10章 中断和动态时钟显示 从本章开始,按照书籍的划分,第10章开始就进入保护模式(Protected Mode)部分了,感觉从这里开始难度突然就增加了。 书中介绍了为什么有中断(Interrupt)的设计,中断的几种方式:外部硬件中断、内部中断和软中断。通过中断做了一个会走的时钟和屏幕上输入字符的程序。 我自己理解中断的一些作用: 为了更好的利用处理器的性能。协同快速和慢速设备一起工作

动态规划---打家劫舍

题目: 你是一个专业的小偷,计划偷窃沿街的房屋。每间房内都藏有一定的现金,影响你偷窃的唯一制约因素就是相邻的房屋装有相互连通的防盗系统,如果两间相邻的房屋在同一晚上被小偷闯入,系统会自动报警。 给定一个代表每个房屋存放金额的非负整数数组,计算你 不触动警报装置的情况下 ,一夜之内能够偷窃到的最高金额。 思路: 动态规划五部曲: 1.确定dp数组及含义 dp数组是一维数组,dp[i]代表

poj2406(连续重复子串)

题意:判断串s是不是str^n,求str的最大长度。 解题思路:kmp可解,后缀数组的倍增算法超时。next[i]表示在第i位匹配失败后,自动跳转到next[i],所以1到next[n]这个串 等于 n-next[n]+1到n这个串。 代码如下; #include<iostream>#include<algorithm>#include<stdio.h>#include<math.

poj3261(可重复k次的最长子串)

题意:可重复k次的最长子串 解题思路:求所有区间[x,x+k-1]中的最小值的最大值。求sa时间复杂度Nlog(N),求最值时间复杂度N*N,但实际复杂度很低。题目数据也比较水,不然估计过不了。 代码入下: #include<iostream>#include<algorithm>#include<stdio.h>#include<math.h>#include<cstring

spoj705( 求不相同的子串个数)

题意:求串s的不同子串的个数 解题思路:任何子串都是某个后缀的前缀,对n个后缀排序,求某个后缀的前缀的个数,减去height[i](第i个后缀与第i-1 个后缀有相同的height[i]个前缀)。 代码如下: #include<iostream>#include<algorithm>#include<stdio.h>#include<math.h>#include<cstrin

poj1330(LCA最近公共祖先)

题意:求最近公共祖先 思路:之前学习了树链剖分,然后我就用树链剖分的一小部分知识就可以解这个题目了,记录每个结点的fa和depth。然后查找时,每次将depth大的结点往上走直到x = y。 代码如下: #include<iostream>#include<algorithm>#include<stdio.h>#include<math.h>#include<cstring>

软考系统规划与管理师考试证书含金量高吗?

2024年软考系统规划与管理师考试报名时间节点: 报名时间:2024年上半年软考将于3月中旬陆续开始报名 考试时间:上半年5月25日到28日,下半年11月9日到12日 分数线:所有科目成绩均须达到45分以上(包括45分)方可通过考试 成绩查询:可在“中国计算机技术职业资格网”上查询软考成绩 出成绩时间:预计在11月左右 证书领取时间:一般在考试成绩公布后3~4个月,各地领取时间有所不同

poj 3723 kruscal,反边取最大生成树。

题意: 需要征募女兵N人,男兵M人。 每征募一个人需要花费10000美元,但是如果已经招募的人中有一些关系亲密的人,那么可以少花一些钱。 给出若干的男女之间的1~9999之间的亲密关系度,征募某个人的费用是10000 - (已经征募的人中和自己的亲密度的最大值)。 要求通过适当的招募顺序使得征募所有人的费用最小。 解析: 先设想无向图,在征募某个人a时,如果使用了a和b之间的关系

poj 2976 分数规划二分贪心(部分对总体的贡献度) poj 3111

poj 2976: 题意: 在n场考试中,每场考试共有b题,答对的题目有a题。 允许去掉k场考试,求能达到的最高正确率是多少。 解析: 假设已知准确率为x,则每场考试对于准确率的贡献值为: a - b * x,将贡献值大的排序排在前面舍弃掉后k个。 然后二分x就行了。 代码: #include <iostream>#include <cstdio>#incl

poj 3258 二分最小值最大

题意: 有一些石头排成一条线,第一个和最后一个不能去掉。 其余的共可以去掉m块,要使去掉后石头间距的最小值最大。 解析: 二分石头,最小值最大。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <c