#include <bits/stdc++.h>
#include <ext/pb_ds/tree_policy.hpp>
#include <ext/pb_ds/assoc_container.hpp>
using namespace std;
using namespace __gnu_pbds;
typedef long long ll;
typedef long double ld;
typedef complex<ld> cd;
typedef pair<int, int> pi;
typedef pair<ll,ll> pl;
typedef pair<ld,ld> pd;
typedef vector<int> vi;
typedef vector<ld> vd;
typedef vector<ll> vl;
typedef vector<pi> vpi;
typedef vector<pl> vpl;
typedef vector<cd> vcd;
template <class T> using Tree = tree<T, null_type, less<T>, rb_tree_tag,tree_order_statistics_node_update>;
#define FOR(i, a, b) for (int i=a; i<(b); i++)
#define F0R(i, a) for (int i=0; i<(a); i++)
#define FORd(i,a,b) for (int i = (b)-1; i >= a; i--)
#define F0Rd(i,a) for (int i = (a)-1; i >= 0; i--)
#define sz(x) (int)(x).size()
#define mp make_pair
#define pb push_back
#define f first
#define s second
#define lb lower_bound
#define ub upper_bound
#define all(x) x.begin(), x.end()
const int MOD = 1000000007;
const ll INF = 1e18;
const int MX = 300001;
int n, r[MX], ans[MX];
pi p[MX];
set<pi> S;
map<int,int> m; vi rm;
vpi tmp[MX];
ll sq(ll x) { return x*x; }
ll dist(pi a, pi b) { a.f -= b.f, a.s -= b.s; return (ll)a.f*a.f+(ll)a.s*a.s; }
bool inter(int a, int b) { return dist(p[a],p[b]) <= sq(r[a]+r[b]); }
void init() {
ios_base::sync_with_stdio(0); cin.tie(0);
cin >> n;
FOR(i,1,n+1) {
cin >> p[i].f >> p[i].s >> r[i];
p[i].f += 1000000000; p[i].s += 1000000000;
S.insert({r[i],-i});
}
}
ll crad = INF;
void del(int x, int e) { ans[x] = e; S.erase({r[x],-x}); }
pi getPos(int x) { return {p[x].f/crad,p[x].s/crad}; }
void recalc() {
F0R(i,sz(m)) tmp[i].clear();
m.clear(), rm.clear();
FOR(i,1,n+1) if (!ans[i]) m[getPos(i).f] = 0;
for (auto& a: m) { a.s = sz(rm); rm.pb(a.f); }
FOR(i,1,n+1) if (!ans[i]) {
pi P = getPos(i);
tmp[m[P.f]].pb({P.s,i});
}
}
void solve() {
pi x = *S.rbegin(); x.s *= -1;
if (2*x.f <= crad) { crad = x.f; recalc(); }
pi pos = getPos(x.s);
for (int t = lb(all(rm),pos.f-2)-rm.begin();
t < sz(rm) && abs(rm[t]-pos.f) <= 2; t ++)
for (int t2 = lb(all(tmp[t]),mp(pos.s-2,-MOD))-tmp[t].begin();
t2 < sz(tmp[t]) && abs(tmp[t][t2].f-pos.s) <= 2; t2++)
if (!ans[tmp[t][t2].s] && inter(x.s,tmp[t][t2].s))
del(tmp[t][t2].s,x.s);
}
int main() {
init();
while (sz(S)) solve();
FOR(i,1,n+1) cout << ans[i] << " ";
}
/* Look for:
* the exact constraints (multiple sets are too slow for n=10^6 :( )
* special cases (n=1?)
* overflow (ll vs int?)
* array bounds
*/
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
8 ms |
7416 KB |
Output is correct |
2 |
Correct |
8 ms |
7528 KB |
Output is correct |
3 |
Correct |
9 ms |
7528 KB |
Output is correct |
4 |
Correct |
8 ms |
7528 KB |
Output is correct |
5 |
Correct |
9 ms |
7528 KB |
Output is correct |
6 |
Correct |
8 ms |
7528 KB |
Output is correct |
7 |
Correct |
9 ms |
7528 KB |
Output is correct |
8 |
Correct |
8 ms |
7528 KB |
Output is correct |
9 |
Correct |
11 ms |
7712 KB |
Output is correct |
10 |
Correct |
9 ms |
7712 KB |
Output is correct |
11 |
Correct |
10 ms |
7724 KB |
Output is correct |
12 |
Correct |
9 ms |
7748 KB |
Output is correct |
13 |
Correct |
8 ms |
7748 KB |
Output is correct |
14 |
Execution timed out |
3041 ms |
7752 KB |
Time limit exceeded |
15 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
782 ms |
37932 KB |
Output is correct |
2 |
Correct |
716 ms |
44764 KB |
Output is correct |
3 |
Correct |
721 ms |
51072 KB |
Output is correct |
4 |
Correct |
676 ms |
58516 KB |
Output is correct |
5 |
Correct |
731 ms |
64080 KB |
Output is correct |
6 |
Correct |
1167 ms |
77680 KB |
Output is correct |
7 |
Correct |
926 ms |
77680 KB |
Output is correct |
8 |
Correct |
1062 ms |
80956 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Execution timed out |
3042 ms |
80956 KB |
Time limit exceeded |
2 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Execution timed out |
3036 ms |
87848 KB |
Time limit exceeded |
2 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
8 ms |
7416 KB |
Output is correct |
2 |
Correct |
8 ms |
7528 KB |
Output is correct |
3 |
Correct |
9 ms |
7528 KB |
Output is correct |
4 |
Correct |
8 ms |
7528 KB |
Output is correct |
5 |
Correct |
9 ms |
7528 KB |
Output is correct |
6 |
Correct |
8 ms |
7528 KB |
Output is correct |
7 |
Correct |
9 ms |
7528 KB |
Output is correct |
8 |
Correct |
8 ms |
7528 KB |
Output is correct |
9 |
Correct |
11 ms |
7712 KB |
Output is correct |
10 |
Correct |
9 ms |
7712 KB |
Output is correct |
11 |
Correct |
10 ms |
7724 KB |
Output is correct |
12 |
Correct |
9 ms |
7748 KB |
Output is correct |
13 |
Correct |
8 ms |
7748 KB |
Output is correct |
14 |
Execution timed out |
3041 ms |
7752 KB |
Time limit exceeded |
15 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
8 ms |
7416 KB |
Output is correct |
2 |
Correct |
8 ms |
7528 KB |
Output is correct |
3 |
Correct |
9 ms |
7528 KB |
Output is correct |
4 |
Correct |
8 ms |
7528 KB |
Output is correct |
5 |
Correct |
9 ms |
7528 KB |
Output is correct |
6 |
Correct |
8 ms |
7528 KB |
Output is correct |
7 |
Correct |
9 ms |
7528 KB |
Output is correct |
8 |
Correct |
8 ms |
7528 KB |
Output is correct |
9 |
Correct |
11 ms |
7712 KB |
Output is correct |
10 |
Correct |
9 ms |
7712 KB |
Output is correct |
11 |
Correct |
10 ms |
7724 KB |
Output is correct |
12 |
Correct |
9 ms |
7748 KB |
Output is correct |
13 |
Correct |
8 ms |
7748 KB |
Output is correct |
14 |
Execution timed out |
3041 ms |
7752 KB |
Time limit exceeded |
15 |
Halted |
0 ms |
0 KB |
- |