Lyndon分解

最新推荐文章于 2024-05-05 16:35:53 发布

selia1078

最新推荐文章于 2024-05-05 16:35:53 发布

阅读量1k

点赞数

分类专栏：字符串

本文链接：https://blog.csdn.net/micaudience/article/details/108110906

版权

字符串专栏收录该内容

7 篇文章 0 订阅

订阅专栏

Lyndon分解

这是什么

Lyndon串：对于一个串 $s_{1…n}$ 的所有后缀， $s_{1…n}$ 的字典序最小，此时，我们称 $s_{1…n}$ 是Lyndon串。

近似Lyndon串：设一个串s是Lyndon串，s’是s的前缀，那么形如ssss’的形式的串称之为近似Lyndon串。

Lyndon分解：将一个串拆成若干个Lyndon串的分解方式。
例如：cuucuucud拆成cuu,cuu,cud

怎么求

设 i 是目前的起始位置，令j = i, k = i + 1,比较 $s_{j}$ 与 $s_{k}$ 的大小：

1. 如果 $s_{j} = s_{k}$ ，那么说明在这一位置无法比较出大小，继续向后比较，j， k同时向后移动。

2. 如果 $s_{j} < s_{k}$ ，那么说明 $suf_{k - (j - i)}$ 比 $suf_{i}$ 字典序大，且其中以 $suf_{k - (j - i)}……suf_{k}$ 都要对应地比 $suf_{i}……suf_{j}$ 大，那么也就必 $suf_{i}$ 大。此时，将 j 移动到 i，从头开始新一轮的比较。

3. 如果 $s_{j}>s_{i}$ ，那么说明 $suf_{k - (j - i)}$ 比 $suf_{i}$ 字典序小，此时意味着当前Lyndon串再延续下去就不合法了。此时，终止内层循环。

在每一轮循环的结束，我们将得到一个近似Lyndon串 $s_{i…k-1}$ ，其中真正的Lyndon串的大小为 k - j ，那么我们每隔 k - j 位做一次分割即可。

分割完成后，将 i 跳到下一个起始位置即可（注意近似Lyndon串的"零头"要重新做）。

洛谷例题

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>
#define LL long long
using namespace std;
const int N = 5e6 + 10;
char s[N];
int main()
{
	int ans = 0, n;
	scanf("%s", s + 1);
	n = strlen(s + 1);
	for(int i = 1; i <= n; )
	{
		int j = i;
		int k = i + 1;
		while(k <= n && s[j] <= s[k])
		{
			if(s[j] == s[k])	j++;
			else	j = i;
			k++;
		}
		while(i <= j)
		{
			i += k-j;
			ans = ans ^ (i - 1);
		}
	}
	cout<<ans;
	return 0;
 }

有什么用

求最小后缀

HDU-6761

题意：对于一个串的每一个前缀，求其最小后缀。

思路：求串的Lyndon分割，如果 $s_{j} < s_{k}$ 那么跳回到当前近似Lyndon串的串首，否则，由 j 位置顺延一个Lyndon串（也就是回到当前近似Lyndon串的最后一个循环节的起始位置）。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>
#define LL long long
using namespace std;
const int N = 1e6 + 10;
const LL mod = 1e9 + 7; 
LL f[N], d[N];
char s[N];
int main()
{
	f[1] = 1;
	for(int i = 2; i <= 1e6; i++)
		f[i] = f[i - 1] * (LL)1112 % mod; 
	int T;
	scanf("%d", &T);
	while(T)
	{
		T--;
		LL ans = 0;
		scanf("%s", s + 1);
		int n = strlen(s + 1);
		d[1] = 1;
		for(int i = 1; i <= n;)
		{
			int j = i;
			int k = i + 1;
			while(k <= n && s[j] <= s[k])
			{
				if(s[j] == s[k])	d[k] = d[j] + k - j, j++;	//顺延一个串的位置 
				else	d[k] = i, j = i;	//回串首 
				k++;
			}
			d[k] = k; //下一个串的开始 
			while(i <= j)	i += k - j;
		}
		for(int i = 1; i <= n; i++)	ans = (ans + d[i] * f[i] % mod) % mod;
		cout<< ans <<endl;
	}
	return 0;
}