oj.uz

Submission #618982

# Submission time Handle Problem Language Result Execution time Memory
618982 2022-08-02T08:45:23 Z CSQ31 Team Contest (JOI22_team) C++17
37 / 100
 2000 ms 27352 KB 
#include <bits/stdc++.h>
using namespace std;
#define pb push_back
#define fi first
#define se second
#define sz(a) (int)(a.size())
#define all(a) a.begin(),a.end()
#define owo ios_base::sync_with_stdio(0);cin.tie(0);
typedef pair<int,int> pii;
typedef long long int ll;
const int MAXN = 5e5;
array<ll,3>pts1[MAXN],pts2[MAXN];
vector<array<ll,3>>a;
vector<int>crd[3];
//i have rectangle queries but i cant do it offline?
//i can but its double log on 2s
//ok try n^2 first
int sz0,sz1,sz2;
struct seg{
	vector<int>t;
	int n;
	seg(){}
	seg(int _n){
		n = _n;
		t.assign(4*n+1,-2e9);
	}
	void upd(int v,int pos,int l,int r,int val){
		if(l==r){
			t[v] = max(t[v],val);
			return;
		}
		int tm = (l+r)/2;
		if(pos<=tm)upd(2*v,pos,l,tm,val);
		else upd(2*v+1,pos,tm+1,r,val);
		t[v] = max(t[2*v],t[2*v+1]);
	}
	void upd2(int v,int pos,int l,int r,int val){
		if(l==r){
			t[v] = val;
			return;
		}
		int tm = (l+r)/2;
		if(pos<=tm)upd(2*v,pos,l,tm,val);
		else upd(2*v+1,pos,tm+1,r,val);
		t[v] = max(t[2*v],t[2*v+1]);
	}
	int query(int v,int l,int r,int tl,int tr){
		if(l>r)return -2e9;
		if(l==tl && r==tr)return t[v];
		int tm = (tl+tr)/2;
		if(r<=tm)return query(2*v,l,r,tl,tm);
		else if(l>tm)return query(2*v+1,l,r,tm+1,tr);
		else return max(query(2*v,l,tm,tl,tm),query(2*v+1,tm+1,r,tm+1,tr));
	}
}amx;
ll ans = -1;
vector<int>pts[MAXN];
 
void solve(int l,int r){
	if(l==r){
		return;
	}
	//all same a
	int mid = (l+r)/2;
	solve(l,mid);
	solve(mid+1,r);
	vector<pii>que;
	vector<array<ll,3>>add;
	for(int i=l;i<=mid;i++){
		for(int x:pts[i]){
			if(pts1[x][0] != -2e9){
				add.pb({pts1[x][1],pts1[x][2],pts1[x][0]});
			}
			if(pts2[x][0] != -2e9){
				add.pb({pts2[x][1],pts2[x][2],pts2[x][0]});
			}
		}
	}
	for(int i=mid+1;i<=r;i++){
		for(int x:pts[i])que.push_back({a[x][1],x});
	}
	for(auto x:add){
		for(auto y:que){
			int id = y.se;
			if(y.fi  < x[0] && a[id][2] < x[1]){
				ans = max(ans,crd[0][a[id][0]] + x[2]);
				
			}
			
		}	
	}
	/*
	sort(all(add),greater<array<ll,3>>());
	sort(all(que),greater<pii>());
	int ptr = -1;
	for(pii qq : que){
		while(ptr+1 < sz(add) && add[ptr+1][0] > qq.fi){ //add still has bigger b
			ptr++;
			array<ll,3>v = add[ptr];
			amx.upd(1,v[1],0,sz2-1,v[2]);
		} 
		int id = qq.se;
		if(a[id][2]+1 < sz2){
			ll res = amx.query(1,a[id][2]+1,sz2-1,0,sz2-1); //we need a bigger ccout<<"query"<<'\n';
			ans = max(ans , crd[0][a[id][0]] + res);
		}
	}
	for(auto x:add)amx.upd(1,x[1],0,sz2-1,-2e9);*/
}
int main()
{
	owo
	int n;
	cin>>n;
	a.resize(n);
	for(int i=0;i<n;i++){
		for(int j=0;j<3;j++){
			cin>>a[i][j];
			crd[j].pb(a[i][j]);
		}
	}
	sort(all(a),[&](array<ll,3>i,array<ll,3>j){return j[0] > i[0];});
	for(int j=0;j<3;j++){
		sort(all(crd[j]));
		crd[j].resize(unique(all(crd[j])) - crd[j].begin());
	}
	//cout<<"sorted"<<'\n';
	//for(int i=0;i<n;i++)cout<<a[i][0]<<" "<<a[i][1]<<" "<<a[i][2]<<'\n';
	sz0 = sz(crd[0]);
	sz1 = sz(crd[1]);
	sz2 = sz(crd[2]);
	seg bmx(sz1),cmx(sz2);
	//cout<<"points"<<'\n';
	for(int i=0;i<n;i++){
		pts1[i][0] = pts2[i][0] = -2e9;
		pts1[i][1] = a[i][1]; //i has bigger b
		pts2[i][2] = a[i][2]; //i has bigger c
		int s1 = lower_bound(all(crd[1]),a[i][1]) - crd[1].begin();
		int s2 = lower_bound(all(crd[2]),a[i][2]) - crd[2].begin();
		if(s1>0)pts1[i][2] = bmx.query(1,0,s1-1,0,sz1-1); //want the max of c among smaller bs
		if(s2>0)pts2[i][1] = cmx.query(1,0,s2-1,0,sz2-1); //want the max of b among smaller cs
		if(pts1[i][2] > a[i][2])pts1[i][0] = pts1[i][1] + pts1[i][2]; 
		if(pts2[i][1] > a[i][1])pts2[i][0] = pts2[i][1] + pts2[i][2];
		bmx.upd(1,s1,0,sz1-1,a[i][2]);
		cmx.upd(1,s2,0,sz2-1,a[i][1]);
		
		if(pts1[i][0] != -2e9){
			pts1[i][1] = lower_bound(all(crd[1]),pts1[i][1]) - crd[1].begin();
			pts1[i][2] = lower_bound(all(crd[2]),pts1[i][2]) - crd[2].begin();
		}
		if(pts2[i][0] != -2e9){
			pts2[i][1] = lower_bound(all(crd[1]),pts2[i][1]) - crd[1].begin();
			pts2[i][2] = lower_bound(all(crd[2]),pts2[i][2]) - crd[2].begin();
		}
		//this is correct
	}
	for(int i=0;i<n;i++){
		a[i][0] = lower_bound(all(crd[0]),a[i][0]) - crd[0].begin();
		a[i][1] = lower_bound(all(crd[1]),a[i][1]) - crd[1].begin();
		a[i][2] = lower_bound(all(crd[2]),a[i][2]) - crd[2].begin();
		pts[a[i][0]].pb(i);
	}
	/*
	for(int i=0;i<n;){
		int k = i;
		while(k+1<n && a[k][0] == a[k+1][0])k++;
		for(int x=i;x<=k;x++){
			for(int j=0;j<x;j++){
				if(a[x][0] == a[j][0])break;
				if(pts1[j][1] > a[x][1] && pts1[j][2] >a[x][2])ans = max(ans,crd[0][a[x][0]] + pts1[j][0]);
				if(pts2[j][1] > a[x][1] && pts2[j][2] >a[x][2])ans = max(ans,crd[0][a[x][0]] + pts2[j][0]);
			}
		}
		i = k+1;
	}*/
	amx = seg(sz2);
	solve(0,sz0-1);
	cout<<ans;
	
	
	
}

Subtask #1
8.0 / 8.0

# Verdict Execution time Memory Grader output
1 Correct 7 ms 11988 KB Output is correct
2 Correct 6 ms 11988 KB Output is correct
3 Correct 9 ms 12016 KB Output is correct
4 Correct 6 ms 12024 KB Output is correct
5 Correct 6 ms 11988 KB Output is correct
6 Correct 6 ms 11988 KB Output is correct
7 Correct 6 ms 12076 KB Output is correct
8 Correct 6 ms 11988 KB Output is correct
9 Correct 6 ms 11988 KB Output is correct
10 Correct 6 ms 11988 KB Output is correct
11 Correct 6 ms 11988 KB Output is correct
12 Correct 7 ms 11988 KB Output is correct
13 Correct 6 ms 11988 KB Output is correct
14 Correct 7 ms 12116 KB Output is correct
15 Correct 8 ms 12128 KB Output is correct
16 Correct 7 ms 12116 KB Output is correct
17 Correct 8 ms 12100 KB Output is correct
18 Correct 7 ms 12116 KB Output is correct
19 Correct 7 ms 12140 KB Output is correct
20 Correct 8 ms 12124 KB Output is correct
21 Correct 9 ms 12208 KB Output is correct
22 Correct 8 ms 12116 KB Output is correct
23 Correct 7 ms 12112 KB Output is correct
24 Correct 8 ms 12116 KB Output is correct
25 Correct 8 ms 12116 KB Output is correct

Subtask #2
29.0 / 29.0

# Verdict Execution time Memory Grader output
1 Correct 7 ms 11988 KB Output is correct
2 Correct 6 ms 11988 KB Output is correct
3 Correct 9 ms 12016 KB Output is correct
4 Correct 6 ms 12024 KB Output is correct
5 Correct 6 ms 11988 KB Output is correct
6 Correct 6 ms 11988 KB Output is correct
7 Correct 6 ms 12076 KB Output is correct
8 Correct 6 ms 11988 KB Output is correct
9 Correct 6 ms 11988 KB Output is correct
10 Correct 6 ms 11988 KB Output is correct
11 Correct 6 ms 11988 KB Output is correct
12 Correct 7 ms 11988 KB Output is correct
13 Correct 6 ms 11988 KB Output is correct
14 Correct 7 ms 12116 KB Output is correct
15 Correct 8 ms 12128 KB Output is correct
16 Correct 7 ms 12116 KB Output is correct
17 Correct 8 ms 12100 KB Output is correct
18 Correct 7 ms 12116 KB Output is correct
19 Correct 7 ms 12140 KB Output is correct
20 Correct 8 ms 12124 KB Output is correct
21 Correct 9 ms 12208 KB Output is correct
22 Correct 8 ms 12116 KB Output is correct
23 Correct 7 ms 12112 KB Output is correct
24 Correct 8 ms 12116 KB Output is correct
25 Correct 8 ms 12116 KB Output is correct
26 Correct 68 ms 12908 KB Output is correct
27 Correct 56 ms 12756 KB Output is correct
28 Correct 48 ms 12756 KB Output is correct
29 Correct 18 ms 12588 KB Output is correct
30 Correct 20 ms 12488 KB Output is correct
31 Correct 72 ms 12732 KB Output is correct
32 Correct 22 ms 12500 KB Output is correct
33 Correct 11 ms 12500 KB Output is correct
34 Correct 15 ms 12644 KB Output is correct
35 Correct 8 ms 12188 KB Output is correct
36 Correct 8 ms 12096 KB Output is correct
37 Correct 10 ms 12500 KB Output is correct
38 Correct 11 ms 12500 KB Output is correct
39 Correct 21 ms 12244 KB Output is correct

Subtask #3
0 / 9.0

# Verdict Execution time Memory Grader output
1 Correct 7 ms 12032 KB Output is correct
2 Correct 6 ms 11988 KB Output is correct
3 Correct 6 ms 11988 KB Output is correct
4 Correct 7 ms 11980 KB Output is correct
5 Correct 8 ms 11988 KB Output is correct
6 Correct 6 ms 11988 KB Output is correct
7 Correct 7 ms 11988 KB Output is correct
8 Correct 6 ms 11988 KB Output is correct
9 Correct 6 ms 11988 KB Output is correct
10 Correct 6 ms 11988 KB Output is correct
11 Execution timed out 2091 ms 27352 KB Time limit exceeded
12 Halted 0 ms 0 KB -

Subtask #4
0 / 9.0

# Verdict Execution time Memory Grader output
1 Correct 7 ms 12032 KB Output is correct
2 Correct 6 ms 11988 KB Output is correct
3 Correct 6 ms 11988 KB Output is correct
4 Correct 7 ms 11980 KB Output is correct
5 Correct 8 ms 11988 KB Output is correct
6 Correct 6 ms 11988 KB Output is correct
7 Correct 7 ms 11988 KB Output is correct
8 Correct 6 ms 11988 KB Output is correct
9 Correct 6 ms 11988 KB Output is correct
10 Correct 6 ms 11988 KB Output is correct
11 Execution timed out 2091 ms 27352 KB Time limit exceeded
12 Halted 0 ms 0 KB -

Subtask #5
0 / 9.0

# Verdict Execution time Memory Grader output
1 Correct 7 ms 12032 KB Output is correct
2 Correct 6 ms 11988 KB Output is correct
3 Correct 6 ms 11988 KB Output is correct
4 Correct 7 ms 11980 KB Output is correct
5 Correct 8 ms 11988 KB Output is correct
6 Correct 6 ms 11988 KB Output is correct
7 Correct 7 ms 11988 KB Output is correct
8 Correct 6 ms 11988 KB Output is correct
9 Correct 6 ms 11988 KB Output is correct
10 Correct 6 ms 11988 KB Output is correct
11 Execution timed out 2091 ms 27352 KB Time limit exceeded
12 Halted 0 ms 0 KB -

Subtask #6
0 / 9.0

# Verdict Execution time Memory Grader output
1 Correct 7 ms 12032 KB Output is correct
2 Correct 6 ms 11988 KB Output is correct
3 Correct 6 ms 11988 KB Output is correct
4 Correct 7 ms 11980 KB Output is correct
5 Correct 8 ms 11988 KB Output is correct
6 Correct 6 ms 11988 KB Output is correct
7 Correct 7 ms 11988 KB Output is correct
8 Correct 6 ms 11988 KB Output is correct
9 Correct 6 ms 11988 KB Output is correct
10 Correct 6 ms 11988 KB Output is correct
11 Execution timed out 2091 ms 27352 KB Time limit exceeded
12 Halted 0 ms 0 KB -

Subtask #7
0 / 27.0

# Verdict Execution time Memory Grader output
1 Correct 7 ms 11988 KB Output is correct
2 Correct 6 ms 11988 KB Output is correct
3 Correct 9 ms 12016 KB Output is correct
4 Correct 6 ms 12024 KB Output is correct
5 Correct 6 ms 11988 KB Output is correct
6 Correct 6 ms 11988 KB Output is correct
7 Correct 6 ms 12076 KB Output is correct
8 Correct 6 ms 11988 KB Output is correct
9 Correct 6 ms 11988 KB Output is correct
10 Correct 6 ms 11988 KB Output is correct
11 Correct 6 ms 11988 KB Output is correct
12 Correct 7 ms 11988 KB Output is correct
13 Correct 6 ms 11988 KB Output is correct
14 Correct 7 ms 12116 KB Output is correct
15 Correct 8 ms 12128 KB Output is correct
16 Correct 7 ms 12116 KB Output is correct
17 Correct 8 ms 12100 KB Output is correct
18 Correct 7 ms 12116 KB Output is correct
19 Correct 7 ms 12140 KB Output is correct
20 Correct 8 ms 12124 KB Output is correct
21 Correct 9 ms 12208 KB Output is correct
22 Correct 8 ms 12116 KB Output is correct
23 Correct 7 ms 12112 KB Output is correct
24 Correct 8 ms 12116 KB Output is correct
25 Correct 8 ms 12116 KB Output is correct
26 Correct 68 ms 12908 KB Output is correct
27 Correct 56 ms 12756 KB Output is correct
28 Correct 48 ms 12756 KB Output is correct
29 Correct 18 ms 12588 KB Output is correct
30 Correct 20 ms 12488 KB Output is correct
31 Correct 72 ms 12732 KB Output is correct
32 Correct 22 ms 12500 KB Output is correct
33 Correct 11 ms 12500 KB Output is correct
34 Correct 15 ms 12644 KB Output is correct
35 Correct 8 ms 12188 KB Output is correct
36 Correct 8 ms 12096 KB Output is correct
37 Correct 10 ms 12500 KB Output is correct
38 Correct 11 ms 12500 KB Output is correct
39 Correct 21 ms 12244 KB Output is correct
40 Correct 7 ms 12032 KB Output is correct
41 Correct 6 ms 11988 KB Output is correct
42 Correct 6 ms 11988 KB Output is correct
43 Correct 7 ms 11980 KB Output is correct
44 Correct 8 ms 11988 KB Output is correct
45 Correct 6 ms 11988 KB Output is correct
46 Correct 7 ms 11988 KB Output is correct
47 Correct 6 ms 11988 KB Output is correct
48 Correct 6 ms 11988 KB Output is correct
49 Correct 6 ms 11988 KB Output is correct
50 Execution timed out 2091 ms 27352 KB Time limit exceeded
51 Halted 0 ms 0 KB -