oj.uz

Submission #66129

# Submission time Handle Problem Language Result Execution time Memory
66129 2018-08-09T18:44:17 Z Benq Circle selection (APIO18_circle_selection) C++11
12 / 100
 3000 ms 87848 KB 
#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
*/

Subtask #1
0 / 7.0

# Verdict Execution time Memory 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 -

Subtask #2
12.0 / 12.0

# Verdict Execution time Memory 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

Subtask #3
0 / 15.0

# Verdict Execution time Memory Grader output
1 Execution timed out 3042 ms 80956 KB Time limit exceeded
2 Halted 0 ms 0 KB -

Subtask #4
0 / 23.0

# Verdict Execution time Memory Grader output
1 Execution timed out 3036 ms 87848 KB Time limit exceeded
2 Halted 0 ms 0 KB -

Subtask #5
0 / 30.0

# Verdict Execution time Memory 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 -

Subtask #6
0 / 13.0

# Verdict Execution time Memory 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 -