oj.uz

Submission #850818

# Submission time Handle Problem Language Result Execution time Memory
850818 2023-09-17T13:28:33 Z EJIC_B_KEDAX Sprinkler (JOI22_sprinkler) C++14
41 / 100
 4000 ms 142008 KB 
#include <bits/stdc++.h>
#include <random>

#ifndef LOCAL
    #pragma GCC optimize("O3")
    //#pragma GCC optimize("Ofast")
    #pragma GCC optimize("unroll-loops")
    #pragma GCC target("avx,avx2,bmi,bmi2,popcnt,lzcnt")
#endif
using namespace std;
typedef long long ll;
typedef double dd;
typedef long double ld;
typedef unsigned int uii;
#define x first
#define y second
#define all(x) x.begin(), x.end()
#define rall(x) x.rbegin(), x.rend()



void solve();

mt19937_64 mt(1);

int32_t main() {
#ifdef LOCAL
    freopen("input.in", "r", stdin);
    freopen("output.out", "w", stdout);
#else
    ios::sync_with_stdio(false);
    cin.tie(nullptr); cout.tie(nullptr);
    //freopen("amusing.in", "r", stdin);
    //freopen("amusing.out", "w", stdout);
#endif
    cout << fixed << setprecision(30);
    int tests = 1;
    //cin >> tests;
    while (tests--) {
        solve();
    }
}

ll mod;

class SegmentTree {
public:
    void build(const vector<int>& b) {
        int sz = 1;
        while (sz < b.size()) {
            sz *= 2;
        }
        x.resize(2 * sz - 1);
        vn.resize(2 * sz - 1, 1);
        size = 2 * sz - 1;
        for (int i = sz - 1; i < 2 * sz - 1; i++) {
            if (i - sz + 1 < b.size()) {
                x[i] = b[i - sz + 1];
            } else {
                x[i] = 0;
            }
        }
        for (int i = sz - 2; i >= 0; i--) {
            //update(i);
            x[i] = 0;
        }
    }

    int get(int l, int r, int i = 0, int lx = 0, int rx = -1) {
        if (rx == -1) {
            rx = size / 2;
        }
        if (lx > r || rx < l) {
            return 0;
        }

        if (lx >= l && rx <= r) {
            return x[i];
        }

        updateVN(i);
        return get(l, r, 2 * i + 1, lx, (rx + lx) / 2) + get(l, r, 2 * i + 2, (rx + lx) / 2 + 1, rx);
    }

    void change(int ind, int c, int i = 0, int lx = 0, int rx = -1) {
        if (rx == -1) {
            rx = size / 2;
        }
        if (lx == rx) {
            x[i] = c;
            return;
        }
        updateVN(i);
        if ((rx + lx) / 2 >= ind) {
            change(ind, c, 2 * i + 1, lx, (rx + lx) / 2);
        } else {
            change(ind, c, 2 * i + 2, (rx + lx) / 2 + 1, rx);
        }
        //update(i);
    }
    void plusOtr(int l, int r, ll p, int i = 0, int lx = 0, int rx = -1) {
        if (r < l) {
            return;
        }
        if (rx == -1) {
            rx = size / 2;
        }
        if (lx >= l && rx <= r) {
            x[i] = (x[i] * p) % mod;
            vn[i] = (vn[i] * p) % mod;
            return;
        }
        if (lx > r || rx < l) {
            return;
        }
        updateVN(i);
        plusOtr(l, r, p, 2 * i + 1, lx, (rx + lx) / 2);
        plusOtr(l, r, p, 2 * i + 2, (rx + lx) / 2 + 1, rx);
        //update(i);
    }
private:
    vector<ll> x, vn;
    int size;
    void updateVN(int i) {
        if (2 * i + 2 < size) {
            if (vn[i] != 1) {
                if (x[2 * i + 1]) {
                    x[2 * i + 1] = (x[2 * i + 1] * vn[i]) % mod;
                }
                vn[2 * i + 1] = (vn[2 * i + 1] * vn[i]) % mod;
                if (x[2 * i + 2]) {
                    x[2 * i + 2] = (x[2 * i + 2] * vn[i]) % mod;
                }
                vn[2 * i + 2] = (vn[2 * i + 2] * vn[i]) % mod;
                vn[i] = 1;
            }
        }
    }
    void update(int i) {
        if (2 * i + 2 < size) {
            x[i] = x[2 * i + 1] + x[2 * i + 2];
        }
    }
};

pair<int, int> seg[200200][41];

void dfs(int s, int p, vector<vector<int>>& g, int dist[], vector<vector<int>>& d, int pr[]) {
    if (p != -1) {
        dist[s] = dist[p] + 1;
    } else {
        dist[s] = 0;
    }
    pr[s] = p;
    d[dist[s]].push_back(s);
    seg[s][0] = {d[dist[s]].size() - 1, d[dist[s]].size() - 1};
    for (int i = 1; i < 41; i++) {
        seg[s][i] = {INT32_MAX, INT32_MIN};
    }
    bool f = true;
    for (int i : g[s]) {
        if (i != p) {
            dfs(i, s, g, dist, d, pr);
            for (int j = 1; j < 41; j++) {
                seg[s][j].x = min(seg[s][j].x, seg[i][j - 1].x);
            }
            for (int j = 1; j < 41; j++) {
                seg[s][j].y = max(seg[s][j].y, seg[i][j - 1].y);
            }
        }
    }
}

void solve() {
    int n;
    cin >> n >> mod;
    vector<vector<int>> g(n), d(n);
    int deg[200200], dist[200200], h[200200], p[200200];
    for (int i = 0; i < n; i++) {
        deg[i] = 0;
    }
    for (int i = 0; i < n - 1; i++) {
        int u, v;
        cin >> u >> v; u--; v--;
        g[u].push_back(v);
        g[v].push_back(u);
        deg[u]++;
        deg[v]++;
    }
    int root = 0;
    for (int i = 0; i < n; i++) {
        if (deg[i] == 1) {
            root = i;
            break;
        }
    }
    dfs(root, -1, g, dist, d, p);
    for (int i = 0; i < n; i++) {
        cin >> h[i];
    }
    for (int i = 0; i < n; i++) {
        for (int j = 0; j < d[i].size(); j++) {
            d[i][j] = h[d[i][j]];
        }
    }
    SegmentTree dec[200200];
    for (int i = 0; i < n; i++) {
        dec[i].build(d[i]);
    }
    int q;
    cin >> q;
    while (q--) {
        int t;
        cin >> t;
        if (t == 1) {
            int x, dk, wk;
            cin >> x >> dk >> wk; x--;
            int nw = x, lev = 0, nlev = 0;
            for (int l = -dk; l <= dk; l++) {
                if (dist[x] - l < 0) {
                    break;
                }
                if (dist[x] - l >= n) {
                    continue;
                }
                int c = (dk + l) / 2;
                while (c > nlev) {
                    nlev++;
                    if (p[nw] != -1) {
                        lev++;
                        nw = p[nw];
                    }
                }
                if (c - l >= 0 && c - l <= dk) {
                    int l1, r1;
                    l1 = seg[nw][lev - l].x;
                    r1 = seg[nw][lev - l].y;
                    dec[dist[x] - l].plusOtr(l1, r1, wk);
                }
            }
        } else {
            int x;
            cin >> x; x--;
            cout << (dec[dist[x]].get(seg[x][0].x, seg[x][0].y) % mod + mod) % mod << '\n';
        }
    }
}

Compilation message

sprinkler.cpp: In member function 'void SegmentTree::build(const std::vector<int>&)':
sprinkler.cpp:50:19: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<int>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
   50 |         while (sz < b.size()) {
      |                ~~~^~~~~~~~~~
sprinkler.cpp:57:28: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<int>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
   57 |             if (i - sz + 1 < b.size()) {
      |                 ~~~~~~~~~~~^~~~~~~~~~
sprinkler.cpp: In function 'void dfs(int, int, std::vector<std::vector<int> >&, int*, std::vector<std::vector<int> >&, int*)':
sprinkler.cpp:160:10: warning: unused variable 'f' [-Wunused-variable]
  160 |     bool f = true;
      |          ^
sprinkler.cpp: In function 'void solve()':
sprinkler.cpp:202:27: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<int>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
  202 |         for (int j = 0; j < d[i].size(); j++) {
      |                         ~~^~~~~~~~~~~~~

Subtask #1
3.0 / 3.0

# Verdict Execution time Memory Grader output
1 Correct 10 ms 14428 KB Output is correct
2 Correct 8 ms 14428 KB Output is correct
3 Correct 10 ms 14424 KB Output is correct
4 Correct 13 ms 16984 KB Output is correct
5 Correct 11 ms 16732 KB Output is correct
6 Correct 13 ms 16984 KB Output is correct
7 Correct 13 ms 16732 KB Output is correct
8 Correct 10 ms 16732 KB Output is correct
9 Correct 9 ms 14484 KB Output is correct
10 Correct 12 ms 14324 KB Output is correct
11 Correct 10 ms 14428 KB Output is correct
12 Correct 11 ms 14532 KB Output is correct
13 Correct 13 ms 14424 KB Output is correct
14 Correct 11 ms 14420 KB Output is correct
15 Correct 9 ms 14564 KB Output is correct
16 Correct 9 ms 14428 KB Output is correct
17 Correct 9 ms 14424 KB Output is correct
18 Correct 10 ms 14428 KB Output is correct
19 Correct 10 ms 14424 KB Output is correct
20 Correct 9 ms 14560 KB Output is correct
21 Correct 9 ms 14428 KB Output is correct
22 Correct 9 ms 14440 KB Output is correct
23 Correct 11 ms 14428 KB Output is correct
24 Correct 9 ms 14560 KB Output is correct
25 Correct 9 ms 14428 KB Output is correct
26 Correct 10 ms 14428 KB Output is correct
27 Correct 9 ms 14428 KB Output is correct
28 Correct 11 ms 14428 KB Output is correct

Subtask #2
9.0 / 9.0

# Verdict Execution time Memory Grader output
1 Correct 9 ms 14424 KB Output is correct
2 Correct 415 ms 120548 KB Output is correct
3 Correct 542 ms 117824 KB Output is correct
4 Correct 465 ms 140156 KB Output is correct
5 Correct 443 ms 119036 KB Output is correct
6 Correct 491 ms 119704 KB Output is correct
7 Correct 524 ms 118764 KB Output is correct
8 Correct 416 ms 118492 KB Output is correct
9 Correct 495 ms 141476 KB Output is correct
10 Correct 567 ms 138532 KB Output is correct
11 Correct 444 ms 120472 KB Output is correct
12 Correct 520 ms 117856 KB Output is correct
13 Correct 381 ms 118020 KB Output is correct
14 Correct 472 ms 119388 KB Output is correct
15 Correct 477 ms 118292 KB Output is correct
16 Correct 467 ms 118404 KB Output is correct
17 Correct 483 ms 118216 KB Output is correct
18 Correct 12 ms 14500 KB Output is correct
19 Correct 9 ms 14428 KB Output is correct
20 Correct 12 ms 14428 KB Output is correct
21 Correct 12 ms 14428 KB Output is correct
22 Correct 9 ms 14428 KB Output is correct

Subtask #3
29.0 / 29.0

# Verdict Execution time Memory Grader output
1 Correct 9 ms 14424 KB Output is correct
2 Correct 415 ms 120548 KB Output is correct
3 Correct 542 ms 117824 KB Output is correct
4 Correct 465 ms 140156 KB Output is correct
5 Correct 443 ms 119036 KB Output is correct
6 Correct 491 ms 119704 KB Output is correct
7 Correct 524 ms 118764 KB Output is correct
8 Correct 416 ms 118492 KB Output is correct
9 Correct 495 ms 141476 KB Output is correct
10 Correct 567 ms 138532 KB Output is correct
11 Correct 444 ms 120472 KB Output is correct
12 Correct 520 ms 117856 KB Output is correct
13 Correct 381 ms 118020 KB Output is correct
14 Correct 472 ms 119388 KB Output is correct
15 Correct 477 ms 118292 KB Output is correct
16 Correct 467 ms 118404 KB Output is correct
17 Correct 483 ms 118216 KB Output is correct
18 Correct 12 ms 14500 KB Output is correct
19 Correct 9 ms 14428 KB Output is correct
20 Correct 12 ms 14428 KB Output is correct
21 Correct 12 ms 14428 KB Output is correct
22 Correct 9 ms 14428 KB Output is correct
23 Correct 12 ms 14428 KB Output is correct
24 Correct 450 ms 120744 KB Output is correct
25 Correct 698 ms 118188 KB Output is correct
26 Correct 525 ms 140160 KB Output is correct
27 Correct 500 ms 119060 KB Output is correct
28 Correct 544 ms 119264 KB Output is correct
29 Correct 493 ms 119836 KB Output is correct
30 Correct 412 ms 118600 KB Output is correct
31 Correct 426 ms 142008 KB Output is correct
32 Correct 568 ms 138620 KB Output is correct
33 Correct 427 ms 120656 KB Output is correct
34 Correct 640 ms 117216 KB Output is correct
35 Correct 12 ms 14428 KB Output is correct
36 Correct 12 ms 14424 KB Output is correct
37 Correct 10 ms 14428 KB Output is correct
38 Correct 10 ms 14680 KB Output is correct
39 Correct 9 ms 14432 KB Output is correct
40 Correct 11 ms 14428 KB Output is correct
41 Correct 9 ms 14444 KB Output is correct
42 Correct 10 ms 14428 KB Output is correct
43 Correct 11 ms 14428 KB Output is correct
44 Correct 9 ms 14424 KB Output is correct
45 Correct 10 ms 14496 KB Output is correct
46 Correct 9 ms 14428 KB Output is correct
47 Correct 9 ms 14428 KB Output is correct
48 Correct 10 ms 14428 KB Output is correct
49 Correct 9 ms 14428 KB Output is correct

Subtask #4
0 / 12.0

# Verdict Execution time Memory Grader output
1 Correct 9 ms 14428 KB Output is correct
2 Correct 706 ms 139068 KB Output is correct
3 Correct 2691 ms 138324 KB Output is correct
4 Correct 1095 ms 138580 KB Output is correct
5 Correct 1962 ms 117304 KB Output is correct
6 Correct 1606 ms 118008 KB Output is correct
7 Correct 1111 ms 118220 KB Output is correct
8 Correct 428 ms 117052 KB Output is correct
9 Correct 733 ms 139264 KB Output is correct
10 Correct 2708 ms 138532 KB Output is correct
11 Correct 1018 ms 118508 KB Output is correct
12 Execution timed out 4054 ms 117160 KB Time limit exceeded
13 Halted 0 ms 0 KB -

Subtask #5
0 / 30.0

# Verdict Execution time Memory Grader output
1 Correct 9 ms 14428 KB Output is correct
2 Correct 699 ms 141672 KB Output is correct
3 Correct 2618 ms 138668 KB Output is correct
4 Correct 1097 ms 140120 KB Output is correct
5 Correct 1895 ms 119100 KB Output is correct
6 Correct 1704 ms 119096 KB Output is correct
7 Correct 1048 ms 118064 KB Output is correct
8 Correct 380 ms 118556 KB Output is correct
9 Correct 701 ms 141600 KB Output is correct
10 Correct 2745 ms 138872 KB Output is correct
11 Correct 1079 ms 120284 KB Output is correct
12 Execution timed out 4061 ms 117556 KB Time limit exceeded
13 Halted 0 ms 0 KB -

Subtask #6
0 / 17.0

# Verdict Execution time Memory Grader output
1 Correct 10 ms 14428 KB Output is correct
2 Correct 8 ms 14428 KB Output is correct
3 Correct 10 ms 14424 KB Output is correct
4 Correct 13 ms 16984 KB Output is correct
5 Correct 11 ms 16732 KB Output is correct
6 Correct 13 ms 16984 KB Output is correct
7 Correct 13 ms 16732 KB Output is correct
8 Correct 10 ms 16732 KB Output is correct
9 Correct 9 ms 14484 KB Output is correct
10 Correct 12 ms 14324 KB Output is correct
11 Correct 10 ms 14428 KB Output is correct
12 Correct 11 ms 14532 KB Output is correct
13 Correct 13 ms 14424 KB Output is correct
14 Correct 11 ms 14420 KB Output is correct
15 Correct 9 ms 14564 KB Output is correct
16 Correct 9 ms 14428 KB Output is correct
17 Correct 9 ms 14424 KB Output is correct
18 Correct 10 ms 14428 KB Output is correct
19 Correct 10 ms 14424 KB Output is correct
20 Correct 9 ms 14560 KB Output is correct
21 Correct 9 ms 14428 KB Output is correct
22 Correct 9 ms 14440 KB Output is correct
23 Correct 11 ms 14428 KB Output is correct
24 Correct 9 ms 14560 KB Output is correct
25 Correct 9 ms 14428 KB Output is correct
26 Correct 10 ms 14428 KB Output is correct
27 Correct 9 ms 14428 KB Output is correct
28 Correct 11 ms 14428 KB Output is correct
29 Correct 9 ms 14424 KB Output is correct
30 Correct 415 ms 120548 KB Output is correct
31 Correct 542 ms 117824 KB Output is correct
32 Correct 465 ms 140156 KB Output is correct
33 Correct 443 ms 119036 KB Output is correct
34 Correct 491 ms 119704 KB Output is correct
35 Correct 524 ms 118764 KB Output is correct
36 Correct 416 ms 118492 KB Output is correct
37 Correct 495 ms 141476 KB Output is correct
38 Correct 567 ms 138532 KB Output is correct
39 Correct 444 ms 120472 KB Output is correct
40 Correct 520 ms 117856 KB Output is correct
41 Correct 381 ms 118020 KB Output is correct
42 Correct 472 ms 119388 KB Output is correct
43 Correct 477 ms 118292 KB Output is correct
44 Correct 467 ms 118404 KB Output is correct
45 Correct 483 ms 118216 KB Output is correct
46 Correct 12 ms 14500 KB Output is correct
47 Correct 9 ms 14428 KB Output is correct
48 Correct 12 ms 14428 KB Output is correct
49 Correct 12 ms 14428 KB Output is correct
50 Correct 9 ms 14428 KB Output is correct
51 Correct 12 ms 14428 KB Output is correct
52 Correct 450 ms 120744 KB Output is correct
53 Correct 698 ms 118188 KB Output is correct
54 Correct 525 ms 140160 KB Output is correct
55 Correct 500 ms 119060 KB Output is correct
56 Correct 544 ms 119264 KB Output is correct
57 Correct 493 ms 119836 KB Output is correct
58 Correct 412 ms 118600 KB Output is correct
59 Correct 426 ms 142008 KB Output is correct
60 Correct 568 ms 138620 KB Output is correct
61 Correct 427 ms 120656 KB Output is correct
62 Correct 640 ms 117216 KB Output is correct
63 Correct 12 ms 14428 KB Output is correct
64 Correct 12 ms 14424 KB Output is correct
65 Correct 10 ms 14428 KB Output is correct
66 Correct 10 ms 14680 KB Output is correct
67 Correct 9 ms 14432 KB Output is correct
68 Correct 11 ms 14428 KB Output is correct
69 Correct 9 ms 14444 KB Output is correct
70 Correct 10 ms 14428 KB Output is correct
71 Correct 11 ms 14428 KB Output is correct
72 Correct 9 ms 14424 KB Output is correct
73 Correct 10 ms 14496 KB Output is correct
74 Correct 9 ms 14428 KB Output is correct
75 Correct 9 ms 14428 KB Output is correct
76 Correct 10 ms 14428 KB Output is correct
77 Correct 9 ms 14428 KB Output is correct
78 Correct 9 ms 14428 KB Output is correct
79 Correct 706 ms 139068 KB Output is correct
80 Correct 2691 ms 138324 KB Output is correct
81 Correct 1095 ms 138580 KB Output is correct
82 Correct 1962 ms 117304 KB Output is correct
83 Correct 1606 ms 118008 KB Output is correct
84 Correct 1111 ms 118220 KB Output is correct
85 Correct 428 ms 117052 KB Output is correct
86 Correct 733 ms 139264 KB Output is correct
87 Correct 2708 ms 138532 KB Output is correct
88 Correct 1018 ms 118508 KB Output is correct
89 Execution timed out 4054 ms 117160 KB Time limit exceeded
90 Halted 0 ms 0 KB -