注册 登录  
 加关注
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

fudq's AC Road

何以解忧,唯有AC!

 
 
 

日志

 
 

hdu 4681 String  

2013-08-16 16:05:35|  分类: ACM-hdu |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |
http://acm.hdu.edu.cn/showproblem.php?pid=4681
给出a,b,c三个字符串,找一个最长的字符串d,满足d是a和b的子序列,c是d的字串。
找到c在a和b中所有符合条件的且是最短的区间,枚举区间,d的长度即是a和b的前半段LCS+lenc+a和b的后半段LCS。
防超时,可以先预处理a和b从左到右的LCS和从又到左的LCS

/*
* pro.cpp
*
* Created on: 2013-08-16
* Author: fudq
*/
#include <functional>
#include <algorithm>
#include <iostream>
#include <fstream>
#include <sstream>
#include <iomanip>
#include <numeric>
#include <cstring>
#include <cassert>
#include <cstdio>
#include <string>
#include <vector>
#include <bitset>
#include <queue>
#include <stack>
#include <cmath>
#include <ctime>
#include <list>
#include <set>
#include <map>
using namespace std;
//#pragma comment(linker,"/STACK:102400000,102400000")

#define FOR(i,a) for((i)=0;i<(a);(i)++)
#define MEM(a) (memset((a),0,sizeof(a)))
#define LL __int64

const int N=1010;
const int M=32010;
const int MOD=1000000007ll;
const int INF=0x7fffffff;
const int dir[4][2]={{-1,0},{1,0},{0,-1},{0,1}};
const double eps=1e-7;
const double PI=acos(-1.0);

inline int sign(double x){return (x>eps)-(x<-eps);}
template<class T> T gcd(T a,T b){return b?gcd(b,a%b):a;}
template<class T> T lcm(T a,T b){return a/gcd(a,b)*b;}
template<class T> inline T lcm(T a,T b,T d){return a/d*b;}
template<class T> inline T Min(T a,T b){return a<b?a:b;}
template<class T> inline T Max(T a,T b){return a>b?a:b;}
template<class T> inline T Min(T a,T b,T c){return min(min(a, b),c);}
template<class T> inline T Max(T a,T b,T c){return max(max(a, b),c);}
template<class T> inline T Min(T a,T b,T c,T d){return min(min(a, b),min(c,d));}
template<class T> inline T Max(T a,T b,T c,T d){return max(max(a, b),max(c,d));}
/*************************/

struct Point{
int x,y;
}pa[N],pb[N];


char a[N],b[N],c[N];
int ka,kb,lena,lenb,lenc,f1[N][N],f2[N][N];

void LCS()
{
int i,j;
MEM(f1);MEM(f2);
for(i=1;i<=lena;i++)
for(j=1;j<=lenb;j++)
{
if(a[i-1]==b[j-1])
f1[i][j]=f1[i-1][j-1]+1;
else
f1[i][j]=max(f1[i-1][j],f1[i][j-1]);
}
for(i=lena;i>=1;i--)
for(j=lenb;j>=1;j--)
{
if(a[i-1]==b[j-1])
f2[i][j]=f2[i+1][j+1]+1;
else
f2[i][j]=max(f2[i+1][j],f2[i][j+1]);
}
}

void Find(int &k,char a[],int len,Point p[])
{
int t,i,j;
k=0;
for(i=0;i<len;i++)
if(a[i] == c[0])
{
t=1;
for(j=i+1;j<len;j++)
{
if(a[j] == c[t])
t++;
if(t == lenc)
break;
}
if(t == lenc)
{
p[k].x=i;
p[k++].y=j;
}
}
}


void init()
{
lena=(int)strlen(a);
lenb=(int)strlen(b);
lenc=(int)strlen(c);
Find(ka,a,lena,pa);
Find(kb,b,lenb,pb);
LCS();
}

void solve()
{
init();
int ans=-1;
for(int i=0;i<ka;i++)
for(int j=0;j<kb;j++)
{
ans=max(ans,f1[pa[i].x][pb[j].x]+f2[pa[i].y+2][pb[j].y+2]);
}
printf("%d\n",ans+lenc);
}

int main()
{
#ifndef ONLINE_JUDGE
freopen("testin.txt","r",stdin);
//freopen("testout.txt","w",stdout);
#endif
int t,cas=1;
scanf("%d",&t);
while(t--)
{
scanf("%s%s%s",a,b,c);
printf("Case #%d: ",cas++);
solve();
}
return 0;
}



  评论这张
 
阅读(95)| 评论(0)
推荐 转载

历史上的今天

在LOFTER的更多文章

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2017