SPOJ - NSUBSTR Substrings

You are given a string S which consists of 250000 lowercase latin letters at most. We define F(x) as the maximal number of times that some string with length x appears in S. For example for string 'ababa' F(3) will be 2 because there is a string 'aba' that occurs twice. Your task is to output F(i) for every i so that 1<=i<=|S|.

Input

String S consists of at most 250000 lowercase latin letters.

Output

Output |S| lines. On the i-th line output F(i).

Example

Input:
ababa

Output:
3
2
2
1
1

题解:SAM模板题。

代码:

#include<cstdio>
#include<vector>
#include<cstring>
#include<algorithm>
 
using namespace std;
 
typedef long long LL;
 
const int MAXN = 1000005;
 
typedef struct Node
{
	int len, pre;//pre:对应图片中的绿线(一个节点显然最多只有一条绿线)
	int Next[26];//Next[]:对应图片中的蓝线
}Node;
 
Node tre[MAXN*2];
 
vector<int> G[MAXN*2];//用于建Suffix Links树 
int cnt, last, siz[MAXN*2], ans[MAXN*2];
char str[MAXN];
 
void Insert(char ch)
{
	int p, q, now, rev;
	p = last, now = cnt++;
	tre[now].len = tre[last].len+1;
	siz[now]++;		//如果节点u包含子串S[1..i],那么满足|siz(u)| = ∑|siz(son(u))|+1,这里先加上那个1,这样的话对于SL树就可以直接求和了
	while(p!=-1 && tre[p].Next[ch-'a']==0)	//每次跳到当前最长且siz集合与当前不同的后缀上,对应图片中的绿线回退(例如7→8→5→S)
	{
		tre[p].Next[ch-'a'] = now;		//tran(st[p], ch)=now,对应图片中的蓝线连接
		p = tre[p].pre;
	}
	if(p==-1)
		tre[now].pre = 0;		//情况①,递归到了初始节点S(空子串)(例如图片中的9号节点pre[9]=0)
	else			//如果中途某个子串tran(st[p], ch)已经存在
	{
		q = tre[p].Next[ch-'a'];
		if(tre[q].len==tre[p].len+1)		//情况②:节点q的最长子串刚好就是节点p的最长子串+S[i],也就是len[q] = len[p]+1
			tre[now].pre = q;
		else						//情况③
		{
			rev = cnt++;
			tre[rev] = tre[q];
			tre[rev].len = tre[p].len+1;	//这三行就是对Suffix Links内向树的插点操作
			tre[q].pre = tre[now].pre = rev;
			while(p!=-1 && tre[p].Next[ch-'a']==q)
			{
				tre[p].Next[ch-'a'] = rev;
				p = tre[p].pre;
			}
		}
	}
	last = now;
}
 
void SechSL(int u)		//求出所有节点的|endpos()|
{
	int v;
	for(int i=0 ; i<G[u].size() ; ++i)
	{
		v = G[u][i];
		SechSL(v);
		siz[u] += siz[v];
	}
	ans[tre[u].len] = max(ans[tre[u].len], siz[u]);
}
 
inline void Init()
{
	cnt = last = 0;
	memset(tre, 0, sizeof(tre));
	for(int i=0 ; i<MAXN*2 ; ++i)G[i].clear();
	tre[cnt++].pre = -1;
}
 
int main()
{
	int n;
	scanf("%s", str+1);
	n = strlen(str+1);
	Init();
	for(int i=1 ; i<=n ; ++i)Insert(str[i]);
	for(int i=1 ; i<=cnt-1 ; ++i)G[tre[i].pre].push_back(i);//建树 
	SechSL(0);
	/*---------------------------
	LL ans = 0;
	for(i=1;i<=cnt-1;i++)
		ans += tre[i].len-tre[tre[i].pre].len;	//求出有多少个本质不同的子串
	-----------------------------*/
	for(int i=n ; i>=1 ; i--)
		ans[i] = max(ans[i], ans[i+1]);
	for(int i=1 ; i<=n ; ++i)
		printf("%d\n", ans[i]);
		
	return 0;
}

 

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
洛谷的SPOJ需要注册一个SPOJ账号并进行绑定才能进行交题。您可以按照以下步骤进行注册: 1. 打开洛谷网站(https://www.luogu.com.cn/)并登录您的洛谷账号。 2. 在网站顶部导航栏中找到“题库”选项,将鼠标悬停在上面,然后选择“SPOJ”。 3. 在SPOJ页面上,您会看到一个提示,要求您注册SPOJ账号并进行绑定。点击提示中的链接,将会跳转到SPOJ注册页面。 4. 在SPOJ注册页面上,按照要求填写您的用户名、密码和邮箱等信息,并完成注册。 5. 注册完成后,返回洛谷网站,再次进入SPOJ页面。您会看到一个输入框,要求您输入刚刚注册的SPOJ用户名。输入用户名后,点击“绑定”按钮即可完成绑定。 现在您已经成功注册并绑定了SPOJ账号,可以开始在洛谷的SPOJ题库上刷题了。祝您顺利完成编程练习!\[1\]\[2\] #### 引用[.reference_title] - *1* *3* [(洛谷入门系列,适合洛谷新用户)洛谷功能全解](https://blog.csdn.net/rrc12345/article/details/122500057)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^control_2,239^v3^insert_chatgpt"}} ] [.reference_item] - *2* [luogu p7492 序列](https://blog.csdn.net/zhu_yin233/article/details/122051384)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^control_2,239^v3^insert_chatgpt"}} ] [.reference_item] [ .reference_list ]
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值