oj.uz

Submission #128392

# Submission time Handle Problem Language Result Execution time Memory
128392 2019-07-10T21:11:50 Z qkxwsm New Home (APIO18_new_home) C++14
47 / 100
 5000 ms 408160 KB 
#include <bits/stdc++.h>

using namespace std;

template<class T, class U>
void ckmin(T &a, U b)
{
	if (a > b) a = b;
}
template<class T, class U>
void ckmax(T &a, U b)
{
	if (a < b) a = b;
}

#define MP make_pair
#define PB push_back
#define LB lower_bound
#define UB upper_bound
#define fi first
#define se second
#define FOR(i, a, b) for (auto i = (a); i < (b); i++)
#define FORD(i, a, b) for (auto i = (a) - 1; i >= (b); i--)
#define SZ(x) ((int) ((x).size()))
#define ALL(x) (x).begin(), (x).end()
#define LIM 100000007
#define INF 1000000007
#define LLINF 2696969696969696969ll
#define MAXN 2100013

typedef long long ll;
typedef long double ld;
typedef pair<int, int> pii;
typedef pair<ll, ll> pll;
typedef vector<int> vi;
typedef vector<ll> vl;
typedef vector<pii> vpi;
typedef vector<pll> vpl;
typedef pair<pii, pii> ppp;
typedef pair<pii, int> ppi;

int N, K, Q, S;
vector<ppp> arr;
vector<ppi> quer;
int ans[MAXN];
multiset<int> pts[MAXN];
vector<array<int, 4> > segments;
vector<array<int, 4> > events;
map<pii, vi> cancer;
vpi seg[MAXN];

void ins(int l, int r, int t)
{
	// cerr << "ins " << l << ' ' << r << " time " << t << endl;
	l = (l + r + 1) >> 1;
	cancer[{l, r}].PB(SZ(segments));
	segments.PB({l, r, t, SZ(events)});
}
void del(int l, int r, int t)
{
	// cerr << "del " << l << ' ' << r << " time " << t << endl;
	l = (l + r + 1) >> 1;
	// cerr << "rem " << (cancer[{l, r}]).size() << endl;
	segments[cancer[{l, r}].back()][3] = t;
	cancer[{l, r}].pop_back();
}

void inseg(int w, int L, int R, int a, int b, pii p)
{
	if (a <= L && R <= b)
	{
		seg[w].PB(p);
		return;
	}
	int mid = (L + R) >> 1;
	if (a <= mid) inseg(w << 1, L, mid, a, b, p);
	if (mid < b) inseg(w << 1 | 1, mid + 1, R, a, b, p);
	return;
}
void build(int w, int L, int R)
{
	if (L != R)
	{
		int mid = (L + R) >> 1;
		build(w << 1, L, mid);
		build(w << 1 | 1, mid + 1, R);
	}
	sort(ALL(seg[w]));
	FOR(i, 1, SZ(seg[w]))
	{
		ckmax(seg[w][i].se, seg[w][i - 1].se);
	}
}
int query(int w, int L, int R, int a, int x)
{
	int res = -INF;
	int ind = UB(ALL(seg[w]), MP(x, INF)) - seg[w].begin() - 1;
	if (ind != -1) res = -x + seg[w][ind].se;
	if (L == R) return res;
	int mid = (L + R) >> 1;
	if (a <= mid) ckmax(res, query(w << 1, L, mid, a, x));
	else ckmax(res, query(w << 1 | 1, mid + 1, R, a, x));
	return res;
}

void solve()
{
	//consider distnace functions with negative slope ONLY.
	//arr: something of type t and position p in times (t1, t2)
	//quer: query of position p at time t, qid i
	FOR(i, 0, (1 << (33 - __builtin_clz(SZ(events)))))
	{
		seg[i].clear();
	}
	events.clear();
	cancer.clear();
	segments.clear();
	FOR(i, 0, K) pts[i].clear();
	FOR(i, 0, SZ(arr))
	{
		int t1 = arr[i].se.fi, t2 = arr[i].se.se, ty = arr[i].fi.fi, x = arr[i].fi.se;
		events.PB({t1, -1, ty, x});
		events.PB({t2 + 1, -2, ty, x});
	}
	FOR(i, 0, SZ(quer))
	{
		int t = quer[i].fi.se, x = quer[i].fi.fi, qid = quer[i].se;
		events.PB({t, 0, qid, x});
	}
	sort(ALL(events));
	FOR(i, 0, K)
	{
		pts[i].insert(2 * LIM);
		pts[i].insert(-2 * LIM);
		ins(-2 * LIM, 2 * LIM, 0);
	}
	FOR(i, 0, SZ(events))
	{
		int typ = events[i][2], pos = events[i][3];
		if (events[i][1] == -1)
		{
			auto it = pts[typ].UB(pos);
			int rt = *it, lt = *prev(it);
			pts[typ].insert(pos);
			//then between lt and rt you generate some new segments!
			//at time i you have to erase a segment! and then insert two new segments!
			//WARNING: IF NOT DISTINCT, PLEASE CHECK THIS PLACE
			ins(lt, pos, i + 1);
			ins(pos, rt, i + 1);
			del(lt, rt, i + 1);
		}
		if (events[i][1] == -2)
		{
			pts[typ].erase(pts[typ].find(pos));
			auto it = pts[typ].UB(pos);
			int rt = *it, lt = *prev(it);
			del(lt, pos, i + 1);
			del(pos, rt, i + 1);
			ins(lt, rt, i + 1);
		}
	}
	FOR(i, 0, SZ(segments))
	{
		int l = segments[i][2], r = segments[i][3];
		inseg(1, 0, SZ(events), l, r, {segments[i][0], segments[i][1]});
	}
	//ok cool we have the segments
	build(1, 0, SZ(events));
	FOR(i, 0, SZ(events))
	{
		if (events[i][1] == 0)
		{
			int qid = events[i][2], pos = events[i][3];
			ckmax(ans[qid], query(1, 0, SZ(events), i + 1, pos));
			// cerr << "ans " << qid << ' ' << ans[qid] << endl;
		}
	}
}
bool cmp(ppi a, ppi b)
{
	return a.fi.se < b.fi.se;
}

int32_t main()
{
	if (fopen("file.in", "r"))
	{
		freopen("file.in", "r", stdin);
	}
	ios_base::sync_with_stdio(0); cin.tie(0);
	cin >> N >> K >> Q;
	arr.resize(N); quer.resize(Q);
	FOR(i, 0, N)
	{
		cin >> arr[i].fi.se >> arr[i].fi.fi >> arr[i].se.fi >> arr[i].se.se;
		arr[i].fi.fi--;
	}
	sort(ALL(arr));
	FOR(i, 0, Q)
	{
		cin >> quer[i].fi.fi >> quer[i].fi.se; quer[i].se = i;
	}
	sort(ALL(quer), cmp);
	solve();
	FOR(i, 0, N)
	{
		arr[i].fi.se = 100000001 - arr[i].fi.se;
	}
	FOR(i, 0, Q)
	{
		quer[i].fi.fi = 100000001 - quer[i].fi.fi;
	}
	solve();
	//events: OPEN/CLOSE a store at time t, position p, type x, ASK at time t, position p
	FOR(i, 0, Q)
	{
		if (ans[i] >= LIM)
		{
			cout << "-1\n";
			continue;
		}
		// cerr << "ans " << ans[i] << endl;
		cout << ans[i] << '\n';
	}
	return 0;
}

Compilation message

new_home.cpp: In function 'int32_t main()':
new_home.cpp:188:10: warning: ignoring return value of 'FILE* freopen(const char*, const char*, FILE*)', declared with attribute warn_unused_result [-Wunused-result]
   freopen("file.in", "r", stdin);
   ~~~~~~~^~~~~~~~~~~~~~~~~~~~~~~

Subtask #1
5.0 / 5.0

# Verdict Execution time Memory Grader output
1 Correct 141 ms 148260 KB Output is correct
2 Correct 140 ms 148396 KB Output is correct
3 Correct 141 ms 148216 KB Output is correct
4 Correct 140 ms 148344 KB Output is correct
5 Correct 158 ms 148584 KB Output is correct
6 Correct 148 ms 148728 KB Output is correct
7 Correct 152 ms 148756 KB Output is correct
8 Correct 156 ms 148728 KB Output is correct
9 Correct 164 ms 148600 KB Output is correct
10 Correct 146 ms 148608 KB Output is correct
11 Correct 142 ms 148600 KB Output is correct
12 Correct 143 ms 148600 KB Output is correct
13 Correct 161 ms 148600 KB Output is correct
14 Correct 143 ms 148600 KB Output is correct
15 Correct 144 ms 148600 KB Output is correct
16 Correct 143 ms 148608 KB Output is correct
17 Correct 142 ms 148680 KB Output is correct
18 Correct 144 ms 148728 KB Output is correct
19 Correct 159 ms 148600 KB Output is correct
20 Correct 159 ms 148600 KB Output is correct
21 Correct 141 ms 148624 KB Output is correct
22 Correct 145 ms 148728 KB Output is correct
23 Correct 146 ms 148728 KB Output is correct
24 Correct 144 ms 148728 KB Output is correct
25 Correct 142 ms 148600 KB Output is correct
26 Correct 142 ms 148524 KB Output is correct
27 Correct 163 ms 148444 KB Output is correct
28 Correct 164 ms 148572 KB Output is correct
29 Correct 142 ms 148600 KB Output is correct
30 Correct 145 ms 148504 KB Output is correct

Subtask #2
7.0 / 7.0

# Verdict Execution time Memory Grader output
1 Correct 141 ms 148260 KB Output is correct
2 Correct 140 ms 148396 KB Output is correct
3 Correct 141 ms 148216 KB Output is correct
4 Correct 140 ms 148344 KB Output is correct
5 Correct 158 ms 148584 KB Output is correct
6 Correct 148 ms 148728 KB Output is correct
7 Correct 152 ms 148756 KB Output is correct
8 Correct 156 ms 148728 KB Output is correct
9 Correct 164 ms 148600 KB Output is correct
10 Correct 146 ms 148608 KB Output is correct
11 Correct 142 ms 148600 KB Output is correct
12 Correct 143 ms 148600 KB Output is correct
13 Correct 161 ms 148600 KB Output is correct
14 Correct 143 ms 148600 KB Output is correct
15 Correct 144 ms 148600 KB Output is correct
16 Correct 143 ms 148608 KB Output is correct
17 Correct 142 ms 148680 KB Output is correct
18 Correct 144 ms 148728 KB Output is correct
19 Correct 159 ms 148600 KB Output is correct
20 Correct 159 ms 148600 KB Output is correct
21 Correct 141 ms 148624 KB Output is correct
22 Correct 145 ms 148728 KB Output is correct
23 Correct 146 ms 148728 KB Output is correct
24 Correct 144 ms 148728 KB Output is correct
25 Correct 142 ms 148600 KB Output is correct
26 Correct 142 ms 148524 KB Output is correct
27 Correct 163 ms 148444 KB Output is correct
28 Correct 164 ms 148572 KB Output is correct
29 Correct 142 ms 148600 KB Output is correct
30 Correct 145 ms 148504 KB Output is correct
31 Correct 1860 ms 205944 KB Output is correct
32 Correct 510 ms 171236 KB Output is correct
33 Correct 1749 ms 204384 KB Output is correct
34 Correct 1849 ms 204524 KB Output is correct
35 Correct 1883 ms 205744 KB Output is correct
36 Correct 1955 ms 205068 KB Output is correct
37 Correct 1249 ms 200360 KB Output is correct
38 Correct 1237 ms 200184 KB Output is correct
39 Correct 915 ms 193252 KB Output is correct
40 Correct 992 ms 194896 KB Output is correct
41 Correct 1389 ms 198796 KB Output is correct
42 Correct 1350 ms 199268 KB Output is correct
43 Correct 473 ms 175328 KB Output is correct
44 Correct 1254 ms 198432 KB Output is correct
45 Correct 1140 ms 194540 KB Output is correct
46 Correct 1022 ms 188560 KB Output is correct
47 Correct 683 ms 186040 KB Output is correct
48 Correct 669 ms 184656 KB Output is correct
49 Correct 858 ms 189540 KB Output is correct
50 Correct 1071 ms 196544 KB Output is correct
51 Correct 784 ms 187872 KB Output is correct

Subtask #3
0 / 10.0

# Verdict Execution time Memory Grader output
1 Execution timed out 5126 ms 405816 KB Time limit exceeded
2 Halted 0 ms 0 KB -

Subtask #4
0 / 23.0

# Verdict Execution time Memory Grader output
1 Execution timed out 5134 ms 408160 KB Time limit exceeded
2 Halted 0 ms 0 KB -

Subtask #5
35.0 / 35.0

# Verdict Execution time Memory Grader output
1 Correct 141 ms 148260 KB Output is correct
2 Correct 140 ms 148396 KB Output is correct
3 Correct 141 ms 148216 KB Output is correct
4 Correct 140 ms 148344 KB Output is correct
5 Correct 158 ms 148584 KB Output is correct
6 Correct 148 ms 148728 KB Output is correct
7 Correct 152 ms 148756 KB Output is correct
8 Correct 156 ms 148728 KB Output is correct
9 Correct 164 ms 148600 KB Output is correct
10 Correct 146 ms 148608 KB Output is correct
11 Correct 142 ms 148600 KB Output is correct
12 Correct 143 ms 148600 KB Output is correct
13 Correct 161 ms 148600 KB Output is correct
14 Correct 143 ms 148600 KB Output is correct
15 Correct 144 ms 148600 KB Output is correct
16 Correct 143 ms 148608 KB Output is correct
17 Correct 142 ms 148680 KB Output is correct
18 Correct 144 ms 148728 KB Output is correct
19 Correct 159 ms 148600 KB Output is correct
20 Correct 159 ms 148600 KB Output is correct
21 Correct 141 ms 148624 KB Output is correct
22 Correct 145 ms 148728 KB Output is correct
23 Correct 146 ms 148728 KB Output is correct
24 Correct 144 ms 148728 KB Output is correct
25 Correct 142 ms 148600 KB Output is correct
26 Correct 142 ms 148524 KB Output is correct
27 Correct 163 ms 148444 KB Output is correct
28 Correct 164 ms 148572 KB Output is correct
29 Correct 142 ms 148600 KB Output is correct
30 Correct 145 ms 148504 KB Output is correct
31 Correct 1860 ms 205944 KB Output is correct
32 Correct 510 ms 171236 KB Output is correct
33 Correct 1749 ms 204384 KB Output is correct
34 Correct 1849 ms 204524 KB Output is correct
35 Correct 1883 ms 205744 KB Output is correct
36 Correct 1955 ms 205068 KB Output is correct
37 Correct 1249 ms 200360 KB Output is correct
38 Correct 1237 ms 200184 KB Output is correct
39 Correct 915 ms 193252 KB Output is correct
40 Correct 992 ms 194896 KB Output is correct
41 Correct 1389 ms 198796 KB Output is correct
42 Correct 1350 ms 199268 KB Output is correct
43 Correct 473 ms 175328 KB Output is correct
44 Correct 1254 ms 198432 KB Output is correct
45 Correct 1140 ms 194540 KB Output is correct
46 Correct 1022 ms 188560 KB Output is correct
47 Correct 683 ms 186040 KB Output is correct
48 Correct 669 ms 184656 KB Output is correct
49 Correct 858 ms 189540 KB Output is correct
50 Correct 1071 ms 196544 KB Output is correct
51 Correct 784 ms 187872 KB Output is correct
52 Correct 1515 ms 211484 KB Output is correct
53 Correct 1376 ms 204484 KB Output is correct
54 Correct 1800 ms 206608 KB Output is correct
55 Correct 1363 ms 204712 KB Output is correct
56 Correct 1499 ms 207096 KB Output is correct
57 Correct 1398 ms 200692 KB Output is correct
58 Correct 1500 ms 204016 KB Output is correct
59 Correct 1482 ms 206792 KB Output is correct
60 Correct 1360 ms 200656 KB Output is correct
61 Correct 640 ms 201480 KB Output is correct
62 Correct 1517 ms 216732 KB Output is correct
63 Correct 1757 ms 209336 KB Output is correct
64 Correct 1707 ms 208868 KB Output is correct
65 Correct 1649 ms 206708 KB Output is correct
66 Correct 1421 ms 201200 KB Output is correct
67 Correct 873 ms 183140 KB Output is correct

Subtask #6
0 / 20.0

# Verdict Execution time Memory Grader output
1 Correct 141 ms 148260 KB Output is correct
2 Correct 140 ms 148396 KB Output is correct
3 Correct 141 ms 148216 KB Output is correct
4 Correct 140 ms 148344 KB Output is correct
5 Correct 158 ms 148584 KB Output is correct
6 Correct 148 ms 148728 KB Output is correct
7 Correct 152 ms 148756 KB Output is correct
8 Correct 156 ms 148728 KB Output is correct
9 Correct 164 ms 148600 KB Output is correct
10 Correct 146 ms 148608 KB Output is correct
11 Correct 142 ms 148600 KB Output is correct
12 Correct 143 ms 148600 KB Output is correct
13 Correct 161 ms 148600 KB Output is correct
14 Correct 143 ms 148600 KB Output is correct
15 Correct 144 ms 148600 KB Output is correct
16 Correct 143 ms 148608 KB Output is correct
17 Correct 142 ms 148680 KB Output is correct
18 Correct 144 ms 148728 KB Output is correct
19 Correct 159 ms 148600 KB Output is correct
20 Correct 159 ms 148600 KB Output is correct
21 Correct 141 ms 148624 KB Output is correct
22 Correct 145 ms 148728 KB Output is correct
23 Correct 146 ms 148728 KB Output is correct
24 Correct 144 ms 148728 KB Output is correct
25 Correct 142 ms 148600 KB Output is correct
26 Correct 142 ms 148524 KB Output is correct
27 Correct 163 ms 148444 KB Output is correct
28 Correct 164 ms 148572 KB Output is correct
29 Correct 142 ms 148600 KB Output is correct
30 Correct 145 ms 148504 KB Output is correct
31 Correct 1860 ms 205944 KB Output is correct
32 Correct 510 ms 171236 KB Output is correct
33 Correct 1749 ms 204384 KB Output is correct
34 Correct 1849 ms 204524 KB Output is correct
35 Correct 1883 ms 205744 KB Output is correct
36 Correct 1955 ms 205068 KB Output is correct
37 Correct 1249 ms 200360 KB Output is correct
38 Correct 1237 ms 200184 KB Output is correct
39 Correct 915 ms 193252 KB Output is correct
40 Correct 992 ms 194896 KB Output is correct
41 Correct 1389 ms 198796 KB Output is correct
42 Correct 1350 ms 199268 KB Output is correct
43 Correct 473 ms 175328 KB Output is correct
44 Correct 1254 ms 198432 KB Output is correct
45 Correct 1140 ms 194540 KB Output is correct
46 Correct 1022 ms 188560 KB Output is correct
47 Correct 683 ms 186040 KB Output is correct
48 Correct 669 ms 184656 KB Output is correct
49 Correct 858 ms 189540 KB Output is correct
50 Correct 1071 ms 196544 KB Output is correct
51 Correct 784 ms 187872 KB Output is correct
52 Execution timed out 5126 ms 405816 KB Time limit exceeded
53 Halted 0 ms 0 KB -