Submission #502008

# Submission time Handle Problem Language Result Execution time Memory
502008 2022-01-05T05:53:40 Z 600Mihnea Fish (IOI08_fish) C++17
95 / 100
577 ms 65540 KB
#include <iostream>
#include <vector>
#include <algorithm>
#include <cstdio>

using namespace std;


struct Fish {
  int sz;
  int tp;
};

bool operator < (Fish a, Fish b) {
  return a.sz < b.sz;
}

const int N = (int) 5e5 + 7;
int n, k, mod;
Fish a[N];
vector<int> aux[N], oftp[N];

int caneat(int f, int s) {
  int low = 0, high = (int) oftp[s].size() - 1, sol = 0;
  while (low <= high) {
    int mid = (low + high) / 2;
    if (a[oftp[f].back()].sz >= 2 * a[oftp[s][mid]].sz) {
      sol = mid + 1;
      low = mid + 1;
    } else {
      high = mid - 1;
    }
  }
  return sol;
}

int aux_order[N], aux_inv_order[N];

bool cmp(int i, int j) {
  return aux[i].back() < aux[j].back();
}
const int TM = 1048576 + 7;

int tree[TM];

void upd(int v, int tl, int tr, int i, int x) {
  if (tr < i || i < tl) {
    return;
  }
  if (tl == tr) {
    tree[v] = x % mod;
    return;
  }
  int tm = (tl + tr) / 2;
  upd(2 * v, tl, tm, i, x);
  upd(2 * v + 1, tm + 1, tr, i, x);
  tree[v] = tree[2 * v] * tree[2 * v + 1] % mod;
}

int get(int v, int tl, int tr, int l, int r) {
  if (tr < l || r < tl) {
    return 1;
  }
  if (l <= tl && tr <= r) {
    return tree[v];
  }
  int tm = (tl + tr) / 2;
  return get(2 * v, tl, tm, l, r) * get(2 * v + 1, tm + 1, tr, l, r) % mod;
}

int nxt;

void upd(int i, int x) {
  upd(1, 1, nxt, i, x);
}

int get(int l, int r) {
  return get(1, 1, nxt, l, r);
}

int thesol[N], frq[N], lc[N], rc[N], full[N], semi[N];

int main() {
  ios::sync_with_stdio(0); cin.tie(0);

  for (int i = 1; i < TM; i++) {
    tree[i] = 1;
  }

 /// freopen ("input", "r", stdin);

  cin >> n >> k >> mod;
  int init_n = n;
  for (int i = 1; i <= n; i++) {
    cin >> a[i].sz >> a[i].tp;
  }
  bool debug = (a[1].tp == 374545);

  sort(a + 1, a + n + 1);


  for (int i = 1; i <= n; i++) {
    aux[a[i].tp].push_back(i);
  }
  n = k;
  for (int i = 1; i <= n; i++) {
    aux_order[i] = i;
  }
  nxt = 1;
  while (nxt < n) nxt *= 2;

  sort(aux_order + 1, aux_order + n + 1, cmp);

  for (int i = 1; i <= n; i++) aux_inv_order[aux_order[i]] = i;
  for (int i = 1; i <= n; i++) oftp[i] = aux[aux_order[i]];
  if (debug) exit(0);
  for (int i = 1; i <= n; i++) aux[i] = oftp[i];

  for (int i = 1; i <= init_n; i++) a[i].tp = aux_inv_order[a[i].tp];

  int go = 1, print = 0;

  if (debug) exit(0);
  for (int i = 1; i <= n; i++) {
    while (2 * a[go].sz <= a[oftp[i].back()].sz) {
      frq[a[go].tp]++;
      upd(a[go].tp, frq[a[go].tp] + 1);
      go++;
    }
    thesol[i] = get(1, i - 1);
    full[i] = 1;
    semi[i] = caneat(i, i) % mod;
    int last = i, low = i + 1, high = n;
    while (low <= high) {
      int j = (low + high) / 2;
      if (caneat(j, i) == caneat(i, i)) {
        last = j;
        low = j + 1;
      } else {
        high = j - 1;
      }
    }
    full[i] = get(i + 1, last);
    full[i] = full[i] * thesol[i] % mod;
    semi[i] = semi[i] * thesol[i] % mod;

    print = (print + full[i]) % mod;
    print = (print + semi[i]) % mod;
  }
  cout << print << "\n";

  return 0;
}

/**
**/
# Verdict Execution time Memory Grader output
1 Correct 14 ms 27852 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 18 ms 27920 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 13 ms 27852 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 13 ms 27852 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 15 ms 27904 KB Output is correct
2 Correct 13 ms 27980 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 14 ms 27916 KB Output is correct
2 Correct 171 ms 36004 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 17 ms 27984 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 15 ms 28044 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 70 ms 31652 KB Output is correct
2 Correct 104 ms 32608 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 18 ms 28060 KB Output is correct
2 Correct 17 ms 28108 KB Output is correct
3 Correct 19 ms 28100 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 119 ms 33612 KB Output is correct
2 Correct 196 ms 36688 KB Output is correct
3 Correct 185 ms 37464 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 176 ms 36852 KB Output is correct
2 Correct 194 ms 37172 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 123 ms 33612 KB Output is correct
2 Correct 193 ms 37300 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 192 ms 36672 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 225 ms 37848 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 184 ms 37260 KB Output is correct
2 Correct 274 ms 41932 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 441 ms 44764 KB Output is correct
2 Correct 332 ms 43876 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 368 ms 47752 KB Output is correct
2 Correct 383 ms 42412 KB Output is correct
3 Correct 410 ms 52468 KB Output is correct
4 Correct 389 ms 42464 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 577 ms 49284 KB Output is correct
2 Runtime error 360 ms 65540 KB Execution killed with signal 9
3 Halted 0 ms 0 KB -