Submission #520483

# Submission time Handle Problem Language Result Execution time Memory
520483 2022-01-30T07:06:27 Z KoD Diversity (CEOI21_diversity) C++17
64 / 100
619 ms 6272 KB
#include <bits/stdc++.h>

using std::vector;
using std::array;
using std::pair;
using std::tuple;

using i64 = std::int64_t;

constexpr int B = 600;
constexpr int MAX = 300000;

constexpr i64 sum_k(const int n) {
    return (i64)n * (n + 1) / 2;
}
 
constexpr i64 sum_k2(const int n) {
    return (i64)n * (n + 1) * (2 * n + 1) / 6;
}

int main() {
    std::ios_base::sync_with_stdio(false);
    std::cin.tie(nullptr);
    int N, Q;
    std::cin >> N >> Q;
    vector<int> A(N);
    for (auto& x : A) {
        std::cin >> x;
        x -= 1;
    }
    vector<vector<int>> query((N + B - 1) / B);
    vector<int> L(Q), R(Q);
    for (int i = 0; i < Q; ++i) {
        std::cin >> L[i] >> R[i];
        L[i] -= 1;
        query[L[i] / B].push_back(i);
    }
    vector<i64> ans(Q);
    vector<int> count(MAX), block(N + 1), nonzero;
    const auto change = [&](const int x, const int c) {
        block[count[x]] -= 1;
        count[x] += c;
        if (count[x] != 0 and block[count[x]] == 0) {
            nonzero.push_back(count[x]);
        }
        block[count[x]] += 1;
    };
    for (auto& qidx : query) {
        if (qidx.empty()) {
            continue;
        }
        std::sort(qidx.begin(), qidx.end(), [&](const int i, const int j) {
            return R[i] < R[j]; 
        });
        int l = 0, r = 0;
        for (const int q : qidx) {
            while (l < L[q]) {
                change(A[l++], -1);
            }
            while (l > L[q]) {
                change(A[--l], 1);
            }
            while (r < R[q]) {
                change(A[r++], 1);
            }
            std::sort(nonzero.begin(), nonzero.end());
            nonzero.erase(std::unique(nonzero.begin(), nonzero.end()), nonzero.end());
            const int len = R[q] - L[q];
            ans[q] = sum_k(len);
            const auto add_to = [&](const i64 s, const int a, const int n) {
                ans[q] += (n * s + sum_k(n) * a) * len;
                ans[q] -= n * s * s + 2 * s * sum_k(n) * a + sum_k2(n) * a * a;
            };
            int left = 0, right = 0;
            vector<int> next;
            for (const int a : nonzero) {
                const int b = block[a];
                if (b == 0) {
                    continue;
                }
                next.push_back(a);
                if (left > right) {
                    std::swap(left, right);
                }
                const int need = std::min(b, (right - left + a - 1) / a);
                const int to_left = need + (b - need) / 2;
                const int to_right = b - to_left;
                add_to(left, a, to_left);
                add_to(right, a, to_right);
                left += a * to_left;
                right += a * to_right;
            }
            ans[q] -= (i64)left * (len - left);
            nonzero = std::move(next);
        }
    }
    for (const i64 x : ans) {
        std::cout << x << '\n';
    }
    return 0;
}
# Verdict Execution time Memory Grader output
1 Correct 1 ms 1484 KB Output is correct
2 Correct 1 ms 1484 KB Output is correct
3 Correct 1 ms 1484 KB Output is correct
4 Correct 1 ms 1484 KB Output is correct
5 Correct 2 ms 1484 KB Output is correct
6 Correct 1 ms 1472 KB Output is correct
7 Correct 1 ms 1468 KB Output is correct
8 Correct 1 ms 1484 KB Output is correct
9 Correct 1 ms 1484 KB Output is correct
10 Correct 1 ms 1484 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 1 ms 1484 KB Output is correct
2 Correct 1 ms 1484 KB Output is correct
3 Correct 2 ms 1612 KB Output is correct
4 Correct 13 ms 3028 KB Output is correct
5 Correct 26 ms 4352 KB Output is correct
6 Correct 41 ms 6168 KB Output is correct
7 Correct 37 ms 6120 KB Output is correct
8 Correct 42 ms 6080 KB Output is correct
9 Correct 38 ms 6080 KB Output is correct
10 Correct 39 ms 6080 KB Output is correct
11 Correct 39 ms 6148 KB Output is correct
12 Correct 40 ms 6104 KB Output is correct
13 Correct 40 ms 6080 KB Output is correct
14 Correct 41 ms 6088 KB Output is correct
15 Correct 39 ms 6068 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 1 ms 1484 KB Output is correct
2 Correct 1 ms 1484 KB Output is correct
3 Correct 2 ms 1612 KB Output is correct
4 Correct 13 ms 3028 KB Output is correct
5 Correct 26 ms 4352 KB Output is correct
6 Correct 41 ms 6168 KB Output is correct
7 Correct 37 ms 6120 KB Output is correct
8 Correct 42 ms 6080 KB Output is correct
9 Correct 38 ms 6080 KB Output is correct
10 Correct 39 ms 6080 KB Output is correct
11 Correct 39 ms 6148 KB Output is correct
12 Correct 40 ms 6104 KB Output is correct
13 Correct 40 ms 6080 KB Output is correct
14 Correct 41 ms 6088 KB Output is correct
15 Correct 39 ms 6068 KB Output is correct
16 Correct 1 ms 1484 KB Output is correct
17 Correct 1 ms 1484 KB Output is correct
18 Correct 2 ms 1612 KB Output is correct
19 Correct 16 ms 3032 KB Output is correct
20 Correct 26 ms 4300 KB Output is correct
21 Correct 41 ms 6228 KB Output is correct
22 Correct 40 ms 6204 KB Output is correct
23 Correct 41 ms 6200 KB Output is correct
24 Correct 41 ms 6268 KB Output is correct
25 Correct 39 ms 6208 KB Output is correct
26 Correct 40 ms 6184 KB Output is correct
27 Correct 40 ms 6272 KB Output is correct
28 Correct 44 ms 6244 KB Output is correct
29 Correct 41 ms 6184 KB Output is correct
30 Correct 40 ms 6240 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 1 ms 1484 KB Output is correct
2 Correct 1 ms 1484 KB Output is correct
3 Correct 2 ms 1612 KB Output is correct
4 Correct 13 ms 3028 KB Output is correct
5 Correct 26 ms 4352 KB Output is correct
6 Correct 41 ms 6168 KB Output is correct
7 Correct 37 ms 6120 KB Output is correct
8 Correct 42 ms 6080 KB Output is correct
9 Correct 38 ms 6080 KB Output is correct
10 Correct 39 ms 6080 KB Output is correct
11 Correct 39 ms 6148 KB Output is correct
12 Correct 40 ms 6104 KB Output is correct
13 Correct 40 ms 6080 KB Output is correct
14 Correct 41 ms 6088 KB Output is correct
15 Correct 39 ms 6068 KB Output is correct
16 Correct 1 ms 1484 KB Output is correct
17 Correct 1 ms 1484 KB Output is correct
18 Correct 2 ms 1612 KB Output is correct
19 Correct 16 ms 3032 KB Output is correct
20 Correct 26 ms 4300 KB Output is correct
21 Correct 41 ms 6228 KB Output is correct
22 Correct 40 ms 6204 KB Output is correct
23 Correct 41 ms 6200 KB Output is correct
24 Correct 41 ms 6268 KB Output is correct
25 Correct 39 ms 6208 KB Output is correct
26 Correct 40 ms 6184 KB Output is correct
27 Correct 40 ms 6272 KB Output is correct
28 Correct 44 ms 6244 KB Output is correct
29 Correct 41 ms 6184 KB Output is correct
30 Correct 40 ms 6240 KB Output is correct
31 Correct 1 ms 1484 KB Output is correct
32 Correct 1 ms 1484 KB Output is correct
33 Correct 1 ms 1464 KB Output is correct
34 Correct 2 ms 1612 KB Output is correct
35 Correct 2 ms 1612 KB Output is correct
36 Correct 2 ms 1532 KB Output is correct
37 Correct 5 ms 2140 KB Output is correct
38 Correct 5 ms 1996 KB Output is correct
39 Correct 10 ms 2508 KB Output is correct
40 Correct 23 ms 3656 KB Output is correct
41 Correct 28 ms 4728 KB Output is correct
42 Correct 28 ms 4668 KB Output is correct
43 Correct 29 ms 4620 KB Output is correct
44 Correct 28 ms 4668 KB Output is correct
45 Correct 31 ms 4680 KB Output is correct
46 Correct 29 ms 4676 KB Output is correct
47 Correct 28 ms 4716 KB Output is correct
48 Correct 28 ms 4664 KB Output is correct
49 Correct 28 ms 4704 KB Output is correct
50 Correct 31 ms 4684 KB Output is correct
51 Correct 28 ms 4620 KB Output is correct
52 Correct 28 ms 4708 KB Output is correct
53 Correct 28 ms 4708 KB Output is correct
54 Correct 30 ms 4816 KB Output is correct
55 Correct 28 ms 4680 KB Output is correct
56 Correct 30 ms 4708 KB Output is correct
57 Correct 30 ms 4808 KB Output is correct
58 Correct 28 ms 4692 KB Output is correct
59 Correct 31 ms 4680 KB Output is correct
60 Correct 31 ms 4696 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 1 ms 1484 KB Output is correct
2 Correct 1 ms 1484 KB Output is correct
3 Correct 1 ms 1484 KB Output is correct
4 Correct 1 ms 1484 KB Output is correct
5 Correct 2 ms 1484 KB Output is correct
6 Correct 1 ms 1472 KB Output is correct
7 Correct 1 ms 1468 KB Output is correct
8 Correct 1 ms 1484 KB Output is correct
9 Correct 1 ms 1484 KB Output is correct
10 Correct 1 ms 1484 KB Output is correct
11 Correct 1 ms 1484 KB Output is correct
12 Correct 1 ms 1484 KB Output is correct
13 Correct 2 ms 1612 KB Output is correct
14 Correct 13 ms 3028 KB Output is correct
15 Correct 26 ms 4352 KB Output is correct
16 Correct 41 ms 6168 KB Output is correct
17 Correct 37 ms 6120 KB Output is correct
18 Correct 42 ms 6080 KB Output is correct
19 Correct 38 ms 6080 KB Output is correct
20 Correct 39 ms 6080 KB Output is correct
21 Correct 39 ms 6148 KB Output is correct
22 Correct 40 ms 6104 KB Output is correct
23 Correct 40 ms 6080 KB Output is correct
24 Correct 41 ms 6088 KB Output is correct
25 Correct 39 ms 6068 KB Output is correct
26 Correct 1 ms 1484 KB Output is correct
27 Correct 1 ms 1484 KB Output is correct
28 Correct 2 ms 1612 KB Output is correct
29 Correct 16 ms 3032 KB Output is correct
30 Correct 26 ms 4300 KB Output is correct
31 Correct 41 ms 6228 KB Output is correct
32 Correct 40 ms 6204 KB Output is correct
33 Correct 41 ms 6200 KB Output is correct
34 Correct 41 ms 6268 KB Output is correct
35 Correct 39 ms 6208 KB Output is correct
36 Correct 40 ms 6184 KB Output is correct
37 Correct 40 ms 6272 KB Output is correct
38 Correct 44 ms 6244 KB Output is correct
39 Correct 41 ms 6184 KB Output is correct
40 Correct 40 ms 6240 KB Output is correct
41 Correct 1 ms 1484 KB Output is correct
42 Correct 1 ms 1484 KB Output is correct
43 Correct 1 ms 1464 KB Output is correct
44 Correct 2 ms 1612 KB Output is correct
45 Correct 2 ms 1612 KB Output is correct
46 Correct 2 ms 1532 KB Output is correct
47 Correct 5 ms 2140 KB Output is correct
48 Correct 5 ms 1996 KB Output is correct
49 Correct 10 ms 2508 KB Output is correct
50 Correct 23 ms 3656 KB Output is correct
51 Correct 28 ms 4728 KB Output is correct
52 Correct 28 ms 4668 KB Output is correct
53 Correct 29 ms 4620 KB Output is correct
54 Correct 28 ms 4668 KB Output is correct
55 Correct 31 ms 4680 KB Output is correct
56 Correct 29 ms 4676 KB Output is correct
57 Correct 28 ms 4716 KB Output is correct
58 Correct 28 ms 4664 KB Output is correct
59 Correct 28 ms 4704 KB Output is correct
60 Correct 31 ms 4684 KB Output is correct
61 Correct 28 ms 4620 KB Output is correct
62 Correct 28 ms 4708 KB Output is correct
63 Correct 28 ms 4708 KB Output is correct
64 Correct 30 ms 4816 KB Output is correct
65 Correct 28 ms 4680 KB Output is correct
66 Correct 30 ms 4708 KB Output is correct
67 Correct 30 ms 4808 KB Output is correct
68 Correct 28 ms 4692 KB Output is correct
69 Correct 31 ms 4680 KB Output is correct
70 Correct 31 ms 4696 KB Output is correct
71 Correct 5 ms 1936 KB Output is correct
72 Correct 5 ms 1996 KB Output is correct
73 Correct 5 ms 1996 KB Output is correct
74 Correct 5 ms 1984 KB Output is correct
75 Correct 5 ms 1996 KB Output is correct
76 Correct 9 ms 2380 KB Output is correct
77 Correct 9 ms 2368 KB Output is correct
78 Correct 8 ms 2368 KB Output is correct
79 Correct 9 ms 2380 KB Output is correct
80 Correct 9 ms 2396 KB Output is correct
81 Correct 17 ms 3304 KB Output is correct
82 Correct 16 ms 3312 KB Output is correct
83 Correct 16 ms 3268 KB Output is correct
84 Correct 17 ms 3292 KB Output is correct
85 Correct 17 ms 3268 KB Output is correct
86 Correct 19 ms 3188 KB Output is correct
87 Correct 19 ms 3268 KB Output is correct
88 Correct 18 ms 3288 KB Output is correct
89 Correct 18 ms 3252 KB Output is correct
90 Correct 19 ms 3268 KB Output is correct
91 Correct 28 ms 4052 KB Output is correct
92 Correct 27 ms 4076 KB Output is correct
93 Correct 28 ms 4084 KB Output is correct
94 Correct 28 ms 3988 KB Output is correct
95 Correct 27 ms 3972 KB Output is correct
96 Correct 32 ms 4008 KB Output is correct
97 Correct 30 ms 3984 KB Output is correct
98 Correct 31 ms 4000 KB Output is correct
99 Correct 31 ms 4088 KB Output is correct
100 Correct 30 ms 3820 KB Output is correct
101 Correct 30 ms 3828 KB Output is correct
102 Correct 36 ms 3780 KB Output is correct
103 Correct 32 ms 3800 KB Output is correct
104 Correct 32 ms 3796 KB Output is correct
105 Correct 30 ms 3812 KB Output is correct
106 Correct 30 ms 3832 KB Output is correct
107 Correct 30 ms 3808 KB Output is correct
108 Correct 31 ms 3788 KB Output is correct
109 Correct 31 ms 3804 KB Output is correct
110 Correct 31 ms 3788 KB Output is correct
111 Correct 32 ms 3796 KB Output is correct
112 Correct 30 ms 3844 KB Output is correct
113 Correct 32 ms 3812 KB Output is correct
114 Correct 32 ms 3912 KB Output is correct
115 Correct 30 ms 3812 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 1 ms 1484 KB Output is correct
2 Correct 1 ms 1484 KB Output is correct
3 Correct 1 ms 1484 KB Output is correct
4 Correct 1 ms 1484 KB Output is correct
5 Correct 2 ms 1484 KB Output is correct
6 Correct 1 ms 1472 KB Output is correct
7 Correct 1 ms 1468 KB Output is correct
8 Correct 1 ms 1484 KB Output is correct
9 Correct 1 ms 1484 KB Output is correct
10 Correct 1 ms 1484 KB Output is correct
11 Correct 1 ms 1484 KB Output is correct
12 Correct 1 ms 1484 KB Output is correct
13 Correct 2 ms 1612 KB Output is correct
14 Correct 13 ms 3028 KB Output is correct
15 Correct 26 ms 4352 KB Output is correct
16 Correct 41 ms 6168 KB Output is correct
17 Correct 37 ms 6120 KB Output is correct
18 Correct 42 ms 6080 KB Output is correct
19 Correct 38 ms 6080 KB Output is correct
20 Correct 39 ms 6080 KB Output is correct
21 Correct 39 ms 6148 KB Output is correct
22 Correct 40 ms 6104 KB Output is correct
23 Correct 40 ms 6080 KB Output is correct
24 Correct 41 ms 6088 KB Output is correct
25 Correct 39 ms 6068 KB Output is correct
26 Correct 1 ms 1484 KB Output is correct
27 Correct 1 ms 1484 KB Output is correct
28 Correct 2 ms 1612 KB Output is correct
29 Correct 16 ms 3032 KB Output is correct
30 Correct 26 ms 4300 KB Output is correct
31 Correct 41 ms 6228 KB Output is correct
32 Correct 40 ms 6204 KB Output is correct
33 Correct 41 ms 6200 KB Output is correct
34 Correct 41 ms 6268 KB Output is correct
35 Correct 39 ms 6208 KB Output is correct
36 Correct 40 ms 6184 KB Output is correct
37 Correct 40 ms 6272 KB Output is correct
38 Correct 44 ms 6244 KB Output is correct
39 Correct 41 ms 6184 KB Output is correct
40 Correct 40 ms 6240 KB Output is correct
41 Correct 1 ms 1484 KB Output is correct
42 Correct 1 ms 1484 KB Output is correct
43 Correct 1 ms 1464 KB Output is correct
44 Correct 2 ms 1612 KB Output is correct
45 Correct 2 ms 1612 KB Output is correct
46 Correct 2 ms 1532 KB Output is correct
47 Correct 5 ms 2140 KB Output is correct
48 Correct 5 ms 1996 KB Output is correct
49 Correct 10 ms 2508 KB Output is correct
50 Correct 23 ms 3656 KB Output is correct
51 Correct 28 ms 4728 KB Output is correct
52 Correct 28 ms 4668 KB Output is correct
53 Correct 29 ms 4620 KB Output is correct
54 Correct 28 ms 4668 KB Output is correct
55 Correct 31 ms 4680 KB Output is correct
56 Correct 29 ms 4676 KB Output is correct
57 Correct 28 ms 4716 KB Output is correct
58 Correct 28 ms 4664 KB Output is correct
59 Correct 28 ms 4704 KB Output is correct
60 Correct 31 ms 4684 KB Output is correct
61 Correct 28 ms 4620 KB Output is correct
62 Correct 28 ms 4708 KB Output is correct
63 Correct 28 ms 4708 KB Output is correct
64 Correct 30 ms 4816 KB Output is correct
65 Correct 28 ms 4680 KB Output is correct
66 Correct 30 ms 4708 KB Output is correct
67 Correct 30 ms 4808 KB Output is correct
68 Correct 28 ms 4692 KB Output is correct
69 Correct 31 ms 4680 KB Output is correct
70 Correct 31 ms 4696 KB Output is correct
71 Correct 5 ms 1936 KB Output is correct
72 Correct 5 ms 1996 KB Output is correct
73 Correct 5 ms 1996 KB Output is correct
74 Correct 5 ms 1984 KB Output is correct
75 Correct 5 ms 1996 KB Output is correct
76 Correct 9 ms 2380 KB Output is correct
77 Correct 9 ms 2368 KB Output is correct
78 Correct 8 ms 2368 KB Output is correct
79 Correct 9 ms 2380 KB Output is correct
80 Correct 9 ms 2396 KB Output is correct
81 Correct 17 ms 3304 KB Output is correct
82 Correct 16 ms 3312 KB Output is correct
83 Correct 16 ms 3268 KB Output is correct
84 Correct 17 ms 3292 KB Output is correct
85 Correct 17 ms 3268 KB Output is correct
86 Correct 19 ms 3188 KB Output is correct
87 Correct 19 ms 3268 KB Output is correct
88 Correct 18 ms 3288 KB Output is correct
89 Correct 18 ms 3252 KB Output is correct
90 Correct 19 ms 3268 KB Output is correct
91 Correct 28 ms 4052 KB Output is correct
92 Correct 27 ms 4076 KB Output is correct
93 Correct 28 ms 4084 KB Output is correct
94 Correct 28 ms 3988 KB Output is correct
95 Correct 27 ms 3972 KB Output is correct
96 Correct 32 ms 4008 KB Output is correct
97 Correct 30 ms 3984 KB Output is correct
98 Correct 31 ms 4000 KB Output is correct
99 Correct 31 ms 4088 KB Output is correct
100 Correct 30 ms 3820 KB Output is correct
101 Correct 30 ms 3828 KB Output is correct
102 Correct 36 ms 3780 KB Output is correct
103 Correct 32 ms 3800 KB Output is correct
104 Correct 32 ms 3796 KB Output is correct
105 Correct 30 ms 3812 KB Output is correct
106 Correct 30 ms 3832 KB Output is correct
107 Correct 30 ms 3808 KB Output is correct
108 Correct 31 ms 3788 KB Output is correct
109 Correct 31 ms 3804 KB Output is correct
110 Correct 31 ms 3788 KB Output is correct
111 Correct 32 ms 3796 KB Output is correct
112 Correct 30 ms 3844 KB Output is correct
113 Correct 32 ms 3812 KB Output is correct
114 Correct 32 ms 3912 KB Output is correct
115 Correct 30 ms 3812 KB Output is correct
116 Incorrect 619 ms 3296 KB Output isn't correct
117 Halted 0 ms 0 KB -