Submission #970945

#TimeUsernameProblemLanguageResultExecution timeMemory
97094512345678밀림 점프 (APIO21_jumps)C++17
48 / 100
876 ms73332 KiB
#include "jumps.h"
#include <bits/stdc++.h>

using namespace std;

const int nx=2e5+5, kx=18;

int l[nx], r[nx], deg[nx], lvl[nx], p[kx][nx], pa[nx][kx], rt, mn[nx], mx[nx];
vector<int> d[nx], g[nx];

void dfs(int u, int p)
{
    //cout<<"dfs "<<u<<'\n';
    mn[u]=mx[u]=u;
    pa[u][0]=p;
    for (int i=1; i<kx; i++) pa[u][i]=pa[pa[u][i-1]][i-1];
    for (auto v:g[u]) dfs(v, u), mn[u]=min(mn[u], mn[v]), mx[u]=max(mx[u], mx[v]);
}

int lca(int u, int v)
{
    if (lvl[u]>lvl[v]) swap(u, v);
    for (int i=kx-1; i>=0; i--) if (lvl[pa[v][i]]>=lvl[u]) v=pa[v][i];
    if (u==v) return v;
    for (int i=kx-1; i>=0; i--) if (pa[u][i]!=pa[v][i]) u=pa[u][i], v=pa[v][i];
    return pa[u][0];
}

void init(int N, std::vector<int> H) {
    stack<int> s;
    for (int i=0; i<N; i++) l[i]=r[i]=N;
    for (int i=0; i<N; i++)
    {
        while (!s.empty()&&H[s.top()]<H[i]) s.pop();
        if (!s.empty()) l[i]=s.top(), d[s.top()].push_back(i), deg[i]++;
        s.push(i);
    }
    while (!s.empty()) s.pop();
    for (int i=N-1; i>=0; i--)
    {
        while (!s.empty()&&H[s.top()]<H[i]) s.pop();
        if (!s.empty()) r[i]=s.top(), d[s.top()].push_back(i), deg[i]++;
        s.push(i);
    }
    queue<int> q;
    for (int i=0; i<N; i++) if (l[i]==N&&r[i]==N) q.push(i), rt=i;
    while (!q.empty())
    {
        auto u=q.front();
        q.pop();
        for (auto x:d[u]) if (--deg[x]==0) q.push(x);
        lvl[u]=max(lvl[l[u]], lvl[r[u]])+1;
    }
    for (int i=0; i<N; i++)
    {
        if (lvl[l[i]]>lvl[r[i]]) g[l[i]].push_back(i);
        else g[r[i]].push_back(i);
    }
    for (int i=0; i<kx; i++) p[i][N]=N;
    for (int i=0; i<N; i++) p[0][i]=(lvl[l[i]]<lvl[r[i]])?l[i]:r[i];
    for (int i=1; i<kx; i++) for (int j=0; j<N; j++) p[i][j]=p[i-1][p[i-1][j]];
    dfs(rt, rt);
    for (int i=0; i<N; i++) sort(g[i].begin(), g[i].end());
}

int minimum_jumps(int A, int B, int C, int D) {
    int rt=lca(C, D);
    if (B<mn[rt]) return -1;
    for (int i=C; i<=D; i++)
    {
        if (B>=mn[i])
        {
            C=i;
            break;
        }
    }
    A=max(A, mn[C]);
    A=lca(A, B);
    int res=0;
    for (int i=kx-1; i>=0; i--) if (lvl[p[i][A]]>=lvl[C]) res+=(1<<i), A=p[i][A];
    return res+lvl[A]-lvl[C];
}
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...