Submission #25904

# Submission time Handle Problem Language Result Execution time Memory
25904 2017-06-25T03:36:11 Z model_code Watering can (POI13_kon) C++11
100 / 100
809 ms 18556 KB
/*************************************************************************
 *                                                                       *
 *                    XX Olimpiada Informatyczna                         *
 *                                                                       *
 *   Zadanie:              Konewka                                       *
 *   Autor:                Michal Adamczyk                               *
 *   Zlozonosc czasowa:    O((n + q) log n)                              *
 *   Zlozonosc pamieciowa: O(n)                                          *
 *   Opis:                 Rozwiazanie wzorcowe                          *
 *                         A - drzewo przedzialowe                       *
 *                         B - drzewo przedzialowe                       *
 *                                                                       *
 *************************************************************************/

#include <algorithm>
using namespace std;

const int INF = 1000*1000*1000;
const int M = 524288; // 2^19

int N, K;
int dd[2*M]; // drzewa dojrzale
int nd_w[2*M]; // drzewa niedojrzale - wartosc w wezle
int nd_max[2*M]; // drzewa niedojrzale - max z synow + wartosc dla wezla
int nd_maxind[2*M]; // drzewa niedojrzale - numer wezla, w którym jest max

/** operacje na drzewach dojrzalych **/

/* wstawia dojrzale drzewo o numerze a */
void wstawDojrzale(int a) {
    int v = M + a;
    dd[v]++;
    while (v != 1) {
        v /= 2;
        dd[v] = dd[2*v] + dd[2*v+1];
    }
}

/* zwraca liczbe dojrzalych drzew na przedziale [a, b] */
int ileDojrzalych(int a, int b) {
    int va = M + a, vb = M + b;
    int wyn = dd[va];
    if (va != vb) wyn += dd[vb];
    while(va / 2 != vb / 2) {
        if (va%2 == 0) wyn += dd[va+1];
        if (vb%2 == 1) wyn += dd[vb-1];
        va /= 2, vb /= 2;
    }
    return wyn;
}

/** operacje na drzewach niedojrzalych */

void przesun_nizej(int v) {
    int ile = nd_w[v];
    nd_w[v] = 0;
    nd_w[2*v] += ile;
    nd_max[2*v] += ile;
    nd_w[2*v+1] += ile;
    nd_max[2*v+1] += ile;
}

/* ustawia wysokosc drzewa a na x */
void ustawWysokosc(int a, int x) {
    int v = M + a;
    nd_w[v] = x;
    nd_max[v] = x;
    nd_maxind[v] = v;
    while (v != 1) {
        v /= 2;
        if (nd_max[2*v] > nd_max[2*v+1])
            nd_max[v] = nd_w[v] + nd_max[2*v], nd_maxind[v] = nd_maxind[2*v];
        else
            nd_max[v] = nd_w[v] + nd_max[2*v+1], nd_maxind[v] = nd_maxind[2*v+1];
    }
}

/* zwieksza o jeden wysokosci drzew na przedziale [p, k] */
void zwieksz(int v, int p, int k, int p_b, int k_b) {
    // v - indeks tablicy odpowiadający przedzialowi bazowemu [p_b, k_b]
    if (k < p_b || k_b < p)
        return;

    if (p <= p_b && k_b <= k) {
        nd_w[v]++;
        nd_max[v]++;
        return;
    }

    if(nd_w[v] > 0)
        przesun_nizej(v);

    int m = (p_b + k_b) / 2;
    zwieksz(2*v, p, k, p_b, m);
    zwieksz(2*v+1, p, k, m+1, k_b);

    if (nd_max[2*v] > nd_max[2*v+1])
        nd_max[v] = nd_max[2*v], nd_maxind[v] = nd_maxind[2*v];
    else
        nd_max[v] = nd_max[2*v+1], nd_maxind[v] = nd_maxind[2*v+1];

}

/* zwraca parę <wysokosc najwyzszego drzewa na przedziale [p, k], jego indeks w tablicy> */
pair<int, int> najwyzsze(int v, int p, int k, int p_b, int k_b) {
    // v - indeks tablicy odpowiadający przedzialowi bazowemu [p_b, k_b]
    if (k < p_b || k_b < p)
        return make_pair(-1,-1);

    if (p <= p_b && k_b <= k)
        return make_pair(nd_max[v], nd_maxind[v]);

    if(nd_w[v] > 0)
        przesun_nizej(v);

    int m = (p_b + k_b) / 2;
    pair<int,int> a = najwyzsze(2*v, p, k, p_b, m);
    pair<int,int> b = najwyzsze(2*v+1, p, k, m+1, k_b);

    return max(a,b);
}

/** funkcje **/

void inicjuj(int n, int k, int *D) {
    N = n, K = k;
    for (int i = 0; i < N; ++i)
        ustawWysokosc(i, D[i]);
}

void podlej(int a, int b) {
    zwieksz(1,a,b,0,M-1);
}

int dojrzale(int a, int b) {

    pair<int, int> n = najwyzsze(1,a,b,0,M-1);
    while(n.first >= K) {
        n.second -= M;
        ustawWysokosc(n.second, -INF);
        wstawDojrzale(n.second);

        n = najwyzsze(1,a,b,0,M-1);
    }

    return ileDojrzalych(a,b);
}
# Verdict Execution time Memory Grader output
1 Correct 0 ms 18556 KB Output is correct
2 Correct 0 ms 18556 KB Output is correct
3 Correct 0 ms 18552 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 0 ms 18548 KB Output is correct
2 Correct 3 ms 18552 KB Output is correct
3 Correct 0 ms 18556 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 83 ms 18556 KB Output is correct
2 Correct 66 ms 18552 KB Output is correct
3 Correct 73 ms 18556 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 129 ms 18556 KB Output is correct
2 Correct 119 ms 18552 KB Output is correct
3 Correct 149 ms 18552 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 193 ms 18556 KB Output is correct
2 Correct 203 ms 18552 KB Output is correct
3 Correct 209 ms 18552 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 366 ms 18556 KB Output is correct
2 Correct 276 ms 18552 KB Output is correct
3 Correct 329 ms 18552 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 466 ms 18556 KB Output is correct
2 Correct 273 ms 18556 KB Output is correct
3 Correct 409 ms 18552 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 666 ms 18552 KB Output is correct
2 Correct 533 ms 18556 KB Output is correct
3 Correct 533 ms 18552 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 706 ms 18552 KB Output is correct
2 Correct 586 ms 18552 KB Output is correct
3 Correct 669 ms 18552 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 809 ms 18556 KB Output is correct
2 Correct 779 ms 18552 KB Output is correct
3 Correct 699 ms 18556 KB Output is correct