manacher——字符串连接（原BZOJ3790）

题目

给你一个目标字符串，你可以：

• 生成所有形式的回文串。
• 把两个字符串连接起来得到一个新字符串，如果一个字符串的后缀和一个字符串的前缀是完全相同的，那么可以将这个重复部分重叠（也可以不重叠）。例如：aba 和 aca 连接起来，可以生成串 abaaca 或 abaca 。

现在给出一个字符串，询问你至少需要连接多少次才能得到这个字符串。

多组数据。

思路

设想一下，如果我们把 aba 和 aca 看做是两条线段，那么当它们连接起来时，像不像两条线段重合了一个单位长度（即一个字符串的后缀和另一个字符串的前缀）后连接成了新的一条更长的线段？

有了这么一个神奇的理解，再回头看看要求：

至少需要连接多少次才能得到这个字符串。

其实就是要求连接出一条与这个字符串等长的线段。

然后显而易见了，这不就是那道贪心算法的“线段覆盖”？

所以我们可以在字符串上跑一遍 $manacher$，以点 $1$ 到 $n$ 为中心可以扩展出 $n$ 个最长回文串，把这 $n$ 个回文串看做是 $n$ 条线段，做一次线段覆盖即可。

细节不多，具体见代码实现：

	#include <bits/stdc++.h>
	using namespace std;
	
	const int N = 5e7 + 10;
	char a[N], ts[N];
	int len[N], n;
	
	struct Edge {
		int l, r;
	} e[N];
	
	inline void manacher() {
		for (int i = 1, maxx = 1, maxn = 1; i <= n; ++ i) {
			if (i <= maxn) len[i] = min(len[(maxx << 1) - i], maxn - i + 1);
			else len[i] = 1;
			while(a[i - len[i]] == a[i + len[i]]) 
				len[i] ++;
			e[i].l = i - len[i] + 1;
			e[i].r = i + len[i] - 1;
			if(i + len[i] - 1 > maxn) {
				maxn = i + len[i] - 1;
				maxx = i;
			}
		}
		return ;
	}
	
	inline bool cmp(Edge a, Edge b) {
		return a.l != b.l ? a.l < b.l : a.r < b.r;
	}
	
	int main() {
		while (scanf("%s", ts + 1) != EOF) {
			a[0] = '$'; a[n = 1] = '#';
			for (int i = 1; i <= strlen(ts + 1); ++ i) {
				a[++ n] = ts[i];
				a[++ n] = '#';
			}
			a[n + 1] = '@';
			manacher();
			sort(e + 1, e + 1 + n, cmp);
			int cnt = 0, mxr = 0, now = 0;
			for (int i = 1; i <= n && now != n; ) {
	            mxr = 0;
	            while (e[i].l <= now + 1)
					mxr = max(mxr, e[i].r), ++ i;
	            -- i;
	            ++ cnt;
	            now = mxr;
	        }
	        printf("%d\n", cnt - 1);
		}
		return 0;
	}