洛谷4287 BZOJ2342 SHOI2011 双倍回文 manacher set

题目链接

题意：
给你一个长度为$n$的字符串，问你长度最长的子串，要求它形如$A{A}^{R}A{A}^R$，其中$A^R$表示把$A$这个串反过来之后的串，并且$A$与$A^R$不能有公用的字符，换句话说，就是最后的答案串$A{A}^{R}A{A}^R$的长度一定是4的倍数。$n<=5e5$。

题解：
一看就是与回文串相关的问题。然而我现在还不会回文树和回文自动机，那么只能考虑用manacher来做。不难看出，要求的串$A{A}^{R}A{A}^R$是一个回文串，并且$A{A}^R$也是一个回文串。

我们可以用manacher做的事情是$O(n)$求出以每一个字符为对称中心的最长回文子串的长度，那么我们先对原串做一遍manacher。我们不难发现，合法的对称轴应该是在两个字符之间的空隙，并且$A$与$A^R$的对称轴也是这样。接下来考虑怎么求答案。

我们设$x$为整个串的对称轴，设$y$为后半个串的对称轴位置，设$r[i]$为$i$这个位置的最长回文半径。我们知道，要让$y$可以作为$x$的后半个串的对称轴，那么要满足$x+r[x]/2>=y$，并且$y-r[y]<=x$。这个应该并不难理解，也就是你可以回文半径能覆盖的更长不要紧，但是要是覆盖不过来就不行。此时我们找到一组$x,y$，那么可以尝试将答案更新为$(y-x)\ast 4$。我们对于每一个$x$，其实只要找到满足$y-r[y]>=x$的最大的一个$y$就可以了。于是我们考虑把$y$按照$y-r[y]$排序，然后枚举$x$，看看在当前$x$应该从小到大插入哪些$y$。我们把这些$y$用一个set维护，这样我们每次更新答案时就只需要在set中询问不超过$x+r[x]/2$的最大的$y$就可以了，这样我们就可以用$(y-x)\ast 4$来尝试更新答案。对于那些已经不满足$y-r[y]>=x$的$y$，我们其实并不用从set中删去，因为很明显它们一定不会是最优的了。

这样就做完这个题了。复杂度$O(nlogn)$。似乎还有复杂度更优秀的并查集做法和线性做法。

代码：

#include <bits/stdc++.h>
using namespace std;

int n,r[2000010],ans,a[2000010];
char t[2000010];
set<int> s;
inline void manacher()
{
	int mx=0,pos=0;
	for(int i=1;i<=n;++i)
	{
		if(mx>i)
		r[i]=min(mx-i,r[pos*2-i]);
		else
		r[i]=0;
		while(r[i]+i+1<=n&&i-r[i]>=1&&t[r[i]+i+1]==t[i-r[i]])
		++r[i];
		if(r[i]+i>mx)
		{
			mx=r[i]+i;
			pos=i;
		}
	}
}
inline int cmp(int x,int y)
{
	return x-r[x]<y-r[y];
}
int main()
{
	scanf("%d",&n);
	scanf("%s",t+1);
	manacher();
	for(int i=1;i<=n;++i)
	a[i]=i;
	sort(a+1,a+n+1,cmp);
	int cur=0;
	for(int i=1;i<=n;++i)
	{
		while(a[cur]-r[a[cur]]<=i&&cur<=n)
		{
			s.insert(a[cur]);
			++cur;
		}
		set<int>::iterator it;
		it=s.upper_bound(i+r[i]/2);
		if(it!=s.begin())
		ans=max(ans,(*(--it)-i)*4);
	}
	printf("%d\n",ans);
	return 0;
}