Submission #340514

#	Time	Username	Problem	Language	Result	Execution time	Memory
340514		AmineWeslati	Pairs (IOI07_pairs)	C++14	70 / 100	932 ms	524288 KiB

pairs

//Never stop trying
/*#pragma GCC target ("avx2")
#pragma GCC optimize ("Ofast")
#pragma GCC optimization ("O3")
#pragma GCC optimization ("unroll-loops")*/
#include "bits/stdc++.h"
using namespace std;
#define boost ios_base::sync_with_stdio(false); cin.tie(0); cout.tie(0)

typedef long long ll;
typedef string str;
typedef double db;
typedef long double ld;
typedef pair<int, int> pi;
#define fi first
#define se second
typedef vector<int> vi;
typedef vector<pi> vpi;
typedef vector<str> vs;
typedef vector<ld> vd;
#define pb push_back
#define eb emplace_back
#define sz(x) (int)x.size()
#define all(x) begin(x), end(x)
#define rall(x) rbegin(x), rend(x)
#define endl "\n"

#define FOR(i,a,b) for (int i = (a); i < (b); ++i)
#define ROF(i,a,b) for (int i = (b)-1; i >= (a); --i)

const int MOD = 1e9 + 7; //998244353
const ll INF = 1e18;
const int MX = 2e5 + 10;
const int nx[4] = {0, 0, 1, -1}, ny[4] = {1, -1, 0, 0}; //right left down up

template<class T> using V = vector<T>;
template<class T> bool ckmin(T& a, const T& b) { return a > b ? a = b, 1 : 0; }
template<class T> bool ckmax(T& a, const T& b) { return a < b ? a = b, 1 : 0; }

ll cdiv(ll a, ll b) { return a / b + ((a ^ b) > 0 && a % b); } // divide a by b rounded up
//constexpr int log2(int x) { return 31 - __builtin_clz(x); } // floor(log2(x))

mt19937 rng(chrono::system_clock::now().time_since_epoch().count());
//mt19937_64 rng(chrono::system_clock::now().time_since_epoch().count());

ll random(ll a, ll b){
    return a + rng() % (b - a + 1);
}

#ifndef LOCAL  
#define cerr if(false) cerr
#endif
#define dbg(x) cerr << #x << " : " << x << endl; 
#define dbgs(x,y) cerr << #x << " : " << x << " / " << #y << " : " << y << endl;
#define dbgv(v) cerr << #v << " : " << "[ "; for(auto it : v) cerr << it << ' '; cerr << ']' << endl;
#define here() cerr << "here" << endl;

void IO() {
#ifdef LOCAL
    freopen("input.txt", "r", stdin);
    freopen("output.txt", "w", stdout);
#endif
}

int B,N,D,M; 

vi bit;
void upd(int i, int v){
	for(;i<MX;i+=i&-i) bit[i]+=v;
}
int sum(int i){
	int ans=0;
	for(;i;i-=i&-i) ans+=bit[i];
	return ans;
}
int get(int l, int r){return sum(r)-sum(l-1); }
	
int cnt;
V<V<vi>> ft;
void update(int ii, int jj, int kk, int v){
	for(int i=ii; i<cnt; i+=i&-i) 
		for(int j=jj ;j<cnt; j+=j&-j) 
			for(int k=kk ;k<cnt; k+=k&-k){
				ft[i][j][k]+=v;
			}
}
int sum2(int ii, int jj, int kk){	
	int ans=0;
	for(int i=ii;i;i-=i&-i)
		for(int j=jj;j;j-=j&-j)
			for(int k=kk;k;k-=k&-k){
				ans+=ft[i][j][k];
			}
	return ans;
}
int get2(int x, int xx, int y, int yy, int z, int zz){
	int ans=0;
	ans+=sum2(xx,yy,zz);
	ans-=sum2(x-1,yy,zz);
	ans-=sum2(xx,y-1,zz);
	ans+=sum2(x-1,y-1,zz);

	ans-=sum2(xx,yy,z-1);
	ans+=sum2(x-1,yy,z-1);
	ans+=sum2(xx,y-1,z-1);
	ans-=sum2(x-1,y-1,z-1);

	return ans;
}

int main() {
    boost; IO();

    cin>>B>>N>>D>>M;
    if(B==1){
    	vi x(N);
    	FOR(i,0,N){
    		cin>>x[i];
    	}
    	sort(all(x));
    	int j=0;
    	ll ans=0;
    	FOR(i,0,N){
    		while(x[i]-x[j]>D) j++;
    		ans+=i-j;
    	}
    	cout << ans << endl;
    }
    else if(B==2){
        bit.assign(MX,0);
    	vpi vec;
    	FOR(i,0,N){
    		int x,y; cin>>x>>y;
    		vec.pb({x-y,x+y});
    	}
    	sort(all(vec));
    		
    	int j=0;
    	ll ans=0;
    	FOR(i,0,N){
    		while(vec[i].fi-vec[j].fi>D){
    			upd(vec[j].se,-1);
    			j++;
			}
			ans+=get(max(vec[i].se-D,1),vec[i].se+D);
    		upd(vec[i].se,1);
    	}
    	cout << ans << endl;
    }
    else{
    	V<pair<pi,pi>> vec;
        vi v;
    	FOR(i,0,N){
    		int x,y,z; cin>>x>>y>>z;
    		vec.pb({{x+y+z,x+y-z},{x-y+z,x-y-z}});
    		v.pb(x+y-z); v.pb(x+y-z-D); v.pb(x+y-z+D);
    		v.pb(x-y+z); v.pb(x-y+z-D); v.pb(x-y+z+D);
    		v.pb(x-y-z); v.pb(x-y-z-D); v.pb(x-y-z+D);
    	} 
        sort(all(vec));
    	sort(all(v));
    	cnt=0;
    	map<int,int> mp;
    	for(auto x: v) if(!mp.count(x)) mp[x]=++cnt;

        assert(cnt<=600);
        cnt++;
        ft.resize(cnt);
        FOR(i,0,cnt){ 
            ft[i].resize(cnt);
            FOR(j,0,cnt) ft[i][j].assign(cnt,0);
        }
    	
    	int j=0;
    	ll ans=0;
    	FOR(i,0,N){
    		pair<pi,pi> &x=vec[i];
    		int a=x.fi.se,b=x.se.fi,c=x.se.se;
    		x.fi.se=mp[x.fi.se];
    		x.se.fi=mp[x.se.fi];
    		x.se.se=mp[x.se.se];

    		
    		while(x.fi.fi-vec[j].fi.fi>D){
    			update(vec[j].fi.se,vec[j].se.fi,vec[j].se.se,-1);
				j++;
	    	}
    		ans+=get2(mp[a-D],mp[a+D], mp[b-D],mp[b+D], mp[c-D],mp[c+D]);
    		update(x.fi.se,x.se.fi,x.se.se,1);
    	}
    	cout << ans << endl;
    }

    return 0;
}

/* Careful!!!
    .Array bounds
    .Infinite loops
    .Uninitialized variables / empty containers
    .Multisets are shit

   Some insights:
    .Binary search
    .Graph representation
    .Write brute force code
    .Change your approach
*/

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