#include "hexagon.h"
#include <stdlib.h>
#include <vector>
using namespace std;
typedef vector<int> vi;
const int N = 200000, MD = 1000000007, V2 = 500000004, V6 = 166666668;
int min(int a, int b) { return a < b ? a : b; }
int max(int a, int b) { return a > b ? a : b; }
unsigned int X = 12345;
int rand_() {
return (X *= 3) >> 1;
}
int dx[] = { 1, 1, 0, -1, -1, 0 };
int dy[] = { 0, 1, 1, 0, -1, -1 };
int choose2(int n) {
return (long long) n * (n - 1) % MD * V2 % MD;
}
int choose3(int n) {
return (long long) n * (n - 1) % MD * (n - 2) % MD * V6 % MD;
}
long long cross(int x1, int y1, int x2, int y2) {
return (long long) x1 * y2 - (long long) x2 * y1;
}
int hh[N], ll[N], ii[N], xx[N][2], yy[N], n, n_;
void add(long long x1, long long y1, long long x2, long long y2) {
if (y1 < y2)
xx[n_][0] = x1, xx[n_][1] = x2, yy[n_] = y1, n_++;
else if (y1 > y2)
xx[n_][0] = x2, xx[n_][1] = x1, yy[n_] = y2, n_++;
}
int compare(int i, int j) {
int c = yy[i] != yy[j] ? yy[i] - yy[j] : (xx[i][0] + xx[i][1]) - (xx[j][0] + xx[j][1]);
return c == 0 ? 0 : (c < 0 ? -1 : 1);
}
void sort(int *ii, int l, int r) {
while (l < r) {
int i = l, j = l, k = r, i_ = ii[l + rand_() % (r - l)], tmp;
while (j < k) {
int c = compare(ii[j], i_);
if (c == 0)
j++;
else if (c < 0) {
tmp = ii[i], ii[i] = ii[j], ii[j] = tmp;
i++, j++;
} else {
k--;
tmp = ii[j], ii[j] = ii[k], ii[k] = tmp;
}
}
sort(ii, l, i);
l = k;
}
}
int *ej[N], eo[N];
void append(int i, int j) {
int o = eo[i]++;
if (o >= 2 && (o & o - 1) == 0)
ej[i] = (int *) realloc(ej[i], o * 2 * sizeof *ej[i]);
ej[i][o] = j;
}
int area(int xi1, int xi2, int xj1, int xj2, int z, int l) {
int d = ((xj2 - xj1) - (xi2 - xi1)) / l, h, ans;
ans = 0;
for (h = 1; h < l; h++)
ans = (ans + (long long) (xi2 - xi1 + 1 + d * h) * (z + h)) % MD;
return ans;
}
int dfs(int p, int i, int d, int l) {
int i_, s, xi1, xi2, o, x;
i_ = i >> 1, s = i & 1, xi1 = xx[ii[i_ << 1 | 0]][s], xi2 = xx[ii[i_ << 1 | 1]][s];
x = (long long) d * (xi2 - xi1 + 1) % MD;
for (o = eo[i]; o--; ) {
int j = ej[i][o];
if (j != p) {
int j_ = j >> 1, t = j & 1, xj1 = xx[ii[j_ << 1 | 0]][t], xj2 = xx[ii[j_ << 1 | 1]][t];
if ((i ^ j) == 1)
x = (x + dfs(i, j, d + l, l) + area(xi1, xi2, xj1, xj2, d, l)) % MD;
else {
x = (x + dfs(i, j, d, l)) % MD;
x = (x - (long long) d * (min(xi2, xj2) - max(xi1, xj1) + 1)) % MD;
}
}
}
return x;
}
int solve() {
int l_, h, i, j, k, l, x, y, ans;
l_ = ll[0];
for (i = 0; i < n; i++)
if (ll[i] != ll[0]) {
l_ = 1;
break;
}
if (l_ == 1) {
n_ = 0, x = 0, y = 0;
for (i = 0; i < n; i++)
for (l = 0; l < ll[i]; l++)
add(x, y, x + dx[hh[i]], y + dy[hh[i]]), x += dx[hh[i]], y += dy[hh[i]];
} else {
n_ = 0, x = 0, y = 0;
for (i = 0; i < n; i++)
add(x, y, x + dx[hh[i]] * l_, y + dy[hh[i]] * l_), x += dx[hh[i]] * l_, y += dy[hh[i]] * l_;
}
for (i = 0; i < n_; i++)
ii[i] = i;
sort(ii, 0, n_);
n_ /= 2;
for (i = 0; i < n_ * 2; i++)
ej[i] = (int *) malloc(2 * sizeof *ej[i]), eo[i] = 0;
for (i = 0; i < n_; i++)
append(i << 1 | 0, i << 1 | 1), append(i << 1 | 1, i << 1 | 0);
for (h = -1, i = 0; i < n_; h = i, i = j) {
int y = yy[ii[i << 1 | 0]];
j = i + 1;
while (j < n_ && yy[ii[j << 1 | 0]] == y)
j++;
if (h != -1) {
k = h, l = i;
while (k < i && l < j)
if (xx[ii[k << 1 | 1]][1] < xx[ii[l << 1 | 1]][0]) {
if (xx[ii[l << 1 | 0]][0] <= xx[ii[k << 1 | 1]][1])
append(k << 1 | 1, l << 1 | 0), append(l << 1 | 0, k << 1 | 1);
k++;
} else {
if (xx[ii[k << 1 | 0]][1] <= xx[ii[l << 1 | 1]][0])
append(k << 1 | 1, l << 1 | 0), append(l << 1 | 0, k << 1 | 1);
l++;
}
}
}
ans = 0;
for (i = 0; i < n_; i++) {
int i0 = ii[i << 1 | 0], i1 = ii[i << 1 | 1];
if (yy[i0] == 0 && xx[i0][0] <= 0 && 0 <= xx[i1][0]) {
ans = dfs(-1, i << 1 | 0, 0, l_);
for (i = 0; i < n_ * 2; i++)
free(ej[i]);
break;
}
if (yy[i1] + l_ == 0 && xx[i0][1] <= 0 && 0 <= xx[i1][1]) {
ans = dfs(-1, i << 1 | 1, 0, l_);
for (i = 0; i < n_ * 2; i++)
free(ej[i]);
break;
}
}
return ans;
}
int draw_territory(int n_, int a, int b, vi hh_, vi ll_) {
int h, i, x, y, ans1, ans2;
long long area2, boundary, internal;
n = n_;
for (i = 0; i < n; i++)
hh[i] = hh_[i] - 1, ll[i] = ll_[i];
x = 0, y = 0, area2 = 0, boundary = 0;
for (i = 0; i < n; i++) {
h = hh[i];
area2 += cross(x, y, x + dx[h] * ll[i], y + dy[h] * ll[i]);
boundary += ll[i];
x += dx[h] * ll[i], y += dy[h] * ll[i];
}
if (area2 < 0)
area2 = -area2;
internal = (area2 - boundary) / 2 + 1;
ans1 = (boundary + internal) % MD, ans2 = 0;
if (n == 3)
ans2 = choose3(ll[0] + 2) * 2 % MD;
else if (b > 0) {
ans2 = 0;
for (h = 0; h < 3; h++) {
ans2 = (ans2 + solve()) % MD;
for (i = 0; i < n; i++)
hh[i] = (hh[i] + 1) % 6;
}
if (ans2 < 0)
ans2 += MD;
ans2 = (long long) ans2 * V2 % MD;
}
return ((long long) ans1 * a + (long long) ans2 * b) % MD;
}
Compilation message
hexagon.cpp: In function 'void append(int, int)':
hexagon.cpp:77:23: warning: suggest parentheses around '-' in operand of '&' [-Wparentheses]
77 | if (o >= 2 && (o & o - 1) == 0)
| ~~^~~
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
204 KB |
Output is correct |
2 |
Correct |
1 ms |
204 KB |
Output is correct |
3 |
Correct |
0 ms |
204 KB |
Output is correct |
4 |
Correct |
0 ms |
204 KB |
Output is correct |
5 |
Correct |
1 ms |
204 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
204 KB |
Output is correct |
2 |
Correct |
0 ms |
204 KB |
Output is correct |
3 |
Correct |
0 ms |
204 KB |
Output is correct |
4 |
Correct |
1 ms |
204 KB |
Output is correct |
5 |
Correct |
0 ms |
204 KB |
Output is correct |
6 |
Correct |
1 ms |
204 KB |
Output is correct |
7 |
Correct |
0 ms |
204 KB |
Output is correct |
8 |
Correct |
0 ms |
204 KB |
Output is correct |
9 |
Correct |
0 ms |
204 KB |
Output is correct |
10 |
Correct |
0 ms |
204 KB |
Output is correct |
11 |
Correct |
0 ms |
204 KB |
Output is correct |
12 |
Correct |
1 ms |
204 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
204 KB |
Output is correct |
2 |
Correct |
1 ms |
332 KB |
Output is correct |
3 |
Correct |
1 ms |
332 KB |
Output is correct |
4 |
Correct |
1 ms |
332 KB |
Output is correct |
5 |
Correct |
2 ms |
460 KB |
Output is correct |
6 |
Correct |
1 ms |
332 KB |
Output is correct |
7 |
Correct |
1 ms |
332 KB |
Output is correct |
8 |
Correct |
1 ms |
332 KB |
Output is correct |
9 |
Correct |
1 ms |
460 KB |
Output is correct |
10 |
Correct |
2 ms |
460 KB |
Output is correct |
11 |
Correct |
1 ms |
460 KB |
Output is correct |
12 |
Correct |
1 ms |
540 KB |
Output is correct |
13 |
Correct |
1 ms |
460 KB |
Output is correct |
14 |
Correct |
1 ms |
460 KB |
Output is correct |
15 |
Correct |
1 ms |
460 KB |
Output is correct |
16 |
Correct |
1 ms |
588 KB |
Output is correct |
17 |
Correct |
1 ms |
588 KB |
Output is correct |
18 |
Correct |
0 ms |
204 KB |
Output is correct |
19 |
Correct |
0 ms |
204 KB |
Output is correct |
20 |
Correct |
0 ms |
204 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
204 KB |
Output is correct |
2 |
Correct |
2 ms |
460 KB |
Output is correct |
3 |
Correct |
1 ms |
332 KB |
Output is correct |
4 |
Correct |
1 ms |
204 KB |
Output is correct |
5 |
Correct |
1 ms |
332 KB |
Output is correct |
6 |
Correct |
2 ms |
460 KB |
Output is correct |
7 |
Correct |
6 ms |
972 KB |
Output is correct |
8 |
Correct |
1 ms |
204 KB |
Output is correct |
9 |
Correct |
1 ms |
204 KB |
Output is correct |
10 |
Correct |
0 ms |
204 KB |
Output is correct |
11 |
Correct |
9 ms |
1448 KB |
Output is correct |
12 |
Correct |
10 ms |
1484 KB |
Output is correct |
13 |
Correct |
9 ms |
1356 KB |
Output is correct |
14 |
Correct |
13 ms |
1532 KB |
Output is correct |
15 |
Correct |
1 ms |
332 KB |
Output is correct |
16 |
Correct |
1 ms |
204 KB |
Output is correct |
17 |
Correct |
0 ms |
204 KB |
Output is correct |
18 |
Correct |
0 ms |
204 KB |
Output is correct |
19 |
Correct |
0 ms |
204 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
204 KB |
Output is correct |
2 |
Correct |
0 ms |
204 KB |
Output is correct |
3 |
Correct |
0 ms |
204 KB |
Output is correct |
4 |
Correct |
2 ms |
460 KB |
Output is correct |
5 |
Correct |
1 ms |
332 KB |
Output is correct |
6 |
Correct |
0 ms |
204 KB |
Output is correct |
7 |
Correct |
1 ms |
332 KB |
Output is correct |
8 |
Correct |
2 ms |
460 KB |
Output is correct |
9 |
Correct |
7 ms |
972 KB |
Output is correct |
10 |
Correct |
1 ms |
204 KB |
Output is correct |
11 |
Correct |
1 ms |
204 KB |
Output is correct |
12 |
Correct |
1 ms |
204 KB |
Output is correct |
13 |
Correct |
12 ms |
1356 KB |
Output is correct |
14 |
Correct |
10 ms |
1484 KB |
Output is correct |
15 |
Correct |
9 ms |
1356 KB |
Output is correct |
16 |
Correct |
10 ms |
1472 KB |
Output is correct |
17 |
Correct |
1 ms |
332 KB |
Output is correct |
18 |
Correct |
0 ms |
204 KB |
Output is correct |
19 |
Correct |
0 ms |
204 KB |
Output is correct |
20 |
Correct |
17 ms |
2076 KB |
Output is correct |
21 |
Correct |
2 ms |
460 KB |
Output is correct |
22 |
Correct |
1 ms |
332 KB |
Output is correct |
23 |
Correct |
21 ms |
2800 KB |
Output is correct |
24 |
Correct |
35 ms |
4332 KB |
Output is correct |
25 |
Correct |
46 ms |
4548 KB |
Output is correct |
26 |
Correct |
17 ms |
2380 KB |
Output is correct |
27 |
Correct |
17 ms |
1692 KB |
Output is correct |
28 |
Correct |
11 ms |
1316 KB |
Output is correct |
29 |
Correct |
39 ms |
4920 KB |
Output is correct |
30 |
Correct |
39 ms |
4936 KB |
Output is correct |
31 |
Correct |
39 ms |
4928 KB |
Output is correct |
32 |
Correct |
41 ms |
4932 KB |
Output is correct |
33 |
Correct |
18 ms |
2508 KB |
Output is correct |
34 |
Correct |
9 ms |
1356 KB |
Output is correct |
35 |
Correct |
1 ms |
204 KB |
Output is correct |
36 |
Correct |
0 ms |
204 KB |
Output is correct |
37 |
Correct |
0 ms |
204 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
204 KB |
Output is correct |
2 |
Correct |
1 ms |
332 KB |
Output is correct |
3 |
Correct |
1 ms |
332 KB |
Output is correct |
4 |
Correct |
1 ms |
332 KB |
Output is correct |
5 |
Correct |
1 ms |
460 KB |
Output is correct |
6 |
Correct |
1 ms |
332 KB |
Output is correct |
7 |
Correct |
1 ms |
332 KB |
Output is correct |
8 |
Correct |
1 ms |
332 KB |
Output is correct |
9 |
Correct |
1 ms |
460 KB |
Output is correct |
10 |
Correct |
1 ms |
460 KB |
Output is correct |
11 |
Correct |
1 ms |
460 KB |
Output is correct |
12 |
Correct |
1 ms |
460 KB |
Output is correct |
13 |
Correct |
1 ms |
460 KB |
Output is correct |
14 |
Correct |
1 ms |
460 KB |
Output is correct |
15 |
Correct |
1 ms |
460 KB |
Output is correct |
16 |
Correct |
1 ms |
588 KB |
Output is correct |
17 |
Correct |
1 ms |
588 KB |
Output is correct |
18 |
Correct |
0 ms |
204 KB |
Output is correct |
19 |
Correct |
2 ms |
460 KB |
Output is correct |
20 |
Correct |
1 ms |
332 KB |
Output is correct |
21 |
Correct |
1 ms |
204 KB |
Output is correct |
22 |
Correct |
1 ms |
332 KB |
Output is correct |
23 |
Correct |
2 ms |
460 KB |
Output is correct |
24 |
Correct |
6 ms |
972 KB |
Output is correct |
25 |
Correct |
1 ms |
332 KB |
Output is correct |
26 |
Correct |
1 ms |
204 KB |
Output is correct |
27 |
Correct |
1 ms |
204 KB |
Output is correct |
28 |
Correct |
9 ms |
1356 KB |
Output is correct |
29 |
Correct |
10 ms |
1484 KB |
Output is correct |
30 |
Correct |
10 ms |
1356 KB |
Output is correct |
31 |
Correct |
10 ms |
1472 KB |
Output is correct |
32 |
Correct |
1 ms |
332 KB |
Output is correct |
33 |
Correct |
1 ms |
204 KB |
Output is correct |
34 |
Correct |
55 ms |
10820 KB |
Output is correct |
35 |
Correct |
53 ms |
5700 KB |
Output is correct |
36 |
Correct |
51 ms |
6080 KB |
Output is correct |
37 |
Correct |
77 ms |
9852 KB |
Output is correct |
38 |
Correct |
82 ms |
9668 KB |
Output is correct |
39 |
Correct |
86 ms |
11548 KB |
Output is correct |
40 |
Correct |
86 ms |
11276 KB |
Output is correct |
41 |
Correct |
138 ms |
22468 KB |
Output is correct |
42 |
Correct |
144 ms |
18576 KB |
Output is correct |
43 |
Correct |
126 ms |
16008 KB |
Output is correct |
44 |
Correct |
85 ms |
20500 KB |
Output is correct |
45 |
Correct |
77 ms |
17308 KB |
Output is correct |
46 |
Correct |
75 ms |
16920 KB |
Output is correct |
47 |
Correct |
128 ms |
22596 KB |
Output is correct |
48 |
Correct |
154 ms |
28612 KB |
Output is correct |
49 |
Correct |
144 ms |
33296 KB |
Output is correct |
50 |
Correct |
0 ms |
204 KB |
Output is correct |
51 |
Correct |
1 ms |
204 KB |
Output is correct |
52 |
Correct |
0 ms |
204 KB |
Output is correct |
53 |
Correct |
0 ms |
204 KB |
Output is correct |
54 |
Correct |
0 ms |
204 KB |
Output is correct |
55 |
Correct |
0 ms |
204 KB |
Output is correct |
56 |
Correct |
1 ms |
204 KB |
Output is correct |
57 |
Correct |
0 ms |
204 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
204 KB |
Output is correct |
2 |
Correct |
0 ms |
204 KB |
Output is correct |
3 |
Correct |
0 ms |
204 KB |
Output is correct |
4 |
Correct |
0 ms |
204 KB |
Output is correct |
5 |
Execution timed out |
2059 ms |
204 KB |
Time limit exceeded |
6 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
204 KB |
Output is correct |
2 |
Correct |
0 ms |
204 KB |
Output is correct |
3 |
Correct |
0 ms |
204 KB |
Output is correct |
4 |
Correct |
1 ms |
204 KB |
Output is correct |
5 |
Correct |
0 ms |
204 KB |
Output is correct |
6 |
Correct |
1 ms |
332 KB |
Output is correct |
7 |
Correct |
1 ms |
332 KB |
Output is correct |
8 |
Correct |
1 ms |
332 KB |
Output is correct |
9 |
Correct |
2 ms |
460 KB |
Output is correct |
10 |
Correct |
1 ms |
332 KB |
Output is correct |
11 |
Correct |
2 ms |
332 KB |
Output is correct |
12 |
Correct |
1 ms |
332 KB |
Output is correct |
13 |
Correct |
1 ms |
460 KB |
Output is correct |
14 |
Correct |
1 ms |
460 KB |
Output is correct |
15 |
Correct |
1 ms |
460 KB |
Output is correct |
16 |
Correct |
1 ms |
460 KB |
Output is correct |
17 |
Correct |
1 ms |
460 KB |
Output is correct |
18 |
Correct |
1 ms |
460 KB |
Output is correct |
19 |
Correct |
1 ms |
460 KB |
Output is correct |
20 |
Correct |
2 ms |
588 KB |
Output is correct |
21 |
Correct |
1 ms |
588 KB |
Output is correct |
22 |
Correct |
1 ms |
204 KB |
Output is correct |
23 |
Correct |
3 ms |
460 KB |
Output is correct |
24 |
Correct |
1 ms |
332 KB |
Output is correct |
25 |
Correct |
1 ms |
204 KB |
Output is correct |
26 |
Correct |
1 ms |
332 KB |
Output is correct |
27 |
Correct |
2 ms |
460 KB |
Output is correct |
28 |
Correct |
6 ms |
1048 KB |
Output is correct |
29 |
Correct |
1 ms |
204 KB |
Output is correct |
30 |
Correct |
1 ms |
204 KB |
Output is correct |
31 |
Correct |
1 ms |
204 KB |
Output is correct |
32 |
Correct |
9 ms |
1448 KB |
Output is correct |
33 |
Correct |
10 ms |
1356 KB |
Output is correct |
34 |
Correct |
9 ms |
1300 KB |
Output is correct |
35 |
Correct |
10 ms |
1484 KB |
Output is correct |
36 |
Correct |
1 ms |
204 KB |
Output is correct |
37 |
Correct |
0 ms |
204 KB |
Output is correct |
38 |
Correct |
0 ms |
204 KB |
Output is correct |
39 |
Correct |
16 ms |
2072 KB |
Output is correct |
40 |
Correct |
2 ms |
480 KB |
Output is correct |
41 |
Correct |
1 ms |
332 KB |
Output is correct |
42 |
Correct |
23 ms |
2792 KB |
Output is correct |
43 |
Correct |
35 ms |
4336 KB |
Output is correct |
44 |
Correct |
35 ms |
4548 KB |
Output is correct |
45 |
Correct |
17 ms |
2404 KB |
Output is correct |
46 |
Correct |
14 ms |
1740 KB |
Output is correct |
47 |
Correct |
9 ms |
1228 KB |
Output is correct |
48 |
Correct |
39 ms |
4932 KB |
Output is correct |
49 |
Correct |
39 ms |
4936 KB |
Output is correct |
50 |
Correct |
40 ms |
4932 KB |
Output is correct |
51 |
Correct |
50 ms |
4932 KB |
Output is correct |
52 |
Correct |
19 ms |
2456 KB |
Output is correct |
53 |
Correct |
9 ms |
1316 KB |
Output is correct |
54 |
Correct |
47 ms |
10772 KB |
Output is correct |
55 |
Correct |
51 ms |
5708 KB |
Output is correct |
56 |
Correct |
48 ms |
6088 KB |
Output is correct |
57 |
Correct |
79 ms |
9792 KB |
Output is correct |
58 |
Correct |
82 ms |
9668 KB |
Output is correct |
59 |
Correct |
85 ms |
11460 KB |
Output is correct |
60 |
Correct |
82 ms |
11276 KB |
Output is correct |
61 |
Correct |
143 ms |
22504 KB |
Output is correct |
62 |
Correct |
132 ms |
18636 KB |
Output is correct |
63 |
Correct |
123 ms |
16048 KB |
Output is correct |
64 |
Correct |
77 ms |
20492 KB |
Output is correct |
65 |
Correct |
78 ms |
17236 KB |
Output is correct |
66 |
Correct |
76 ms |
17168 KB |
Output is correct |
67 |
Correct |
125 ms |
22592 KB |
Output is correct |
68 |
Correct |
138 ms |
28656 KB |
Output is correct |
69 |
Correct |
158 ms |
33356 KB |
Output is correct |
70 |
Execution timed out |
2082 ms |
204 KB |
Time limit exceeded |
71 |
Halted |
0 ms |
0 KB |
- |