Submission #726241

#	Time	Username	Problem	Language	Result	Execution time	Memory
726241		alvingogo	Rectangles (IOI19_rect)	C++14	72 / 100	5095 ms	294388 KiB

rect

#include "rect.h"
#include <bits/stdc++.h>
#pragma GCC optimize("Ofast")
#define fs first
#define sc second
using namespace std;

typedef pair<int,int> pii;

namespace{

struct BIT{
	int n;
	vector<int> bt;
	BIT(int x){
		n=x;
		bt.resize(x+1);
	}
	BIT(){
		n=0;
	}
	void update(int x,int y){
		for(;x<=n;x+=(x&-x)){
			bt[x]+=y;
		}
	}
	int query(int x){
		int ans=0;
		for(;x>0;x-=(x&-x)){
			ans+=bt[x];
		}
		return ans;
	}
};

}
long long count_rectangles(vector<vector<int> > v) {
	int n=v.size();
	int m=v[0].size();
	if(n<=2 || m<=2){
		return 0;
	}
	vector<vector<pii> > az(n),bz(m);
	for(int i=1;i<n-1;i++){
		vector<pair<int,int> > gg;
		vector<int> vis(m),up(m),down(m);
		for(int j=0;j<m;j++){
			gg.push_back({v[i][j],j});
		}
		sort(gg.begin(),gg.end());
		for(auto h:gg){
			int x=h.sc-1,y=h.sc+1;
			int l=h.sc,r=h.sc;
			if(x<0 || !vis[x]){

			}
			else{
				l=up[x];
			}
			if(y>=m || !vis[y]){
				
			}
			else{
				r=down[y];
			}
			vis[h.sc]=1;
			up[h.sc]=l;
			down[h.sc]=r;
			down[l]=r;
			up[r]=l;
			if(l!=0 && r!=m-1 && (v[i][r+1]>h.fs && v[i][l-1]>h.fs)){
				az[i].push_back({l,r});
			}
		}
	}
	for(int i=1;i<m-1;i++){
		vector<pair<int,int> > gg;
		vector<int> vis(n),up(n),down(n);
		for(int j=0;j<n;j++){
			gg.push_back({v[j][i],j});
		}
		sort(gg.begin(),gg.end());
		for(auto h:gg){
			int x=h.sc-1,y=h.sc+1;
			int l=h.sc,r=h.sc;
			if(x<0 || !vis[x]){

			}
			else{
				l=up[x];
			}
			if(y>=n || !vis[y]){
				
			}
			else{
				r=down[y];
			}
			vis[h.sc]=1;
			up[h.sc]=l;
			down[h.sc]=r;
			down[l]=r;
			up[r]=l;
			if(l!=0 && r!=n-1 && (v[r+1][i]>h.fs && v[l-1][i]>h.fs)){
				bz[i].push_back({l,r});
			}
		}
	}
	vector<pair<int,pii> > ev[2505];
	map<pii,int> mp;
	for(int i=1;i<n-1;i++){
		map<pii,int> tmp;
		for(auto h:az[i]){
			tmp[h]=1;
			if(mp.find(h)==mp.end()){
				mp[h]=i;
			}
		}
		if(i>1){
			for(auto h:az[i-1]){
				if(tmp.find(h)==tmp.end()){
					ev[h.fs].push_back({-h.sc,{mp[h],i-1}});
					mp.erase(h);
				}
			}
		}
	}
	for(auto h:mp){
		ev[h.fs.fs].push_back({-h.fs.sc,{h.sc,n-2}});
	}
	mp.clear();
	for(int i=1;i<m-1;i++){
		map<pii,int> tmp;
		for(auto h:bz[i]){
			tmp[h]=1;
			if(mp.find(h)==mp.end()){
				mp[h]=i;
			}
		}
		if(i>1){
			for(auto h:bz[i-1]){
				if(tmp.find(h)==tmp.end()){
					for(int j=mp[h];j<=i-1;j++){
						ev[j].push_back({-(i-1),{-h.fs,-h.sc}});
					}
					mp.erase(h);
				}
			}
		}
	}
	for(auto h:mp){
		for(int j=h.sc;j<=m-2;j++){
			ev[j].push_back({-(m-2),{-h.fs.fs,-h.fs.sc}});
		}
	}
	
	BIT vb[2502];
	for(int i=0;i<n;i++){
		vb[i]=BIT(n);
	}
	//vector<BIT> vb(n,BIT(n));
	long long ans=0;
	for(int i=1;i<m-1;i++){
		sort(ev[i].begin(),ev[i].end());
	}
	for(int i=1;i<m-1;i++){
		//sort(ev[i].begin(),ev[i].end());
		for(auto h:ev[i]){
			if(h.sc.sc>0){
				for(int i=h.sc.fs;i<=h.sc.sc;i++){
					ans+=vb[i].query(h.sc.sc);
				}	
			}
			else{
				vb[-h.sc.fs].update(-h.sc.sc,1);
			}
		}
		for(auto h:ev[i]){
			if(h.sc.sc<0){
				vb[-h.sc.fs].update(-h.sc.sc,-1);
			}
		}
	}
	return ans;
}

• Subtask 1: 8 / 8
• Subtask 2: 7 / 7
• Subtask 3: 12 / 12
• Subtask 4: 22 / 22
• Subtask 5: 10 / 10
• Subtask 6: 13 / 13
• Subtask 7: 0 / 28