| # |
Submission time |
Handle |
Problem |
Language |
Result |
Execution time |
Memory |
| 556691 |
2022-05-03T16:53:27 Z |
Bungmint |
Pairs (IOI07_pairs) |
C++17 |
|
423 ms |
524288 KB |
// Copyright © 2022 Youngmin Park. All rights reserved.
//#pragma GCC optimize("O3")
//#pragma GCC target("avx2")
#include <bits/stdc++.h>
using namespace std;
using ll = long long;
using vi = vector<int>;
using pii = pair<int, int>;
using vpi = vector<pii>;
using pll = pair<ll, ll>;
using vl = vector<ll>;
using vpl = vector<pll>;
using ld = long double;
template <typename T, size_t SZ>
using ar = array<T, SZ>;
template <typename T>
using pqg = priority_queue<T, vector<T>, greater<T>>;
#define all(v) (v).begin(), (v).end()
#define pb push_back
#define sz(x) (int)(x).size()
#define fi first
#define se second
#define lb lower_bound
#define ub upper_bound
constexpr int INF = 1e9;
constexpr ll LINF = 1e18;
const ld PI = acos((ld)-1.0);
constexpr int dx[4] = {1, 0, -1, 0}, dy[4] = {0, 1, 0, -1};
mt19937_64 rng(chrono::steady_clock::now().time_since_epoch().count());
template <typename T>
constexpr bool ckmin(T &a, const T &b) { return b < a ? a = b, 1 : 0; }
template <typename T>
constexpr bool ckmax(T &a, const T &b) { return b > a ? a = b, 1 : 0; }
template <typename A, typename B>
ostream &operator<<(ostream &os, const pair<A, B> &p)
{
return os << '(' << p.first << ", " << p.second << ')';
}
template <typename T_container, typename T = typename enable_if<!is_same<T_container, string>::value, typename T_container::value_type>::type>
ostream &operator<<(ostream &os, const T_container &v)
{
os << '{';
string sep;
for (const T &x : v)
os << sep << x, sep = ", ";
return os << '}';
}
template <typename T>
ostream &operator<<(ostream &os, const deque<T> &v) {
os << vector<T>(all(v));
return os;
}
template <typename T, typename S, typename C>
ostream &operator<<(ostream &os, priority_queue<T, S, C> pq) {
vector<T> v;
while (sz(pq)) {
v.pb(pq.top());
pq.pop();
}
os << v;
return os;
}
void dbg_out()
{
cerr << "\033[0m" << endl;
}
template <typename Head, typename... Tail>
void dbg_out(Head H, Tail... T)
{
cerr << ' ' << H;
dbg_out(T...);
}
#ifdef LOCAL
#define dbg(...) cerr << "\033[1;35m(" << #__VA_ARGS__ << "):\033[33m", dbg_out(__VA_ARGS__)
#else
#define dbg(...) 42
#endif
inline namespace RecursiveLambda
{
template <typename Fun>
struct y_combinator_result
{
Fun fun_;
template <typename T>
explicit y_combinator_result(T &&fun) : fun_(forward<T>(fun)) {}
template <typename... Args>
decltype(auto) operator()(Args &&...args)
{
return fun_(ref(*this), forward<Args>(args)...);
}
};
template <typename Fun>
decltype(auto) y_combinator(Fun &&fun)
{
return y_combinator_result<decay_t<Fun>>(forward<Fun>(fun));
}
};
class Range {
struct Iter
{
int iter;
constexpr Iter(int it) noexcept : iter(it) {}
constexpr void operator++() noexcept { iter++; }
constexpr bool operator!=(const Iter &other) const noexcept { return iter != other.iter; }
constexpr int operator*() const noexcept { return iter; }
};
const Iter first, last;
public:
explicit constexpr Range(const int f, const int l) noexcept : first(f), last(max(f, l)) {}
constexpr Iter begin() const noexcept { return first; }
constexpr Iter end() const noexcept { return last; }
};
constexpr Range rep(const int l, const int r) noexcept { return Range(l, r); }
constexpr Range rep(const int n) noexcept { return Range(0, n); }
class ReversedRange {
struct Iter {
int itr;
constexpr Iter(const int pos) noexcept : itr(pos) {}
constexpr void operator++() noexcept { --itr; }
constexpr bool operator!=(const Iter& other) const noexcept { return itr != other.itr; }
constexpr int operator*() const noexcept { return itr; }
};
const Iter first, last;
public:
explicit constexpr ReversedRange(const int f, const int l) noexcept
: first(l - 1), last(min(f, l) - 1) {}
constexpr Iter begin() const noexcept { return first; }
constexpr Iter end() const noexcept { return last; }
};
constexpr ReversedRange per(const int l, const int r) noexcept { return ReversedRange(l, r); }
constexpr ReversedRange per(const int n) noexcept { return ReversedRange(0, n); }
/**
* Description: 2D Binary Indexed Tree
* Source: Own
* Verification: https://cses.fi/problemset/task/1739
* Time complexity: O(log^2 n) update/query
* Memory complexity: O(n^2)
* O-indexing
*/
template <typename T>
struct BIT2D {
vector<vector<T>> bit;
int N, M;
BIT2D(int n, int m) : N(n), M(m) {
bit.resize(n + 1);
for (auto &v : bit) v.resize(m + 1);
}
void upd(int x, int y, T v) {
for (x++; x <= N; x += x & -x) {
for (int i = y + 1; i <= M; i += i & -i) {
bit[x][i] += v;
}
}
}
T query(int x, int y) {
T res = T();
for (x++; x; x -= x & -x) {
for (int i = y + 1; i; i -= i & -i) {
res += bit[x][i];
}
}
return res;
}
};
template <typename T>
struct BIT3D {
vector<vector<vector<T>>> bit;
int N, M, K;
BIT3D(int n, int m, int k) : N(n), M(m), K(k) {
bit.resize(n + 1);
for (auto &v : bit) {
v.resize(m + 1);
for (auto &vv : v) vv.resize(k + 1);
}
}
void upd(int x, int y, int z, T v) {
for (x++; x <= N; x += x & -x) {
for (int i = y + 1; i <= M; i += i & -i) {
for (int j = z + 1; j <= K; j += j & -j) {
bit[x][i][j] += v;
}
}
}
}
T query(int x, int y, int z) {
T res = T();
for (x++; x; x -= x & -x) {
for (int i = y + 1; i; i -= i & -i) {
for (int j = z + 1; j; j -= j & -j) {
res += bit[x][i][j];
}
}
}
return res;
}
};
vector<ar<int, 3>> a;
vi xs, ys, zs;
int b, n, d, m;
ll solve1() {
sort(all(a));
ll res = 0;
for (int i : rep(n)) {
ar<int, 3> targ = a[i];
targ[0] += d;
int it = ub(all(a), targ) - a.begin();
res += it - i - 1;
}
return res;
}
ll solve2() {
for (int i : rep(n)) {
auto [x, y, _] = a[i];
xs.pb(x + y), ys.pb(x - y);
xs.pb(x + y + d), xs.pb(x + y - d);
ys.pb(x - y + d), ys.pb(x - y - d);
}
sort(all(xs)), sort(all(ys));
xs.resize(unique(all(xs)) - xs.begin());
ys.resize(unique(all(ys)) - ys.begin());
BIT2D<int> bit(sz(xs), sz(ys));
ll res = 0;
auto getX = [&](int x) -> int { return lb(all(xs), x) - xs.begin(); };
auto getY = [&](int y) -> int { return lb(all(ys), y) - ys.begin(); };
for (int i : rep(n)) {
auto [x, y, _] = a[i];
bit.upd(getX(x + y), getY(x - y), 1);
}
for (int i : rep(n)) {
auto [x, y, _] = a[i];
int x1 = (x + y - d), x2 = (x + y + d);
int y1 = (x - y - d), y2 = (x - y + d);
x1 = getX(x1), x2 = getX(x2);
y1 = getY(y1), y2 = getY(y2);
res += bit.query(x2, y2) - bit.query(x1 - 1, y2) - bit.query(x2, y1 - 1) + bit.query(x1 - 1, y1 - 1);
res--;
}
return res / 2LL;
}
ll solve3() {
for (int i : rep(n)) {
auto [x, y, z] = a[i];
xs.pb(x + y + z), ys.pb(x + y - z), zs.pb(x - y - z);
xs.pb(x + y + z + d), xs.pb(x + y + z - d);
ys.pb(x + y - z + d), ys.pb(x + y - z - d);
zs.pb(x - y - z + d), zs.pb(x - y - z - d);
}
sort(all(xs)), sort(all(ys)), sort(all(zs));
xs.resize(unique(all(xs)) - xs.begin());
ys.resize(unique(all(ys)) - ys.begin());
zs.resize(unique(all(zs)) - zs.begin());
BIT3D<int> bit(sz(xs), sz(ys), sz(zs));
dbg(xs, ys, zs);
ll res = 0;
auto getX = [&](int x) -> int { return lb(all(xs), x) - xs.begin(); };
auto getY = [&](int y) -> int { return lb(all(ys), y) - ys.begin(); };
auto getZ = [&](int z) -> int { return lb(all(zs), z) - zs.begin(); };
for (int i : rep(n)) {
auto [x, y, z] = a[i];
dbg(x + y + z, x + y - z, x - y - z);
bit.upd(getX(x + y + z), getY(x + y - z), getZ(x - y - z), 1);
}
for (int i : rep(n)) {
auto [x, y, z] = a[i];
int x1 = (x + y + z - d), x2 = (x + y + z + d);
int y1 = (x + y - z - d), y2 = (x + y - z + d);
int z1 = (x - y - z - d), z2 = (x - y - z + d);
x1 = getX(x1), x2 = getX(x2);
y1 = getY(y1), y2 = getY(y2);
z1 = getZ(z1), z2 = getZ(z2);
res += bit.query(x2, y2, z2);
res -= bit.query(x2, y2, z1 - 1);
res -= bit.query(x2, y1 - 1, z2);
res -= bit.query(x1 - 1, y2, z2);
res += bit.query(x2, y1 - 1, z1 - 1);
res += bit.query(x1 - 1, y2, z1 - 1);
res += bit.query(x1 - 1, y1 - 1, z2);
res -= bit.query(x1 - 1, y1 - 1, z1 - 1);
res--;
dbg(res);
}
dbg(res);
return res / 2LL;
}
void solve()
{
cin >> b >> n >> d >> m;
a.resize(n);
for (int i : rep(n)) for (int j : rep(b)) {
cin >> a[i][j];
}
if (b == 1) cout << solve1() << '\n';
else if (b == 2) cout << solve2() << '\n';
else cout << solve3() << '\n';
}
int main()
{
cin.tie(0)->sync_with_stdio(0);
cin.exceptions(cin.failbit);
int testcase = 1;
// cin >> testcase;
while (testcase--)
{
solve();
}
#ifdef LOCAL
cerr << "Time elapsed: " << 1.0 * (double)clock() / CLOCKS_PER_SEC << " s.\n";
#endif
}
Compilation message
pairs.cpp: In function 'll solve3()':
pairs.cpp:80:18: warning: statement has no effect [-Wunused-value]
80 | #define dbg(...) 42
| ^~
pairs.cpp:269:2: note: in expansion of macro 'dbg'
269 | dbg(xs, ys, zs);
| ^~~
pairs.cpp:80:18: warning: statement has no effect [-Wunused-value]
80 | #define dbg(...) 42
| ^~
pairs.cpp:276:3: note: in expansion of macro 'dbg'
276 | dbg(x + y + z, x + y - z, x - y - z);
| ^~~
pairs.cpp:80:18: warning: statement has no effect [-Wunused-value]
80 | #define dbg(...) 42
| ^~
pairs.cpp:296:3: note: in expansion of macro 'dbg'
296 | dbg(res);
| ^~~
pairs.cpp:80:18: warning: statement has no effect [-Wunused-value]
80 | #define dbg(...) 42
| ^~
pairs.cpp:298:2: note: in expansion of macro 'dbg'
298 | dbg(res);
| ^~~
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Correct |
0 ms |
212 KB |
Output is correct |
| 2 |
Correct |
0 ms |
212 KB |
Output is correct |
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Correct |
1 ms |
332 KB |
Output is correct |
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Correct |
21 ms |
1876 KB |
Output is correct |
| 2 |
Correct |
20 ms |
1876 KB |
Output is correct |
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Correct |
26 ms |
2260 KB |
Output is correct |
| 2 |
Correct |
24 ms |
2264 KB |
Output is correct |
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Correct |
26 ms |
2260 KB |
Output is correct |
| 2 |
Correct |
27 ms |
2256 KB |
Output is correct |
| 3 |
Correct |
27 ms |
2268 KB |
Output is correct |
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Correct |
16 ms |
34628 KB |
Output is correct |
| 2 |
Correct |
19 ms |
34760 KB |
Output is correct |
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Correct |
79 ms |
5312 KB |
Output is correct |
| 2 |
Correct |
83 ms |
5472 KB |
Output is correct |
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Correct |
161 ms |
23156 KB |
Output is correct |
| 2 |
Correct |
174 ms |
47556 KB |
Output is correct |
| 3 |
Correct |
147 ms |
25204 KB |
Output is correct |
| 4 |
Correct |
171 ms |
43152 KB |
Output is correct |
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Runtime error |
250 ms |
524288 KB |
Execution killed with signal 9 |
| 2 |
Halted |
0 ms |
0 KB |
- |
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Incorrect |
29 ms |
46040 KB |
Output isn't correct |
| 2 |
Halted |
0 ms |
0 KB |
- |
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Incorrect |
102 ms |
6308 KB |
Output isn't correct |
| 2 |
Halted |
0 ms |
0 KB |
- |
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Incorrect |
334 ms |
37108 KB |
Output isn't correct |
| 2 |
Halted |
0 ms |
0 KB |
- |
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Incorrect |
423 ms |
95888 KB |
Output isn't correct |
| 2 |
Halted |
0 ms |
0 KB |
- |
