BZOJ 3676: [Apio2014]回文串

昨天BC打得我好伤心QAQ

T1模数为毛是(1e8)+7，坑了我一个小时

T2 T4果断不会

T3看出来是回文树了

发现好像没学过哎

只好现学（抄）了一下

然后由于不是很懂，于是写挂了

今天学了一下顺便复习了马拉车

感觉还是很简单的

（但是依然不会证复杂度，感觉有点慢啊）

于是把这题水了一下

大概是回文树+拓扑排序吧

#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
using namespace std;
const int N=300000+5;
typedef long long ll;
struct Node{
	int len,suflink;
	int ch[26];
	ll sum;
	void clear(){
		memset(ch,0,sizeof(ch));
		len=suflink=0;sum=0;
	}
}tr[N];
int suf,num;
void init(){
	suf=num=2;
	tr[1].clear();tr[1].len=-1;tr[1].suflink=1;
	tr[2].clear();tr[2].len=0;tr[2].suflink=1;
}
char s[N];
bool AddLetter(int pos){
	int cur=suf,curlen=0;
	int c=s[pos]-'a';
	while(true){
		curlen=tr[cur].len;
		if(pos-1-curlen>=0&&s[pos-1-curlen]==s[pos])break;
		cur=tr[cur].suflink;
	}
	if(tr[cur].ch[c]){
		suf=tr[cur].ch[c];
		tr[suf].sum++;
		return false;
	}
	suf=++num;
	tr[suf].clear();
	tr[suf].len=tr[cur].len+2;
	tr[suf].sum=1;
	tr[cur].ch[c]=suf;
	if(tr[suf].len==1){
		tr[suf].suflink=2;
		return true;
	}
	while(true){
		cur=tr[cur].suflink;
		curlen=tr[cur].len;
		if(pos-1-curlen>=0&&s[pos-1-curlen]==s[pos]){
			tr[suf].suflink=tr[cur].ch[c];
			return true;
		}
	}
}
int deg[N];
ll toposort(){
	for(int i=1;i<=num;i++)
	deg[tr[i].suflink]++;
	queue<int>q;
	for(int i=1;i<=num;i++)
	if(!deg[i])q.push(i);
	while(!q.empty()){
		int u=q.front();q.pop();
		int v=tr[u].suflink;
		tr[v].sum+=tr[u].sum;
		deg[v]--;
		if(!deg[v])q.push(v);
	}
	ll ans=0;
	for(int i=1;i<=num;i++)
	ans=max(ans,tr[i].len*tr[i].sum);
	return ans;
}
int main(){
	//freopen("a.in","r",stdin);
	scanf("%s",s);
	int n=strlen(s);
	init();
	for(int i=0;i<n;i++)
	AddLetter(i);
	printf("%lld\n",toposort());
	return 0;
}