Submission #94725

# Submission time Handle Problem Language Result Execution time Memory
94725 2019-01-23T11:32:35 Z SuperJava Pairs (IOI07_pairs) C++17
70 / 100
410 ms 31296 KB
//fold
#ifndef KHALIL
#include <bits/stdc++.h>
#else
#include "header.h"
#endif
#define endl '\n'
#define mp make_pair
#define tostr(x) static_cast<ostringstream&>((ostringstream()<<dec<<x)).str()
#define rep(i,begin,end) for(auto i = begin;i < end;i++)
#define repr(i,begin,end) for(auto i = begin-1;i >= end;i--)
#define pb push_back
#define sz(a) ((int)(a).size())
#define fi first
#define se second
#define abs(a) ((a) < (0) ? (-1)*(a) : (a))
#define SQ(a) ((a)*(a))
#define eqd(a,b) (abs(a-b)<1e-9)
#define X real()
#define Y imag()
using namespace std;
typedef long long ll;
typedef long double ld;
template <typename t> t in(t q){cin >> q;return q;}
template <typename T> ostream& operator<<(ostream& os, const vector<T>& v){os << "[";for (int i = 0; i < sz(v); ++i) { os << v[i]; if (i != sz(v) - 1) os << ",";}os << "]";return os;}
template <typename T, typename S>ostream& operator<<(ostream& os, const map<T, S>& v){for (auto it : v)os << "(" << it.first << ":" << it.second << ")";return os;}
template <typename T, typename S>ostream& operator<<(ostream& os, const pair<T, S>& v){os << "(" << v.first << "," << v.second << ")";return os;}
const long double PI = acosl(-1);
mt19937 rng(chrono::steady_clock::now().time_since_epoch().count());
mt19937_64 rng64(chrono::steady_clock::now().time_since_epoch().count());
inline int rand(int l,int r){return uniform_int_distribution<int>(l, r)(rng);}
inline ll rand(ll l,ll r){return uniform_int_distribution<ll>(l, r)(rng64);}
//endfold
#define  N  (100'005)
#define MOD (1000000000l + 7l)
#define OO (1050000000)
#define OOL (1100000000000000000)
 
vector<pair<ll,ll>> base;
vector<ll> seg[N*4];

void build(int n,int l,int r){
	if(l == r){
		seg[n].push_back(base[l].first-base[l].second);
		return;
	}
	int mid = (l+r)/2;
	build(2*n,l,mid);
	build(2*n+1,mid+1,r);
	seg[n] = seg[2*n];
	seg[n].insert(seg[n].end(),seg[2*n+1].begin(),seg[2*n+1].end());
	sort(seg[n].begin(), seg[n].end());
}

ll get(int n,int l,int r,ll x,ll e,ll q,ll w){
	if(l > e || r < x) return 0;
	if(x <= l && r <= e)
		return upper_bound(seg[n].begin(), seg[n].end(), w) -
				lower_bound(seg[n].begin(), seg[n].end(), q);
	ll mid = (l+r)/2;
	return get(2*n,l,mid,x,e,q,w)+get(2*n+1,mid+1,r,x,e,q,w);
}

int u[80][80][80];

int main(){
	//fold
	ios::sync_with_stdio(0);
	cin.tie(0); cout.tie(0);
	cout << setprecision(10);
	//endfold
	ll b,n,d,m;
	cin >> b >> n >> d >> m;
	if(b == 1){
		vector<ll> v;
		rep(i,0,n){
			ll w;
			cin >> w;
			v.push_back(w);
		}
		sort(v.begin(), v.end());
		ll ans = 0;
		rep(i,0,n){
			ans += upper_bound(v.begin(), v.end(),v[i]+d)-lower_bound(v.begin(), v.end(),v[i]-d)-1;
		}
		cout << ans/2;
	}else if(b == 2){
		vector<ll> z;
		rep(i,0,n){
			ll w,q;
			cin >> q >> w;
			base.push_back({q,w});
			z.push_back(q+w);
		}
		sort(z.begin(), z.end());
		sort(base.begin(), base.end(),
			[](pair<int,int> a,pair<int,int> b){
				return a.first+a.second < b.first+b.second;
			});
		build(1,0,n-1);
		ll ans = 0;
		rep(i,0,n){
			ll l = lower_bound(z.begin(), z.end(),z[i]-d)-z.begin();
			ll r = upper_bound(z.begin(), z.end(),z[i]+d)-z.begin()-1;
			ans += get(1,0,n-1,l,r,base[i].first-d-base[i].second,
									base[i].first+d-base[i].second);
		}
		cout << (ans-n)/2 << endl;
	}else{
		vector<tuple<ll,ll,ll>> v;
		rep(i,0,n){
			ll w,q,e;
			cin >> q >> w >> e;
			v.push_back({q,w,e});
			u[q][w][e]++;
		}
		sort(v.begin(), v.end());
		if(n <= 1000){
			ll ans = 0;
			rep(i,0,n){
				rep(j,i+1,n){
					ll q1,w1,e1;
					ll q2,w2,e2;
					tie(q1,w1,e1) = v[i];
					tie(q2,w2,e2) = v[j];			
					if(abs(q1-q2)+abs(w2-w1)+abs(e2-e1) <= d) ans++;
					if(abs(q1-q2) > d) break;
				}
			}
			cout << ans;
		}else{
			ll ans = 0;
			rep(i,0,n){
				int x,y,z;
				tie(x,y,z) = v[i];
				rep(j,x,x+d+1){
					rep(k,y,y+d-(j-x)+1){
						rep(l,z,z+d-(j-x)-(k-y)+1){
							ans += u[j][k][l];
						}
					}
				}
				ans -= u[x][y][z];
			}
			rep(i,1,76){
				rep(j,1,76){
					rep(k,1,76){
						ans += u[i][j][k]*(u[i][j][k]-1)/2;
					}
				}
			}
			cout << ans;
		}
	}
	return 0;
}
# Verdict Execution time Memory Grader output
1 Correct 8 ms 9720 KB Output is correct
2 Correct 9 ms 9720 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 8 ms 9720 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 29 ms 11072 KB Output is correct
2 Correct 27 ms 10996 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 39 ms 11108 KB Output is correct
2 Correct 31 ms 10992 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 38 ms 10996 KB Output is correct
2 Correct 36 ms 10996 KB Output is correct
3 Correct 33 ms 11124 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 10 ms 9848 KB Output is correct
2 Correct 11 ms 9848 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 209 ms 31056 KB Output is correct
2 Correct 124 ms 31080 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 261 ms 31296 KB Output is correct
2 Correct 260 ms 31084 KB Output is correct
3 Correct 236 ms 31080 KB Output is correct
4 Correct 243 ms 31208 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 290 ms 31112 KB Output is correct
2 Correct 308 ms 31220 KB Output is correct
3 Correct 229 ms 31208 KB Output is correct
4 Correct 213 ms 31084 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 12 ms 11384 KB Output is correct
2 Correct 13 ms 11484 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 40 ms 13160 KB Output is correct
2 Incorrect 60 ms 13672 KB Output isn't correct
3 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Incorrect 68 ms 14572 KB Output isn't correct
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Incorrect 410 ms 14956 KB Output isn't correct
2 Halted 0 ms 0 KB -