BZOJ3676:[Apio2014]回文串(回文自动机)

很水的我刚刚打出了听说半年的回文自动机,就开始吹了

长为n的串不同的回文子串只有n个
回文自动机有两个初始状态分别表示长度为0和-1的空串
其余状态表示一个回文子串

有两个属性
Right集:表示其出现位置
len:子串的长度

三个指针
son[x]:在串的两端加一个字母x所形成的状态
fa:son的逆指针
pre:后缀中的最长回文串的状态

与sam有些差不多的性质
Right要么不交,要么包含,pre形成一个树形结构
某节点的Right集为其所有儿子Right的并集

对于后端插入字符x,同sam类似
找出表示最长后缀回文串的状态last(上次所走到的状态)
通过跳pre的方式找到一个最长的,插入x后仍是回文的后缀
将x设为其儿子
再通过跳pre的方法找到x的pre

由于回文的对称性
pre也能表示最长前缀回文串
将上一段的后缀换成前缀也通顺
故上述方法也适用于前端插入字符

对于可持久化和在trie上拓展我虽然大概知道,但也打不出来

对于本题
题面
Sam的做法
回文自动机:每个状态的len*Right集大小能贡献答案

#include <iostream>
#include <fstream>
#include <algorithm>
#include <cmath>
#include <ctime>
#include <cstdio>
#include <cstdlib>
#include <cstring>

using namespace std;
#define mmst(a, b) memset(a, b, sizeof(a))
#define mmcp(a, b) memcpy(a, b, sizeof(b))

typedef long long LL;

const int N=1001000;

int n;
int len[N],pre[N],son[N][26],r[N],cnt=2,last=1;
int s[N];
LL ans;
char cc[N];

void Insert(int x,int pos)
{
    s[pos]=x;
    while(s[pos]!=s[pos-len[last]-1])
    last=pre[last];
    int cur;
    if(!son[last][x])
    {
        son[last][x]=++cnt;
        len[cnt]=len[last]+2;
        cur=son[last][x];
        last=pre[last];

        while(last&&s[pos]!=s[pos-len[last]-1])
        last=pre[last];

        if(last)
        pre[cur]=son[last][x];
        else
        pre[cur]=2;
    }
    else
    cur=son[last][x];

    last=cur;
    r[last]++;
}

int main()
{
    pre[1]=0,pre[2]=1;
    len[1]=-1,len[2]=0;

    scanf("%s",cc);
    n=strlen(cc);
    for(int i=1;i<=n;i++)
    Insert(cc[i-1],i);

    for(int i=cnt;i>=1;i--)
    r[pre[i]]+=r[i];

    for(int i=3;i<=cnt;i++)
    ans=max(ans,(LL)len[i]*r[i]);

    cout<<ans<<endl;

    return 0;
}

这里写图片描述

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值