Submission #107299

# Submission time Handle Problem Language Result Execution time Memory
107299 2019-04-23T07:42:55 Z leduykhongngu Circle selection (APIO18_circle_selection) C++11
72 / 100
3000 ms 66020 KB
#include <bits/stdc++.h>
#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);
        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);
            int Ly=max(1ll*y-2*r,-1000000000ll);
            int Ry=min(1ll*y+2*r,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:88: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:90: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 2 ms 384 KB Output is correct
2 Correct 2 ms 384 KB Output is correct
3 Correct 3 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 3 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 3 ms 384 KB Output is correct
11 Correct 3 ms 384 KB Output is correct
12 Correct 2 ms 384 KB Output is correct
13 Correct 2 ms 384 KB Output is correct
14 Correct 3 ms 384 KB Output is correct
15 Correct 2 ms 384 KB Output is correct
16 Correct 3 ms 512 KB Output is correct
17 Correct 3 ms 640 KB Output is correct
18 Correct 4 ms 512 KB Output is correct
19 Correct 8 ms 1408 KB Output is correct
20 Correct 9 ms 1408 KB Output is correct
21 Correct 10 ms 1452 KB Output is correct
22 Correct 9 ms 1536 KB Output is correct
23 Correct 15 ms 1408 KB Output is correct
24 Correct 12 ms 1536 KB Output is correct
25 Correct 9 ms 1536 KB Output is correct
26 Correct 8 ms 1536 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 1297 ms 65848 KB Output is correct
2 Correct 1372 ms 65920 KB Output is correct
3 Correct 1348 ms 65612 KB Output is correct
4 Correct 1246 ms 66020 KB Output is correct
5 Correct 1044 ms 60924 KB Output is correct
6 Correct 1083 ms 65764 KB Output is correct
7 Correct 1101 ms 64936 KB Output is correct
8 Correct 1047 ms 65480 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 2 ms 384 KB Output is correct
2 Correct 476 ms 22072 KB Output is correct
3 Correct 1577 ms 65892 KB Output is correct
4 Correct 1566 ms 65860 KB Output is correct
5 Correct 2283 ms 64100 KB Output is correct
6 Correct 854 ms 34536 KB Output is correct
7 Correct 299 ms 18316 KB Output is correct
8 Correct 29 ms 3968 KB Output is correct
9 Correct 1732 ms 65832 KB Output is correct
10 Execution timed out 3026 ms 63944 KB Time limit exceeded
# Verdict Execution time Memory Grader output
1 Correct 813 ms 65760 KB Output is correct
2 Correct 560 ms 65892 KB Output is correct
3 Correct 2451 ms 65512 KB Output is correct
4 Correct 670 ms 65636 KB Output is correct
5 Correct 563 ms 65896 KB Output is correct
6 Correct 2276 ms 65280 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 2 ms 384 KB Output is correct
2 Correct 2 ms 384 KB Output is correct
3 Correct 3 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 3 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 3 ms 384 KB Output is correct
11 Correct 3 ms 384 KB Output is correct
12 Correct 2 ms 384 KB Output is correct
13 Correct 2 ms 384 KB Output is correct
14 Correct 3 ms 384 KB Output is correct
15 Correct 2 ms 384 KB Output is correct
16 Correct 3 ms 512 KB Output is correct
17 Correct 3 ms 640 KB Output is correct
18 Correct 4 ms 512 KB Output is correct
19 Correct 8 ms 1408 KB Output is correct
20 Correct 9 ms 1408 KB Output is correct
21 Correct 10 ms 1452 KB Output is correct
22 Correct 9 ms 1536 KB Output is correct
23 Correct 15 ms 1408 KB Output is correct
24 Correct 12 ms 1536 KB Output is correct
25 Correct 9 ms 1536 KB Output is correct
26 Correct 8 ms 1536 KB Output is correct
27 Correct 16 ms 2560 KB Output is correct
28 Correct 18 ms 2552 KB Output is correct
29 Correct 20 ms 2560 KB Output is correct
30 Correct 21 ms 2560 KB Output is correct
31 Correct 20 ms 2552 KB Output is correct
32 Correct 15 ms 2560 KB Output is correct
33 Correct 348 ms 22064 KB Output is correct
34 Correct 315 ms 22136 KB Output is correct
35 Correct 333 ms 22200 KB Output is correct
36 Correct 328 ms 22132 KB Output is correct
37 Correct 319 ms 22144 KB Output is correct
38 Correct 308 ms 22208 KB Output is correct
39 Correct 1307 ms 17476 KB Output is correct
40 Correct 1008 ms 17392 KB Output is correct
41 Correct 969 ms 17500 KB Output is correct
42 Correct 248 ms 17340 KB Output is correct
43 Correct 170 ms 22128 KB Output is correct
44 Correct 189 ms 22088 KB Output is correct
45 Correct 195 ms 22256 KB Output is correct
46 Correct 183 ms 22156 KB Output is correct
47 Correct 185 ms 22128 KB Output is correct
48 Correct 183 ms 22172 KB Output is correct
49 Correct 204 ms 22132 KB Output is correct
50 Correct 194 ms 22252 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 2 ms 384 KB Output is correct
2 Correct 2 ms 384 KB Output is correct
3 Correct 3 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 3 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 3 ms 384 KB Output is correct
11 Correct 3 ms 384 KB Output is correct
12 Correct 2 ms 384 KB Output is correct
13 Correct 2 ms 384 KB Output is correct
14 Correct 3 ms 384 KB Output is correct
15 Correct 2 ms 384 KB Output is correct
16 Correct 3 ms 512 KB Output is correct
17 Correct 3 ms 640 KB Output is correct
18 Correct 4 ms 512 KB Output is correct
19 Correct 8 ms 1408 KB Output is correct
20 Correct 9 ms 1408 KB Output is correct
21 Correct 10 ms 1452 KB Output is correct
22 Correct 9 ms 1536 KB Output is correct
23 Correct 15 ms 1408 KB Output is correct
24 Correct 12 ms 1536 KB Output is correct
25 Correct 9 ms 1536 KB Output is correct
26 Correct 8 ms 1536 KB Output is correct
27 Correct 1297 ms 65848 KB Output is correct
28 Correct 1372 ms 65920 KB Output is correct
29 Correct 1348 ms 65612 KB Output is correct
30 Correct 1246 ms 66020 KB Output is correct
31 Correct 1044 ms 60924 KB Output is correct
32 Correct 1083 ms 65764 KB Output is correct
33 Correct 1101 ms 64936 KB Output is correct
34 Correct 1047 ms 65480 KB Output is correct
35 Correct 2 ms 384 KB Output is correct
36 Correct 476 ms 22072 KB Output is correct
37 Correct 1577 ms 65892 KB Output is correct
38 Correct 1566 ms 65860 KB Output is correct
39 Correct 2283 ms 64100 KB Output is correct
40 Correct 854 ms 34536 KB Output is correct
41 Correct 299 ms 18316 KB Output is correct
42 Correct 29 ms 3968 KB Output is correct
43 Correct 1732 ms 65832 KB Output is correct
44 Execution timed out 3026 ms 63944 KB Time limit exceeded