后缀自动机 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-06-15matplotlib作图不显示3D图,怎么办?
- 2024-06-1503-Loki 日志监控
- 2024-06-1504-让LLM理解知识 -Prompt
- 2024-06-05做软件测试需要懂代码吗?
- 2024-06-0514-ShardingSphere的分布式主键实现
- 2024-06-03为什么以及如何要进行架构设计权衡?
- 2024-05-31全网首发第二弹!软考2024年5月《软件设计师》真题+解析+答案!(11-20题)
- 2024-05-31全网首发!软考2024年5月《软件设计师》真题+解析+答案!(21-30题)
- 2024-05-30【Java】百万数据excel导出功能如何实现
- 2024-05-30我们小公司,哪像华为一样,用得上IPD(集成产品开发)?