hdu 4681 string
2024-08-25 10:57:46
字符串DP
题意:给你三个字符串a,b,c求字符串d的长度。
字符串d满足的要求:是a和b的公共子序列,c是它的子串。
定义dp1[i][j]表示a的第 i 位与b的第 j 位之前相同的子序列长度(包括 i ,j),dp2[i][j]表示a的第 i 位与b的第 j 位之后相同的子序列长度(包括 i ,j)。
查找c的首尾字符在a 和b 的位置,求出c首字符之前a和b的公共子序列的长度以及c尾字符之后a和b的公共子序列的长度,加上c的长度即为所求。
#include<stdio.h>
#include<string.h>
#define maxn 1010
#define max(a,b) (a)>(b)?(a):(b) char a[maxn],b[maxn],c[maxn];
int dp1[maxn][maxn],dp2[maxn][maxn];
int n,m,len;
int sa[maxn][],sb[maxn][];
void LCS()
{
for(int i=;i<=n;i++)
for(int j=;j<=m;j++)
{
if(a[i]==b[j])
dp1[i][j]=dp1[i-][j-]+;
else
dp1[i][j]=max(dp1[i-][j],dp1[i][j-]);
}
for(int i=n;i>;i--)
for(int j=m;j>;j--)
{
if(a[i]==b[j])
dp2[i][j]=dp2[i+][j+]+;
else
dp2[i][j]=max(dp2[i+][j],dp2[i][j+]);
}
}
void Deal()
{
int k,j,t1,t2;
t1=t2=;
for(int i=;i<=n;i++)
if(a[i]==c[])
{
for(j=i+,k=;j<=n;j++)
{
if(a[j]==c[k+]) k++;
if(k==len) break;
}
if(k==len) sa[t1][]=i,sa[t1++][]=j;
}
for(int i=;i<=m;i++)
if(b[i]==c[])
{
for(j=i+,k=;j<=m;j++)
{
if(b[j]==c[k+]) k++;
if(k==len) break;
}
if(k==len) sb[t2][]=i,sb[t2++][]=j;
}
int ans=;
for(int i=;i<t1;i++)
for(int j=;j<t2;j++)
ans=max(ans,dp1[sa[i][]-][sb[j][]-]+dp2[sa[i][]+][sb[j][]+]);
printf("%d\n",len+ans);
}
int main()
{
int T,cas=;
scanf("%d",&T);
while(T--)
{
scanf("%s%s%s",a+,b+,c+);
n=strlen(a+);
m=strlen(b+);
len=strlen(c+);
memset(dp1,,sizeof(dp1));
memset(dp2,,sizeof(dp2));
printf("Case #%d: ",cas++);
LCS();
Deal();
}
return ;
}
最新文章
- iis7 压缩js文件和启用gzip压缩
- 【Java每日一题】20161125
- 二、快速起步(Mysql镜像)
- WP&;Win10仿微信消息框代码分享
- Linux的服务器初始优化脚本。
- 用c#写一个json的万能解析器
- 【C# 反射泛型】
- Tomcat无法部署项目
- Pattern()和Matcher() 用法
- Unity3D 画线插件 Vectrosity 画一个一直循环的正弦函数曲线
- Scala学习文档-样本类与模式匹配(match,case,Option)
- 大白菜U盘启动制作工具装机维护版V5.0–大白菜U盘下载中心
- Django配置session
- c语言_FILE结构体解释及相关操作
- Struts2中的值栈
- mybatis(3)---传参数的方法
- XSS 漏洞原理及防御方法
- falcon常用参数解析
- 20190319xlVBA_根据考勤数据统计缺勤缺考数据
- ElasticSearch集群介绍二