题意:给你n个数,你可以将这n个数分成任意个长度相等个的区间,每个区间的权值为区间内数的异或和,要求分割后的区间的权值组成回文串,问你回文串最大的长度是多长,可以认为这n个数组成了一个环,第1个数和第n个数是相邻的
题解:因为要将这n个数分成任意个长度相等的区间,所以k是n的因子,假设我们现在将这n个数分成m个区间了,并且已经知道这m个区间的起点的原来n个数中的第z个,那么我们得到了m个权值
由异或的性质我们可以知道,我们可以用一个前缀异或和来很快(O(1))的得到这些权值
我们将得到的权值倍增一次,对于倍增后的数组进行Manacher后,如果得到的p数组中出现了长度为区间长度的p那么我们就可以知道,我们将区间分割为m个是可以得到一个回文串的
所以不断的枚举n的因子,更新答案,即可得到回文串的最大长度
/**
* ┏┓ ┏┓
* ┏┛┗━━━━━━━┛┗━┓
* ┃ ┃
* ┃ ━ ┃
* ┃ > < ┃
* ┃ ┃
* ┃... ⌒ ... ┃
* ┃ ┃
* ┗━┓ ┏━┛
* ┃ ┃ Code is far away from bug with the animal protecting
* ┃ ┃ 神兽保佑,代码无bug
* ┃ ┃
* ┃ ┃
* ┃ ┃
* ┃ ┃
* ┃ ┗━━━┓
* ┃ ┣┓
* ┃ ┏┛
* ┗┓┓┏━┳┓┏┛
* ┃┫┫ ┃┫┫
* ┗┻┛ ┗┻┛
*/
// warm heart, wagging tail,and a smile just for you!
//
// _ooOoo_
// o8888888o
// 88" . "88
// (| -_- |)
// O\ = /O
// ____/`---'\____
// .' \| |// `.
// / \||| : |||// \
// / _||||| -:- |||||- \
// | | \ - /// | |
// | \_| ''\---/'' | |
// \ .-\__ `-` ___/-. /
// ___`. .' /--.--\ `. . __
// ."" '< `.___\_<|>_/___.' >'"".
// | | : `- \`.;`\ _ /`;.`/ - ` : | |
// \ \ `-. \_ __\ /__ _/ .-` / /
// ======`-.____`-.___\_____/___.-`____.-'======
// `=---='
// ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
// 佛祖保佑 永无BUG
#include <set>
#include <map>
#include <stack>
#include <cmath>
#include <queue>
#include <cstdio>
#include <string>
#include <vector>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long LL;
typedef pair<int, int> pii;
typedef unsigned long long uLL;
#define ls rt<<1
#define rs rt<<1|1
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define bug printf("*********\n")
#define FIN freopen("input.txt","r",stdin);
#define FON freopen("output.txt","w+",stdout);
#define IO ios::sync_with_stdio(false),cin.tie(0)
#define debug1(x) cout<<"["<<#x<<" "<<(x)<<"]\n"
#define debug2(x,y) cout<<"["<<#x<<" "<<(x)<<" "<<#y<<" "<<(y)<<"]\n"
#define debug3(x,y,z) cout<<"["<<#x<<" "<<(x)<<" "<<#y<<" "<<(y)<<" "<<#z<<" "<<z<<"]\n"
const int maxn = 1e6 + 5;
const int INF = 0x3f3f3f3f;
const int mod = 1e9 + 7;
const double Pi = acos(-1);
LL gcd(LL a, LL b) {
return b ? gcd(b, a % b) : a;
}
LL lcm(LL a, LL b) {
return a / gcd(a, b) * b;
}
double dpow(double a, LL b) {
double ans = 1.0;
while(b) {
if(b % 2)ans = ans * a;
a = a * a;
b /= 2;
} return ans;
}
LL quick_pow(LL x, LL y) {
LL ans = 1;
while(y) {
if(y & 1) {
ans = ans * x % mod;
} x = x * x % mod;
y >>= 1;
} return ans;
}
int n;
int a[maxn << 1];
int b[maxn << 1];
int p[maxn << 1];
int s[maxn << 1];
int pre[maxn << 1];
void manacher(int *str, int n) {
int m = 0;
s[0] = -1;
for(int i = 0; i < n; i++) {
s[++m] = -1;
s[++m] = str[i];
}
s[++m] = 0;
// debug2(n, m);
// for(int i = 1; i <= m; i++) {
// printf("%d ", s[i]);
// }
// printf("\n");
for(int i = 0; i <= m; i++) p[i] = 0;
int mx = 0, id = 0;
for(int i = 1; i <= m; i++) {
if(mx > i) p[i] = min(p[2 * id - i], mx - i);
else p[i] = 1;
while(i - p[i] > 0 && s[i - p[i]] == s[i + p[i]]) p[i]++;
if(i + p[i] > mx) {
mx = i + p[i];
id = i;
}
}
}
bool check(int len) {
int num = n / len;
for(int i = 0; i < len; i++) {
//枚举起点
for(int j = 0; j < num; j++) {
b[j] = pre[(j + 1) * len + i] ^ pre[j * len + i];//
}
for(int i = 0; i < num; ++i) {
b[i + num] = b[i];
// b[i + num + num] = b[i];
}
manacher(b, num * 2);
for(int i = 1; i <= 4 * num; ++i) {
if(p[i] >= num) return true;
}
}
return false;
}
int main() {
#ifndef ONLINE_JUDGE
FIN
#endif
scanf("%d", &n);
for(int i = 1; i <= n; i++) {
scanf("%d", &a[i]);
pre[i] = pre[i - 1] ^ a[i];
}
for(int i = n + 1; i <= n + n; i++) {
pre[i] = pre[i - 1] ^ a[i - n];
}
int ans = 1;
for(int i = 1; i * i <= n; i++) {
if(n % i) continue;
if(check(n / i)) ans = max(ans, i);
if(check(i)) ans = max(ans, n / i);
}
printf("%d\n", ans);
return 0;
}
京公网安备 11010502036488号