BZOJ2342【manacher】【set】

/* I will wait for you*/

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <ctime>
#include <algorithm>
#include <iostream>
#include <fstream>
#include <vector>
#include <queue>
#include <deque>
#include <map>
#include <set>
#include <string>
#define make make_pair
#define fi first
#define se second

using namespace std;

typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> pii;

const int maxn = 1000010;
const int maxm = 1010;
const int maxs = 26;
const int inf = 0x3f3f3f3f;
const int P = 1000000007;
const double error = 1e-9;

inline int read()
{
	int x = 0, f = 1; char ch = getchar();
	while (ch <= 47 || ch >= 58)
		f = (ch == 45 ? -1 : 1), ch = getchar();
	while (ch >= 48 && ch <= 57)
		x = x * 10 + ch - 48, ch = getchar();
	return x * f;
}

int n, pos, ans, mx[maxn], y[maxn];
char s[maxn]; set<int> m;

void manacher()
{  
	for (int i = 0; i <= n; i++) {
		if (pos + mx[pos] >= i)
			mx[i] = min(mx[2 * pos - i], pos + mx[pos] - i);
		for(;i - mx[i] - 1 >= 0 && i + mx[i] + 1 <= n && 
		     s[i - mx[i] -1 ] == s[i + mx[i] + 1]; mx[i]++);
		if (i + mx[i] > pos + mx[pos])
			pos = i;
	}
}			

bool cmp (int a, int b)
{
	return a - mx[a] < b - mx[b];
}

int main()
{
	n = read(), scanf("%s", s + 1);

	for (int i = n; i >= 1; i--)
		s[2 * i - 1] = s[i];
	for (int i = 0; i <= 2 * n; i += 2)
		s[i] = '#';
	n *= 2, manacher(), n /= 2;

	for (int i = 1; i <= n; i++)
		y[i] = i, mx[i] = mx[2 * i] / 2;
	sort(y + 1, y + 1 + n, cmp);

	for (int i = 1, j = 1; i <= n; i++) {
		for (; j <= n && y[j] - mx[y[j]] <= i; j++)
			m.insert(y[j]);
		set<int> :: iterator pos = m.upper_bound(i + mx[i] / 2);
		if (pos != m.begin())
			ans = max(ans, 4 * (*--pos - i));
	}
	
	printf("%d\n", ans);		

	return 0;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值