bzoj2565 最长双回文串

最长双回文串

题目背景：

BZOJ - 2565

分析：据说用回文自动机很裸，据说回文自动机做这个题很简单，宝宝这种蒟蒻不会，还是就manacher吧，记录下信息即可，详见代码注释。

Source：

#include 
     
     
      
      
#include 
      
      
       
       
#include 
       
       
        
        
#include 
        
        
         
         

using namespace std;

const int N = 1e6 + 5;
int n, m, T, p[N * 2];
char s[N * 2], t[N];
int dp[N][2];

void manacher() {
	memset(dp, 0, sizeof(dp));
 	int id = 0, pos = 0, x = 0;
	for (int i = 1; i <= n; ++i) {
		if (pos > i) x = min(p[id * 2 - i], pos - i);
		else x = 1;
		while (s[i - x] == s[i + x]) ++x;
		if (i + x > pos) {
			/*因为寻找到的id的位置显然是往有递增的，
			所以对于每一个pos所能找到的最左边的id
			一定就是第一个覆盖当前的，也就是原来的
			pos + 1到现在的pos之间*/
			for(int j = pos + 1; j <= i + x; ++j) dp[j][0] = i; 
			pos = i + x, id = i;
		}
		p[i] = x;
	}
	p[n + 1] = 1;
	id = pos = n + 1;
	for(int i = n; i >= 1; --i) {
		if(pos < i) x = min(p[id * 2 - i], i - pos);
		else x = 1;
		while(s[i - x] == s[i + x]) ++x;
		if(i - x < pos) {
		/*同上理解，注意方向*/ 
			for(int j = pos - 1; j >= i - x; --j) dp[j][1] = i;
			pos = i - x, id = i;
		}
		p[i] = x;
	}
}

int main() {
	cin >> t + 1;
	m = strlen(t + 1);
	s[n = 0] = '!';
	for (int i = 1; i <= m; ++i) {
		s[++n] = '#';
		s[++n] = t[i];
	}
	s[++n] = '#';
	s[n + 1] = ' ';
	manacher();
	int ans = 0, len;
	for (int i = 3; i < 2 * m; i += 2)
		len = dp[i][1] - dp[i][0], ans = max(ans, len);
		/*考虑左边的回文串长度是i - dp[i][0] + 1, 
		右边的回文串长度是dp[i][1] - i - 1(分情况讨论可得) 
		两边的长度一共就是dp[i][1] - dp[i][0]*/ 
	cout << ans << '\n';
	return 0;
}