UVA719 Glass Beads 最小表示 / 后缀自动机

最新推荐文章于 2024-09-11 15:28:58 发布

YB Lin

最新推荐文章于 2024-09-11 15:28:58 发布

阅读量71

点赞数

本文链接：https://blog.csdn.net/m0_51068403/article/details/120538773

版权

最小表示法后缀数组字典序重复子串贪心算法

关键词由CSDN通过智能技术生成

传送门

题意：

求 $s$ 的最小表示。

$|s|\le 1e4$

思路：

可以选择直接套最小表示法的板子，但也可以将 $s$ 复制，即 $s s = s + s$ ，对 $s s$ 建 $s a m$ ，让后贪心的找一个长度为 $n$ 的字典序最小的子串，这个可以直接遍历得到，让后其 $l e n - n + 1$ 即为答案。

// Problem: Glass Beads
// Contest: Virtual Judge - UVA
// URL: https://vjudge.net/problem/UVA-719
// Memory Limit: 1024 MB
// Time Limit: 3000 ms
// 
// Powered by CP Editor (https://cpeditor.org)

//#pragma GCC optimize("Ofast,no-stack-protector,unroll-loops,fast-math")
//#pragma GCC target("sse,sse2,sse3,ssse3,sse4.1,sse4.2,avx,avx2,popcnt,tune=native")
//#pragma GCC optimize(2)
#include<cstdio>
#include<iostream>
#include<string>
#include<cstring>
#include<map>
#include<cmath>
#include<cctype>
#include<vector>
#include<set>
#include<queue>
#include<algorithm>
#include<sstream>
#include<ctime>
#include<cstdlib>
#include<random>
#include<cassert>
#define X first
#define Y second
#define L (u<<1)
#define R (u<<1|1)
#define pb push_back
#define mk make_pair
#define Mid ((tr[u].l+tr[u].r)>>1)
#define Len(u) (tr[u].r-tr[u].l+1)
#define random(a,b) ((a)+rand()%((b)-(a)+1))
#define db puts("---")
using namespace std;

//void rd_cre() { freopen("d://dp//data.txt","w",stdout); srand(time(NULL)); }
//void rd_ac() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//AC.txt","w",stdout); }
//void rd_wa() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//WA.txt","w",stdout); }

typedef long long LL;
typedef unsigned long long ULL;
typedef pair<int,int> PII;

const int N=1000010,mod=1e9+7,INF=0x3f3f3f3f;
const double eps=1e-6;

int n;
int last=1,tot=1,len;
int ed[N];
char s[N];
struct Node {
	int ch[26];
	int len,fa;
}sam[N*2];

void init() {
	for(int i=1;i<=tot;i++) {
		for(int j=0;j<26;j++) sam[i].ch[j]=0;
		sam[i].len=sam[i].fa=0;
	}
	last=tot=1; len=0;
}

void add(int c,int id) {
	int p=last,np=last=++tot; 
	sam[np].len=sam[p].len+1;
	while(p&&!sam[p].ch[c]) {
		sam[p].ch[c]=np;
		p=sam[p].fa;
	} 
	if(!p) sam[np].fa=1;
	else {
		int q=sam[p].ch[c];
		if(sam[q].len==sam[p].len+1) sam[np].fa=q;
		else {
			int nq=++tot; sam[nq]=sam[q];
			sam[nq].len=sam[p].len+1;
			sam[q].fa=sam[np].fa=nq;
			while(p&&sam[p].ch[c]==q) {
				sam[p].ch[c]=nq;
				p=sam[p].fa;
			}
		}
	}
}

int main()
{
//	ios::sync_with_stdio(false);
//	cin.tie(0);

	int _; scanf("%d",&_);
	while(_--) {
		init();
		scanf("%s",s+1); n=strlen(s+1);
		for(int i=1;i<=n;i++) add(s[i]-'a',i);
		for(int i=1;i<=n;i++) add(s[i]-'a',i);
		int st=1,nn=n;
		st=sam[st].ch['m'-'a'];
		st=sam[st].ch['a'-'a'];
		printf("%d\n",sam[st].len);
	}


	return 0;
}
/*

*/