UVALive 4975 (LA 4975) Casting Spells Manacher + Set维护

题目大意：

就是对于一个串, 长度不超过3*10^5, 询问其中最长的形似是ww^ww^类型的串, 其中w^是w反转得到的串

大致思路：

好早以前写的题了....当时交UVALive一直返回WA一直找不到错哪里了...现在回去看发现Rejudge成AC了....

就是先Manacher处理出每个位置的回文半径

从左至右扫一遍用set维护当前能覆盖到的右界, 然后枚举后半的ww^的中心应该选取的最左的左中心即可

代码如下：

Result  :  Accepted     Memory  :  ? KB     Time  :  1996 ms

/*
* Author: Gatevin
* Created Time: 2015/5/29 10:02:17
* File Name: Rin_Tohsaka.cpp
*/
#include<iostream>
#include<sstream>
#include<fstream>
#include<vector>
#include<list>
#include<deque>
#include<queue>
#include<stack>
#include<map>
#include<set>
#include<bitset>
#include<algorithm>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cctype>
#include<cmath>
#include<ctime>
#include<iomanip>
using namespace std;
const double eps(1e-8);
typedef long long lint;
#define foreach(e, x) for(__typeof(x.begin()) e = x.begin(); e != x.end(); ++e)
#define SHOW_MEMORY(x) cout<<sizeof(x)/(1024*1024.)<<"MB"<<endl

char in[300010];
char s[600010];
int R[600010];

void Manacher(char *s, int *R, int len)
{
int p = 0, mx = 0;
R[0] = 1;
for(int i = 1; i <= len; i++)
{
if(mx > i) R[i] = min(R[2*p - i], mx - i);
else R[i] = 1;
while(s[i - R[i]] == s[i + R[i]])
R[i]++;
if(i + R[i] > mx)
p = i, mx = i + R[i];
}
return;
}

set<pair<int, int> > S;
set<pair<int, int> > S2;
vector < pair<int, int> > V;
int main()
{
int Z;
scanf("%d", &Z);
while(Z--)
{
scanf("%s", in);
int len = strlen(in);
s[0] = '@';
for(int i = 0; i < len; i++)
s[2*i + 1] = in[i], s[2*i + 2] = '#';
s[2*len] = '$';
Manacher(s, R, 2*len);
/*
* 枚举ww^ww^的中心位置为i,
* 那么对于i是最优解的位置ww^的对称中心j满足覆盖半径R[j]: j + R[j] - 1 >= i - 1, j肯定也是偶数
*/
//for(int i = 0; i <= 2*len; i++)
// cout<<R[i]<<endl;
S.clear();
S2.clear();
int ans = 0;
for(int i = 2; i < 2*len; i += 2)
{
V.clear();
set<pair<int, int> > :: iterator it;
for(it = S.begin(); it != S.end(); it++)
if((*it).first < i - 1) V.push_back(make_pair((*it).first, (*it).second));
else break;
for(int j = 0, sz = V.size(); j < sz; j++)
S.erase(V[j]), S2.erase(make_pair(V[j].second, V[j].first));
//2*y - i + 1 >= i - R[i] + 1
//即 y >= (2*i - R[i]) / 2
it = S2.lower_bound(make_pair(((2*i - R[i]) >> 1) + (R[i] & 1), 0));
if(it != S2.end())
ans = max(ans, (i - (*it).first) << 1);
S.insert(make_pair(i + R[i] - 1, i));
S2.insert(make_pair(i, i + R[i] - 1));
}
printf("%d\n", ans);
}
return 0;
}