https://blog.csdn.net/u013368721/article/details/42100363
看了这一篇博客大概学习了一下,原理细节其实还不是特别明白,但至少是能基本使用了.
然后做了博客推荐的几道题目,算是学学回文自动机的应用吧.
首先最基本就是要知道回文自动机上各个数组记录的东西,然后就能根据这些信息完成一些题目要求的统计吧.
Palindromes and Super Abilities URAL - 1960
就是找前缀的本质不同回文串种类,回文自动机上除了一开始的两个节点,有几个节点就说明有几个不同回文串,所以一边插入一边查询有多少节点就能输出了.
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
typedef pair<double,double> pdd;
typedef pair<ll,ll> pll;
#define fi first
#define se second
#define mp make_pair
#define pb push_back
const int maxn=1e5+15;
const int maxe=maxn*2;
const int maxm=26;
const int mod=1e9+7;
const int inf=0x3f3f3f3f;
int n;
struct palindromic_tree{
int next[maxn][maxm];
int fail[maxn];
int cnt[maxn];
int num[maxn];
int len[maxn];
int S[maxn];
int last;
int n,p;
int newnode(int l){
for(int i=0;i<maxm;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=0;
n=0;
S[n]=-1;
fail[0]=1;
}
int get_fail(int x){
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];
next[cur][c]=now;
num[now]=num[fail[now]]+1;
}
last=next[cur][c];
cnt[last]++;
}
void count(){
for(int i=p-1;i>=0;i--)cnt[fail[i]]+=cnt[i];
}
}pt;
char str[maxn];
int main(){
scanf("%s",str);
n=strlen(str);
pt.init();
for(int i=0;i<n;i++){
pt.add(str[i]);
printf("%d",pt.p-2);
if(i<n-1)printf(" ");
}
return 0;
}
最长双回文串 HYSBZ - 2565
就是求最长的双回文串,回文自动机每次插入字符,就可以得出以当前字符为右端点的最长回文串长度,所以我们正反各做一遍就能简单得到答案了.
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
typedef pair<double,double> pdd;
typedef pair<ll,ll> pll;
#define fi first
#define se second
#define mp make_pair
#define pb push_back
const int maxn=1e5+15;
const int maxe=maxn*2;
const int maxm=26;
const int mod=1e9+7;
const int inf=0x3f3f3f3f;
int n;
struct palindromic_tree{
int next[maxn][maxm];
int fail[maxn];
int cnt[maxn];
int num[maxn];
int len[maxn];
int S[maxn];
int last;
int n,p;
int newnode(int l){
for(int i=0;i<maxm;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=0;
n=0;
S[n]=-1;
fail[0]=1;
}
int get_fail(int x){
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];
next[cur][c]=now;
num[now]=num[fail[now]]+1;
}
last=next[cur][c];
cnt[last]++;
}
void count(){
for(int i=p-1;i>=0;i--)cnt[fail[i]]+=cnt[i];
}
}pt;
char str[maxn];
int anslen[maxn];
int main(){
scanf("%s",str);
n=strlen(str);
pt.init();
for(int i=0;i<n;i++){
pt.add(str[i]);
anslen[i]=pt.len[pt.last];
}
int ans=0;
pt.init();
pt.add(str[n-1]);
for(int i=n-2;i>=1;i--){
pt.add(str[i]);
ans=max(ans,anslen[i-1]+pt.len[pt.last]);
}
printf("%d\n",ans);
return 0;
}
拉拉队排练 HYSBZ - 2160
要将回文串按照长度从大到小选择k个出来,并求出长度乘积.
k在极限数据下很大,也可以知道应该可以构造出回文串特别多的数据.不过肯定是大量重复的回文串,只需要按照长度统计回文传的数量,然后倒着遍历一边,用快速幂就能算答案了.
按照长度统计只需要先把回文自动机的cnt数组处理出来,然后遍历回文自动机的节点根据len数组统计就行了.
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
typedef pair<double,double> pdd;
typedef pair<ll,ll> pll;
#define fi first
#define se second
#define mp make_pair
#define pb push_back
const int maxn=1e6+15;
const int maxe=maxn*2;
const int maxm=26;
const int mod=19930726;
const int inf=0x3f3f3f3f;
int n;
ll k;
struct palindromic_tree{
int next[maxn][maxm];
int fail[maxn];
ll cnt[maxn];
int num[maxn];
int len[maxn];
int S[maxn];
int last;
int n,p;
int newnode(int l){
for(int i=0;i<maxm;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=0;
n=0;
S[n]=-1;
fail[0]=1;
}
int get_fail(int x){
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];
next[cur][c]=now;
num[now]=num[fail[now]]+1;
}
last=next[cur][c];
cnt[last]++;
}
void count(){
for(int i=p-1;i>=0;i--)cnt[fail[i]]+=cnt[i];
}
}pt;
ll pow_mod(ll a,ll b){
ll ans=1;
while(b){
if(b&1)ans=ans*a%mod;
a=a*a%mod;
b/=2;
}
return ans;
}
char str[maxn];
ll lencount[maxn];
int main(){
scanf("%d%lld",&n,&k);
scanf("%s",str);
pt.init();
for(int i=0;i<n;i++){
pt.add(str[i]);
}
pt.count();
for(int i=2;i<pt.p;i++){
lencount[pt.len[i]]+=pt.cnt[i];
}
ll cc=0;
ll ans=1;
for(int i=n;i>0;i--){
if(!(i&1))continue;
if(lencount[i]==0)continue;
if(cc+lencount[i]>=k){
ll kk=k-cc;
ans=ans*pow_mod(i,kk)%mod;
cc=k;
break;
}
ans=ans*pow_mod(i,lencount[i])%mod;
cc+=lencount[i];
}
if(cc<k)printf("-1\n");
else printf("%lld\n",ans);
return 0;
}
Palisection CodeForces - 17E
求相交的回文串数量,并不好求,转化成求不相交的情况数再用总数减去.
这就需要知道每一个位置为右端点的回文串数量,一开始怎么想也不知道怎么统计,后来发现有个num数组已经统计好了这个东西.
所以只需要正反做两遍就能比较轻松的统计出来了.
#include <bits/stdc++.h>
using namespace std;
#define fi first
#define se second
#define pii pair<int,int>
#define pll pair<long long ,long long>
#define mp make_pair
typedef long long ll;
#define pb push_back
const int maxn=2*1e6+15;
const int maxm=30;
const int maxe=4000015;
const int mod=51123987;
const int inf=0x3f3f3f3f;
struct palindromic_tree{
//int next[maxn][maxm];
//vector<int> next[maxn];
vector<pii> next[maxn];
int fail[maxn];
int cnt[maxn];
int num[maxn];
int len[maxn];
int S[maxn];
int last;
int n,p;
int newnode(int l){
next[p].clear();
//for(int i=0;i<maxm;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=0;
n=0;
S[n]=-1;
fail[0]=1;
}
int get_fail(int x){
while(S[n-len[x]-1]!=S[n])x=fail[x];
return x;
}
int find_next(int a,int b){
for(int i=0;i<next[a].size();i++){
if(next[a][i].fi==b)return next[a][i].se;
}
return 0;
}
void add(int c){
c-='a';
S[++n]=c;
int cur=get_fail(last);
if(!find_next(cur,c)){
int now=newnode(len[cur]+2);
fail[now]=find_next(get_fail(fail[cur]),c);
//fail[now]=next[get_fail(fail[cur])][c];
//next[cur][c]=now;
next[cur].pb(mp(c,now));
num[now]=num[fail[now]]+1;
}
last=find_next(cur,c);
cnt[last]++;
}
void count(){
for(int i=p-1;i>=0;--i)cnt[fail[i]]+=cnt[i];
}
}st;
int n;
char str[maxn];
ll pre[maxn];
int main(){
scanf("%d",&n);
scanf("%s",str+1);
st.init();
for(int i=1;i<=n;i++){
st.add(str[i]);
pre[i]=(st.num[st.last]+pre[i-1])%mod;
}
ll ans=(pre[n]*(pre[n]-1)/2)%mod;
st.init();
ll ans1=0;
for(int i=n;i>=1;i--){
st.add(str[i]);
ans1=(ans1+st.num[st.last]*pre[i-1]%mod)%mod;
}
ans=(ans-ans1+mod)%mod;
printf("%lld\n",ans);
return 0;
}
The Problem to Slow Down You Gym - 100548G
找两个字符串的公共子回文串.
一开始还想着使用同一个回文自动机,先插入字符串a,然后插入’#’,然后插入字符串b,就能怎么统计来着.
感觉可行,然而对回文自动机不够了解,不会做.
看别人的做法是建两个回文自动机,然后dfs找.
因为有next[][]数组,就可以保证dfs往下走总是两个自动机上相同的节点.
#include <bits/stdc++.h>
using namespace std;
#define fi first
#define se second
#define pii pair<int,int>
#define pll pair<long long ,long long>
#define mp make_pair
typedef long long ll;
#define pb push_back
const int maxn=2*1e5+15;
const int maxm=26;
const int maxe=4000015;
const int mod=51123987;
const int inf=0x3f3f3f3f;
struct palindromic_tree{
int next[maxn][maxm];
int fail[maxn];
int cnt[maxn];
int num[maxn];
int len[maxn];
int S[maxn];
int last;
int n,p;
int newnode(int l){
for(int i=0;i<maxm;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=0;
n=0;
S[n]=-1;
fail[0]=1;
}
int get_fail(int x){
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];
next[cur][c]=now;
num[now]=num[fail[now]]+1;
}
last=next[cur][c];
cnt[last]++;
}
void count(){
for(int i=p-1;i>=0;--i)cnt[fail[i]]+=cnt[i];
}
}ta,tb;
ll ans;
int n;
char str[maxn];
void dfs(int posa,int posb){
for(int i=0;i<maxm;i++){
if(ta.next[posa][i]&&tb.next[posb][i]){
int x=ta.next[posa][i];
int y=tb.next[posb][i];
ans+=ll(ta.cnt[x])*tb.cnt[y];
dfs(x,y);
}
}
}
int main(){
int T;
scanf("%d",&T);
int cas=0;
while(T--){
cas++;
ans=0;
scanf("%s",str);
n=strlen(str);
ta.init();
for(int i=0;i<n;i++){
ta.add(str[i]);
}
ta.count();
scanf("%s",str);
n=strlen(str);
tb.init();
for(int i=0;i<n;i++){
tb.add(str[i]);
}
tb.count();
dfs(0,0);
dfs(1,1);
printf("Case #%d: %lld\n",cas,ans);
}
return 0;
}
Comment in Markdown
Preview: