#undef _GLIBCXX_DEBUG // disable run-time bound checking, etc
#pragma GCC optimize("Ofast,inline") // Ofast = O3,fast-math,allow-store-data-races,no-protect-parens
#pragma GCC optimize ("unroll-loops")
#pragma GCC target("bmi,bmi2,lzcnt,popcnt") // bit manipulation
#pragma GCC target("movbe") // byte swap
#pragma GCC target("aes,pclmul,rdrnd") // encryption
#pragma GCC target("avx,avx2,f16c,fma,sse3,ssse3,sse4.1,sse4.2") // SIMD
#include <bits/stdc++.h>
#include <bits/extc++.h>
using namespace std;
#define enablell 1
typedef long long ll;
typedef unsigned long long ull;
#if enablell
#define int ll
#define inf int(1e18)
#define float double
#else
const int inf = int(2e9);
#endif
typedef vector<ull> vull;
typedef vector<int> vi;
typedef vector<vi> vvi;
typedef vector<vvi> vvvi;
typedef vector<vvvi> vvvvi;
typedef vector<bool> vb;
typedef vector<vb> vvb;
typedef vector<vvb> vvvb;
typedef pair<int, int> p2;
typedef vector<p2> vp2;
typedef vector<vp2> vvp2;
typedef vector<vvp2> vvvp2;
typedef tuple<int, int, int> p3;
typedef vector<p3> vp3;
typedef vector<vp3> vvp3;
typedef vector<vvp3> vvvp3;
typedef tuple<int, int, int, int> p4;
typedef vector<p4> vp4;
#define PBDS 0
#define _LOCAL _MSC_VER > 0
#if _LOCAL
#define gc() getchar()
#define popcount(x) __popcnt(x)
#define leading_zeros(x) _lzcnt_u32(x)
uint32_t clz(uint32_t x) { return _lzcnt_u32(x); }
uint32_t ctz(uint32_t x) { return _tzcnt_u32(x); }
#define bswap64(x) _byteswap_uint64(x)
#define assert(x) debassert(x)
#else
#define popcount(x) __builtin_popcount(x)
uint32_t clz(uint32_t x) { return __builtin_clz(x); }
uint32_t ctz(uint32_t x) { return __builtin_ctzll(x); }
#define bswap64(x) __builtin_bswap64(x)
#define gc() getchar_unlocked()
#if PBDS
using namespace __gnu_pbds;
// lower_bound is now upper_bound and vice versa (multiset). generally a bit broken
template<typename T> using indexed_multiset = tree<int, null_type, less_equal<T>, rb_tree_tag, tree_order_statistics_node_update>;
template<typename T> using indexed_set = tree<T, null_type, less<T>, rb_tree_tag, tree_order_statistics_node_update>;
struct chash { // large odd number for C
const uint64_t C = ll(4e18 * acos(0)) | 71;
ll operator()(ll x) const { return __builtin_bswap64(x * C); }
};
template<typename T, typename U> using fast_map = __gnu_pbds::gp_hash_table<T, U, chash>;
template<typename T> using fast_set = __gnu_pbds::gp_hash_table<T, null_type, chash>;
template<typename T, typename H> using fast_set_h = __gnu_pbds::gp_hash_table<T, null_type, H>;
#endif
#endif
#define FAST_INPUT 0
#define FILE_TC 0
#if FILE_TC && _LOCAL
//ifstream filein("C:\\users\\matis\\source\\repos\\comp_prog\\x64\\release\\in.txt");
ifstream filein("C:\\users\\matis\\source\\repos\\comp_prog\\x64\\debug\\in.txt");
//ifstream filein("E:\\downloads\\test_data\\test_data\\005-case05.in");
//ifstream filein("E:\\desktop\\po-repos\\swedish-olympiad-2023\\online\\tomtplanering\\data\\secret\\group10\\010-case10.in");
#define cin filein
void fast() {}
#else
inline void fast() { ios::sync_with_stdio(false); cin.tie(NULL); cout.tie(NULL); }
#endif
#if FAST_INPUT && (!FILE_TC || !_LOCAL)
inline void read(int& v) { v = 0; int sign = 1; char c = gc(); if (c == '-') { sign = -1; } else { v += c - '0'; } while ((c = gc()) && c != ' ' && c != '\n') { if (c == EOF) { v = -1; return; } v *= 10; v += c - '0'; } v *= sign; }
inline void read(int& u, int& v) { read(u); read(v); }
inline void read(int& u, int& v, int& k) { read(u); read(v); read(k); }
//inline void read(int& v) { char c; while ((c = getchar()) != EOF && c != ' ' && c != '\n') { v *= 10; v += c - '0'; } }
inline void read(string& s) { char c; while ((c = gc()) != EOF && c != '\n' && c != ' ') { s.push_back(c); } }
inline void readline(string& s) { char c; while ((c = gc()) != EOF && c != '\n') { s.push_back(c); } }
#else
template <typename T> inline void read(T& a) { cin >> a; }
template <typename T> inline void read(T& a, T& b) { cin >> a >> b; }
template <typename T> inline void read(T& a, T& b, T& c) { cin >> a >> b >> c; }
#endif
#define quit cout << endl; _Exit(0);
#define dread(type, a) type a; read(a)
#define dread2(type, a, b) dread(type, a); dread(type, b)
#define dread3(type, a, b, c) dread2(type, a, b); dread(type, c)
#define dread4(type, a, b, c, d) dread3(type, a, b, c); dread(type, d)
#define dread5(type, a, b, c, d, e) dread4(type, a, b, c, d); dread(type, e)
#define readvector(type, name, size) vector<type> name(size); rep(i,size) {dread(type,temp); name[i]=temp;}
#ifdef _DEBUG
#define noop cout << "";
#define deb __debugbreak();
#define debassert(expr) if(!(expr)) deb;
#define debif(expr) if(expr) deb;
#else
#define noop ;
#define deb ;
#define debif(expr) ;
#define debassert(expr) ;
#endif
#define rep(i, high) for (int i = 0; i < high; i++)
#define repp(i, low, high) for (int i = low; i < high; i++)
#define repe(i, container) for (auto& i : container)
#define per(i, high) for (int i = high-1; i >= 0; i--)
#define perr(i, low, high) for (int i = high-1; i >= low; i--)
#define all(a) a.begin(),a.end()
#define rall(a) a.rbegin(),a.rend()
#define setcontains(set, x) (set.find(x) != set.end())
#define within(a, b, c, d) (a >= 0 && a < b && c >= 0 && c < d)
#define sz(container) ((int)container.size())
#define mp(a,b) (make_pair(a,b))
#define ceildiv(x,y) ((x + y - 1) / y)
template <typename T, typename U> inline void operator+=(pair<T, U>& l, const pair<T, U>& r) { l = { l.first + r.first,l.second + r.second }; }
template <typename T, typename U> inline pair<T, U> operator+(const pair<T, U> l, const pair<T, U> r) { return { l.first + r.first, l.second + r.second }; }
template <typename T, typename U> inline pair<T, U> operator-(const pair<T, U> l, const pair<T, U> r) { return { l.first - r.first, l.second - r.second }; }
template <typename T, typename U> inline pair<T, U> operator*(const pair<T, U> l, const int m) { return { l.first * m, l.second * m }; }
template <typename Out> inline void split(const string& s, char delim, Out result) { istringstream iss(s); string item; while (getline(iss, item, delim)) { *result++ = item; } }
inline vector<string> split(const string& s, char delim) { vector<string> elems; split(s, delim, back_inserter(elems)); return elems; }
vector<string> split(string s, string d) { size_t k = 0, n, e = d.length(); string t; vector<string> r; while ((n = s.find(d, k)) != string::npos) { t = s.substr(k, n - k); k = n + e; r.push_back(t); } r.push_back(s.substr(k)); return r; }
ll binpow(ll a, ll b) { ll res = 1; while (b > 0) { if (b & 1) res = res * a; a = a * a; b >>= 1; } return res; }
ll binpow(ll a, ll b, ll m) { a %= m; long long res = 1; while (b > 0) { if (b & 1) res = res * a % m; a = a * a % m; b >>= 1; } return res; } // For a < 2^31
#if 1
auto Start = chrono::high_resolution_clock::now();
void resettimer() { Start = chrono::high_resolution_clock::now(); }
int elapsedmillis() { return chrono::duration_cast<chrono::milliseconds>(chrono::high_resolution_clock::now() - Start).count(); }
random_device rd;
mt19937 rng(rd());
template<typename T, typename U> inline int randint(T lo, U hi) { return uniform_int_distribution<int>((int)lo, (int)hi)(rng); } // [lo,hi]
template<typename T> inline T randel(vector<T>& v) { return v[uniform_int_distribution<int>(int(0), int(v.size()) - int(1))(rng)]; } // [lo,hi]
#endif
const ll mod = 1e9 + 7;
vp2 dirs = { {0,1},{0,-1},{1,0},{-1,0}, {0,0} };
int square(int x)
{
return x * x;
}
int dist(p2 a, p2 b)
{
return square(a.first - b.first) + square(a.second - b.second);
}
struct UF
{
vi par;
vi size;
UF(int n) : par(n), size(n, 1)
{
rep(i, n)par[i] = i;
}
int find(int x) { return x == par[x] ? x : par[x] = find(par[x]); }
void merge(int a, int b)
{
a = find(a); b = find(b);
if (a == b) return;
if (size[a] < size[b]) swap(a, b);
par[b] = a;
size[a] += size[b];
}
};
int32_t main()
{
fast();
dread2(int, n, m);
dread2(int, w, h);
vp3 trees;
rep(i, n)
{
dread3(int, x, y, r);
trees.emplace_back(x, y, r);
}
typedef tuple<double, int, int> dp3;
vector<dp3> intersections;
int bottomwall = n;
int leftwall = n + 1;
int topwall = n + 2;
int rightwall = n + 3;
rep(i, n)
{
int x, y, r;
tie(x, y, r) = trees[i];
intersections.emplace_back(y-r,i, bottomwall);
intersections.emplace_back(x-r,i, leftwall);
intersections.emplace_back(abs(h-y)-r,i, topwall);
intersections.emplace_back(abs(w-x)-r,i, rightwall);
repp(j, i+1, n)
{
if (i == j) continue;
int x2, y2, r2;
tie(x2, y2, r2) = trees[j];
intersections.emplace_back(sqrt(dist(p2(x, y), p2(x2, y2))) - r - r2, i, j);
}
}
sort(all(intersections));
vp3 queries;
rep(i, m)
{
dread2(int, r, c);
queries.emplace_back(r, c-1, i);
}
sort(all(queries));
vi ans(m);
UF uf(n + 4);
int p = 0;
int bottomleft = 1;
int topleft = 8;
int topright = 4;
int bottomright = 2;
repe(q, queries)
{
int r, c, i;
tie(r, c, i) = q;
c = 1 << c;
while (p<intersections.size()&&get<0>(intersections[p])+0.001<r*2)
{
int d, a, b;
tie(d, a, b) = intersections[p];
uf.merge(a, b);
p++;
}
int qans = 1|2|4|8;
if (uf.find(leftwall)==uf.find(topwall))
{
if (c==topleft) qans = topleft;
else qans &= ~topleft;
}
if (uf.find(topwall)==uf.find(rightwall))
{
if (c==topright) qans = topright;
else qans &= ~topright;
}
if (uf.find(rightwall)==uf.find(bottomwall))
{
if (c==bottomright) qans = bottomright;
else qans &= ~bottomright;
}
if (uf.find(bottomwall)==uf.find(leftwall))
{
if (c==bottomleft) qans = bottomleft;
else qans &= ~bottomleft;
}
if (uf.find(leftwall)==uf.find(rightwall))
{
if (c == bottomleft || c == bottomright) qans &= ~(topleft | topright);
else qans &= ~(bottomleft | bottomright);
}
if (uf.find(topwall)==uf.find(bottomwall))
{
if (c == bottomleft || c == topleft) qans &= ~(topright | bottomright);
else qans &= ~(bottomleft | topleft);
}
ans[i] = qans;
}
rep(i, m)
{
rep(j, 4)
{
if (ans[i] & (1 << j)) cout << j + 1;
}
cout << "\n";
}
quit;
}
Compilation message
park.cpp: In function 'int32_t main()':
park.cpp:252:11: warning: comparison of integer expressions of different signedness: 'll' {aka 'long long int'} and 'std::vector<std::tuple<double, long long int, long long int> >::size_type' {aka 'long unsigned int'} [-Wsign-compare]
252 | while (p<intersections.size()&&get<0>(intersections[p])+0.001<r*2)
| ~^~~~~~~~~~~~~~~~~~~~~
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Correct |
392 ms |
49808 KB |
Output is correct |
| 2 |
Correct |
386 ms |
49848 KB |
Output is correct |
| 3 |
Correct |
400 ms |
49824 KB |
Output is correct |
| 4 |
Correct |
384 ms |
49896 KB |
Output is correct |
| 5 |
Correct |
384 ms |
49892 KB |
Output is correct |
| 6 |
Correct |
392 ms |
49860 KB |
Output is correct |
| 7 |
Correct |
404 ms |
49856 KB |
Output is correct |
| 8 |
Correct |
361 ms |
49872 KB |
Output is correct |
| 9 |
Correct |
1 ms |
212 KB |
Output is correct |
| 10 |
Correct |
1 ms |
212 KB |
Output is correct |
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Correct |
64 ms |
4456 KB |
Output is correct |
| 2 |
Correct |
42 ms |
5264 KB |
Output is correct |
| 3 |
Correct |
42 ms |
5308 KB |
Output is correct |
| 4 |
Correct |
41 ms |
5316 KB |
Output is correct |
| 5 |
Correct |
49 ms |
5340 KB |
Output is correct |
| 6 |
Correct |
51 ms |
5380 KB |
Output is correct |
| 7 |
Correct |
38 ms |
4928 KB |
Output is correct |
| 8 |
Correct |
41 ms |
4876 KB |
Output is correct |
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Correct |
392 ms |
49808 KB |
Output is correct |
| 2 |
Correct |
386 ms |
49848 KB |
Output is correct |
| 3 |
Correct |
400 ms |
49824 KB |
Output is correct |
| 4 |
Correct |
384 ms |
49896 KB |
Output is correct |
| 5 |
Correct |
384 ms |
49892 KB |
Output is correct |
| 6 |
Correct |
392 ms |
49860 KB |
Output is correct |
| 7 |
Correct |
404 ms |
49856 KB |
Output is correct |
| 8 |
Correct |
361 ms |
49872 KB |
Output is correct |
| 9 |
Correct |
1 ms |
212 KB |
Output is correct |
| 10 |
Correct |
1 ms |
212 KB |
Output is correct |
| 11 |
Correct |
64 ms |
4456 KB |
Output is correct |
| 12 |
Correct |
42 ms |
5264 KB |
Output is correct |
| 13 |
Correct |
42 ms |
5308 KB |
Output is correct |
| 14 |
Correct |
41 ms |
5316 KB |
Output is correct |
| 15 |
Correct |
49 ms |
5340 KB |
Output is correct |
| 16 |
Correct |
51 ms |
5380 KB |
Output is correct |
| 17 |
Correct |
38 ms |
4928 KB |
Output is correct |
| 18 |
Correct |
41 ms |
4876 KB |
Output is correct |
| 19 |
Correct |
428 ms |
54320 KB |
Output is correct |
| 20 |
Correct |
422 ms |
54196 KB |
Output is correct |
| 21 |
Correct |
426 ms |
54320 KB |
Output is correct |
| 22 |
Correct |
430 ms |
54196 KB |
Output is correct |
| 23 |
Correct |
418 ms |
54264 KB |
Output is correct |
| 24 |
Correct |
404 ms |
54272 KB |
Output is correct |
