#include <iostream>
#include <vector>
#include <set>
#include <unordered_set>
#include <algorithm>
#define debug(x) cerr << #x << " " << x << "\n"
#define debugs(x) cerr << #x << " " << x << " "
#pragma GCC optimize("Ofast")
using namespace std;
typedef long long ll;
typedef pair <int, int> pii;
const int bSize = 11;
const int bits = 23;
const int NMAX = 500001;
int MOD;
int dm;
class Aint {
public:
int aint[NMAX * 4];
void init() {
for(int i = 0; i < NMAX * 4; i++)
aint[i] = 1;
}
void add(int node, int st, int dr, int l, int x) {
if(st == dr) {
aint[node] += x;
aint[node] = max(aint[node], 1);
return;
}
int mid = (st + dr) / 2;
if(l <= mid) {
add(node * 2, st, mid, l, x);
} else {
add(node * 2 + 1, mid + 1, dr, l, x);
}
aint[node] = aint[node * 2] * aint[node * 2 + 1];
aint[node] %= MOD;
}
int query(int node, int st, int dr, int a, int b){
if(a <= st && dr <= b)
return aint[node];
int mid = (st + dr) / 2;
int prod = 1;
if(a <= mid){
prod *= query(node * 2, st, mid, a, b);
prod %= MOD;
}
if(b > mid){
prod *= query(node * 2 + 1, mid + 1, dr, a, b);
prod %= MOD;
}
return prod;
}
void update(int node, int st, int dr, int l, int x) {
if(st == dr) {
dm = aint[node];
aint[node] = x;
//aint[node] = max(aint[node], 1);
return;
}
int mid = (st + dr) / 2;
if(l <= mid) {
update(node * 2, st, mid, l, x);
} else {
update(node * 2 + 1, mid + 1, dr, l, x);
}
aint[node] = (aint[node * 2] % MOD) * (aint[node * 2 + 1] % MOD);
aint[node] %= MOD;
}
} aint, banned;
int imp[NMAX];
int col[NMAX];
int maxi[NMAX];
pii v[NMAX];
int f[NMAX];
pii cine[NMAX];
vector <int> g[NMAX];
int main() {
ios_base::sync_with_stdio(false);
cin.tie(0);
cout.tie(0);
int n, k, i;
cin >> n >> k >> MOD;
for(i = 1; i <= n; i++) {
cin >> v[i].first >> v[i].second;
g[v[i].second].push_back(v[i].first);
}
sort(v + 1, v + n + 1);
for(i = n; i >= 1; i--) {
if(!col[v[i].second]) {
imp[i] = 1;
}
col[v[i].second]++;
}
int cnt = 0;
for(i = 1; i <= n; i++){
if(imp[i]){
cine[++cnt] = v[i];
f[cine[cnt].second] = cnt;
}
}
aint.init();
banned.init();
for(i = 1; i <= n; i++){
sort(g[i].begin(), g[i].end());
maxi[i] = col[i];
if(!f[i]) continue;
aint.update(1, 1, n, f[i], col[i] + 1);
banned.update(1, 1, n, f[i], col[i] + 1);
}
int dr = n, st = n;
int rez = 0;
for(i = n; i >= 1; i--){
while(st >= 1 && v[st].first * 2 > v[i].first){
aint.add(1, 1, n, f[v[st].second], -1);
banned.add(1, 1, n, f[v[st].second], -1);
col[v[st].second]--;
st--;
}
if(i == n){
for(int j = 1; j <= n; j++){
maxi[j] = col[j];
}
}
if(imp[i]){
if(i != n){
aint.update(1, 1, n, f[v[i].second], 1);
int r = f[v[i].second];
int pas = (1 << bits);
int a = 0;
if(col[v[i].second] > 0)
a = g[v[i].second][col[v[i].second] - 1];
int b = g[v[i].second][col[v[i].second]];
while(pas){
if(r + pas <= k && a * 2 <= cine[r + pas].first && b * 2 > cine[r + pas].first)
r += pas;
pas /= 2;
}
rez += aint.query(1, 1, n, 1, r);
aint.update(1, 1, n, f[v[i].second], col[v[i].second] + 1);
banned.add(1, 1, n, f[v[i].second], -1);
rez %= MOD;
}
if(col[v[i].second] > 0 || i == n) /// avem frecventa mai mare ca 1..0
rez += banned.aint[1];
rez %= MOD;
banned.update(1, 1, n, f[v[i].second], 1);
}
}
cout << rez;
return 0;
}
Compilation message
fish.cpp: In function 'int main()':
fish.cpp:119:9: warning: unused variable 'dr' [-Wunused-variable]
119 | int dr = n, st = n;
| ^~
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
13 ms |
27732 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
13 ms |
27728 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
14 ms |
27644 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
15 ms |
27768 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
14 ms |
27700 KB |
Output is correct |
2 |
Correct |
14 ms |
27696 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
14 ms |
27732 KB |
Output is correct |
2 |
Correct |
341 ms |
36364 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
16 ms |
27732 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
15 ms |
27828 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
139 ms |
31648 KB |
Output is correct |
2 |
Correct |
161 ms |
32348 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
16 ms |
27860 KB |
Output is correct |
2 |
Correct |
18 ms |
27956 KB |
Output is correct |
3 |
Correct |
17 ms |
27860 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
199 ms |
33604 KB |
Output is correct |
2 |
Incorrect |
323 ms |
36232 KB |
Output isn't correct |
3 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
297 ms |
37700 KB |
Output is correct |
2 |
Correct |
326 ms |
36944 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
193 ms |
34920 KB |
Output is correct |
2 |
Correct |
339 ms |
36964 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
323 ms |
36712 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
388 ms |
38236 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
395 ms |
36456 KB |
Output is correct |
2 |
Correct |
593 ms |
41968 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
752 ms |
42080 KB |
Output is correct |
2 |
Correct |
583 ms |
39760 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
602 ms |
43172 KB |
Output is correct |
2 |
Correct |
660 ms |
41900 KB |
Output is correct |
3 |
Correct |
704 ms |
46064 KB |
Output is correct |
4 |
Correct |
595 ms |
41976 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
916 ms |
45004 KB |
Output is correct |
2 |
Correct |
1319 ms |
59088 KB |
Output is correct |
3 |
Correct |
1538 ms |
58972 KB |
Output is correct |
4 |
Correct |
1319 ms |
54580 KB |
Output is correct |