P3426 [POI2005]SZA-Template

题目：

给定一个字符串$s$，要你用一个印章印出$s$，其中相同的位置可以多次印上相同的字符，问印章上字符串的最小长度。
$(1\le |s|\le 5\times 10^5)$

题解：

首先可以发现印章必是$s$的一个前缀，也是$s$的一个后缀，所以印章是$s$的一个$\text{Border}$。
设印章的字符串为$b$，印章印出的字符串为$t$，可以发现$t$每次经过扩展以后，$b$都是$t$的一个$\text{Border}$，设每次扩展前的位置为$p_1$，扩展后的位置为$p_2$，那么要满足$p_2-p_1\le |b|$。
所以现在问题就变成找到长度最小的$s$的$\text{Border}$ $|b|$满足任意两个以它为$\text{Border}$的$s$的前缀$p_1,p_2$满足$|p_2-p_1|\le |b|$。我们可以求出$s$的$next$数组，建出$fail$树，找到从根到$n$的唯一路径$P$，这条路径上除根以外的所有点即为$s$的所有$\text{Border}$，那么我们找出路径$P$上深度最小的点$u$满足把以$u$为根的子树中的所有结点代表的前缀位置从小到大排序后相邻的位置之差小于等于$u$表示的前缀位置。可以从根往结点$n$走，然后用链表维护子树中的结点按前缀位置排序后的序列，每次往下走时都把其余子树中的点从链表中删掉，然后更新相邻位置之差的最大值，与当前的结点比较即可。

复杂度：$O(|s|)$

代码：

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<vector>
#include<queue>
#include<stack>
#include<map>
#include<set>
#include<string>
#include<bitset>
#include<sstream>
#include<ctime>
//#include<chrono>
//#include<random>
//#include<unordered_map>
using namespace std;

#define ll long long
#define ls o<<1
#define rs o<<1|1
#define pii pair<int,int>
#define fi first
#define se second
#define pb push_back
#define mp make_pair
#define sz(x) (int)(x).size()
#define all(x) (x).begin(),(x).end()
const double pi=acos(-1.0);
const double eps=1e-6;
const int mod=1e9+7;
const int INF=0x3f3f3f3f;
const int maxn=5e5+5;
ll read(){
	ll x=0,f=1;
	char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
	return x*f;
}
int n,curMax;
char s[maxn];
int nxt[maxn],link[maxn];
vector<int>g[maxn];
pii a[maxn];
void getNext(){
	nxt[1]=0;
	for(int i=2,j=0;i<=n;i++){
		while(j&&s[j+1]!=s[i])j=nxt[j];
		if(s[j+1]==s[i])j++;
		nxt[i]=j;
	}
}
void build(){
	for(int i=1;i<=n;i++){
		g[nxt[i]].pb(i);
	}
}
int dfs(int u){
	if(u==n)return 1;
	for(auto v:g[u]){
		if(dfs(v)){
			link[u]=v;
			return 1;
		}
	}
	return 0;
}
void del(int u){
	if(a[u].fi!=-1){
		a[a[u].fi].se=a[u].se;
	}
	a[a[u].se].fi=a[u].fi;
	if(a[u].fi!=-1){
		curMax=max(curMax,a[u].se-a[u].fi);
	}
	for(auto v:g[u]){
		del(v);
	}
}
int solve(){
	getNext();
	build();
	dfs(0);
	for(int i=1;i<=n;i++){
		if(i==1)a[i].fi=-1;
		else a[i].fi=i-1;
		if(i==n)a[i].se=-1;
		else a[i].se=i+1;
	}
	curMax=1;
	int p=0;
	while(p<n){
		for(auto v:g[p]){
			if(v==link[p])continue;
			del(v);
		}
		p=link[p];
		if(curMax<=p){
			return p;
		}
	}
	return p;
}
int main(void){
	// freopen("in.txt","r",stdin);
	scanf("%s",s+1);
	n=strlen(s+1);
	printf("%d\n",solve());
	return 0;
}