PAM学习小结
PAM
回文自动机html
建议先学习AC自动机：AC自动机讲解超详细c++
回文自动机，顾名思义，用来处理回文串的自动机。数组
功能：
1.求\(S\)串内本质不一样的回文串个数学习
2.求\(S\)串内本质不一样的回文串出现次数spa
3.最小回文划分指针
4.\(S\)串中如下标\(i\)结尾的最长回文串长度code
回文树
看看本身感悟一下。感受特别形象，都不用解释了啊htm
仍是稍微解释一下：blog
1.回文数上每个节点表明了原串上出现过的一个本质不一样回文子串，原串上的每个回文子串都在回文树上有对应。回文树上每个点表明的串都是回文串。排序
2.回文树分两部分，奇和偶，奇树上的点表明的回文串长度为奇数，偶树上的为偶
3.儿子节点表明串长度为父亲节点表明串长度\(+2\)
4.和\(Trie\)类似的其余性质，不说了
Fail指针
学过AC自动机的OIer们应该就很熟悉啦QwQ
\(Fail\)指针含义：这个节点所表明的回文串的最长回文后缀
Trans指针
通常作许多PAM题目经常使用的东西
\(Trans\)指针含义：小于等于当前节点长度一半的最长回文后缀
构建PAM
咱们要维护如下信息
char s[maxn]; //原串 int fail[maxn]; //fail指针 int len[maxn]; //该节点表示的字符串长度 int tree[maxn][26]; //同Trie，指向儿子 int trans[maxn]; //trans指针 int tot,pre; //tot表明节点数，pre表明上次插入字符后指向的回文树位置
其中\(fail,len,tree,trans\)为PAM上的信息
构建PAM的方法为增量，即一个一个加入字符构建PAM
奇树和偶树的根长度\(len\)分别为\(-1\)和\(0\)
设当前咱们插入原串中\(i\)位置的字符\(u\)
那么以\(i\)为结尾的最长回文串应该为(以\(i-1\)为结尾的最长回文串\(+u\))，而且那个回文串要知足前一个字符等于\(u\)（否则就不是回文串了啊）
要找到那个点很是简单，不断从\(pre\)开始跳\(fail\)，直到找到一个知足\(s[i-len[x]-1]==u\) 的节点\(Fail\) ，那么从\(Fail\)建一个\(u\)儿子便可以表示新的回文串。
新点的\(fail\)怎么求呢。
明显为从\(pre\)开始跳\(fail\)，找到{ [第二个(知足\(s[i-len[x]-1]==u\)) 的节点\(x\) ]的\(u\)儿子 }
也就是从\(Fail\)开始跳\(fail\)，找到{ [第一个(知足\(s[i-len[x]-1]==u\)) 的节点\(x\) ]的\(u\)儿子 }
跳到根记得判断
特别提醒:节点\(1\)为奇根，节点\(0\)为偶根，\(fail[0]=1\) , \(len[1]=-1\)
时间复杂度证实参考OIwiki：OIwiki-PAM
放代码理解：
int getfail(int x,int i){ //从x开始跳fail，知足字符s[i]的节点
while(i-len[x]-1<0||s[i-len[x]-1]!=s[i])x=fail[x];
return x;
}
void insert(int u,int i){
int Fail=getfail(pre,i); //找到符合要求的点
if(!tree[Fail][u]){ //没建过就新建节点
len[++tot]=len[Fail]+2; //长度天然是父亲长度+2
fail[tot]=tree[getfail(fail[Fail],i)][u]; //fail为知足条件的次短回文串+u
tree[Fail][u]=tot; //认儿子
}
pre=tree[Fail][u]; //更新pre
}
至于\(trans\)维护也和\(fail\)差很少
根据\(trans\)的定义去推一下怎么搞吧
放一下完整代码：
char s[maxn]; //原串
int fail[maxn]; //fail指针
int len[maxn]; //该节点表示的字符串长度
int tree[maxn][26]; //同Trie，指向儿子
int trans[maxn]; //trans指针
int tot,pre; //tot表明节点数，pre表明上次插入字符后指向的回文树位置
int getfail(int x,int i){ //从x开始跳fail，知足字符s[i]的节点
while(i-len[x]-1<0||s[i-len[x]-1]!=s[i])x=fail[x];
return x;
}
int gettrans(int x,int i){
while(((len[x]+2)<<1)>len[tot]||s[i-len[x]-1]!=s[i])x=fail[x];
return x;
}
void insert(int u,int i){
int Fail=getfail(pre,i); //找到符合要求的点
if(!tree[Fail][u]){ //没建过就新建节点
len[++tot]=len[Fail]+2; //长度天然是父亲长度+2
fail[tot]=tree[getfail(fail[Fail],i)][u]; //fail为知足条件的次短回文串+u
tree[Fail][u]=tot; //指儿子
if(len[tot]<=2)trans[tot]=fail[tot]; //特殊trans
else{
int Trans=gettrans(trans[Fail],i); //求trans
trans[tot]=tree[Trans][u];
}
}
pre=tree[Fail][u]; //更新pre
}
应用
P5496【模板】回文自动机（PAM）
求第 i 个整数表示原串以第 i 个字符结尾的回文子串个数，强制在线
明显：一个回文串的答案等于其最长回文后缀的答案\(+1\) （这超好理解的吧
那就在多维护一个信息\(ans\)表示答案，新建节点时更新便可
ans[tot]=ans[fail[tot]]+1;
答案为lastans=ans[pre];
代码：
#include<bits/stdc++.h>
#define maxn 510001
using namespace std;
char s[maxn];
int fail[maxn],len[maxn],ans[maxn],trie[maxn][26];
int pre,slen,lastans,tot;
int getfail(int x,int i){
while(i-len[x]-1<0||s[i-len[x]-1]!=s[i])x=fail[x];
return x;
}
int main(){
scanf("%s",s);slen=strlen(s);
fail[0]=1;len[1]=-1;tot=1;
for(int i=0;i<slen;i++){
if(i>=1)s[i]=(s[i]-97+lastans)%26+97;
int u=s[i]-'a';
int Fail=getfail(pre,i);
if(!trie[Fail][u]){
fail[++tot]=trie[getfail(fail[Fail],i)][u];
trie[Fail][u]=tot;
len[tot]=len[Fail]+2;
ans[tot]=ans[fail[tot]]+1;
}
pre=trie[Fail][u];
lastans=ans[pre];
printf("%d ",lastans);
}
return 0;
}
P4287[SHOI2011]双倍回文
学好\(trans\)指针，秒切此题
明显：当存在\(i\)知足\(len[trans[i]]*2==len[i]\)而且知足题意中\(len[trans[i]]%2==0\)即为符合题意的串，取最长便可。
代码：真·模板
#include<bits/stdc++.h>
#define maxn 510001
using namespace std;
char s[maxn]; //原串
int fail[maxn]; //fail指针
int len[maxn]; //该节点表示的字符串长度
int tree[maxn][26]; //同Trie，指向儿子
int trans[maxn]; //trans指针
int tot,pre; //tot表明节点数，pre表明上次插入字符后指向的回文树位置
int getfail(int x,int i){ //从x开始跳fail，知足字符s[i]的节点
while(i-len[x]-1<0||s[i-len[x]-1]!=s[i])x=fail[x];
return x;
}
int gettrans(int x,int i){
while(((len[x]+2)<<1)>len[tot]||s[i-len[x]-1]!=s[i])x=fail[x];
return x;
}
void insert(int u,int i){
int Fail=getfail(pre,i); //找到符合要求的点
if(!tree[Fail][u]){ //没建过就新建节点
len[++tot]=len[Fail]+2; //长度天然是父亲长度+2
fail[tot]=tree[getfail(fail[Fail],i)][u]; //fail为知足条件的次短回文串+u
tree[Fail][u]=tot; //指儿子
if(len[tot]<=2)trans[tot]=fail[tot]; //特殊trans
else{
int Trans=gettrans(trans[Fail],i); //求trans
trans[tot]=tree[Trans][u];
}
}
pre=tree[Fail][u]; //更新pre
}
int slen,ans;
int main(){
scanf("%d",&slen);
scanf("%s",s);
fail[0]=1;len[1]=-1;tot=1;
for(int i=0;i<slen;i++)insert(s[i]-'a',i);
for(int i=2;i<=tot;i++){
if(len[trans[i]]*2==len[i]&&len[trans[i]]%2==0)
ans=max(ans,len[i]);
}
printf("%d\n",ans);
return 0;
}
P4555[国家集训队]最长双回文串
题目描述： 顺序和逆序读起来彻底同样的串叫作回文串。好比`acbca`是回文串，而`abc`不是（`abc`的顺序为`abc`，逆序为`cba`，不相同）。 输入长度为n的串S，求S的最长双回文子串T,便可将T分为两部分X，Y，（|X|,|Y|≥1）且X和Y都是回文串。
简单PAM题
题解：
正着建一棵PAM，\(a[i]\)记录当前位置\(i\)结尾最长回文串长度
反着建一棵PAM，\(b[i]\)记录当前位置\(i\)结尾最长回文串长度
\(a[i]+b[i+1]\)即为以\(i\)为分界的双回文串
取\(a[i]+b[i+1]\)最大值即为答案
代码：封装版PAM
#include<bits/stdc++.h>
#define maxn 510001
using namespace std;
char s[maxn];
int slen,a[maxn],b[maxn],ans;
struct PAM{
int fail[maxn],len[maxn],trie[maxn][26];
int tot,pre;
void init(){fail[0]=1;len[1]=-1;tot=1;pre=0;}
int getfail(int x,int i){
while(i-len[x]-1<0||s[i-len[x]-1]!=s[i])x=fail[x];
return x;
}
void insert(int u,int i){
int Fail=getfail(pre,i);
if(!trie[Fail][u]){
fail[++tot]=trie[getfail(fail[Fail],i)][u];
trie[Fail][u]=tot;
len[tot]=len[Fail]+2;
}
pre=trie[Fail][u];
}
}A,B;
int main(){
scanf("%s",s);slen=strlen(s);A.init();B.init();
for(int i=0;i<slen;i++)A.insert(s[i]-'a',i),a[i]=A.len[A.pre];
reverse(s,s+slen); //翻转
for(int i=0;i<slen;i++)B.insert(s[i]-'a',i),b[slen-i-1]=B.len[B.pre];
for(int i=0;i<slen-1;i++)ans=max(ans,a[i]+b[i+1]);
printf("%d\n",ans);
return 0;
}
P4762[CERC2014]Virus synthesis
PAM好题，请好好思考
题意：
初始有一个空串，利用下面的操做构造给定串 \(S\) 。
一、串开头或末尾加一个字符
二、串开头或末尾加一个该串的逆串
求最小化操做数， \(∣\ S∣≤10^5\) 。
题解：
PAM上dp
P1659[国家集训队]拉拉队排练
一眼可得PAM
咱们令PAM上多记录一个信息\(sum\)，表示该节点表示串在原串上出现了多少次。
当咱们处理完了\(sum\)，对于长度\(len\)为奇数的节点的信息\(sum\)计入数组\(a[i]\).
\(a[i]\)为长度为\(i\)的回文子串出现次数。
\(a[i]\)降序排序后累加答案快速幂处理一下便可，不需太多点拨
重点来了
讲一下怎么处理\(sum\)
咱们能够发现当一个节点\(u\)的\(sum+1\)，那么\(fail[u]\)的\(sum\)也要\(+1\)
熟悉AC自动机的OIer能够敏锐的察觉到能够用拓扑排序了（例如我
建PAM的时候打个标记，最后统一一个拓扑排序向\(fail\)去更新\(sum\)便可
queue<int >q; //in数组为fail入边数量
void tuopu(){
for(int i=0;i<=tot;i++)if(in[i]==0)q.push(i);
while(!q.empty()){
int u=q.front();q.pop();
sum[fail[u]]+=sum[u];in[fail[u]]--;
if(in[fail[u]]==0)q.push(fail[u]);
}
}
好像没什么问题，多一个拓扑排序就好了
但真的如此吗？
咱们观察PAM和AC自动机的区别
AC自动机是建好\(Trie\)后再进行\(getFail\)的，\(fail\)的节点编号是会大于自身节点编号
而PAM不会出现这种状况，PAM\(fail\)定义不一样于AC自动机，构建使用增量法，保证了\(fail\)的节点编号必定小于自身节点编号。
因此就能够不用拓扑排序了，直接一个\(for\)从后到前更新便可
for(int i=tot;i>=0;i--)sum[fail[i]]+=sum[i];
总代码：
#include<bits/stdc++.h>
#define maxn 1010001
#define ll long long
#define mod 19930726
using namespace std;
char s[maxn];
int fail[maxn],len[maxn],trie[maxn][26],trans[maxn];
long long sum[maxn];
int per,slen,tot;
long long a[maxn],K,ans=1;
int getfail(int x,int i){
while(i-len[x]-1<0||s[i-len[x]-1]!=s[i])x=fail[x];
return x;
}
int gettrans(int x,int i){
while(((len[x]+2)<<1)>len[tot]||s[i-len[x]-1]!=s[i])x=fail[x];
return x;
}
void insert(int u,int i){
int Fail=getfail(per,i);
if(!trie[Fail][u]){
len[++tot]=len[Fail]+2;
fail[tot]=trie[getfail(fail[Fail],i)][u];
trie[Fail][u]=tot;
if(len[tot]<=2)trans[tot]=fail[tot];
else{
int Trans=gettrans(trans[Fail],i);
trans[tot]=trie[Trans][u];
}
}
per=trie[Fail][u];
sum[per]++; //记录sum
}
ll qpow(ll n,ll m){
ll ans=1ll;
while(m){
if(m&1){ans=ans*n;ans%=mod;}
n=n*n;n%=mod;m>>=1;
}return ans%mod;
}
int main(){
scanf("%d%lld",&slen,&K);
scanf("%s",s);
fail[0]=1;len[1]=-1;tot=1;
for(int i=0;i<slen;i++)insert(s[i]-'a',i);
for(int i=tot;i>=1;i--)sum[fail[i]]+=sum[i]; //更新sum
for(int i=2;i<=tot;i++)a[len[i]]+=sum[i],a[len[i]]%=mod; //长度处理
for(int i=slen;i>=1;i--){ //答案处理
if(i%2==1){
if(K>=a[i]){
ans*=qpow(i,a[i]);ans%=mod;
K-=a[i];
}else{
ans*=qpow(i,K);ans%=mod;
K-=K;
break;
}
}
}
if(K==0) //判-1
printf("%lld\n",ans%mod);
else
printf("-1\n");
return 0;
}
CF17E Palisection
卡空间PAM，2010没有PAM，因此都是马拉车
众所周知，PAM拥有十分优秀的时间复杂度，但空间复杂度lj得不行
但这题卡空间，因此得用到邻接链表PAM
先讲思路
题目要求相交的回文子串对，这很难作
因而咱们求补集，求不相交的回文子串对，再用总数减便可
求法和上文的最长双回文子串 相似
正反建一次PAM，存该位置结尾的回文子串个数，而后加法改乘法
本身领悟一下，挺简单的。
如今讲一下邻接链表PAM
注意：邻接链表PAM不是使空间变小了，而是用时间换空间
咱们记边结构体\(line\)
存\(3\)个信息：\(nx,to,w\) 分别表示上一条边，这条边通向的节点编号，这条边是表明哪一个字符
数组\(fir[i]\)表示\(i\)伸出的最后一条边的编号（头插式
当咱们要寻找\(u\)的\(v\)儿子
咱们就像邻接链表同样找，直到有一条边的\(w==v\)为止
找不到记得指根
int getson(int u,int v){
for(int i=u;i!=-1;i=l[i].nx)
if(l[i].w==v)return l[i].to;
return -1;
}
建点的时候把边建上
void insert(int u,int i){
int Fail=getfail(pre,i),ls=getfail(fail[Fail],i);
if(getson(fir[Fail],u)==-1){
if(getson(fir[ls],u)==-1)fail[++tot]=0; //找不到指根
else fail[++tot]=getson(fir[ls],u); //找到了
l[++cnt]=(line){fir[Fail],tot,u};fir[Fail]=cnt; //加边
len[tot]=len[Fail]+2;
ans[tot]=ans[fail[tot]]+1; //结尾回文子串个数
pre=tot;
}else
pre=getson(fir[Fail],u);
}
然鹅事实上你仍然过不了，你还要继续压空间，省掉一堆数组就能够过啦！
总代码：
#include<bits/stdc++.h>
#define maxn 2000005
#define mod 51123987
using namespace std;
char s[maxn];
int slen,b[maxn];
long long res;
int fail[maxn],len[maxn],ans[maxn],fir[maxn];
struct line{int nx,to,w;}l[maxn];
int tot,pre,cnt;
void init(){
memset(fir,-1,sizeof(fir));cnt=0;
fail[0]=1;len[1]=-1;tot=1;pre=0;
}
int getfail(int x,int i){
while(i-len[x]-1<0||s[i-len[x]-1]!=s[i])x=fail[x];
return x;
}
int getson(int u,int v){
for(int i=u;i!=-1;i=l[i].nx)
if(l[i].w==v)return l[i].to;
return -1;
}
void insert(int u,int i){
int Fail=getfail(pre,i),ls=getfail(fail[Fail],i);
if(getson(fir[Fail],u)==-1){
if(getson(fir[ls],u)==-1)fail[++tot]=0;
else fail[++tot]=getson(fir[ls],u);
l[++cnt]=(line){fir[Fail],tot,u};fir[Fail]=cnt;
len[tot]=len[Fail]+2;
ans[tot]=ans[fail[tot]]+1;
pre=tot;
}else
pre=getson(fir[Fail],u);
}
int main(){
int n;
scanf("%d",&n);
scanf("%s",s);slen=strlen(s);init();
reverse(s,s+slen);
for(int i=0;i<slen;i++)insert(s[i]-'a',i),b[slen-i-1]=ans[pre];
for(int i=slen-1;i>=0;i--)b[i]+=b[i+1],b[i]%=mod;
reverse(s,s+slen);init();
for(int i=0;i<slen-1;i++){
insert(s[i]-'a',i);int x=ans[pre];
res+=(1ll*x*b[i+1])%mod,res%=mod;
}
printf("%lld\n",((1ll*b[0]*(b[0]-1)/2ll)%mod-res+mod)%mod);
return 0;
}
To be continue……
- 1. Linux下PAM模块学习总结
- 2. Linux-PAM & PAM-MySQL的总结
- 3. 学习小结
- 4. Linux学习笔记：sudo、tcp_wrapper、pam
- 5. 回文树 PAM 学习笔记
- 6. PAM
- 7. MySQL 学习小结
- 8. vue 学习小结
- 9. kdump 学习小结
- 10. JavaScript学习小结
- 更多相关文章...
- • 您已经学习了 XML Schema，下一步学习什么呢？ - XML Schema 教程
- • 我们已经学习了 SQL，下一步学习什么呢？ - SQL 教程
- • Tomcat学习笔记（史上最全tomcat学习笔记）
- • 适用于PHP初学者的学习线路和建议
-
每一个你不满意的现在，都有一个你没有努力的曾经。
- 1. 「插件」Runner更新Pro版，帮助设计师远离996
- 2. 错误 707 Could not load file or assembly ‘Newtonsoft.Json, Version=12.0.0.0, Culture=neutral, PublicKe
- 3. Jenkins 2018 报告速览，Kubernetes使用率跃升235%！
- 4. TVI-Android技术篇之注解Annotation
- 5. android studio启动项目
- 6. Android的ADIL
- 7. Android卡顿的检测及优化方法汇总（线下+线上）
- 8. 登录注册的业务逻辑流程梳理
- 9. NDK(1)创建自己的C/C++文件
- 10. 小菜的系统框架界面设计-你的评估是我的决策
- 1. Linux下PAM模块学习总结
- 2. Linux-PAM & PAM-MySQL的总结
- 3. 学习小结
- 4. Linux学习笔记：sudo、tcp_wrapper、pam
- 5. 回文树 PAM 学习笔记
- 6. PAM
- 7. MySQL 学习小结
- 8. vue 学习小结
- 9. kdump 学习小结
- 10. JavaScript学习小结