This submission is migrated from previous version of oj.uz, which used different machine for grading. This submission may have different result if resubmitted.
// #pragma GCC optimize("O3")
// #pragma GCC optimize("unroll-loops")
#include <bits/stdc++.h>
using namespace std;
#define ll long long
#ifdef _DEBUG
int find_best(int n);
std::vector<int> ask(int i);
#else
#include "prize.h"
#endif
int N;
vector<int> a, b;
int segSet(int l, int r, int i, int k, int x)
{
if (l > k || r <= k)
return b[i];
if (l + 1 == r)
return a[l] = b[i] = x;
int m = (l + r) / 2;
return b[i] = segSet(l, m, i * 2 + 1, k, x) + segSet(m, r, i * 2 + 2, k, x);
}
int segGet(int l, int r, int i, int ql, int qr)
{
if (ql <= l && qr >= r)
return b[i];
if (ql >= r || qr <= l)
return 0;
int m = (l + r) / 2;
return segGet(l, m, i * 2 + 1, ql, qr) + segGet(m, r, i * 2 + 2, ql, qr);
}
int asks = 0;
int getRecursive(int l, int r, int cnt_here, int cnt_l, int cnt_r)
{
if (l >= r)
return -1;
if (cnt_here - segGet(0, N, 0, l, r) <= 0)
return -1;
int m = (l + r) / 2;
int i = m;
vector<int> answer;
bool is_special = false;
int diff_l = 0;
int diff_r = 0;
do
{
if (i == r)
break;
is_special = false;
asks++;
assert(asks < 10000);
assert(i >= 0 && i < N);
answer = ask(i);
if (answer[0] + answer[1] == 0)
return i;
answer[0] -= cnt_l;
answer[1] -= cnt_r;
if (answer[0] + answer[1] != cnt_here)
{
is_special = true;
i++;
}
else
{
answer[0] -= i - m;
diff_l += i - m;
diff_r += i - m;
}
} while (is_special);
int left = -1;
if (i == r)
left = getRecursive(l, m, cnt_here - (i - m), cnt_l, cnt_r + (i - m));
else
left = getRecursive(l, m, answer[0], cnt_l, cnt_r + answer[1] + diff_r);
if (left != -1)
return left;
return getRecursive(i + 1, r, answer[1], cnt_l + answer[0] + diff_l, cnt_r);
}
int find_best(int n)
{
N = n;
a.resize(n);
b.resize(4 * n);
vector<pair<int, int>> cnt_here;
set<int> values;
for (int k = 0; k < 750; k++)
{
int i = rand() % n;
vector<int> answer = ask(i);
if (answer[0] + answer[1] == 0)
return i;
values.insert(answer[0] + answer[1]);
cnt_here.push_back({answer[0] + answer[1], i});
if (values.size() >= 6)
break;
}
sort(cnt_here.begin(), cnt_here.end());
for (auto [x, i] : cnt_here) {
if (x != cnt_here.back().first) {
segSet(0, N, 0, i, 1);
}
}
int result = getRecursive(0, n, cnt_here.back().first, 0, 0);
assert(result != -1);
vector<int> answer = ask(result);
assert(answer[0] + answer[1] == 0);
return result;
}
#ifdef _DEBUG
#include <iostream>
#include <vector>
#include <algorithm>
using namespace std;
static const int max_q = 10000;
static int n;
static int query_count = 0;
static vector<int> g;
static vector<vector<int>> rank_count;
vector<int> ask(int i)
{
query_count++;
if (query_count > max_q)
{
cerr << "Query limit exceeded" << endl;
exit(0);
}
if (i < 0 || i >= n)
{
cerr << "Bad index: " << i << endl;
exit(0);
}
vector<int> res(2);
res[0] = rank_count[g[i] - 1][i + 1];
res[1] = rank_count[g[i] - 1][n] - res[0];
return res;
}
int main()
{
for (int t = 10; t <= 97; t++)
{
string path = "/home/elias/Downloads/ioi2017tests/prize/tests/2-" + to_string(t) + ".out";
freopen(path.c_str(), "r", stdin);
string blub;
cin >> blub;
cin >> n;
cin >> blub >> blub;
query_count = 0;
g.clear();
rank_count.clear();
asks = 0;
g.resize(n);
for (int i = 0; i < n; i++)
{
cin >> g[i];
if (g[i] < 1)
{
cerr << "Invalid rank " << g[i] << " at index " << i << endl;
exit(0);
}
}
int max_rank = *max_element(g.begin(), g.end());
rank_count.resize(max_rank + 1, vector<int>(n + 1, 0));
for (int r = 0; r <= max_rank; r++)
{
for (int i = 1; i <= n; i++)
{
rank_count[r][i] = rank_count[r][i - 1];
if (g[i - 1] == r)
rank_count[r][i]++;
}
}
for (int i = 0; i <= n; i++)
for (int r = 1; r <= max_rank; r++)
rank_count[r][i] += rank_count[r - 1][i];
int res = find_best(n);
assert(g[res] == 1);
cout << res << " " << g[res] << endl
<< "Query count: " << query_count << endl;
}
return 0;
}
#endif
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |