# | Submission time | Handle | Problem | Language | Result | Execution time | Memory |
---|---|---|---|---|---|---|---|
22096 | 2017-04-29T08:27:35 Z | 삼*전자 그린픽스(#1005, pichulia) | None (KRIII5P_3) | C++ | 129 ms | 9312 KB |
#include<stdio.h> #include<algorithm> #define M 1000000007 #define I 262144 using namespace std; typedef pair<int, int> pii; int n; int x[2*I]; int xl; pii a[I]; long long int tw[I]; void input() { int i, j, k, l; xl = 0; x[xl++] = 0; x[xl++] = M; scanf("%d", &n); for (i = 0; i < n; i++) { scanf("%d %d", &a[i].second, &a[i].first); if (a[i].second >= a[i].first) { n--; i--; continue; } x[xl++] = a[i].first; x[xl++] = a[i].second; } sort(a, a + n); sort(x, x + xl); xl = unique(x, x + xl) - x; for (i = 0; i < n; i++) { //change first and second k = a[i].first; a[i].first = lower_bound(x, x + xl, a[i].second) - x; a[i].second = lower_bound(x, x + xl, k) - x; } tw[0] = 1; for (i = 1; i < I; i++) tw[i] = (tw[i - 1] * 2) % M; } long long int resa, resb; int ms[2 * I]; int get_ms(int dep, int ql, int qr, int ll, int rr) { if (qr < ll || rr < ql) return 0; if (ql <= ll && rr <= qr) return ms[dep]; return get_ms(dep * 2, ql, qr, ll, (ll + rr) / 2) + get_ms(dep * 2 + 1, ql, qr, (ll + rr) / 2 + 1, rr); } void update_ms(int i) { i += I; ms[i]++; i /= 2; while (i > 0) { ms[i] = ms[i * 2] + ms[i * 2 + 1]; i /= 2; } } long long int func(int xi) { int i, j, k; // return 2^s0 + 2^s1 + .. +2^sx[i] - 1 long long int res = 0; //naive solution if (0) { for (i = 1; i < xi; i++) { long long int cur; long long int dist = x[i + 1] - x[i]; long long int val = get_ms(1, 0, i, 0, I - 1); cur = (dist * tw[val]) % M; res = (res + cur); } } return res; } void process() { int i, j, k, l; /* resa += (2^sq - 1) * q - { (2^sp - 1) * p (2^sp+1 - 2^sp) * (p+1) ... (2^sq-1 - 2^sq-2) * (q-1) (2^sq - 2^sq-1) * q } = p - q + 2^sp + 2^sp+! + ... + 2^sq-2 + 2^sq-1 resb += (2^sq-1 - 1) + 1 */ resa = 0; resb = 0; long long int cura, curb; for (i = n - 1; i >= 0; i--) { cura = 0; cura = (cura + func(a[i].second)) % M; cura = (cura - func(a[i].first) + M) % M; k = get_ms(1, 0, a[i].second - 1, 0, I - 1); curb = tw[k]; resa = (resa + cura) % M; resb = (resb + curb) % M; // printf("%d %d --> %lld %lld\n", x[a[i].first], x[a[i].second], cura, curb); update_ms(a[i].first); } printf("%lld %lld\n", resa, resb); } int main() { int i, j, k, l; input(); process(); return 0; }
Compilation message
# | Verdict | Execution time | Memory | Grader output |
---|---|---|---|---|
1 | Partially correct | 3 ms | 9312 KB | Output is partially correct |
2 | Partially correct | 3 ms | 9312 KB | Output is partially correct |
3 | Partially correct | 3 ms | 9312 KB | Output is partially correct |
4 | Partially correct | 0 ms | 9312 KB | Output is partially correct |
5 | Partially correct | 0 ms | 9312 KB | Output is partially correct |
6 | Partially correct | 9 ms | 9312 KB | Output is partially correct |
7 | Partially correct | 3 ms | 9312 KB | Output is partially correct |
8 | Partially correct | 6 ms | 9312 KB | Output is partially correct |
9 | Partially correct | 6 ms | 9312 KB | Output is partially correct |
10 | Partially correct | 6 ms | 9312 KB | Output is partially correct |
11 | Partially correct | 83 ms | 9312 KB | Output is partially correct |
12 | Partially correct | 89 ms | 9312 KB | Output is partially correct |
13 | Partially correct | 66 ms | 9312 KB | Output is partially correct |
14 | Partially correct | 73 ms | 9312 KB | Output is partially correct |
15 | Partially correct | 69 ms | 9312 KB | Output is partially correct |
16 | Partially correct | 66 ms | 9312 KB | Output is partially correct |
17 | Partially correct | 79 ms | 9312 KB | Output is partially correct |
18 | Partially correct | 66 ms | 9312 KB | Output is partially correct |
19 | Partially correct | 66 ms | 9312 KB | Output is partially correct |
20 | Partially correct | 49 ms | 9312 KB | Output is partially correct |
21 | Partially correct | 106 ms | 9312 KB | Output is partially correct |
22 | Partially correct | 126 ms | 9312 KB | Output is partially correct |
23 | Partially correct | 103 ms | 9312 KB | Output is partially correct |
24 | Partially correct | 106 ms | 9312 KB | Output is partially correct |
25 | Partially correct | 83 ms | 9312 KB | Output is partially correct |
26 | Partially correct | 113 ms | 9312 KB | Output is partially correct |
27 | Partially correct | 113 ms | 9312 KB | Output is partially correct |
28 | Partially correct | 99 ms | 9312 KB | Output is partially correct |
29 | Partially correct | 89 ms | 9312 KB | Output is partially correct |
30 | Partially correct | 96 ms | 9312 KB | Output is partially correct |
31 | Partially correct | 106 ms | 9312 KB | Output is partially correct |
32 | Partially correct | 79 ms | 9312 KB | Output is partially correct |
33 | Partially correct | 89 ms | 9312 KB | Output is partially correct |
34 | Partially correct | 83 ms | 9312 KB | Output is partially correct |
35 | Partially correct | 86 ms | 9312 KB | Output is partially correct |
36 | Partially correct | 109 ms | 9312 KB | Output is partially correct |
37 | Partially correct | 116 ms | 9312 KB | Output is partially correct |
38 | Partially correct | 129 ms | 9312 KB | Output is partially correct |
39 | Partially correct | 103 ms | 9312 KB | Output is partially correct |
40 | Partially correct | 106 ms | 9312 KB | Output is partially correct |
41 | Partially correct | 0 ms | 9312 KB | Output is partially correct |
42 | Partially correct | 3 ms | 9312 KB | Output is partially correct |
43 | Partially correct | 0 ms | 9312 KB | Output is partially correct |
44 | Partially correct | 6 ms | 9312 KB | Output is partially correct |
45 | Partially correct | 36 ms | 9312 KB | Output is partially correct |
46 | Partially correct | 36 ms | 9312 KB | Output is partially correct |
47 | Partially correct | 36 ms | 9312 KB | Output is partially correct |
48 | Partially correct | 66 ms | 9312 KB | Output is partially correct |
49 | Partially correct | 76 ms | 9312 KB | Output is partially correct |
50 | Partially correct | 96 ms | 9312 KB | Output is partially correct |