RMQ即求区间(i,j)的最值。通过O(nlogn)处理,O(1)给出答案。

RMQ主要是动态规划来做。dp[i][j]表示从i开始的长为2^j的区间最值。

那么可以得到dp[i][j]=max(dp[i][j-1],dp[i+(1<<(j-1))][j-1]);

dp[i][j],这个区间可以分为2段(可以重叠),那最值就是这两段的最值。

查询时要找到那个j,那j=(int)(log((y-x+1)*1.0)/log(2.0));

对于求回文 可以转变为当前的位子进行枚举 求当前的位置的后缀和当前位置的前面部分的公共长度,又前面一部分就是在后面添加的2*n-i的位置
所以只要求出height[i+1.....2*n-i]的最小值,这里就用到RMQ来做;

#include<stdio.h>
#include<math.h>
#include<string.h>
#include<stdlib.h>
//#include<Windows.h>
#define maxn 2100
#define LL long long
using namespace std;
int wa[maxn],wb[maxn],wv[maxn],WS[maxn],n;
int dp[maxn][];
int cmp(int *r,int a,int b,int l)
{return r[a]==r[b]&&r[a+l]==r[b+l];}
int min(int x,int y)
{return x<y?x:y;}
void da(int *r,int *sa,int n,int m)
{
int i,j,p,*x=wa,*y=wb,*t;
for(i=;i<m;i++) WS[i]=;
for(i=;i<n;i++) WS[x[i]=r[i]]++;
for(i=;i<m;i++) WS[i]+=WS[i-];
for(i=n-;i>=;i--) sa[--WS[x[i]]]=i;
for(j=,p=;p<n;j*=,m=p)
{
for(p=,i=n-j;i<n;i++) y[p++]=i;
for(i=;i<n;i++) if(sa[i]>=j) y[p++]=sa[i]-j;
for(i=;i<n;i++) wv[i]=x[y[i]];
for(i=;i<m;i++) WS[i]=;
for(i=;i<n;i++) WS[wv[i]]++;
for(i=;i<m;i++) WS[i]+=WS[i-];
for(i=n-;i>=;i--) sa[--WS[wv[i]]]=y[i];
for(t=x,x=y,y=t,p=,x[sa[]]=,i=;i<n;i++)
x[sa[i]]=cmp(y,sa[i-],sa[i],j)?p-:p++;
}
return;
}
int Rank[maxn],height[maxn];
void calheight(int *r,int *sa,int n)
{
int i,j,k=;
for(i=;i<=n;i++) Rank[sa[i]]=i;
for(i=;i<n;height[Rank[i++]]=k)
for(k?k--:,j=sa[Rank[i]-];r[i+k]==r[j+k];k++);
return;
}
int r[maxn],sa[maxn];
void RMQ()
{
int i,j;
memset(dp,,sizeof(dp));
for(i=;i<=*n+;i++)
dp[i][]=height[i];
for(j=;j<=;j++)
for(i=;i+(<<j)-<=*n+;i++)
{
dp[i][j]=min(dp[i][j-],dp[i+(<<(j-))][j-]);
}
}
int lcp(int left,int right)
{
int a=Rank[left];
int b=Rank[right];
if(a>b)
{
int t=a;
a=b;
b=t;
}
a++;
int k=(int)(log((b-a+)*1.0)/log(2.0));
return min(dp[a][k],dp[b-(<<k)+][k]);
}
char s[maxn];
int main()
{
int i,j;
scanf("%s",s);
n=strlen(s);
s[n]='#';
int len=n+;
for(i=n-;i>=;i--)
s[len++]=s[i];
//printf("%s\n",s);
for(i=;i<*n+;i++)
r[i]=s[i];
r[*n+]=;
da(r,sa,n*+,);
calheight(r,sa,n*+);
RMQ();
int ans=-;
int set=;
int res;
//对于求回文 可以转变为当前的位子进行枚举 求当前的位置的后缀和当前位置的前面部分的公共长度,又前面一部分就是在后面添加的2*n-i的位置
//所以只要求出height[i+1.....2*n-i]的最小值
for(i=;i<n;i++)
{
res=lcp(i,*n-i)*-;//对于奇数
if(res>ans)
{
ans=res;
set=i;
}
res=lcp(i,*n-i+)*;//对于偶数
if(res>ans)
{
ans=res;
set=i;
}
}
if(ans%)
{
for(i=set-ans/;i<=set+ans/;i++)
{
printf("%c",s[i]);
}
}
else
{
for(i=set-ans/;i<=set+ans/-;i++)
{
printf("%c",s[i]);
}
}
printf("\n");
//system("pause");
}

最新文章

  1. jq倒计时(代码)
  2. js模版引擎handlebars.js实用教程——关于HTML编码
  3. SQL中使用WITH AS提高性能,使用公用表表达式(CTE)简化嵌套SQL
  4. NPOI 读写Excel
  5. Drawable和Bitmap的区别
  6. ORA-3136报错
  7. iOS Layer CABasicAnimation
  8. C#入门经典中的SelectionFont属性为null
  9. Sping中的IOC四种注解的简单记录
  10. jQuery 入门
  11. VS编译代码未通过,常见问题。
  12. AngularJS学习之旅—AngularJS Table(十一)
  13. Codeforces.612E.Square Root of Permutation(构造)
  14. mormot支持TCP/IP
  15. Makefile shell subst $(1)
  16. C# 多线程学习系列二
  17. 基于tcpdump的Android智能移动终端数据包捕获完整解决方案
  18. Gitlab Jenkins WebHook 持续集成配置踩坑记
  19. python中键值叫唤例子
  20. ADB故障时的一些命令

热门文章

  1. 忘记sql server 2008 sa的密码的解决方案
  2. 2019-8-8-WPF-非客户区的触摸和鼠标点击响应
  3. LA2965 Jurassic Remains
  4. MySQL系列(七)--SQL优化的步骤
  5. 通过游戏学python 3.6 第一季 第八章 实例项目 猜数字游戏--核心代码--猜测次数--随机函数和屏蔽错误代码--优化代码及注释--简单账号密码登陆--账号的注册查询和密码的找回修改--锁定账号--锁定次数
  6. tr的display属性出现td的colspan无效问题
  7. 关闭浏览器或者关闭使用window.open打开的页面时添加监听事件
  8. 20190716-T1-礼物
  9. python类相关总结(持续更新)
  10. 3D hover文字特效