P3804 【模板】后缀自动机(SAM)

news2024/9/23 15:26:24

题目描述

给定一个只包含小写字母的字符串 S。

请你求出 S 的所有出现次数不为 11 的子串的出现次数乘上该子串长度的最大值。

输入格式

一行一个仅包含小写字母的字符串 S。

输出格式

一个整数,为所求答案。

题解:这里就不讲后缀自动机的模板相关问题了,网有很多关于后缀自动机的详解(推荐)。我在看完了后缀自动机之后做这题,卡住我的不是后缀自动机,而是以下操作

 

	//Part 2 按len从大到小排序(和SA好像啊)后计算答案
	for(int i=1;i<=node;i++) t[len[i]]++;
	for(int i=1;i<=node;i++) t[i]+=t[i-1];
	for(int i=1;i<=node;i++) A[t[len[i]]--]=i;
	for(int i=node;i>=1;i--)
	{//从小到大枚举,实际上在模拟parent树的DFS
		int now=A[i];siz[fa[now]]+=siz[now];
		if(siz[now]>1) ans=max(ans,1ll*siz[now]*len[now]);
	}

 这段代码的要做的就是把节点 按节点代表最大字符串长度 排序,但我一直看不懂这里的代码逻辑,于是转换方式,用容易看懂的逻辑写代码。

方式一

struct edge
{
	int pos, len;
	bool operator<(edge e1)
	{
		return len > e1.len;
	}
}e[N];//存节点及其代表的最大长度


for (int i = 1; i <= tot; i++)
{
	e[i] = { i,nd[i].len };
}

sort(e + 1, e + tot + 1);//按长度从大到小排序
//因为节点所代表的字符长度越长,这个节点就越接近于叶子节点
//所以排好序后,从小到大遍历也就是从叶子节点到根节点的遍历,把endpos的个数往上累加

for (int i = 1; i <= tot; i++)
{

	siz[nd[e[i].pos].fa] += siz[e[i].pos];
}

for (int i = 1; i <= tot; i++)
{
	if (siz[i] > 1)
		ans = max(ans, (long long)(siz[i] * nd[i].len));
}

 方式二

int deg[N << 2];
queue<int> q;
void solve()
{
	while (!q.empty()) q.pop();
	memset(deg, 0, sizeof deg);
	for (int i = 1; i <= num; i++)
	{
		++deg[a[i].fa];
		siz[i] = 0;
	}

	for (int i = 1; i <= num; i++)
	{
		
		if (!deg[i]) q.push(i);//不是任何节点的父节点,即叶子节点入队
	}
	while (!q.empty())
	{
		int x = q.front();
		q.pop();
		if (pre[x])
			++siz[x];
		if (a[x].fa == 1) continue;
		siz[a[x].fa] += siz[x];
		if ((--deg[a[x].fa]) == 0)//该节点的子节点计算过一次,减一,当减为0后,即子节点全部计            
                                  //算完成,视为叶子节点入队
		{
			q.push(a[x].fa);
		}
	}
	long long ans = 0;
	for (int i = 1; i <= num; i++)
	{
		if (siz[i] > 1)
		{
			ans = max(ans, (long long)(siz[i] * a[i].len));
		}
	}

	cout << ans;
}

现在贴上这三种方式的完整代码

一,

#include<iostream>
#include<map>
#include<algorithm>	
#include<cmath> 
#include<cstring>
#include<stack>
#include<vector>
#include<queue>
using namespace std;
const long long N = 2e6 + 115;
struct node
{
	int fa, len;
	int son[26];
	long long size;
	node()
	{
		fa = len = 0;
		size = 0;
		memset(son, 0, sizeof son);
	}
}d[N];

int las, tot,n;

void add(int now)
{
	int p = las;
	int np  = ++tot;
	las = np;
	d[np].size = 1;//作为实节点,原始endpos=1
	d[np].len = d[p].len + 1;
	for (; p && !d[p].son[now]; p  = d[p].fa) d[p].son[now] = np;
	if (!p) d[np].fa = 1;
	else
	{
		int q = d[p].son[now];
		if (d[q].len == d[p].len + 1) d[np].fa = q;
		else
		{
			int nq = ++tot;
			d[nq] = d[q];
			d[nq].size = 0;//作为虚节点,原始endpos数=0
			d[nq].len = d[p].len + 1;
			d[q].fa = nq;
			d[np].fa = nq;
			for (; p && d[p].son[now] == q; p = d[p].fa) d[p].son[now] = nq;
		}
	}
}

int tmp[N], tp[N];

void so()
{
	for (int i = 0; i <= n; i++)
		tmp[i] = 0;
	for (int i = 1; i <= tot; i++)
		tmp[d[i].len]++;
	for (int i = 1; i <=n ; i++)
		tmp[i] += tmp[i - 1];
	for (int i = 1; i <=tot; i++)
		tp[tmp[d[i].len]--] = i;
	for (int i = tot; i >= 1; i--)
	{
		int now = tp[i];
		d[d[now].fa].size += d[now].size;
	}
}


void getans()
{
	long long ans = 0;
	for (int i = 1; i <= tot;i++)
	{
		
		if (d[i].size > 1)
			ans = max(ans, d[i].size * d[i].len);
	}
	cout << ans;
}

int main()
{
	string s;
	cin >> s;
	n = s.size();
	las = tot = 1;
	for (int i = 0; i < n; i++)
	{
		add(s[i] - 'a');
	}
	so();
	getans();
}

 

 二,

#include<iostream>
#include<map>
#include<algorithm>	
#include<cmath> 
#include<cstring>
#include<stack>
#include<vector>
using namespace std;
const long long N = 2e6 + 115;

int siz[N];

struct node
{
	int fa, len;
	int ch[26];
	node() { memset(ch, 0, sizeof(ch)); len = 0; }
}nd[N ];
int las = 1, tot = 1;
//las最近的一个实节点,tot节点总数 
struct edge
{
	int pos, len;
	bool operator<(edge e1)
	{
		return len > e1.len;
	}
}e[N];
void add(int c)
{
	int p = las, np;
	las = np = ++tot;
	siz[tot] = 1;
	nd[np].len = nd[p].len + 1;
	for (; p && !nd[p].ch[c]; p = nd[p].fa) nd[p].ch[c] = np;
	if (!p) nd[np].fa = 1;
	else
	{
		int q = nd[p].ch[c];
		if (nd[p].len + 1 == nd[q].len)
		{
			
			nd[np].fa = q;
		}
		else
		{
			int nq = ++tot;
			nd[nq] = nd[q];
			nd[nq].len = nd[p].len + 1;
			nd[q].fa = nd[np].fa = nq;
			for (; p && nd[p].ch[c] == q; p = nd[p].fa) nd[p].ch[c] = nq;
		}
	}
}

long long ans = 0;

int main()
{
	string s;
	cin >> s;
	int le = s.size();

	for (int i = 0; i < le; i++)
	{
		add(s[i] - 'a');
	}
	for (int i = 1; i <= tot; i++)
	{
		e[i] = { i,nd[i].len };
	}
	sort(e + 1, e + tot + 1);
	for (int i = 1; i <= tot; i++)
	{

			siz[nd[e[i].pos].fa] += siz[e[i].pos];
	}

	for (int i = 1; i <= tot; i++)
	{
		if (siz[i] > 1)
			ans = max(ans, (long long)(siz[i] * nd[i].len));
	}

	cout << ans;

}

三,

#include<iostream>
#include<map>
#include<algorithm>	
#include<cmath> 
#include<cstring>
#include<stack>
#include<vector>
#include<queue>
using namespace std;
const long long N = 1e6 + 115;

int siz[N << 2];

struct node
{
	int fa, len;
	int ch[26];
	node() { memset(ch, 0, sizeof(ch)); len = 0; }
}a[N << 2];

int cnt[N << 2];

int num = 1, np = 1;
bool pre[N << 2];
void add(int c)
{
	int p = np; np = ++num;
	pre[np] = true;//标记实节点
	a[np].len = a[p].len + 1;
	for (; p && !a[p].ch[c]; p = a[p].fa) a[p].ch[c] = np;
	if (!p) a[np].fa = 1;
	else
	{
		int q = a[p].ch[c];
		if (a[q].len == a[p].len + 1) a[np].fa = q;
		else
		{
			int nq = ++num;
			a[nq] = a[q];
			a[nq].len = a[p].len + 1;
			a[q].fa = a[np].fa = nq;
			for (; p && a[p].ch[c] == q; p = a[p].fa) a[p].ch[c] = nq;
		}
	}
}

int deg[N << 2];
queue<int> q;
void solve()
{
	while (!q.empty()) q.pop();
	memset(deg, 0, sizeof deg);
	for (int i = 1; i <= num; i++)
	{
		++deg[a[i].fa];
		siz[i] = 0;
	}

	for (int i = 1; i <= num; i++)
	{
		
		if (!deg[i]) q.push(i);//不是任何节点的父节点,即叶子节点入队
	}
	while (!q.empty())
	{
		int x = q.front();
		q.pop();
		if (pre[x])
			++siz[x];
		if (a[x].fa == 1) continue;
		siz[a[x].fa] += siz[x];
		if ((--deg[a[x].fa]) == 0)//该节点的子节点计算过一次,减一,当减为0后,即子节点全部计            
                                  //算完成,视为叶子节点入队
		{
			q.push(a[x].fa);
		}
	}
	long long ans = 0;
	for (int i = 1; i <= num; i++)
	{
		if (siz[i] > 1)
		{
			ans = max(ans, (long long)(siz[i] * a[i].len));
		}
	}

	cout << ans;
}

int main()
{
	string s;
	cin >> s;
	int le = s.size();

	for (int i = 0; i < le; i++)
	{
		add(s[i] - 'a');
	}

	solve();

}

 

特别注意:节点里的数组是整数数组,而不是字符数组,int ch[26] 而不是char ch[26] ,你猜我为什么着重提这个.....

 

 

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/732417.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

zigbee学习之DHT11温湿度传感器+zigbee无线通信

开发环境&#xff1a;IAR烧录器串口调试助手CC2530DHT11 两个模块&#xff1a;一个作为协调器&#xff0c;负责接收数据&#xff0c;一个作为终端&#xff0c;负责发送数据 步骤&#xff1a; 1、SampleApp.c里配引脚P0_6(查看硬件上的标识) 2、DTH11.c里配引脚 3、修改PANID和信…

超纯水抛光混床树脂的选择及工艺流程

一、什么是超纯水&#xff1f; 既将水中的导电介质几乎完全去除&#xff0c;又将水中不离解的胶体物质、气体及有机物均去除至很低程度的水。电阻率大于18MΩ*cm&#xff0c;或接近18.3MΩ*cm极限值。 超纯水是科技界为了研制超纯材料&#xff08;半导体原件材料、纳米精细陶…

服务器搭建oracle,并远程连接教程

下载两个压缩包&#xff0c;然后上传到服务器&#xff0c; 软件安装09&#xff1a;CentOS安装Oracle - 虚拟机 - 5997CK - 欢迎您! (hezhilin.online) 这里有全部步骤&#xff0c;反正过了几天我也会忘记&#xff0c;不赘述了。 直接上拆的坑&#xff1a; 开启服务器端口后…

借书问题-2022年全国青少年信息素养大赛Python国赛第3题

[导读]&#xff1a;超平老师计划推出《全国青少年信息素养大赛Python编程真题解析》50讲&#xff0c;这是超平老师解读Python编程挑战赛真题系列的第5讲。 全国青少年信息素养大赛&#xff08;原全国青少年电子信息智能创新大赛&#xff09;是“世界机器人大会青少年机器人设计…

9.33UEC++、容器

1.定义&#xff1a; 2.TArray&#xff1a;快&#xff0c;小&#xff0c;高 &#xff1a;同质容器 &#xff1a;创建栈对象&#xff0c;不能创建堆对象 如何构建TArray&#xff1a; 获取方式&#xff1a; 实现方式&#xff1a;定义一个array容器数组&#xff0c;两种遍历方式…

跳跃游戏(力扣)贪心 JAVA

给定一个非负整数数组 nums &#xff0c;你最初位于数组的 第一个下标 。 数组中的每个元素代表你在该位置可以跳跃的最大长度。 判断你是否能够到达最后一个下标。 示例 1&#xff1a; 输入&#xff1a;nums [2,3,1,1,4] 输出&#xff1a;true 解释&#xff1a;可以先跳 1 步…

【实用工具】(BeyondCompare+CppCheck对某次更新的文件进行静态检查)按文件清单复制文件到指定文件夹的辅助工具【2023.07.07】

摘要 开发一个功能提交了多次代码&#xff0c;现在需要比较最终的提交和某次提交的差异并进行静态代码检查。用到了BeyondCompare和CppCheck软件。但是有个问题&#xff0c;BeyondCompare不能把差异项导出&#xff0c;于是乎花了两个小时写了个辅助工具。 使用说明 第一步&…

校服选购定制小程序开发制作功能介绍

本次为各大校服服装定制生产老板介绍一下校服定制商城小程序系统的功能&#xff0c;为大家做参考。 校服定制小程序主要功能有&#xff1a; 1、每个学校校服定制信息独立&#xff0c;不同的学校打开小程序里面信息可以不同。 2、学校致家长的一封信。 3、学生信息录入和下单信息…

【VUE】项目设置超时6分钟,实际1分钟就超时了

一、背景&#xff1a;VUE项目中&#xff0c;前端接口请求设置6分钟&#xff0c;但实际在浏览器中1分钟就超时了 timeout: 6 * 60 * 1000二、经查&#xff1a;不是前端代码设置的问题&#xff0c;也不是浏览器的问题&#xff0c;而是nginx超时设置默认为60S&#xff0c;导致1分…

hexo个人博客搭建(二)butterfly主题配置

Butterfly主题安装文档(二)之主题配置 一、回顾安装butterfly主题 1、在hexo项目根目录下执行操作clone主题 git clone -b master https://github.com/jerryc127/hexo-theme-butterfly.git themes/butterfly2、如果沒有 pug 以及 stylus 的渲染器&#xff0c;还需要下载&…

matlab语言的前世今生

matlab语言的算法库是用什么语言实现的? Matlab的算法库是用C和C++语言实现的。 Matlab本身是一个解释型的高级编程语言,但它的核心算法库是使用更底层的C和C++语言编写的。这些库提供了许多数值计算、统计分析、信号处理、图像处理、优化等常用算法函数,并以mex文件的形式集…

投票评选活动小程序v2-搜索功能和最新排序功能实现

投票评选活动小程序-搜索功能和最新排序功能实现 优秀人物评选活动&#xff0c;五四奖章评选、优秀工作者人物评选、劳动最光荣评选。通常需要一个搜索功能&#xff0c;用户通过输入名称或编号搜索对应的作品项目或者人物&#xff0c;进行投票。或者通过最新排序功能查看列表情…

pytorch超详细安装教程,Anaconda、PyTorch和PyCharm整套安装流程

本文介绍基于Anaconda环境以及PyCharm软件结合&#xff0c;安装PyTorch深度学习框架。 PyTorch深度学习框架详细安装教程 一、anaconda安装&#xff08;一&#xff09;下载&#xff08;二&#xff09;安装&#xff08;三&#xff09;配置环境变量&#xff08;四&#xff09;检查…

虚假人脸检测实验

虚假人脸检测实验 虚假人脸识别 数据集链接 链接&#xff1a;https://pan.baidu.com/s/1hDyJ91dAwI5j5GTR0hD_cA?pwd4cki 原理 ResNet-18是一种经典的CNN网络&#xff0c;是 Deep Residual Learning 团队在 2017 年提出的。它是为了解决 ImageNet 数据集上的图像分类任务而…

会话机制【Cookie 和 Session】,登陆页面的模拟实现

前言 小亭子正在努力的学习编程&#xff0c;接下来将开启JavaEE的学习~~ 分享的文章都是学习的笔记和感悟&#xff0c;如有不妥之处希望大佬们批评指正~~ 同时如果本文对你有帮助的话&#xff0c;烦请点赞关注支持一波, 感激不尽~~ 目录 前言 Cookie 和 Session 是什么 Cookie…

【计算机网络】计算机网络概述

&#x1f525; 本文由 程序喵正在路上 原创&#xff0c;CSDN首发&#xff01; &#x1f496; 系列专栏&#xff1a;计算机网络 &#x1f320; 首发时间&#xff1a;2023年7月8日 &#x1f98b; 欢迎关注&#x1f5b1;点赞&#x1f44d;收藏&#x1f31f;留言&#x1f43e; &…

基于Spring Boot的社区适龄青年征兵系统设计与实现(Java+spring boot+MySQL)

获取源码或者论文请私信博主 演示视频&#xff1a; 基于Spring Boot的社区适龄青年征兵系统设计与实现&#xff08;Javaspring bootMySQL&#xff09; 使用技术&#xff1a; 前端&#xff1a;html css javascript jQuery ajax thymeleaf 后端&#xff1a;Java springboot框架 …

BOM操作

JavaScript组成 BOM 浏览器对象模型 window对象 是一个全局对象&#xff0c;也就是JavaScript中的顶级对象 像document&#xff0c;alert() console.log() 都是window对象的属性&#xff0c; 基本的BOM的属性和方法都属于window对象 通过var定义在全局作用域中的变量&#x…

Ruoyi单体项目与Echarts4.2.1地图集成时的思路及解决办法

目录 背景 一、相关数据 1、湖南省2021全省地区生产总值 2、湖南Geojson数据 二、Ruoyi集成设计与实现 1、自定义地图注册 2、湖南地图引用 3、图表初始化及数据绑定 4、实际效果 总结 背景 在之前博客中&#xff0c;介绍了Echarts和Ruoyi集成的一些博客基于Echarts2.X的…

时间序列预测 | Matlab麻雀算法(SSA)优化极限梯度提升树XGBoost时间序列预测,SSA-XGBoost时间序列预测模型,单列数据输入模型

文章目录 效果一览文章概述部分源码参考资料效果一览 文章概述 麻雀算法(SSA)优化极限梯度提升树XGBoost时间序列预测,SSA-XGBoost时间序列预测模型,单列数据输入模型 评价指标包括:MAPE、RMSE等,代码质量极高,方便学习和替换数据。要求2018版本及以上。 部分源码