后缀自动机 SAM
2022/6/12 23:23:53
本文主要是介绍后缀自动机 SAM,对大家解决编程问题具有一定的参考价值,需要的程序猿们随着小编来一起学习吧!
struct NODE { int ch[26]; int len,fa; NODE(){memset(ch,0,sizeof(ch));len=0;} }dian[MAXN<<1]; int las=1,tot=1; void add(int c) { int p=las;int np=las=++tot; dian[np].len=dian[p].len+1; for(;p&&!dian[p].ch[c];p=dian[p].fa)dian[p].ch[c]=np; if(!p)dian[np].fa=1;//以上为case 1 else { int q=dian[p].ch[c]; if(dian[q].len==dian[p].len+1)dian[np].fa=q;//以上为case 2 else { int nq=++tot;dian[nq]=dian[q]; dian[nq].len=dian[p].len+1; dian[q].fa=dian[np].fa=nq; for(;p&&dian[p].ch[c]==q;p=dian[p].fa)dian[p].ch[c]=nq;//以上为case 3 } } } char s[MAXN];int len; int main() { scanf("%s",s);len=strlen(s); for(int i=0;i<len;i++)add(s[i]-'a'); }
void add(int c) { int p=las;int np=las=++tot; dian[np].len=dian[p].len+1; ```![](https://www.www.weizhi.cc/i/l/?n=22&i=blog/2747001/202206/2747001-20220612212002420-484956904.png) ![](https://www.www.weizhi.cc/i/l/?n=22&i=blog/2747001/202206/2747001-20220612211757935-1270651735.png) ```cpp for(;p&&!dian[p].ch[c];p=dian[p].fa)dian[p].ch[c]=np;
if(!p)dian[np].fa=1;//以上为case 1`
else { int q=dian[p].ch[c]; if(dian[q].len==dian[p].len+1)dian[np].fa=q;
else { int nq=++tot;dian[nq]=dian[q]; dian[nq].len=dian[p].len+1; dian[q].fa=dian[np].fa=nq; for(;p&&dian[p].ch[c]==q;p=dian[p].fa)dian[p].ch[c]=nq;//以上为case 3
strcpy(dian[nq].ch,dian[q].ch); dian[nq].fa=dian[q].fa;
for(;p&&!dian[p].ch[c];p=dian[p].fa)dian[p].ch[c]=np;
for(;p&&dian[p].ch[c]==q;p=dian[p].fa)dian[p].ch[c]=nq;
P3804 【模板】后缀自动机 (SAM)
#include<cstdio> #include<cstring> #include<string> #include<cmath> #include<algorithm> #include<iostream> using namespace std; int tot=1,las=1; struct NODE { int ch[26]; int len,fa; NODE(){memset(ch,0,sizeof(ch));len=fa=0;} }dian[2000010]; struct Edge { int t,nexty; }edge[2000010]; int head[2000010],cnt=0; void jia(int a,int b) { cnt++; edge[cnt].t=b; edge[cnt].nexty=head[a]; head[a]=cnt; } long long zhi[2000010]; inline void add(int c) { register int p=las,np=las=++tot;zhi[tot]=1; dian[np].len=dian[p].len+1; for(;p&&!dian[p].ch[c];p=dian[p].fa)dian[p].ch[c]=np; if(!p)dian[np].fa=1; else { register int q=dian[p].ch[c]; if(dian[q].len==dian[p].len+1)dian[np].fa=q; else { register int nq=++tot; dian[nq]=dian[q];dian[nq].len=dian[p].len+1; dian[q].fa=dian[np].fa=nq; for(;p&&dian[p].ch[c]==q;p=dian[p].fa)dian[p].ch[c]=nq; } } } char s[2000010]; int cd; long long ans=0; void dfs(int node) { for(register int i=head[node];i;i=edge[i].nexty) { dfs(edge[i].t); zhi[node]+=zhi[edge[i].t]; } if(zhi[node]!=1)ans=max(ans,zhi[node]*dian[node].len); } int main() { scanf("%s",s);cd=strlen(s); for(register int i=0;i<cd;i++)add(s[i]-'a'); for(register int i=2;i<=tot;i++)jia(dian[i].fa,i); dfs(1); printf("%lld\n",ans); return 0; }
这篇关于后缀自动机 SAM的文章就介绍到这儿,希望我们推荐的文章对大家有所帮助,也希望大家多多支持为之网!
- 2024-05-13TiDB + ES:转转业财系统亿级数据存储优化实践
- 2024-05-09“2024鸿蒙零基础快速实战-仿抖音App开发(ArkTS版)”实战课程已上线
- 2024-05-09聊聊如何通过arthas-tunnel-server来远程管理所有需要arthas监控的应用
- 2024-05-09log4j2这么配就对了
- 2024-05-09nginx修改Content-Type
- 2024-05-09Redis多数据源,看这篇就够了
- 2024-05-09Google Chrome驱动程序 124.0.6367.62(正式版本)去哪下载?
- 2024-05-09有没有大佬知道这种数据应该怎么抓取呀?
- 2024-05-09这种运行结果里的10.100000001,怎么能最快改成10.1?
- 2024-05-09企业src漏洞挖掘-有意思的命令执行