oj.uz

Submission #340505

# Submission time Handle Problem Language Result Execution time Memory
340505 2020-12-27T19:22:10 Z AmineWeslati Pairs (IOI07_pairs) C++14
70 / 100
 1099 ms 524288 KB 
//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(MX,0);
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); }
	
const int MX2=500;
int ft[MX2][MX2][MX2];
void update(int ii, int jj, int kk, int v){
	for(int i=ii; i<MX2; i+=i&-i) 
		for(int j=jj ;j<MX2; j+=j&-j) 
			for(int k=kk ;k<MX2; 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){
    	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{
    	memset(ft,0,sizeof(ft));
    	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(v));
    	int cnt=0;
    	unordered_map<int,int> mp;
    	for(auto x: v) if(!mp.count(x)) mp[x]=++cnt;
        assert(cnt<MX2);
    	sort(all(vec));
    	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.0 / 6.0

# Verdict Execution time Memory Grader output
1 Correct 1 ms 1132 KB Output is correct
2 Correct 1 ms 1132 KB Output is correct

Subtask #2
6.0 / 6.0

# Verdict Execution time Memory Grader output
1 Correct 1 ms 1132 KB Output is correct

Subtask #3
6.0 / 6.0

# Verdict Execution time Memory Grader output
1 Correct 15 ms 1516 KB Output is correct
2 Correct 14 ms 1516 KB Output is correct

Subtask #4
6.0 / 6.0

# Verdict Execution time Memory Grader output
1 Correct 23 ms 1644 KB Output is correct
2 Correct 19 ms 1516 KB Output is correct

Subtask #5
6.0 / 6.0

# Verdict Execution time Memory Grader output
1 Correct 21 ms 1516 KB Output is correct
2 Correct 20 ms 1516 KB Output is correct
3 Correct 21 ms 1516 KB Output is correct

Subtask #6
7.0 / 7.0

# Verdict Execution time Memory Grader output
1 Correct 1 ms 1132 KB Output is correct
2 Correct 1 ms 1132 KB Output is correct

Subtask #7
7.0 / 7.0

# Verdict Execution time Memory Grader output
1 Correct 32 ms 2284 KB Output is correct
2 Correct 30 ms 2284 KB Output is correct

Subtask #8
8.0 / 8.0

# Verdict Execution time Memory Grader output
1 Correct 36 ms 2284 KB Output is correct
2 Correct 34 ms 2284 KB Output is correct
3 Correct 34 ms 2284 KB Output is correct
4 Correct 34 ms 2284 KB Output is correct

Subtask #9
8.0 / 8.0

# Verdict Execution time Memory Grader output
1 Correct 37 ms 2284 KB Output is correct
2 Correct 37 ms 2284 KB Output is correct
3 Correct 36 ms 2284 KB Output is correct
4 Correct 36 ms 2284 KB Output is correct

Subtask #10
0 / 10.0

# Verdict Execution time Memory Grader output
1 Correct 334 ms 490464 KB Output is correct
2 Runtime error 1037 ms 524288 KB Execution killed with signal 6 (could be triggered by violating memory limits)

Subtask #11
10.0 / 10.0

# Verdict Execution time Memory Grader output
1 Correct 457 ms 497628 KB Output is correct
2 Correct 447 ms 497628 KB Output is correct
3 Correct 419 ms 497628 KB Output is correct

Subtask #12
0 / 10.0

# Verdict Execution time Memory Grader output
1 Correct 625 ms 497636 KB Output is correct
2 Correct 719 ms 497884 KB Output is correct
3 Runtime error 1099 ms 524288 KB Execution killed with signal 6 (could be triggered by violating memory limits)

Subtask #13
0 / 10.0

# Verdict Execution time Memory Grader output
1 Correct 783 ms 497628 KB Output is correct
2 Correct 846 ms 498396 KB Output is correct
3 Runtime error 1092 ms 524288 KB Execution killed with signal 6 (could be triggered by violating memory limits)