Submission #107302

# Submission time Handle Problem Language Result Execution time Memory
107302 2019-04-23T07:47:46 Z leduykhongngu Circle selection (APIO18_circle_selection) C++11
72 / 100
3000 ms 66064 KB
#include <bits/stdc++.h>
#pragma GCC optimize ("Ofast")
#define FOR(i,a,b) for (int i=(a),_b=(b); i<=_b; ++i)
#define FORR(i,a,b) for (int i=(a),_b=(b); i>=_b; --i)
#define REP(i,b) for (int i=0,_b=(b); i<_b; ++i)
#define endl '\n'
#define sz(x) (int) x.size()
#define mod %
#define fillchar(x,y,z) memset(x,z,y)
#define pii pair<int,int>
#define fi first
#define se second
#define mp make_pair
typedef long long int64;
typedef unsigned long long qword;
using namespace std;
int ReadInt()
{
    char c;
    for (c = getchar(); c!='-'&&(c < '0' || c > '9'); c = getchar());
    int nega=0;
    int ans = c - '0';
    if (c=='-') nega=1,ans=0;
    for (c = getchar(); c >= '0' && c <= '9'; c = getchar())
        ans = ans * 10 + c - '0';
    if (nega) return -ans;
    return ans;
}
void WriteInt(int x)
{
    if (x > 9) WriteInt(x / 10);
    putchar(x % 10 + '0');
}
const int maxn=3e5+5;
inline int64 sqr(int64 x)
{
    return x*x;
}
struct Circle{
    int64 x,y,r;
    int id;

    bool IsIntersect(const Circle &c)
    {
        return sqr(x-c.x)+sqr(y-c.y)<=sqr(r+c.r);
    }
}cir[maxn];
bool operator <(const Circle &b,const Circle &c)
    {
        return mp(-b.r,b.id)<mp(-c.r,c.id);
    }
vector<set<pii> > Point;
vector<int> X;
int n;
int ans[maxn];
pii tmp[maxn];
int main()
{
    #ifdef meomeomeooooo
        freopen("input.txt","r",stdin);
        //freopen(".txt","w",stdout);
    #endif // meomeomeooooo
    iostream::sync_with_stdio(false);
    cin.tie(0);
    n=ReadInt();
    set<Circle> myset;
    FOR(i,1,n) {
        //cin >> cir[i].x >> cir[i].y >> cir[i].r;
        cir[i].x=ReadInt(); cir[i].y=ReadInt(); cir[i].r=ReadInt();
        X.push_back(cir[i].x);
        cir[i].id=i;
        myset.insert(cir[i]);
    }
    sort(X.begin(),X.end());
    X.resize(unique(X.begin(),X.end())-X.begin());
    #define GetDec(x) (lower_bound(X.begin(),X.end(),x)-X.begin())
    Point.resize(X.size());
    FOR(i,1,n) {
        int x=GetDec(cir[i].x);
        Point[x].insert({cir[i].y,cir[i].id});
    }
    while (myset.size()) {
        Circle cur=*myset.begin(); myset.erase(cur);
        ans[cur.id]=cur.id;
        int x=cur.x,y=cur.y,r=cur.r;
        int Lx=max(1ll*x-2*r,-1000000000ll);
        int Rx=min(1ll*x+2*r,1000000000ll);
        Lx=GetDec(Lx);
        int Ly=max(1ll*y-2*r,-1000000000ll);
        int Ry=min(1ll*y+2*r,1000000000ll);

        for (int i=Lx; i<X.size()&&X[i]<=Rx; ++i) {
            int top=0;
            int64 meomeo=sqrt(4ll*r*r-sqr(X[i]-x));
//            int Ly=max(1ll*y-meomeo,-1000000000ll);
//            int Ry=min(1ll*y+meomeo,1000000000ll);
            for (auto j=Point[i].lower_bound(mp(Ly,0)); j!=Point[i].end()&&(*j).fi<=Ry; ++j) {
                int id=(*j).se;
                if (cur.IsIntersect(cir[id])) {
                    assert(ans[id]==0||ans[id]==cur.id);
                    ans[id]=cur.id;
                    myset.erase(cir[id]);
                    tmp[++top]=*j;
                }
            }
            FOR(j,1,top)
                Point[i].erase(tmp[j]);
        }
    }
    FOR(i,1,n) {
        WriteInt(ans[i]); putchar(' ');
    }
    return 0;
}

Compilation message

circle_selection.cpp: In function 'int main()':
circle_selection.cpp:92:25: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
         for (int i=Lx; i<X.size()&&X[i]<=Rx; ++i) {
                        ~^~~~~~~~~
circle_selection.cpp:94:19: warning: unused variable 'meomeo' [-Wunused-variable]
             int64 meomeo=sqrt(4ll*r*r-sqr(X[i]-x));
                   ^~~~~~
# Verdict Execution time Memory Grader output
1 Correct 3 ms 384 KB Output is correct
2 Correct 2 ms 384 KB Output is correct
3 Correct 2 ms 384 KB Output is correct
4 Correct 2 ms 384 KB Output is correct
5 Correct 2 ms 384 KB Output is correct
6 Correct 2 ms 384 KB Output is correct
7 Correct 2 ms 384 KB Output is correct
8 Correct 2 ms 384 KB Output is correct
9 Correct 2 ms 384 KB Output is correct
10 Correct 2 ms 384 KB Output is correct
11 Correct 2 ms 384 KB Output is correct
12 Correct 2 ms 384 KB Output is correct
13 Correct 3 ms 384 KB Output is correct
14 Correct 2 ms 380 KB Output is correct
15 Correct 2 ms 384 KB Output is correct
16 Correct 4 ms 512 KB Output is correct
17 Correct 3 ms 612 KB Output is correct
18 Correct 3 ms 512 KB Output is correct
19 Correct 12 ms 1536 KB Output is correct
20 Correct 9 ms 1380 KB Output is correct
21 Correct 8 ms 1408 KB Output is correct
22 Correct 10 ms 1408 KB Output is correct
23 Correct 12 ms 1408 KB Output is correct
24 Correct 8 ms 1408 KB Output is correct
25 Correct 9 ms 1408 KB Output is correct
26 Correct 10 ms 1408 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 1247 ms 66064 KB Output is correct
2 Correct 1304 ms 66024 KB Output is correct
3 Correct 1261 ms 65844 KB Output is correct
4 Correct 1231 ms 66000 KB Output is correct
5 Correct 1149 ms 61000 KB Output is correct
6 Correct 1062 ms 65752 KB Output is correct
7 Correct 1184 ms 64852 KB Output is correct
8 Correct 1068 ms 65380 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 3 ms 512 KB Output is correct
2 Correct 437 ms 22096 KB Output is correct
3 Correct 1610 ms 65864 KB Output is correct
4 Correct 1662 ms 65860 KB Output is correct
5 Correct 2307 ms 64072 KB Output is correct
6 Correct 779 ms 34412 KB Output is correct
7 Correct 294 ms 18284 KB Output is correct
8 Correct 39 ms 4088 KB Output is correct
9 Correct 1848 ms 65864 KB Output is correct
10 Execution timed out 3026 ms 63920 KB Time limit exceeded
# Verdict Execution time Memory Grader output
1 Correct 806 ms 65752 KB Output is correct
2 Correct 576 ms 65740 KB Output is correct
3 Correct 2271 ms 65344 KB Output is correct
4 Correct 605 ms 65788 KB Output is correct
5 Correct 599 ms 65864 KB Output is correct
6 Correct 2268 ms 65124 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 3 ms 384 KB Output is correct
2 Correct 2 ms 384 KB Output is correct
3 Correct 2 ms 384 KB Output is correct
4 Correct 2 ms 384 KB Output is correct
5 Correct 2 ms 384 KB Output is correct
6 Correct 2 ms 384 KB Output is correct
7 Correct 2 ms 384 KB Output is correct
8 Correct 2 ms 384 KB Output is correct
9 Correct 2 ms 384 KB Output is correct
10 Correct 2 ms 384 KB Output is correct
11 Correct 2 ms 384 KB Output is correct
12 Correct 2 ms 384 KB Output is correct
13 Correct 3 ms 384 KB Output is correct
14 Correct 2 ms 380 KB Output is correct
15 Correct 2 ms 384 KB Output is correct
16 Correct 4 ms 512 KB Output is correct
17 Correct 3 ms 612 KB Output is correct
18 Correct 3 ms 512 KB Output is correct
19 Correct 12 ms 1536 KB Output is correct
20 Correct 9 ms 1380 KB Output is correct
21 Correct 8 ms 1408 KB Output is correct
22 Correct 10 ms 1408 KB Output is correct
23 Correct 12 ms 1408 KB Output is correct
24 Correct 8 ms 1408 KB Output is correct
25 Correct 9 ms 1408 KB Output is correct
26 Correct 10 ms 1408 KB Output is correct
27 Correct 19 ms 2532 KB Output is correct
28 Correct 17 ms 2560 KB Output is correct
29 Correct 19 ms 2688 KB Output is correct
30 Correct 17 ms 2560 KB Output is correct
31 Correct 22 ms 2560 KB Output is correct
32 Correct 16 ms 2560 KB Output is correct
33 Correct 325 ms 22128 KB Output is correct
34 Correct 307 ms 22128 KB Output is correct
35 Correct 325 ms 22128 KB Output is correct
36 Correct 392 ms 22128 KB Output is correct
37 Correct 323 ms 22304 KB Output is correct
38 Correct 343 ms 22308 KB Output is correct
39 Correct 1265 ms 17512 KB Output is correct
40 Correct 1290 ms 17520 KB Output is correct
41 Correct 1327 ms 17520 KB Output is correct
42 Correct 309 ms 17448 KB Output is correct
43 Correct 199 ms 22128 KB Output is correct
44 Correct 187 ms 22128 KB Output is correct
45 Correct 182 ms 22128 KB Output is correct
46 Correct 183 ms 22080 KB Output is correct
47 Correct 178 ms 22128 KB Output is correct
48 Correct 175 ms 22100 KB Output is correct
49 Correct 182 ms 22128 KB Output is correct
50 Correct 215 ms 22172 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 3 ms 384 KB Output is correct
2 Correct 2 ms 384 KB Output is correct
3 Correct 2 ms 384 KB Output is correct
4 Correct 2 ms 384 KB Output is correct
5 Correct 2 ms 384 KB Output is correct
6 Correct 2 ms 384 KB Output is correct
7 Correct 2 ms 384 KB Output is correct
8 Correct 2 ms 384 KB Output is correct
9 Correct 2 ms 384 KB Output is correct
10 Correct 2 ms 384 KB Output is correct
11 Correct 2 ms 384 KB Output is correct
12 Correct 2 ms 384 KB Output is correct
13 Correct 3 ms 384 KB Output is correct
14 Correct 2 ms 380 KB Output is correct
15 Correct 2 ms 384 KB Output is correct
16 Correct 4 ms 512 KB Output is correct
17 Correct 3 ms 612 KB Output is correct
18 Correct 3 ms 512 KB Output is correct
19 Correct 12 ms 1536 KB Output is correct
20 Correct 9 ms 1380 KB Output is correct
21 Correct 8 ms 1408 KB Output is correct
22 Correct 10 ms 1408 KB Output is correct
23 Correct 12 ms 1408 KB Output is correct
24 Correct 8 ms 1408 KB Output is correct
25 Correct 9 ms 1408 KB Output is correct
26 Correct 10 ms 1408 KB Output is correct
27 Correct 1247 ms 66064 KB Output is correct
28 Correct 1304 ms 66024 KB Output is correct
29 Correct 1261 ms 65844 KB Output is correct
30 Correct 1231 ms 66000 KB Output is correct
31 Correct 1149 ms 61000 KB Output is correct
32 Correct 1062 ms 65752 KB Output is correct
33 Correct 1184 ms 64852 KB Output is correct
34 Correct 1068 ms 65380 KB Output is correct
35 Correct 3 ms 512 KB Output is correct
36 Correct 437 ms 22096 KB Output is correct
37 Correct 1610 ms 65864 KB Output is correct
38 Correct 1662 ms 65860 KB Output is correct
39 Correct 2307 ms 64072 KB Output is correct
40 Correct 779 ms 34412 KB Output is correct
41 Correct 294 ms 18284 KB Output is correct
42 Correct 39 ms 4088 KB Output is correct
43 Correct 1848 ms 65864 KB Output is correct
44 Execution timed out 3026 ms 63920 KB Time limit exceeded