Submission #446364

# Submission time Handle Problem Language Result Execution time Memory
446364 2021-07-21T16:38:44 Z ics0503 Constellation 3 (JOI20_constellation3) C++17
100 / 100
1192 ms 431088 KB
#pragma GCC optimize("O3")
#pragma GCC optimize("Ofast")
#pragma GCC optimize("unroll-loops")
#include<stdio.h>
#include<vector>
#include<algorithm>
#include<set>
using namespace std; 
int a[212121];
struct xy {
	int t,x,y,c;
}all[414141]; int an;
bool sort_y(xy a, xy b) {
	if (a.y != b.y)return a.y > b.y;
	if (a.t != b.t)return a.t < b.t;
	return a.x < b.x;
}
set<pair<pair<int, int>, int>>S; int tn;
vector<int>edge[414141];
vector<pair<int, int>>qr[414141];
int S_start[414141], S_end[414141];
vector<int>S_ER[414141]; int point[414141];
void insert_S(int s, int e,int parent) {
	tn++;
	S.insert({ { -s,e },tn });
	S_start[tn] = s; S_end[tn] = e;
	if (parent) edge[parent].push_back(tn);
}
int pp[414141][20];
long long sparse[414141][20];;
vector<int>need_to_fill[414141][20];
void init_sparse(int w, int bef) {
	pp[w][0] = bef;
	for (int i = 0; pp[w][i] != -1; i++) {
		need_to_fill[pp[w][i]][i].push_back(w);
		pp[w][i + 1] = pp[pp[w][i]][i];
	}
	for (int nxt : edge[w]) init_sparse(nxt, w);
}
long long D[414141], depth[414141], CS[414141];
void dfs(int w,int bef,int dep) {
	long long sum = 0;
	depth[w] = dep; 
	for (int nxt : edge[w]) {
		dfs(nxt, w, dep + 1);
		sum += D[nxt];
	}
	CS[w] = sum;
	for (int nxt : edge[w]) {
		long long v = sum - D[nxt];
		sparse[nxt][0] = v;
		for (int i = 1; i < 18; i++) {
			for (int who : need_to_fill[nxt][i]) {
				int p = pp[who][i - 1];
				sparse[who][i] = sparse[p][i - 1] + sparse[who][i - 1];
			}
		}
	}
	long long qmax = 0;
	for (auto q : qr[w]) {
		long long res = q.second;
		int now = point[q.first];
		res += CS[now];
		for (int i = 18; i >= 0; i--) {
			int nxt = pp[now][i];
			if (nxt == -1 || depth[nxt] <= depth[w])continue;
			res += sparse[now][i];
			now = nxt;
		}
		if(now!=w)res += sparse[now][0];
		qmax = max(qmax, res);
	}
	D[w] = max(sum, qmax);
}
int ck[414141];
vector<int>erL;
int main() {
	int n; scanf("%d", &n);
	int i, j;
	for (i = 1; i <= n; i++) {
		scanf("%d", &a[i]);
		all[an++] = { 1,i,a[i],0 };
	}
	int m; scanf("%d", &m);
	long long allS = 0;
	for (i = 1; i <= m; i++) {
		all[an].t = 2;
		scanf("%d%d%d", &all[an].x, &all[an].y, &all[an].c);
		allS += all[an].c;
		an++;
	}
	sort(all, all + an, sort_y);
	int st, en;
	insert_S(1, n, 0);
	int e9 = 1e9;
	for (st = 0; st < an; st = en) {
		for (en = st; all[st].y == all[en].y; en++);
		// [st, en)
		for (i = st; i < en; i++) {
			int t = all[i].t, x = all[i].x, y = all[i].y, c = all[i].c;
			if (t != 1) continue;			
			auto v = S.lower_bound({ { -x,-e9 },-e9 });
			int who = v->second;
			S_ER[who].push_back(x);
			if(ck[who]==0) erL.push_back(who);
			ck[who] = 1;
		}
		for (auto who : erL) {
			int s = S_start[who], e = S_end[who];
			for (i = 0; i < S_ER[who].size(); i++) {
				int x = S_ER[who][i];
				point[x] = who;
				int new_st = s, new_en = x - 1;
				if (new_st <= new_en) insert_S(new_st, new_en, who);
				s = x + 1;
				if (i + 1 == S_ER[who].size()) {
					if (s <= e) insert_S(s, e, who);
				}
			}
			ck[who] = 0;
			S.erase({ {-S_start[who],S_end[who] },who });
		}
		erL.clear();
		for (i = st; i < en; i++) {
			int t = all[i].t, x = all[i].x, y = all[i].y, c = all[i].c;
			if (t != 2) continue;
			auto v = S.lower_bound({ { -x,-e9 },-e9 });
			if (v == S.end())continue;
			int who = v->second;
			if(S_start[who] <= x && x <= S_end[who]) qr[who].push_back({ x, c });
		}
	}
	for (i = 1; i <= tn; i++)for (j = 0; j < 20; j++)pp[i][j] = -1;
	init_sparse(1, -1);
	dfs(1, -1, 0);
	printf("%lld\n", allS - D[1]);
	return 0;
}

Compilation message

constellation3.cpp: In function 'int main()':
constellation3.cpp:100:36: warning: unused variable 'y' [-Wunused-variable]
  100 |    int t = all[i].t, x = all[i].x, y = all[i].y, c = all[i].c;
      |                                    ^
constellation3.cpp:100:50: warning: unused variable 'c' [-Wunused-variable]
  100 |    int t = all[i].t, x = all[i].x, y = all[i].y, c = all[i].c;
      |                                                  ^
constellation3.cpp:110:18: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<int>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
  110 |    for (i = 0; i < S_ER[who].size(); i++) {
      |                ~~^~~~~~~~~~~~~~~~~~
constellation3.cpp:116:15: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<int>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
  116 |     if (i + 1 == S_ER[who].size()) {
      |         ~~~~~~^~~~~~~~~~~~~~~~~~~
constellation3.cpp:125:36: warning: unused variable 'y' [-Wunused-variable]
  125 |    int t = all[i].t, x = all[i].x, y = all[i].y, c = all[i].c;
      |                                    ^
constellation3.cpp:78:14: warning: ignoring return value of 'int scanf(const char*, ...)' declared with attribute 'warn_unused_result' [-Wunused-result]
   78 |  int n; scanf("%d", &n);
      |         ~~~~~^~~~~~~~~~
constellation3.cpp:81:8: warning: ignoring return value of 'int scanf(const char*, ...)' declared with attribute 'warn_unused_result' [-Wunused-result]
   81 |   scanf("%d", &a[i]);
      |   ~~~~~^~~~~~~~~~~~~
constellation3.cpp:84:14: warning: ignoring return value of 'int scanf(const char*, ...)' declared with attribute 'warn_unused_result' [-Wunused-result]
   84 |  int m; scanf("%d", &m);
      |         ~~~~~^~~~~~~~~~
constellation3.cpp:88:8: warning: ignoring return value of 'int scanf(const char*, ...)' declared with attribute 'warn_unused_result' [-Wunused-result]
   88 |   scanf("%d%d%d", &all[an].x, &all[an].y, &all[an].c);
      |   ~~~~~^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# Verdict Execution time Memory Grader output
1 Correct 160 ms 224096 KB Output is correct
2 Correct 156 ms 224032 KB Output is correct
3 Correct 153 ms 224188 KB Output is correct
4 Correct 157 ms 224132 KB Output is correct
5 Correct 158 ms 224112 KB Output is correct
6 Correct 152 ms 224064 KB Output is correct
7 Correct 154 ms 224052 KB Output is correct
8 Correct 151 ms 224064 KB Output is correct
9 Correct 157 ms 224036 KB Output is correct
10 Correct 153 ms 224080 KB Output is correct
11 Correct 154 ms 224328 KB Output is correct
12 Correct 159 ms 224200 KB Output is correct
13 Correct 155 ms 224068 KB Output is correct
14 Correct 158 ms 224068 KB Output is correct
15 Correct 154 ms 224076 KB Output is correct
16 Correct 154 ms 224144 KB Output is correct
17 Correct 158 ms 224096 KB Output is correct
18 Correct 152 ms 224144 KB Output is correct
19 Correct 158 ms 224068 KB Output is correct
20 Correct 156 ms 224068 KB Output is correct
21 Correct 156 ms 224064 KB Output is correct
22 Correct 157 ms 224100 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 160 ms 224096 KB Output is correct
2 Correct 156 ms 224032 KB Output is correct
3 Correct 153 ms 224188 KB Output is correct
4 Correct 157 ms 224132 KB Output is correct
5 Correct 158 ms 224112 KB Output is correct
6 Correct 152 ms 224064 KB Output is correct
7 Correct 154 ms 224052 KB Output is correct
8 Correct 151 ms 224064 KB Output is correct
9 Correct 157 ms 224036 KB Output is correct
10 Correct 153 ms 224080 KB Output is correct
11 Correct 154 ms 224328 KB Output is correct
12 Correct 159 ms 224200 KB Output is correct
13 Correct 155 ms 224068 KB Output is correct
14 Correct 158 ms 224068 KB Output is correct
15 Correct 154 ms 224076 KB Output is correct
16 Correct 154 ms 224144 KB Output is correct
17 Correct 158 ms 224096 KB Output is correct
18 Correct 152 ms 224144 KB Output is correct
19 Correct 158 ms 224068 KB Output is correct
20 Correct 156 ms 224068 KB Output is correct
21 Correct 156 ms 224064 KB Output is correct
22 Correct 157 ms 224100 KB Output is correct
23 Correct 157 ms 224848 KB Output is correct
24 Correct 159 ms 224732 KB Output is correct
25 Correct 156 ms 224752 KB Output is correct
26 Correct 157 ms 224776 KB Output is correct
27 Correct 163 ms 224864 KB Output is correct
28 Correct 160 ms 224888 KB Output is correct
29 Correct 158 ms 224788 KB Output is correct
30 Correct 162 ms 224780 KB Output is correct
31 Correct 156 ms 224900 KB Output is correct
32 Correct 156 ms 225268 KB Output is correct
33 Correct 156 ms 225176 KB Output is correct
34 Correct 159 ms 225280 KB Output is correct
35 Correct 155 ms 224984 KB Output is correct
36 Correct 156 ms 224116 KB Output is correct
37 Correct 155 ms 224144 KB Output is correct
38 Correct 160 ms 225604 KB Output is correct
39 Correct 160 ms 224656 KB Output is correct
40 Correct 182 ms 225092 KB Output is correct
41 Correct 160 ms 224476 KB Output is correct
42 Correct 160 ms 224524 KB Output is correct
43 Correct 179 ms 225088 KB Output is correct
44 Correct 159 ms 224452 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 160 ms 224096 KB Output is correct
2 Correct 156 ms 224032 KB Output is correct
3 Correct 153 ms 224188 KB Output is correct
4 Correct 157 ms 224132 KB Output is correct
5 Correct 158 ms 224112 KB Output is correct
6 Correct 152 ms 224064 KB Output is correct
7 Correct 154 ms 224052 KB Output is correct
8 Correct 151 ms 224064 KB Output is correct
9 Correct 157 ms 224036 KB Output is correct
10 Correct 153 ms 224080 KB Output is correct
11 Correct 154 ms 224328 KB Output is correct
12 Correct 159 ms 224200 KB Output is correct
13 Correct 155 ms 224068 KB Output is correct
14 Correct 158 ms 224068 KB Output is correct
15 Correct 154 ms 224076 KB Output is correct
16 Correct 154 ms 224144 KB Output is correct
17 Correct 158 ms 224096 KB Output is correct
18 Correct 152 ms 224144 KB Output is correct
19 Correct 158 ms 224068 KB Output is correct
20 Correct 156 ms 224068 KB Output is correct
21 Correct 156 ms 224064 KB Output is correct
22 Correct 157 ms 224100 KB Output is correct
23 Correct 157 ms 224848 KB Output is correct
24 Correct 159 ms 224732 KB Output is correct
25 Correct 156 ms 224752 KB Output is correct
26 Correct 157 ms 224776 KB Output is correct
27 Correct 163 ms 224864 KB Output is correct
28 Correct 160 ms 224888 KB Output is correct
29 Correct 158 ms 224788 KB Output is correct
30 Correct 162 ms 224780 KB Output is correct
31 Correct 156 ms 224900 KB Output is correct
32 Correct 156 ms 225268 KB Output is correct
33 Correct 156 ms 225176 KB Output is correct
34 Correct 159 ms 225280 KB Output is correct
35 Correct 155 ms 224984 KB Output is correct
36 Correct 156 ms 224116 KB Output is correct
37 Correct 155 ms 224144 KB Output is correct
38 Correct 160 ms 225604 KB Output is correct
39 Correct 160 ms 224656 KB Output is correct
40 Correct 182 ms 225092 KB Output is correct
41 Correct 160 ms 224476 KB Output is correct
42 Correct 160 ms 224524 KB Output is correct
43 Correct 179 ms 225088 KB Output is correct
44 Correct 159 ms 224452 KB Output is correct
45 Correct 1051 ms 310828 KB Output is correct
46 Correct 1018 ms 309536 KB Output is correct
47 Correct 1006 ms 308312 KB Output is correct
48 Correct 1012 ms 310984 KB Output is correct
49 Correct 995 ms 307732 KB Output is correct
50 Correct 1004 ms 308048 KB Output is correct
51 Correct 1032 ms 308124 KB Output is correct
52 Correct 1040 ms 310020 KB Output is correct
53 Correct 1067 ms 310024 KB Output is correct
54 Correct 1153 ms 382408 KB Output is correct
55 Correct 1077 ms 357928 KB Output is correct
56 Correct 980 ms 345304 KB Output is correct
57 Correct 919 ms 336996 KB Output is correct
58 Correct 329 ms 236356 KB Output is correct
59 Correct 332 ms 236152 KB Output is correct
60 Correct 1192 ms 431088 KB Output is correct
61 Correct 529 ms 285984 KB Output is correct
62 Correct 1033 ms 366660 KB Output is correct
63 Correct 464 ms 272816 KB Output is correct
64 Correct 501 ms 283716 KB Output is correct
65 Correct 1048 ms 368632 KB Output is correct
66 Correct 458 ms 272324 KB Output is correct