牛客Manacher专题~小G的项链

小G的项链

https://ac.nowcoder.com/acm/problem/14943

题意:给你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;
}
全部评论

相关推荐

评论
点赞
收藏
分享
牛客网
牛客企业服务