Submission #979219

# Submission time Handle Problem Language Result Execution time Memory
979219 2024-05-10T11:54:05 Z sleepntsheep Rainforest Jumps (APIO21_jumps) C++17
23 / 100
780 ms 22012 KB
#pragma GCC optimize("O2,unroll-loops")
#include "jumps.h"

#include <vector>
#include <algorithm>

int subtask = -1;

#define MAX_N 200000
#define MAX_M (2*MAX_N)

int N, head[MAX_N], nxt[MAX_M], vv[MAX_M], ng[MAX_N], pg[MAX_N], HI[MAX_N], H[MAX_N];
int dd[MAX_N], pp[MAX_N], jj[MAX_N], dd_[MAX_N], pp_[MAX_N], jj_[MAX_N];
int dd1[MAX_N],pp1[MAX_N],jj1[MAX_N],ee1[MAX_N];
int dd2[MAX_N],pp2[MAX_N],jj2[MAX_N],ee2[MAX_N];

/*
 * dd - high link
 * dd_ - low link
 * dd1 - left link
 * dd2 - right link */

void link(int u, int v)
{
    static int i = 1;
    nxt[i] = head[u];
    vv[i] = v;
    head[u] = i++;
}

void binlift_connect(int*dd,int*pp,int*jj,int p,int u)
{
    dd[u] = dd[p] + 1;
    pp[u] = p;
    jj[u] = dd[p] - dd[jj[p]] == dd[jj[p]] - dd[jj[jj[p]]] ? jj[jj[p]] : p;
}

void init(int N, std::vector<int> H) {
    for(auto&x:H)--x;
    ::N = N;
    if (N <= 200) subtask = 2;
    else if (N <= 2000) subtask = 3;
    int flag = 1;
    for (int i = 1; i < N; ++i) if (H[i] != H[i-1] + 1) flag = 0;
    if(flag) subtask = 1;

    static int stk[MAX_N], so;
    so = 0;

    for (int i = 0; i < N; stk[so++] = i++) while (so && H[stk[so-1]] < H[i]) ng[stk[so-1]]=i, link(stk[--so], i);
    while(so)ng[stk[--so]]=-1;
    for (int i = N - 1; i >= 0; stk[so++] = i--) while (so && H[stk[so-1]] < H[i]) pg[stk[so-1]]=i, link(stk[--so], i);
    while(so)pg[stk[--so]]=-1;

    for (int i = 0; i < N; ++i) HI[::H[i]=H[i]] = i;

    for (int j = N; j--;)
    {
        int i = HI[j];
        int high = (pg[i] == -1) ? ng[i] : (ng[i] == -1 ? pg[i] : (H[ng[i]] > H[pg[i]] ? ng[i] : pg[i]));
        int low = pg[i] + ng[i] - high;
        if (high != -1) binlift_connect(dd, pp, jj, high, i);
        else pp[i] = jj[i] = i;
        if (low != -1) binlift_connect(dd_, pp_, jj_, low, i);
        else pp_[i] = jj_[i] = i;
    }

    if (subtask == -1) subtask = 4;

    for(int i=0;i<N;++i)
    {
        /* left */
        if(pg[i]==-1) { pp1[i]=jj1[i]=i,dd1[i]=0; }
        else{
            binlift_connect(dd1,pp1,jj1,pg[i],i);
            if(jj1[i]==pp1[i])ee1[i]=H[pg[i]];
            else ee1[i]=std::max(ee1[jj1[pp1[i]]],ee1[pp1[i]]);
        }

        /* right */
        if(ng[i]==-1) { pp2[i]=jj2[i]=i,dd2[i]=0; }
        else {
            binlift_connect(dd2,pp2,jj2,ng[i],i);
            if(jj2[i]==pp2[i])ee2[i]=H[ng[i]];
            else ee2[i]=std::max(ee2[jj2[pp2[i]]],ee2[pp2[i]]);
        }
    }
    subtask=7;
}

int Q_COUNTER = 0;

int minimum_jumps(int A, int B, int C, int D) {
    ++Q_COUNTER;

    if (subtask == 1) return (A>D)?-1:(B>=C?0:C-B);

    if (subtask == 4 && Q_COUNTER > 5) subtask = 5;

    if(subtask <= 4)
    {
        static int d[MAX_N], qu[MAX_N], qh=0, qt=0,z=1e9;
        qh=qt=0,z=1e9;
        for(int j=0;j<N;++j)d[j]=1e9;
        for (int i = A; i <= B; ++i) d[qu[qh++]=i]=0;
        while(qh-qt) for(int u=qu[qt++],j=head[u];j;j=nxt[j]) if(d[vv[j]]>d[u]+1)d[vv[j]]=d[u]+1,qu[qh++]=vv[j];
        for (int j = C; j <= D; ++j)if(d[j]!=1e9&&d[j]<z)z=d[j];
        return(z<1e9?z:-1);
    }

    if (subtask == 5 && A != B) subtask = 6;
    if (subtask == 6 && C != D) subtask = 7;

    if (subtask == 6)
    {
        if(H[B]>H[C])return -1;
        int u=B;
        for(;pp1[u]!=u;)
        {
            if(jj1[u]>=A&&ee1[u]<H[C])u=jj1[u];
            else if(pp1[u]>=A&&H[pp1[u]]<H[C])u=pp1[u];
            else break;
        }
        A=u;
    }

    if (subtask == 5 || subtask == 6)
    {
        if(A==-1)return -1;
        int u=A,u1;
        for(;pp[u]!=u;)
        {
            if (H[jj[u]]<=H[C])u=jj[u];
            else if(H[pp[u]]<=H[C])u=pp[u];
            else break;
        }
        u1=u;
        for(;pp_[u]!=u;)
        {
            if(H[jj_[u]]<=H[C])u=jj_[u];
            else if(H[pp_[u]]<=H[C])u=pp_[u];
            else break;
        }
        if(u!=C)return -1;
        return dd_[u1]-dd_[u]+dd[A]-dd[u1];
    }

    if (subtask == 7)
    {
        int highest_left;
        {
            int u=B;
            for(;pp1[u] != u;)
            {
                if(jj1[u]>=A)u=jj1[u];
                else if(pp1[u]>=A)u=pp1[u];
                else break;
            }
            highest_left=u;
        }


        /* the right endpoint is a position that is at lower_bound(highest_left) */
        {
            int u=C;
            for(;pp2[u]!=u;)
            {
                if(jj2[u]<=D&&ee2[u]<H[highest_left])u=jj2[u];
                else if(pp2[u]<=D&&H[pp2[u]]<H[highest_left])u=pp2[u];
                else break;
            }
            if (u!=pp2[u]&&pp2[u]<=D)u=pp2[u];
            C=u;
        }

        /* select highest left endpoint possible */
        {
            int u=B;
            for(;pp1[u]!=u;)
            {
                if(jj1[u]>=A&&ee1[u]<H[C])u=jj1[u];
                else if(pp1[u]>=A&&H[pp1[u]]<H[C])u=pp1[u];
                else break;
            }
        }

        {
            if(A==-1)return -1;
            int u=A,u1;
            for(;pp[u]!=u;)
            {
                if (H[jj[u]]<=H[C])u=jj[u];
                else if(H[pp[u]]<=H[C])u=pp[u];
                else break;
            }
            u1=u;
            for(;pp_[u]!=u;)
            {
                if(H[jj_[u]]<=H[C])u=jj_[u];
                else if(H[pp_[u]]<=H[C])u=pp_[u];
                else break;
            }
            if(u!=C)return -1;
            return dd_[u1]-dd_[u]+dd[A]-dd[u1];
        }
    }

    return -1;
}

# Verdict Execution time Memory Grader output
1 Correct 2 ms 14680 KB Output is correct
2 Correct 2 ms 14680 KB Output is correct
3 Incorrect 85 ms 20008 KB Output isn't correct
4 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Incorrect 2 ms 14680 KB Output isn't correct
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Incorrect 2 ms 14680 KB Output isn't correct
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Incorrect 2 ms 14680 KB Output isn't correct
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 2 ms 14676 KB Output is correct
2 Correct 2 ms 16728 KB Output is correct
3 Correct 2 ms 16728 KB Output is correct
4 Correct 179 ms 17688 KB Output is correct
5 Correct 766 ms 21440 KB Output is correct
6 Correct 430 ms 16212 KB Output is correct
7 Correct 540 ms 21528 KB Output is correct
8 Correct 326 ms 16984 KB Output is correct
9 Correct 674 ms 21300 KB Output is correct
10 Correct 780 ms 21180 KB Output is correct
11 Correct 730 ms 21696 KB Output is correct
12 Correct 734 ms 21188 KB Output is correct
13 Correct 733 ms 21584 KB Output is correct
14 Correct 730 ms 21692 KB Output is correct
15 Correct 671 ms 21592 KB Output is correct
16 Correct 672 ms 21608 KB Output is correct
17 Correct 2 ms 14680 KB Output is correct
18 Correct 2 ms 14932 KB Output is correct
19 Correct 3 ms 14680 KB Output is correct
20 Correct 2 ms 14680 KB Output is correct
21 Correct 4 ms 14680 KB Output is correct
22 Correct 3 ms 14680 KB Output is correct
23 Correct 3 ms 14680 KB Output is correct
24 Correct 3 ms 16728 KB Output is correct
25 Correct 2 ms 14680 KB Output is correct
26 Correct 2 ms 14680 KB Output is correct
27 Correct 12 ms 14680 KB Output is correct
28 Correct 15 ms 14680 KB Output is correct
29 Correct 12 ms 14680 KB Output is correct
30 Correct 11 ms 14680 KB Output is correct
31 Correct 13 ms 16728 KB Output is correct
32 Correct 2 ms 14680 KB Output is correct
33 Correct 17 ms 18512 KB Output is correct
34 Correct 43 ms 21232 KB Output is correct
35 Correct 23 ms 21184 KB Output is correct
36 Correct 32 ms 21424 KB Output is correct
37 Correct 25 ms 22012 KB Output is correct
38 Correct 26 ms 21672 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 2 ms 14676 KB Output is correct
2 Correct 2 ms 16728 KB Output is correct
3 Correct 2 ms 16728 KB Output is correct
4 Correct 179 ms 17688 KB Output is correct
5 Correct 766 ms 21440 KB Output is correct
6 Correct 430 ms 16212 KB Output is correct
7 Correct 540 ms 21528 KB Output is correct
8 Correct 326 ms 16984 KB Output is correct
9 Correct 674 ms 21300 KB Output is correct
10 Correct 780 ms 21180 KB Output is correct
11 Correct 730 ms 21696 KB Output is correct
12 Correct 734 ms 21188 KB Output is correct
13 Correct 733 ms 21584 KB Output is correct
14 Correct 730 ms 21692 KB Output is correct
15 Correct 671 ms 21592 KB Output is correct
16 Correct 672 ms 21608 KB Output is correct
17 Correct 2 ms 14680 KB Output is correct
18 Correct 2 ms 14932 KB Output is correct
19 Correct 3 ms 14680 KB Output is correct
20 Correct 2 ms 14680 KB Output is correct
21 Correct 4 ms 14680 KB Output is correct
22 Correct 3 ms 14680 KB Output is correct
23 Correct 3 ms 14680 KB Output is correct
24 Correct 3 ms 16728 KB Output is correct
25 Correct 2 ms 14680 KB Output is correct
26 Correct 2 ms 14680 KB Output is correct
27 Correct 12 ms 14680 KB Output is correct
28 Correct 15 ms 14680 KB Output is correct
29 Correct 12 ms 14680 KB Output is correct
30 Correct 11 ms 14680 KB Output is correct
31 Correct 13 ms 16728 KB Output is correct
32 Correct 2 ms 14680 KB Output is correct
33 Correct 17 ms 18512 KB Output is correct
34 Correct 43 ms 21232 KB Output is correct
35 Correct 23 ms 21184 KB Output is correct
36 Correct 32 ms 21424 KB Output is correct
37 Correct 25 ms 22012 KB Output is correct
38 Correct 26 ms 21672 KB Output is correct
39 Correct 2 ms 14680 KB Output is correct
40 Correct 2 ms 16728 KB Output is correct
41 Correct 2 ms 16728 KB Output is correct
42 Correct 178 ms 17732 KB Output is correct
43 Correct 698 ms 21444 KB Output is correct
44 Correct 494 ms 15960 KB Output is correct
45 Correct 738 ms 21240 KB Output is correct
46 Correct 411 ms 16984 KB Output is correct
47 Correct 708 ms 21384 KB Output is correct
48 Correct 762 ms 21184 KB Output is correct
49 Correct 739 ms 21696 KB Output is correct
50 Correct 753 ms 21180 KB Output is correct
51 Correct 612 ms 21380 KB Output is correct
52 Correct 645 ms 21860 KB Output is correct
53 Correct 650 ms 21436 KB Output is correct
54 Correct 642 ms 21420 KB Output is correct
55 Correct 3 ms 16884 KB Output is correct
56 Incorrect 64 ms 21200 KB Output isn't correct
57 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 2 ms 14680 KB Output is correct
2 Correct 2 ms 14680 KB Output is correct
3 Incorrect 85 ms 20008 KB Output isn't correct
4 Halted 0 ms 0 KB -