【回文自動機】bzoj3676 [Apio2014]回文串
阿新 • • 發佈:2017-06-03
typedef 最大 clas 位置 log tails while bzoj3 -1
回文自動機講解!http://blog.csdn.net/u013368721/article/details/42100363
pam上每個點代表本質不同的回文子串。len(i)代表長度,cnt(i)代表個數(要最後在fail樹上dp一遍方可)。
答案直接枚舉一遍結點,然後用len(i)*cnt(i),取最大者即可。
回文自動機是非常優越的數據結構,可惜比manacher多一個字符集的空間……
#include<cstdio> #include<cstring> #include<algorithm> using namespace std; #define MAXN 300010 #define MAXC 26 struct PAM{ int next[MAXN][MAXC];//next指針,next指針和字典樹類似,指向的串為當前串兩端加上同一個字符構成 int fail[MAXN];//fail指針,失配後跳轉到fail指針指向的節點 int cnt[MAXN]; int num[MAXN]; int len[MAXN];//len[i]表示節點i表示的回文串的長度 int S[MAXN];//存放添加的字符 int last;//指向上一個字符所在的節點,方便下一次add int n;//字符數組指針 int p;//節點指針 int newnode(int l){//新建節點 for(int i=0;i<MAXC;++i){ next[p][i]=0; } cnt[p]=0; num[p]=0; len[p]=l; return p++; } void init(){//初始化 p=0; newnode(0); newnode(-1); last=n=0; S[n]=-1;//開頭放一個字符集中沒有的字符,減少特判 fail[0]=1; } int get_fail(int x){//和KMP一樣,失配後找一個盡量最長的 while (S[n-len[x]-1]!=S[n]){ x=fail[x]; } return x ; } void add(int c){ c-=‘a‘; S[++n]=c; int cur=get_fail(last);//通過上一個回文串找這個回文串的匹配位置 if (!next[cur][c]){//如果這個回文串沒有出現過,說明出現了一個新的本質不同的回文串 int now=newnode(len[cur]+2);//新建節點 fail[now]=next[get_fail(fail[cur])][c];//和AC自動機一樣建立fail指針,以便失配後跳轉 next[cur][c]=now; num[now]=num[fail[now]]+1; } last=next[cur][c]; cnt[last]++; } void count(){//父親累加兒子的cnt,因為如果fail[v]=u,則u一定是v的子回文串! for(int i=p-1;i>=0;--i){ cnt[fail[i]]+=cnt[i]; } } }pam; typedef long long ll; char s[MAXN]; int main(){ int len; pam.init(); scanf("%s",s+1); len=strlen(s+1); for(int i=1;i<=len;++i){ pam.add(s[i]); } pam.count(); ll ans=0; for(int i=2;i<pam.p;++i){ ans=max(ans,(ll)pam.cnt[i]*pam.len[i]); } printf("%lld\n",ans); return 0; }
【回文自動機】bzoj3676 [Apio2014]回文串