#include <bits/stdc++.h>
using namespace std;
void just_do_it();
int main() {
#ifdef KAMIRULEZ
freopen("kamirulez.inp", "r", stdin);
freopen("kamirulez.out", "w", stdout);
#endif
ios_base::sync_with_stdio(0);
cin.tie(0);
just_do_it();
return 0;
}
const int maxn = 5e5 + 20;
pair<int, int> a[maxn];
vector<int> pos[maxn];
int cnt[maxn];
int tree[maxn];
int id[maxn];
int max_id[maxn];
int mod;
void build(int id, int lt, int rt) {
if (lt == rt) {
tree[id] = cnt[lt] % mod;
return;
}
int mt = (lt + rt) / 2;
build(id * 2, lt, mt);
build(id * 2 + 1, mt + 1, rt);
tree[id] = tree[id * 2] * tree[id * 2 + 1] % mod;
}
void update(int id, int lt, int rt, int pos, int val) {
if (lt == rt) {
tree[id] = val % mod;
return;
}
int mt = (lt + rt) / 2;
if (pos <= mt) {
update(id * 2, lt, mt, pos, val);
}
else {
update(id * 2 + 1, mt + 1, rt, pos, val);
}
tree[id] = tree[id * 2] * tree[id * 2 + 1] % mod;
}
int get(int id, int lt, int rt, int ql, int qr) {
if (ql > qr) {
return 1;
}
if (lt == ql && rt == qr) {
return tree[id];
}
int mt = (lt + rt) / 2;
if (qr <= mt) {
return get(id * 2, lt, mt, ql, qr);
}
else if (ql >= mt + 1) {
return get(id * 2 + 1, mt + 1, rt, ql, qr);
}
else {
return get(id * 2, lt, mt, ql, mt) * get(id * 2 + 1, mt + 1, rt, mt + 1, qr) % mod;
}
}
void just_do_it() {
int n, k;
cin >> n >> k >> mod;
for (int i = 1; i <= n; i++) {
cin >> a[i].first >> a[i].second;
}
sort(a + 1, a + n + 1);
k = 0;
for (int i = n; i >= 1; i--) {
if (!id[a[i].second]) {
id[a[i].second] = ++k;
}
a[i].second = id[a[i].second];
}
for (int i = 1; i <= n; i++) {
cnt[a[i].second]++;
pos[a[i].second].push_back(i);
}
for (int i = 1; i <= k; i++) {
cnt[i]++;
}
for (int i = n; i >= 1; i--) {
max_id[i] = max(max_id[i + 1], a[i].second);
}
build(1, 1, k);
int lt = n;
int res = 0;
for (int i = n; i >= 1; i--) {
while (lt > 0 && a[lt].first * 2 > a[i].first) {
cnt[a[lt].second]--;
update(1, 1, k, a[lt].second, cnt[a[lt].second]);
lt--;
}
if (max_id[i] != max_id[i + 1]) {
int x = max_id[i];
res += (cnt[x] - 1) % mod * get(1, 1, k, x + 1, k) % mod;
res %= mod;
int low = max_id[lower_bound(a + 1, a + n + 1, make_pair(a[pos[x][cnt[x] - 1]].first * 2, -1)) - a] + 1;
res += get(1, 1, k, low, x - 1) * get(1, 1, k, x + 1, k) % mod;
res %= mod;
}
}
cout << res;
}
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
7 ms |
11988 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
8 ms |
11988 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
7 ms |
11988 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
7 ms |
12104 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
8 ms |
12120 KB |
Output is correct |
2 |
Correct |
8 ms |
12084 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
8 ms |
12044 KB |
Output is correct |
2 |
Correct |
201 ms |
20468 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
7 ms |
12116 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
8 ms |
12192 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
83 ms |
15660 KB |
Output is correct |
2 |
Correct |
98 ms |
17160 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
8 ms |
12116 KB |
Output is correct |
2 |
Correct |
9 ms |
12248 KB |
Output is correct |
3 |
Correct |
10 ms |
12116 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
148 ms |
17664 KB |
Output is correct |
2 |
Correct |
217 ms |
21180 KB |
Output is correct |
3 |
Correct |
247 ms |
21792 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
190 ms |
20752 KB |
Output is correct |
2 |
Correct |
268 ms |
21936 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
123 ms |
17508 KB |
Output is correct |
2 |
Correct |
251 ms |
20684 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
221 ms |
20488 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
258 ms |
21868 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
239 ms |
20020 KB |
Output is correct |
2 |
Correct |
365 ms |
32416 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
445 ms |
25736 KB |
Output isn't correct |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
311 ms |
26588 KB |
Output isn't correct |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
491 ms |
27920 KB |
Output isn't correct |
2 |
Halted |
0 ms |
0 KB |
- |