| # | Time | Username | Problem | Language | Result | Execution time | Memory |
|---|---|---|---|---|---|---|---|
| 955775 | emad234 | Horses (IOI15_horses) | C++17 | 0 ms | 0 KiB |
This submission is migrated from previous version of oj.uz, which used different machine for grading. This submission may have different result if resubmitted.
#include "horses.h"
#include <bits/stdc++.h>
#define ll long long
#define F first
#define S second
#define pii pair<ll, ll>
const ll mod = 1e9 + 7;
const ll mxN = 5e5 + 5;
using namespace std;
ll NN;
struct tree
{
vector<pii> seg;
ll s, e;
void init(ll n)
{
NN = exp2(ceil(log2(n)));
seg.resize(NN * 4);
}
ll queryMax(ll l = 1, ll r = NN, ll ind = 1)
{
if (l > e || r < s)
return 0;
if (l >= s && r <= e)
return seg[ind].F;
ll md = (l + r) / 2;
return max(queryMax(l, md, ind * 2), queryMax(md + 1, r, ind * 2 + 1));
}
ll queryProd(ll l = 1, ll r = NN, ll ind = 1)
{
if (l > e || r < s)
return 1;
if (l >= s && r <= e)
return seg[ind].S;
ll md = (l + r) / 2;
return (queryProd(l, md, ind * 2) * queryProd(md + 1, r, ind * 2 + 1)) % mod;
}
void update(pii val, ll ind)
{
ind += NN;
seg[ind] = val;
while (ind /= 2)
{
seg[ind].F = max(seg[ind * 2].F, seg[ind * 2 + 1].F);
seg[ind].S = (seg[ind * 2].S * seg[ind * 2 + 1].S) % mod;
}
}
} seg;
set<int> s;
ll st[mxN], ed[mxN];
ll x[mxN], y[mxN];
ll vals[mxN];
ll n;
ll solve()
{
ll id = 0;
int i = n - 1;
while (i > 0 && id <= 70)
{
if (x[i] == 1)
{
i = st[i];
id--;
}
i--;
id++;
}
if (i == -1)
i = 0;
ll mx = -1, mxid = -1;
ll prd = 1;
for (i; i < n; i++)
{
if (x[i] == 1)
{
seg.s = i + 1, seg.e = ed[i] + 1;
if (vals[i] * prd > mx)
{
mx = vals[i];
mxid = i;
prd = 1;
}
i = ed[i];
continue;
}
if (prd * x[i] > mx)
{
mx = y[i];
mxid = i;
prd = 1;
continue;
}
prd *= x[i];
if (prd * y[i] > mx)
{
mx = y[i];
mxid = i;
prd = 1;
continue;
}
}
seg.s = 1;
seg.e = mxid + 1;
return (mx * seg.queryProd()) % mod;
}
ll NN;
struct tree
{
vector<pii> seg;
ll s, e;
void init(ll n)
{
NN = exp2(ceil(log2(n)));
seg.resize(NN * 4);
}
ll queryMax(ll l = 1, ll r = NN, ll ind = 1)
{
if (l > e || r < s)
return 0;
if (l >= s && r <= e)
return seg[ind].F;
ll md = (l + r) / 2;
return max(queryMax(l, md, ind * 2), queryMax(md + 1, r, ind * 2 + 1));
}
ll queryProd(ll l = 1, ll r = NN, ll ind = 1)
{
if (l > e || r < s)
return 1;
if (l >= s && r <= e)
return seg[ind].S;
ll md = (l + r) / 2;
return (queryProd(l, md, ind * 2) * queryProd(md + 1, r, ind * 2 + 1)) % mod;
}
void update(pii val, ll ind)
{
ind += NN;
seg[ind] = val;
while (ind /= 2)
{
seg[ind].F = max(seg[ind * 2].F, seg[ind * 2 + 1].F);
seg[ind].S = (seg[ind * 2].S * seg[ind * 2 + 1].S) % mod;
}
}
} seg;
set<int> s;
ll st[mxN], ed[mxN];
ll x[mxN], y[mxN];
ll vals[mxN];
ll n;
ll solve()
{
ll id = 0;
int i = n - 1;
while (i > 0 && id <= 70)
{
if (x[i] == 1)
{
i = st[i];
id--;
}
i--;
id++;
}
if (i == -1)
i = 0;
ll mx = -1, mxid = -1;
ll prd = 1;
for (i; i < n; i++)
{
// fprintf(_outputFile, "%d ", i);
if (x[i] == 1)
{
seg.s = i + 1, seg.e = ed[i] + 1;
if (vals[i] * prd > mx)
{
mx = vals[i];
mxid = i;
prd = 1;
}
i = ed[i];
continue;
}
if (prd * x[i] > mx)
{
mx = y[i];
mxid = i;
prd = 1;
continue;
}
prd *= x[i];
if (prd * y[i] > mx)
{
mx = y[i];
mxid = i;
prd = 1;
continue;
}
}
// fprintf(_outputFile, "%d\n", 0);
seg.s = 1;
seg.e = mxid + 1;
return (mx * seg.queryProd()) % mod;
}
int init(int N, int X[], int Y[])
{
seg.init(N);
n = N;
for (ll i = 0; i < N; i++)
{
y[i] = Y[i];
x[i] = X[i];
seg.update({Y[i], X[i]}, i);
}
ll l = -1, r = -1;
for (ll i = 0; i < N; i++)
{
if (X[i] == 1)
{
if (l == -1)
l = i;
r = i;
}
else
{
if (l != -1)
{
s.insert(l);
st[r] = l;
ed[l] = r;
l = r = -1;
}
}
}
if (l != -1)
{
s.insert(l);
st[r] = l;
ed[l] = r;
l = r = -1;
}
// for (auto el : s)
// {
// fprintf(_outputFile, "%d ", el);
// fprintf(_outputFile, "%d\n", ed[el]);
// }
return solve();
}
int updateX(int pos, int val)
{
if (x[pos] != 1 && val == 1)
{
auto it = s.upper_bound(pos);
ll l = pos, r = pos;
if (it != s.end() && (*it) == pos + 1)
{
l = pos, r = ed[(*it)];
s.erase(it);
}
s.insert(l);
it = s.lower_bound(l);
if (it != s.begin())
{
it--;
if (ed[(*it)] == l - 1)
{
l = (*it);
s.erase(it);
s.erase(s.lower_bound(pos));
s.insert(l);
}
}
st[r] = l;
ed[l] = r;
seg.s = l + 1, seg.e = r + 1;
vals[l] = seg.queryMax();
}
if (x[pos] == 1 && val != 1)
{
auto it = s.upper_bound(pos);
it--;
ll l = (*it), r = pos - 1, l1 = pos + 1, r1 = ed[(*it)];
s.erase(it);
if (l <= r)
{
s.insert(l);
st[r] = l;
ed[l] = r;
seg.s = l + 1, seg.e = r + 1;
vals[l] = seg.queryMax();
}
if (l1 <= r1)
{
s.insert(l1);
st[r1] = l1;
ed[l1] = r1;
seg.s = l1 + 1, seg.e = r1 + 1;
vals[l1] = seg.queryMax();
}
}
x[pos] = val;
seg.update({y[pos], x[pos]}, pos);
// for (auto el : s)
// {
// fprintf(_outputFile, "%d ", el);
// fprintf(_outputFile, "%d\n", ed[el]);
// }
return solve();
}
int updateY(int pos, int val)
{
y[pos] = val;
seg.update({y[pos], x[pos]}, pos);
return solve();
}