Description
考虑一个只包含小写拉丁字母的字符串s。我们定义s的一个子串t的“出
现值”为t在s中的出现次数乘以t的长度。请你求出s的所有回文子串中的最
大出现值。
Input
输入只有一行,为一个只包含小写字母(a -z)的非空字符串s。
Output
输出一个整数,为逝查回文子串的最大出现值。
Sample Input
【样例输入l】
abacaba
【样例输入2]
www
Sample Output
【样例输出l】
7
【样例输出2]
4
HINT
一个串是回文的,当且仅当它从左到右读和从右到左读完全一样。
在第一个样例中,回文子串有7个:a,b,c,aba,aca,bacab,abacaba,其中:
● a出现4次,其出现值为4:1:1=4
● b出现2次,其出现值为2:1:1=2
● c出现1次,其出现值为l:1:l=l
● aba出现2次,其出现值为2:1:3=6
● aca出现1次,其出现值为1=1:3=3
●bacab出现1次,其出现值为1:1:5=5
● abacaba出现1次,其出现值为1:1:7=7
故最大回文子串出现值为7。
【数据规模与评分】
数据满足1≤字符串长度≤300000。
Source
为了做这题特地学了Manacher…
用的是黄学长的做法所以我就不再赘述了..
黄学长题解传送门
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define MAXN 300010
#define LL long long
using namespace std;
int n;
LL ans;
int p[MAXN],c[MAXN],tmp[MAXN<<1],Pow[20]={1};
char ch[MAXN];
struct sam
{
int last,cnt,p,q,np,nq;
int a[MAXN<<1][26],fa[MAXN<<1],len[MAXN<<1],right[MAXN<<1],pos[MAXN<<1],deep[MAXN<<1];
int Fa[MAXN<<1][20];
sam() {last=++cnt;}
void insert(int c,int id)
{
p=last;np=last=++cnt;len[np]=len[p]+1;right[np]=1;pos[id]=last;
while (!a[p][c]&&p) a[p][c]=np,p=fa[p];
if (!p) fa[np]=1;
else
{
q=a[p][c];
if (len[q]==len[p]+1) fa[np]=q;
else
{
nq=++cnt;len[nq]=len[p]+1;memcpy(a[nq],a[q],sizeof(a[q]));
fa[nq]=fa[q];fa[q]=fa[np]=nq;
while (a[p][c]==q) a[p][c]=nq,p=fa[p];
}
}
}
void init()
{
for (int i=1;i<=cnt;i++) c[len[i]]++;
for (int i=1;i<=n;i++) c[i]+=c[i-1];
for (int i=cnt;i;i--) tmp[c[len[i]]--]=i;
for (int i=cnt;i;i--) right[fa[tmp[i]]]+=right[tmp[i]];
for (int i=1;i<=cnt;i++)
{
int x=tmp[i];
deep[x]=deep[fa[x]]+1;Fa[x][0]=fa[x];
for (int j=1;Pow[j]<=deep[x];j++) Fa[x][j]=Fa[Fa[x][j-1]][j-1];
}
}
void update(int l,int r)
{
int mid=pos[r];
for (int i=19;i>=0;i--)
{
int x=Fa[mid][i];
if (len[x]>=r-l+1) mid=x;
}
ans=max(ans,(LL)(right[mid])*(r-l+1));
}
}sam;
void manacher()
{
int maxn=0,id=0;
for (int i=1;i<=n;i++)
{
if (maxn>i) p[i]=min(p[2*id-i-1],maxn-i);
else p[i]=0;
while (ch[i+p[i]+1]==ch[i-p[i]]) ++p[i],sam.update(i-p[i]+1,i+p[i]);
if (p[i]+i>maxn) maxn=p[i]+i,id=i;
}
maxn=0;
for (int i=1;i<=n;i++)
{
if (maxn>i) p[i]=min(p[2*id-i],maxn-i-1);
else p[i]=1,sam.update(i-p[i]+1,i+p[i]-1);
while (ch[i+p[i]]==ch[i-p[i]]) ++p[i],sam.update(i-p[i]+1,i+p[i]-1);
if (p[i]+i>maxn) maxn=p[i]+i,id=i;
}
}
int main()
{
for (int i=1;i<=19;i++) Pow[i]=Pow[i-1]<<1;
scanf("%s",ch+1);n=strlen(ch+1);
for (int i=1;i<=n;i++) sam.insert(ch[i]-'a',i);
sam.init();ch[0]='+';ch[n+1]='*';manacher();
cout<<ans<<endl;
}