Submission #51087

# Submission time Handle Problem Language Result Execution time Memory
51087 2018-06-16T03:29:58 Z edisonhello Circle selection (APIO18_circle_selection) C++11
100 / 100
1793 ms 42740 KB
#pragma GCC optimize("Ofast,unroll-loops")
#pragma GCC target("sse,sse2,sse3,ssse3,sse4,avx,tune=native")
#include<cstdio>
#include<unordered_map>
#include<vector>
#include<utility>
#include<algorithm>
using namespace std;
 
int n,pos[300005],ans[300005];
struct circle{
    int x,y,r,i;
} c[300005];
 
namespace std {
    template<> struct hash<pair<int,int>> {
        long long operator()(const pair<int,int> &a) const {
            return 1ll*a.first*1000000009+a.second;
        }
    };
}
 
unordered_map<pair<int,int>,int> mp;
vector<int> idxs[300005];
 
int rit(){
    char c; int x=0; bool mi=0;
    while((c=getchar_unlocked())){ if(c=='-' || (c>='0' && c<='9'))break; }
    if(c=='-')mi=1;
    else x=c-'0';
    while((c=getchar_unlocked())){ if(c<'0' || c>'9')break; x=x*10+(c&15); }
    return mi?-x:x;
}
 
int main(){
    n=rit();
    for(int i=0;i<n;++i){
        c[i].x=rit(),c[i].y=rit(),c[i].r=rit();
        c[i].i=i+1;
    }
    sort(c,c+n,[](const circle &a,const circle &b){ return a.r==b.r?a.i<b.i:a.r>b.r; });
    // for(int i=0;i<n;++i)pos[c[i].i]=i;
    int masterPtr=0;
    int D=c[0].r,killcount=0;
 
    auto buildGraph=[&](){
        for(int i=1;i<=int(mp.size());++i)idxs[i].clear();
        mp.clear();
        D=c[masterPtr].r;
        killcount=0;
        for(int i=masterPtr;i<n;++i){
            if(c[i].i==0)continue;
            pair<int,int> pos(c[i].x/D,c[i].y/D);
            int &it=mp[pos];
            if(it)idxs[it].push_back(i);
            else idxs[(it=mp.size())].push_back(i);
            // cout<<"insert "<<i<<" into "<<it<<" , pos: "<<pos.first<<" , "<<pos.second<<endl;
        }
    };
 
    buildGraph();
    while(masterPtr<n){
        while(masterPtr<n && c[masterPtr].i==0)++masterPtr;
        if(masterPtr>=n)break;
        // cout<<"now masterPtr: "<<masterPtr<<" circle is "<<c[masterPtr].i<<endl;
        // cout<<"circle position: "<<c[masterPtr].x<<" "<<c[masterPtr].y<<endl;
        if(D/c[masterPtr].r>=5 && killcount>=1000){
            buildGraph();
            // cout<<"map contains: "<<endl;
            // for(auto p:mp)cout<<"[("<<p.first.first<<" , "<<p.first.second<<") : "<<p.second<<"] ";
            // cout<<endl;
        }
        pair<int,int> pos(c[masterPtr].x/D,c[masterPtr].y/D);
        // cout<<"base position: "<<pos.first<<" "<<pos.second<<endl;
        for(int dx=-2;dx<=2;++dx){
            for(int dy=-2;dy<=2;++dy){
                auto it=mp.find(make_pair(pos.first+dx,pos.second+dy));
                if(it==mp.end())continue;
                // cout<<"it's not the end: "<<it->second<<endl;
                for(int i:idxs[it->second]){
                    if(c[i].i==0)continue;
                    if(1ll*(c[i].x-c[masterPtr].x)*(c[i].x-c[masterPtr].x)+1ll*(c[i].y-c[masterPtr].y)*(c[i].y-c[masterPtr].y)>1ll*(c[masterPtr].r+c[i].r)*(c[masterPtr].r+c[i].r))continue;
                    ans[c[i].i]=c[masterPtr].i;
                    // cout<<"killed "<<c[i].i<<", killer: "<<c[masterPtr].i<<endl;
                    
                    if(c[i].i!=c[masterPtr].i)c[i].i=0;
                    ++killcount;
                }
            }
        }
        c[masterPtr].i=0;
    }
    for(int i=1;i<=n;++i){
        printf("%d ",ans[i]);
    }
}
# Verdict Execution time Memory Grader output
1 Correct 9 ms 7292 KB Output is correct
2 Correct 7 ms 7396 KB Output is correct
3 Correct 8 ms 7528 KB Output is correct
4 Correct 8 ms 7588 KB Output is correct
5 Correct 8 ms 7588 KB Output is correct
6 Correct 8 ms 7588 KB Output is correct
7 Correct 8 ms 7592 KB Output is correct
8 Correct 9 ms 7592 KB Output is correct
9 Correct 9 ms 7592 KB Output is correct
10 Correct 7 ms 7592 KB Output is correct
11 Correct 9 ms 7608 KB Output is correct
12 Correct 8 ms 7608 KB Output is correct
13 Correct 8 ms 7780 KB Output is correct
14 Correct 9 ms 7780 KB Output is correct
15 Correct 9 ms 7780 KB Output is correct
16 Correct 9 ms 7780 KB Output is correct
17 Correct 9 ms 7780 KB Output is correct
18 Correct 8 ms 7780 KB Output is correct
19 Correct 11 ms 7808 KB Output is correct
20 Correct 10 ms 7808 KB Output is correct
21 Correct 11 ms 7820 KB Output is correct
22 Correct 17 ms 8204 KB Output is correct
23 Correct 18 ms 8220 KB Output is correct
24 Correct 18 ms 8220 KB Output is correct
25 Correct 21 ms 8236 KB Output is correct
26 Correct 17 ms 8252 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 167 ms 16936 KB Output is correct
2 Correct 177 ms 17320 KB Output is correct
3 Correct 149 ms 17320 KB Output is correct
4 Correct 132 ms 17320 KB Output is correct
5 Correct 173 ms 17980 KB Output is correct
6 Correct 810 ms 28872 KB Output is correct
7 Correct 248 ms 28872 KB Output is correct
8 Correct 446 ms 28872 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 8 ms 28872 KB Output is correct
2 Correct 494 ms 28872 KB Output is correct
3 Correct 1793 ms 42516 KB Output is correct
4 Correct 1520 ms 42516 KB Output is correct
5 Correct 1282 ms 42516 KB Output is correct
6 Correct 380 ms 42516 KB Output is correct
7 Correct 249 ms 42516 KB Output is correct
8 Correct 37 ms 42516 KB Output is correct
9 Correct 1618 ms 42516 KB Output is correct
10 Correct 1310 ms 42516 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 1711 ms 42604 KB Output is correct
2 Correct 1730 ms 42604 KB Output is correct
3 Correct 512 ms 42604 KB Output is correct
4 Correct 1572 ms 42604 KB Output is correct
5 Correct 1573 ms 42604 KB Output is correct
6 Correct 213 ms 42604 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 9 ms 7292 KB Output is correct
2 Correct 7 ms 7396 KB Output is correct
3 Correct 8 ms 7528 KB Output is correct
4 Correct 8 ms 7588 KB Output is correct
5 Correct 8 ms 7588 KB Output is correct
6 Correct 8 ms 7588 KB Output is correct
7 Correct 8 ms 7592 KB Output is correct
8 Correct 9 ms 7592 KB Output is correct
9 Correct 9 ms 7592 KB Output is correct
10 Correct 7 ms 7592 KB Output is correct
11 Correct 9 ms 7608 KB Output is correct
12 Correct 8 ms 7608 KB Output is correct
13 Correct 8 ms 7780 KB Output is correct
14 Correct 9 ms 7780 KB Output is correct
15 Correct 9 ms 7780 KB Output is correct
16 Correct 9 ms 7780 KB Output is correct
17 Correct 9 ms 7780 KB Output is correct
18 Correct 8 ms 7780 KB Output is correct
19 Correct 11 ms 7808 KB Output is correct
20 Correct 10 ms 7808 KB Output is correct
21 Correct 11 ms 7820 KB Output is correct
22 Correct 17 ms 8204 KB Output is correct
23 Correct 18 ms 8220 KB Output is correct
24 Correct 18 ms 8220 KB Output is correct
25 Correct 21 ms 8236 KB Output is correct
26 Correct 17 ms 8252 KB Output is correct
27 Correct 18 ms 42604 KB Output is correct
28 Correct 12 ms 42604 KB Output is correct
29 Correct 12 ms 42604 KB Output is correct
30 Correct 36 ms 42604 KB Output is correct
31 Correct 27 ms 42604 KB Output is correct
32 Correct 27 ms 42604 KB Output is correct
33 Correct 57 ms 42604 KB Output is correct
34 Correct 54 ms 42604 KB Output is correct
35 Correct 58 ms 42604 KB Output is correct
36 Correct 379 ms 42604 KB Output is correct
37 Correct 412 ms 42604 KB Output is correct
38 Correct 425 ms 42604 KB Output is correct
39 Correct 301 ms 42604 KB Output is correct
40 Correct 302 ms 42604 KB Output is correct
41 Correct 308 ms 42604 KB Output is correct
42 Correct 309 ms 42604 KB Output is correct
43 Correct 375 ms 42604 KB Output is correct
44 Correct 391 ms 42604 KB Output is correct
45 Correct 413 ms 42604 KB Output is correct
46 Correct 471 ms 42604 KB Output is correct
47 Correct 445 ms 42604 KB Output is correct
48 Correct 357 ms 42604 KB Output is correct
49 Correct 438 ms 42604 KB Output is correct
50 Correct 400 ms 42604 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 9 ms 7292 KB Output is correct
2 Correct 7 ms 7396 KB Output is correct
3 Correct 8 ms 7528 KB Output is correct
4 Correct 8 ms 7588 KB Output is correct
5 Correct 8 ms 7588 KB Output is correct
6 Correct 8 ms 7588 KB Output is correct
7 Correct 8 ms 7592 KB Output is correct
8 Correct 9 ms 7592 KB Output is correct
9 Correct 9 ms 7592 KB Output is correct
10 Correct 7 ms 7592 KB Output is correct
11 Correct 9 ms 7608 KB Output is correct
12 Correct 8 ms 7608 KB Output is correct
13 Correct 8 ms 7780 KB Output is correct
14 Correct 9 ms 7780 KB Output is correct
15 Correct 9 ms 7780 KB Output is correct
16 Correct 9 ms 7780 KB Output is correct
17 Correct 9 ms 7780 KB Output is correct
18 Correct 8 ms 7780 KB Output is correct
19 Correct 11 ms 7808 KB Output is correct
20 Correct 10 ms 7808 KB Output is correct
21 Correct 11 ms 7820 KB Output is correct
22 Correct 17 ms 8204 KB Output is correct
23 Correct 18 ms 8220 KB Output is correct
24 Correct 18 ms 8220 KB Output is correct
25 Correct 21 ms 8236 KB Output is correct
26 Correct 17 ms 8252 KB Output is correct
27 Correct 167 ms 16936 KB Output is correct
28 Correct 177 ms 17320 KB Output is correct
29 Correct 149 ms 17320 KB Output is correct
30 Correct 132 ms 17320 KB Output is correct
31 Correct 173 ms 17980 KB Output is correct
32 Correct 810 ms 28872 KB Output is correct
33 Correct 248 ms 28872 KB Output is correct
34 Correct 446 ms 28872 KB Output is correct
35 Correct 8 ms 28872 KB Output is correct
36 Correct 494 ms 28872 KB Output is correct
37 Correct 1793 ms 42516 KB Output is correct
38 Correct 1520 ms 42516 KB Output is correct
39 Correct 1282 ms 42516 KB Output is correct
40 Correct 380 ms 42516 KB Output is correct
41 Correct 249 ms 42516 KB Output is correct
42 Correct 37 ms 42516 KB Output is correct
43 Correct 1618 ms 42516 KB Output is correct
44 Correct 1310 ms 42516 KB Output is correct
45 Correct 1711 ms 42604 KB Output is correct
46 Correct 1730 ms 42604 KB Output is correct
47 Correct 512 ms 42604 KB Output is correct
48 Correct 1572 ms 42604 KB Output is correct
49 Correct 1573 ms 42604 KB Output is correct
50 Correct 213 ms 42604 KB Output is correct
51 Correct 18 ms 42604 KB Output is correct
52 Correct 12 ms 42604 KB Output is correct
53 Correct 12 ms 42604 KB Output is correct
54 Correct 36 ms 42604 KB Output is correct
55 Correct 27 ms 42604 KB Output is correct
56 Correct 27 ms 42604 KB Output is correct
57 Correct 57 ms 42604 KB Output is correct
58 Correct 54 ms 42604 KB Output is correct
59 Correct 58 ms 42604 KB Output is correct
60 Correct 379 ms 42604 KB Output is correct
61 Correct 412 ms 42604 KB Output is correct
62 Correct 425 ms 42604 KB Output is correct
63 Correct 301 ms 42604 KB Output is correct
64 Correct 302 ms 42604 KB Output is correct
65 Correct 308 ms 42604 KB Output is correct
66 Correct 309 ms 42604 KB Output is correct
67 Correct 375 ms 42604 KB Output is correct
68 Correct 391 ms 42604 KB Output is correct
69 Correct 413 ms 42604 KB Output is correct
70 Correct 471 ms 42604 KB Output is correct
71 Correct 445 ms 42604 KB Output is correct
72 Correct 357 ms 42604 KB Output is correct
73 Correct 438 ms 42604 KB Output is correct
74 Correct 400 ms 42604 KB Output is correct
75 Correct 155 ms 42604 KB Output is correct
76 Correct 142 ms 42604 KB Output is correct
77 Correct 144 ms 42604 KB Output is correct
78 Correct 137 ms 42604 KB Output is correct
79 Correct 223 ms 42604 KB Output is correct
80 Correct 169 ms 42604 KB Output is correct
81 Correct 1531 ms 42604 KB Output is correct
82 Correct 1585 ms 42604 KB Output is correct
83 Correct 1662 ms 42604 KB Output is correct
84 Correct 1554 ms 42604 KB Output is correct
85 Correct 1710 ms 42604 KB Output is correct
86 Correct 1657 ms 42604 KB Output is correct
87 Correct 1422 ms 42740 KB Output is correct
88 Correct 1220 ms 42740 KB Output is correct
89 Correct 1198 ms 42740 KB Output is correct
90 Correct 1283 ms 42740 KB Output is correct
91 Correct 1124 ms 42740 KB Output is correct
92 Correct 1082 ms 42740 KB Output is correct
93 Correct 1461 ms 42740 KB Output is correct
94 Correct 1641 ms 42740 KB Output is correct
95 Correct 1436 ms 42740 KB Output is correct
96 Correct 1598 ms 42740 KB Output is correct
97 Correct 1523 ms 42740 KB Output is correct
98 Correct 952 ms 42740 KB Output is correct
99 Correct 1746 ms 42740 KB Output is correct
100 Correct 1783 ms 42740 KB Output is correct
101 Correct 1145 ms 42740 KB Output is correct
102 Correct 1499 ms 42740 KB Output is correct
103 Correct 1642 ms 42740 KB Output is correct
104 Correct 1664 ms 42740 KB Output is correct
105 Correct 928 ms 42740 KB Output is correct
106 Correct 1290 ms 42740 KB Output is correct
107 Correct 1322 ms 42740 KB Output is correct
108 Correct 1269 ms 42740 KB Output is correct
109 Correct 1335 ms 42740 KB Output is correct
110 Correct 1380 ms 42740 KB Output is correct
111 Correct 1322 ms 42740 KB Output is correct
112 Correct 1260 ms 42740 KB Output is correct
113 Correct 1306 ms 42740 KB Output is correct
114 Correct 1330 ms 42740 KB Output is correct
115 Correct 1405 ms 42740 KB Output is correct
116 Correct 1274 ms 42740 KB Output is correct