oj.uz

Submission #887826

# Submission time Handle Problem Language Result Execution time Memory
887826 2023-12-15T09:25:58 Z phong Palindromes (APIO14_palindrome) C++17
73 / 100
 265 ms 131072 KB 
#pragma GCC optimize("Ofast")
#pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,avx2,fma")
#pragma GCC optimize("unroll-loops")
#include<bits/stdc++.h>

#define ll long long
const int nmax = 3e5 + 5;
const ll oo = 1e18;
const int lg = 18;
const ll base = 311;
const ll mod = 734065025437;
#define pii pair<int, int>
#define fi first
#define se second
#define debug(a, n) for(int i = 1; i <= n; ++i) cout << a[i] << ' ';
using namespace std;


int n, cur = 0;
string s;
ll pw[nmax];
ll rmq[2][nmax][lg + 1];
int D_odd[nmax];
int D_even[nmax];

void Calc_D_odd() {
    int L = 1;
    int R = 0;
    for(int i = 1 ; i <= n; i++) {
        if(i > R) D_odd[i] = 0;
        else D_odd[i] = min(R - i, D_odd[L + (R - i)]);
        while(i - D_odd[i] - 1 > 0 && i + D_odd[i] + 1 <= n && s[i - D_odd[i] - 1] == s[i + D_odd[i] + 1]) {
            D_odd[i]++;
        }

        if(i + D_odd[i] > R) {
            R = i + D_odd[i];
            L = i - D_odd[i];
        }
    }
}

void Calc_D_even() {
    int L = 1;
    int R = 0;
    for(int i = 1 ; i < n; i++) {
        int j = i + 1;
        if(j > R) D_even[i] = 0;
        else D_even[i] = min(R - j + 1, D_even[L + (R - j)]);
        while(i - D_even[i] > 0 && j + D_even[i] <= n && s[i - D_even[i]] == s[j + D_even[i]]) {
            D_even[i]++;
        }
        if(i + D_even[i] > R) {
            R = i + D_even[i];
            L = j - D_even[i];
        }
    }
}

ll combine(ll a, ll b, int x){
    return (a * pw[x] + b) % mod;
}
void Init(){
    pw[0] = 1;
    for(int i = 1, j; i <= n; ++i){
        pw[i] = pw[i - 1] * base % mod;
        rmq[0][i][0] = s[i] - 'a' + 3;
        rmq[1][i][0] = s[n - i + 1] - 'a' + 3;
    }
    for(int j = 1; j <= lg; ++j){
        for(int i = 1; i + (1 << j) - 1 <= n; ++i){
            for(int t = 0; t < 2; ++t){
                rmq[t][i][j] = combine(rmq[t][i][j - 1], rmq[t][i + (1 << (j - 1))][j - 1], 1 << (j - 1));
            }
        }
    }
}
ll get(int l, int r, int t){
    int k = r - l + 1;
    ll res = 0;
    for(int j = 0; j <= __lg(k); ++j){
        if(k >> j & 1){
            res = combine(res, rmq[t][l][j], 1 << j);
            l += 1 << j;
        }
    }
    return res;
}

bool check(int l, int r){
    return get(l, r, 0) == get(n - r + 1, n - l + 1, 1);
}

map<ll, bool> mp;
map<ll, int>  ind;
int cnt[nmax];
ll f[nmax];
vector<int> tmp;
vector<int> adj[nmax];

void update(int l, int r){
    bool ok = 1;
    while(l <= r){
        ll x = get(l, r, 0);
        if(!mp[x]){
            ind[x] = ++cur;
            tmp.push_back(cur);
            if(ok){
                f[cur]++;
                ok = 0;
            }
            cnt[cur] = r - l + 1;
            mp[x] = 1;

        }
        else{
            int id = ind[x];
            tmp.push_back(id);
            if(ok) f[id]++;

            break;
        }
        ++l, --r;
    }
    if(l > r) tmp.push_back(0);
    int x = (int)tmp.size() - 1;
    for(int i = 0; i < x; ++i){
        int x = tmp[i], y = tmp[i + 1];
        adj[y].push_back(x);
    }
    tmp.clear();
}
void dfs(int u){
    for(auto v : adj[u]){
        dfs(v);
        f[u] += f[v];
    }
}
main(){
    ios_base::sync_with_stdio(0);
    cin.tie(0);cout.tie(0);
//    freopen("code.inp", "r", stdin);
//    freopen("code.out", "w", stdout);
    cin >> s;
    n = s.size(); s = ' ' + s + ' ';
    Init();
    Calc_D_odd();
    Calc_D_even();
    for(int i = 1; i <= n; ++i){
        update(i - D_odd[i], i + D_odd[i]);
        if(i < n) update(i - D_even[i] + 1, i + D_even[i]);
    }
    dfs(0);

    ll ans = 0;
    for(int i = 1; i <= cur; ++i){
        ans = max(ans, 1ll * f[i] * cnt[i]);
    }
    cout << ans;
}
/*

*/

Compilation message

palindrome.cpp: In function 'void Init()':
palindrome.cpp:65:20: warning: unused variable 'j' [-Wunused-variable]
   65 |     for(int i = 1, j; i <= n; ++i){
      |                    ^
palindrome.cpp: At global scope:
palindrome.cpp:139:1: warning: ISO C++ forbids declaration of 'main' with no type [-Wreturn-type]
  139 | main(){
      | ^~~~

Subtask #1
8.0 / 8.0

# Verdict Execution time Memory Grader output
1 Correct 5 ms 19036 KB Output is correct
2 Correct 4 ms 19036 KB Output is correct
3 Correct 4 ms 19036 KB Output is correct
4 Correct 4 ms 19088 KB Output is correct
5 Correct 3 ms 18844 KB Output is correct
6 Correct 4 ms 19036 KB Output is correct
7 Correct 4 ms 19036 KB Output is correct
8 Correct 4 ms 19036 KB Output is correct
9 Correct 4 ms 19036 KB Output is correct
10 Correct 4 ms 19036 KB Output is correct
11 Correct 3 ms 19036 KB Output is correct
12 Correct 4 ms 19036 KB Output is correct
13 Correct 4 ms 19040 KB Output is correct
14 Correct 3 ms 19036 KB Output is correct
15 Correct 3 ms 19036 KB Output is correct
16 Correct 4 ms 19036 KB Output is correct
17 Correct 4 ms 19036 KB Output is correct
18 Correct 4 ms 19036 KB Output is correct
19 Correct 4 ms 19036 KB Output is correct
20 Correct 3 ms 19076 KB Output is correct
21 Correct 4 ms 19036 KB Output is correct
22 Correct 3 ms 19036 KB Output is correct
23 Correct 4 ms 19104 KB Output is correct
24 Correct 4 ms 19036 KB Output is correct
25 Correct 4 ms 19036 KB Output is correct
26 Correct 4 ms 19096 KB Output is correct
27 Correct 4 ms 19036 KB Output is correct
28 Correct 3 ms 18864 KB Output is correct
29 Correct 3 ms 19032 KB Output is correct
30 Correct 4 ms 19036 KB Output is correct
31 Correct 4 ms 19036 KB Output is correct
32 Correct 4 ms 19032 KB Output is correct

Subtask #2
15.0 / 15.0

# Verdict Execution time Memory Grader output
1 Correct 4 ms 21084 KB Output is correct
2 Correct 5 ms 21084 KB Output is correct
3 Correct 4 ms 21140 KB Output is correct
4 Correct 5 ms 21084 KB Output is correct
5 Correct 5 ms 21084 KB Output is correct
6 Correct 4 ms 21080 KB Output is correct
7 Correct 5 ms 21336 KB Output is correct
8 Correct 5 ms 21084 KB Output is correct
9 Correct 5 ms 21336 KB Output is correct
10 Correct 4 ms 21244 KB Output is correct
11 Correct 6 ms 21084 KB Output is correct
12 Correct 5 ms 21084 KB Output is correct

Subtask #3
24.0 / 24.0

# Verdict Execution time Memory Grader output
1 Correct 16 ms 22620 KB Output is correct
2 Correct 14 ms 22820 KB Output is correct
3 Correct 15 ms 22620 KB Output is correct
4 Correct 18 ms 22620 KB Output is correct
5 Correct 11 ms 22620 KB Output is correct
6 Correct 14 ms 22620 KB Output is correct
7 Correct 15 ms 22620 KB Output is correct
8 Correct 6 ms 21088 KB Output is correct
9 Correct 6 ms 21084 KB Output is correct
10 Correct 10 ms 22108 KB Output is correct

Subtask #4
26.0 / 26.0

# Verdict Execution time Memory Grader output
1 Correct 200 ms 67252 KB Output is correct
2 Correct 178 ms 66412 KB Output is correct
3 Correct 265 ms 67136 KB Output is correct
4 Correct 219 ms 66644 KB Output is correct
5 Correct 135 ms 66784 KB Output is correct
6 Correct 109 ms 61992 KB Output is correct
7 Correct 162 ms 63480 KB Output is correct
8 Correct 29 ms 50168 KB Output is correct
9 Correct 63 ms 54008 KB Output is correct
10 Correct 118 ms 64424 KB Output is correct

Subtask #5
0 / 27.0

# Verdict Execution time Memory Grader output
1 Runtime error 259 ms 131072 KB Execution killed with signal 9
2 Halted 0 ms 0 KB -