1. 程式人生 > >BZOJ3676 [Apio2014]迴文串

BZOJ3676 [Apio2014]迴文串

傳送門

迴文自動機板子題~

迴文自動機和ACA以及SAM都是很類似的[畢竟都是自動機嗎233]

迴文自動機的樹形結構是 fail指標構成的

用增量法 構造即可

(其實我也沒完全學懂呢T^T)

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define inf 20021225
#define ll long long
#define mxn 310000
using namespace std;

struct node{int ch[26],len,fa,s;}t[mxn];
// ch add character fa fail s times
char ch[mxn]; int lt,poi,s[mxn],n; ll ans;
void init()
{
	s[0]=-1;
	t[0].len=0; t[1].len=-1;
	t[0].fa=1; t[1].fa=1;
	poi=lt=1;
}
int id(char c){return c-'a';}
void extend(int c)
{
	int p=lt; s[++n]=c;
	while(s[n-1-t[p].len] != s[n]) p=t[p].fa;
	if(!t[p].ch[c])
	{
		int q=++poi; t[q].len=t[p].len+2;int np = t[p].fa;
		while(s[n-1-t[np].len] != s[n]) np=t[np].fa;
		t[q].fa=t[np].ch[c]; t[p].ch[c]=q; 
	}
	lt=t[p].ch[c]; t[lt].s++;
}
void calc()
{
	for(int i=poi;~i;i--)	t[t[i].fa].s += t[i].s;
	for(int i=0;i<=poi;i++)	ans=max(ans,1ll*t[i].len*t[i].s);
}
int main()
{
	scanf("%s",ch); init(); int n = strlen(ch);
	for(int i=0;i<n;i++)	extend(id(ch[i]));
	calc(); printf("%lld\n",ans);
	return 0;
}