#include <bits/stdc++.h>
using namespace std;
const int INF = 1e9;
struct segment_tree {
int size;
vector<array<int, 2>> tree;
vector<int> push;
array<int, 2> merge(array<int, 2> a, array<int, 2> b) {
return a[0] < b[0] ? a : (a[0] > b[0] ? b : array<int, 2> {a[0], a[1] + b[1]});
}
void apply(int x, int value) {
tree[x][0] += value;
push[x] += value;
}
void propagate(int x) {
if (push[x] != 0) {
apply(x * 2, push[x]);
apply(x * 2 + 1, push[x]);
push[x] = 0;
}
}
segment_tree(int n) : size(2 << __lg(n)), tree(2 * size, {INF, 0}), push(2 * size) {
fill(tree.begin() + size, tree.begin() + size + n, array<int, 2> {0, 1});
for (int i = size - 1; i > 0; i--) {
tree[i] = merge(tree[i * 2], tree[i * 2 + 1]);
}
}
void update(int l, int r, int value) {
if (l < r) {
update(l, r, value, 1, 0, size);
}
}
void update(int l, int r, int value, int x, int lx, int rx) {
if (lx >= r || rx <= l) {
return;
} else if (lx >= l && rx <= r) {
apply(x, value);
} else {
propagate(x);
int mid = (lx + rx) / 2;
update(l, r, value, x * 2, lx, mid);
update(l, r, value, x * 2 + 1, mid, rx);
tree[x] = merge(tree[x * 2], tree[x * 2 + 1]);
}
}
int query() {
return (tree[1][0] > 0 ? 0 : tree[1][1]);
}
};
int main() {
cin.tie(nullptr)->sync_with_stdio(false);
int N, K;
cin >> N >> K;
vector<vector<array<int, 3>>> rect(N + 1);
for (int i = 0; i < K; i++) {
int x1, y1, x2, y2;
cin >> x1 >> y1 >> x2 >> y2;
x1--, y1--, y2--;
rect[x1].push_back({y1, y2, 1});
rect[x2].push_back({y1, y2, -1});
}
auto Even = [&](int L, int R, int d) {
int l = L % (2 * d) >= d ? L + d - L % d : L;
int r = R % (2 * d) >= d ? R - R % d - 1 : R;
return l > r ? pair {0, -1} : pair {l / (2 * d) * d + l % (2 * d), r / (2 * d) * d + r % (2 * d)};
};
auto Odd = [&](int L, int R, int d) {
int l = L % (2 * d) < d ? L + d - L % d : L;
int r = R % (2 * d) < d ? R - R % d - 1 : R;
return l > r ? pair {0, -1} : pair {l / (2 * d) * d + (l - d) % (2 * d), r / (2 * d) * d + (r - d) % (2 * d)};
};
long long answer = numeric_limits<long long>::max();
for (int d = 1; d < N; d++) {
if (N % d == 0) {
int N_even = N / (2 * d) * d + min(d, N % (2 * d)), N_odd = N / (2 * d) * d + max(0, N % (2 * d) - d);
segment_tree st_even(N_even), st_odd(N_odd);
long long answer_e = 0, answer_o = 0;
for (int i = 0; i < N; i++) {
for (auto [j1, j2, t] : rect[i]) {
auto [j1_e, j2_e] = Even(j1, j2, d);
auto [j1_o, j2_o] = Odd(j1, j2, d);
st_even.update(j1_e, j2_e + 1, t);
st_odd.update(j1_o, j2_o + 1, t);
}
int cnt_e = st_even.query(), cnt_o = st_odd.query();
int r = (i / d) % 2;
answer_o += (r ? cnt_e : N_even - cnt_e) + (!r ? cnt_o : N_odd - cnt_o);
answer_e += (!r ? cnt_e : N_even - cnt_e) + (r ? cnt_o : N_odd - cnt_o);
}
answer = min({answer, answer_e, answer_o});
}
}
cout << answer << "\n";
}
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
344 KB |
Output is correct |
2 |
Correct |
0 ms |
348 KB |
Output is correct |
3 |
Correct |
0 ms |
348 KB |
Output is correct |
4 |
Correct |
0 ms |
348 KB |
Output is correct |
5 |
Correct |
0 ms |
348 KB |
Output is correct |
6 |
Correct |
0 ms |
348 KB |
Output is correct |
7 |
Correct |
0 ms |
416 KB |
Output is correct |
8 |
Correct |
0 ms |
348 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
46 ms |
5976 KB |
Output is correct |
2 |
Correct |
11 ms |
2140 KB |
Output is correct |
3 |
Correct |
34 ms |
8180 KB |
Output is correct |
4 |
Correct |
24 ms |
3416 KB |
Output is correct |
5 |
Correct |
41 ms |
6744 KB |
Output is correct |
6 |
Correct |
31 ms |
5212 KB |
Output is correct |
7 |
Correct |
7 ms |
2140 KB |
Output is correct |
8 |
Correct |
28 ms |
5468 KB |
Output is correct |
9 |
Correct |
63 ms |
7684 KB |
Output is correct |
10 |
Correct |
40 ms |
5976 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
344 KB |
Output is correct |
2 |
Correct |
0 ms |
348 KB |
Output is correct |
3 |
Correct |
1 ms |
348 KB |
Output is correct |
4 |
Correct |
2 ms |
348 KB |
Output is correct |
5 |
Correct |
1 ms |
348 KB |
Output is correct |
6 |
Correct |
1 ms |
348 KB |
Output is correct |
7 |
Correct |
1 ms |
348 KB |
Output is correct |
8 |
Correct |
1 ms |
488 KB |
Output is correct |
9 |
Correct |
1 ms |
348 KB |
Output is correct |
10 |
Correct |
0 ms |
348 KB |
Output is correct |
11 |
Correct |
1 ms |
348 KB |
Output is correct |
12 |
Correct |
1 ms |
348 KB |
Output is correct |
13 |
Correct |
1 ms |
348 KB |
Output is correct |
14 |
Correct |
1 ms |
408 KB |
Output is correct |
15 |
Correct |
1 ms |
348 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
344 KB |
Output is correct |
2 |
Correct |
0 ms |
348 KB |
Output is correct |
3 |
Correct |
1 ms |
348 KB |
Output is correct |
4 |
Correct |
2 ms |
348 KB |
Output is correct |
5 |
Correct |
1 ms |
348 KB |
Output is correct |
6 |
Correct |
1 ms |
348 KB |
Output is correct |
7 |
Correct |
1 ms |
348 KB |
Output is correct |
8 |
Correct |
1 ms |
488 KB |
Output is correct |
9 |
Correct |
1 ms |
348 KB |
Output is correct |
10 |
Correct |
0 ms |
348 KB |
Output is correct |
11 |
Correct |
1 ms |
348 KB |
Output is correct |
12 |
Correct |
1 ms |
348 KB |
Output is correct |
13 |
Correct |
1 ms |
348 KB |
Output is correct |
14 |
Correct |
1 ms |
408 KB |
Output is correct |
15 |
Correct |
1 ms |
348 KB |
Output is correct |
16 |
Correct |
26 ms |
1484 KB |
Output is correct |
17 |
Correct |
39 ms |
3624 KB |
Output is correct |
18 |
Correct |
95 ms |
3708 KB |
Output is correct |
19 |
Correct |
624 ms |
3728 KB |
Output is correct |
20 |
Correct |
643 ms |
3732 KB |
Output is correct |
21 |
Correct |
36 ms |
3928 KB |
Output is correct |
22 |
Correct |
2 ms |
348 KB |
Output is correct |
23 |
Correct |
79 ms |
2156 KB |
Output is correct |
24 |
Correct |
82 ms |
3676 KB |
Output is correct |
25 |
Correct |
19 ms |
864 KB |
Output is correct |
26 |
Correct |
74 ms |
2396 KB |
Output is correct |
27 |
Correct |
102 ms |
3416 KB |
Output is correct |
28 |
Correct |
86 ms |
3700 KB |
Output is correct |
29 |
Correct |
26 ms |
1884 KB |
Output is correct |
30 |
Correct |
3 ms |
604 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
46 ms |
5976 KB |
Output is correct |
2 |
Correct |
11 ms |
2140 KB |
Output is correct |
3 |
Correct |
34 ms |
8180 KB |
Output is correct |
4 |
Correct |
24 ms |
3416 KB |
Output is correct |
5 |
Correct |
41 ms |
6744 KB |
Output is correct |
6 |
Correct |
31 ms |
5212 KB |
Output is correct |
7 |
Correct |
7 ms |
2140 KB |
Output is correct |
8 |
Correct |
28 ms |
5468 KB |
Output is correct |
9 |
Correct |
63 ms |
7684 KB |
Output is correct |
10 |
Correct |
40 ms |
5976 KB |
Output is correct |
11 |
Correct |
1 ms |
344 KB |
Output is correct |
12 |
Correct |
0 ms |
348 KB |
Output is correct |
13 |
Correct |
1 ms |
348 KB |
Output is correct |
14 |
Correct |
2 ms |
348 KB |
Output is correct |
15 |
Correct |
1 ms |
348 KB |
Output is correct |
16 |
Correct |
1 ms |
348 KB |
Output is correct |
17 |
Correct |
1 ms |
348 KB |
Output is correct |
18 |
Correct |
1 ms |
488 KB |
Output is correct |
19 |
Correct |
1 ms |
348 KB |
Output is correct |
20 |
Correct |
0 ms |
348 KB |
Output is correct |
21 |
Correct |
1 ms |
348 KB |
Output is correct |
22 |
Correct |
1 ms |
348 KB |
Output is correct |
23 |
Correct |
1 ms |
348 KB |
Output is correct |
24 |
Correct |
1 ms |
408 KB |
Output is correct |
25 |
Correct |
1 ms |
348 KB |
Output is correct |
26 |
Correct |
26 ms |
1484 KB |
Output is correct |
27 |
Correct |
39 ms |
3624 KB |
Output is correct |
28 |
Correct |
95 ms |
3708 KB |
Output is correct |
29 |
Correct |
624 ms |
3728 KB |
Output is correct |
30 |
Correct |
643 ms |
3732 KB |
Output is correct |
31 |
Correct |
36 ms |
3928 KB |
Output is correct |
32 |
Correct |
2 ms |
348 KB |
Output is correct |
33 |
Correct |
79 ms |
2156 KB |
Output is correct |
34 |
Correct |
82 ms |
3676 KB |
Output is correct |
35 |
Correct |
19 ms |
864 KB |
Output is correct |
36 |
Correct |
74 ms |
2396 KB |
Output is correct |
37 |
Correct |
102 ms |
3416 KB |
Output is correct |
38 |
Correct |
86 ms |
3700 KB |
Output is correct |
39 |
Correct |
26 ms |
1884 KB |
Output is correct |
40 |
Correct |
3 ms |
604 KB |
Output is correct |
41 |
Correct |
976 ms |
12784 KB |
Output is correct |
42 |
Correct |
176 ms |
11968 KB |
Output is correct |
43 |
Correct |
491 ms |
12708 KB |
Output is correct |
44 |
Correct |
190 ms |
13204 KB |
Output is correct |
45 |
Correct |
99 ms |
12244 KB |
Output is correct |
46 |
Correct |
1097 ms |
11924 KB |
Output is correct |
47 |
Correct |
76 ms |
11424 KB |
Output is correct |
48 |
Correct |
334 ms |
11552 KB |
Output is correct |
49 |
Correct |
149 ms |
11260 KB |
Output is correct |
50 |
Execution timed out |
2041 ms |
11372 KB |
Time limit exceeded |
51 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
344 KB |
Output is correct |
2 |
Correct |
0 ms |
348 KB |
Output is correct |
3 |
Correct |
0 ms |
348 KB |
Output is correct |
4 |
Correct |
0 ms |
348 KB |
Output is correct |
5 |
Correct |
0 ms |
348 KB |
Output is correct |
6 |
Correct |
0 ms |
348 KB |
Output is correct |
7 |
Correct |
0 ms |
416 KB |
Output is correct |
8 |
Correct |
0 ms |
348 KB |
Output is correct |
9 |
Correct |
46 ms |
5976 KB |
Output is correct |
10 |
Correct |
11 ms |
2140 KB |
Output is correct |
11 |
Correct |
34 ms |
8180 KB |
Output is correct |
12 |
Correct |
24 ms |
3416 KB |
Output is correct |
13 |
Correct |
41 ms |
6744 KB |
Output is correct |
14 |
Correct |
31 ms |
5212 KB |
Output is correct |
15 |
Correct |
7 ms |
2140 KB |
Output is correct |
16 |
Correct |
28 ms |
5468 KB |
Output is correct |
17 |
Correct |
63 ms |
7684 KB |
Output is correct |
18 |
Correct |
40 ms |
5976 KB |
Output is correct |
19 |
Correct |
1 ms |
344 KB |
Output is correct |
20 |
Correct |
0 ms |
348 KB |
Output is correct |
21 |
Correct |
1 ms |
348 KB |
Output is correct |
22 |
Correct |
2 ms |
348 KB |
Output is correct |
23 |
Correct |
1 ms |
348 KB |
Output is correct |
24 |
Correct |
1 ms |
348 KB |
Output is correct |
25 |
Correct |
1 ms |
348 KB |
Output is correct |
26 |
Correct |
1 ms |
488 KB |
Output is correct |
27 |
Correct |
1 ms |
348 KB |
Output is correct |
28 |
Correct |
0 ms |
348 KB |
Output is correct |
29 |
Correct |
1 ms |
348 KB |
Output is correct |
30 |
Correct |
1 ms |
348 KB |
Output is correct |
31 |
Correct |
1 ms |
348 KB |
Output is correct |
32 |
Correct |
1 ms |
408 KB |
Output is correct |
33 |
Correct |
1 ms |
348 KB |
Output is correct |
34 |
Correct |
26 ms |
1484 KB |
Output is correct |
35 |
Correct |
39 ms |
3624 KB |
Output is correct |
36 |
Correct |
95 ms |
3708 KB |
Output is correct |
37 |
Correct |
624 ms |
3728 KB |
Output is correct |
38 |
Correct |
643 ms |
3732 KB |
Output is correct |
39 |
Correct |
36 ms |
3928 KB |
Output is correct |
40 |
Correct |
2 ms |
348 KB |
Output is correct |
41 |
Correct |
79 ms |
2156 KB |
Output is correct |
42 |
Correct |
82 ms |
3676 KB |
Output is correct |
43 |
Correct |
19 ms |
864 KB |
Output is correct |
44 |
Correct |
74 ms |
2396 KB |
Output is correct |
45 |
Correct |
102 ms |
3416 KB |
Output is correct |
46 |
Correct |
86 ms |
3700 KB |
Output is correct |
47 |
Correct |
26 ms |
1884 KB |
Output is correct |
48 |
Correct |
3 ms |
604 KB |
Output is correct |
49 |
Correct |
976 ms |
12784 KB |
Output is correct |
50 |
Correct |
176 ms |
11968 KB |
Output is correct |
51 |
Correct |
491 ms |
12708 KB |
Output is correct |
52 |
Correct |
190 ms |
13204 KB |
Output is correct |
53 |
Correct |
99 ms |
12244 KB |
Output is correct |
54 |
Correct |
1097 ms |
11924 KB |
Output is correct |
55 |
Correct |
76 ms |
11424 KB |
Output is correct |
56 |
Correct |
334 ms |
11552 KB |
Output is correct |
57 |
Correct |
149 ms |
11260 KB |
Output is correct |
58 |
Execution timed out |
2041 ms |
11372 KB |
Time limit exceeded |
59 |
Halted |
0 ms |
0 KB |
- |