oj.uz

Submission #272729

# Submission time Handle Problem Language Result Execution time Memory
272729 2020-08-18T14:15:41 Z MKopchev Circle selection (APIO18_circle_selection) C++14
87 / 100
 3000 ms 377948 KB 
#include<bits/stdc++.h>
using namespace std;
const int nmax=3e5+42,inf=1e9+42,MX=(1<<20)+42;

struct circle
{
    int x,y,r,id;
};

int n;
circle inp[nmax],sorted_inp[nmax],other_sorted[nmax];

bool touch(circle a,circle b)
{
    long long d=1LL*(a.x-b.x)*(a.x-b.x)+1LL*(a.y-b.y)*(a.y-b.y);

    long long mx=a.r+b.r;

    mx=mx*mx;

    return d<=mx;
}

bool cmp(circle a,circle b)
{
    if(a.r!=b.r)return a.r>b.r;
    return a.id<b.id;
}

bool cmp_y(circle a,circle b)
{
    if(a.y!=b.y)return a.y<b.y;
    return a.id<b.id;
}
int outp[nmax];

vector< pair<int/*y*/,int/*id*/> > nodes[MX];
vector<int> prv[MX];
vector<int> nxt[MX];
vector<int> tree[MX];

void update(int node,int l,int r,int pos,circle cur)
{
    //cout<<node<<" -> "<<cur.y<<" "<<cur.id<<endl;
    nodes[node].push_back({cur.y,cur.id});

    if(l==r)return;

    int av=(l+r)/2;

    if(pos<=av)update(node*2,l,av,pos,cur);
    else update(node*2+1,av+1,r,pos,cur);
}
void build(int where,int node,int l,int r)
{
    while(tree[where].size()<=node)tree[where].push_back(0);
    tree[where][node]=r-l+1;

    if(l==r)return;

    int av=(l+r)/2;

    build(where,node*2,l,av);
    build(where,node*2+1,av+1,r);
}
void sub_tree(int where,int node,int l,int r,int pos)
{
    while(1)
    {
        tree[where][node]--;

        if(l==r)return;

        int av=(l+r)/2;

        if(pos<=av)node=node*2,r=av;
        else node=node*2+1,l=av+1;
    }
}

int faster(int where,int node,int l,int r)
{
    while(l!=r)
    {
        int av=(l+r)/2;
        if(tree[where][node*2])
        {
            node=node*2;
            r=av;
        }
        else
        {
            node=node*2+1;
            l=av+1;
        }
    }
    return l;
}
int query(int where,int node,int l,int r,int lq,int rq)
{
    if(tree[where][node]==0)return -1;

    if(l==lq&&r==rq)return faster(where,node,l,r);

    int av=(l+r)/2;

    if(lq<=av)
    {
        int mem=query(where,node*2,l,av,lq,min(av,rq));
        if(mem!=-1)return mem;
    }
    if(av<rq)
    {
        int mem=query(where,node*2+1,av+1,r,max(av+1,lq),rq);
        if(mem!=-1)return mem;
    }
    return -1;
}
int find_first(int where,int start)
{
    return query(where,1,0,nodes[where].size()-1,start,nodes[where].size()-1);
}

int pick(int a,int b)
{
    long long sum=a;
    sum=sum+b;

    if(inf<sum)sum=inf;
    if(-inf>sum)sum=-inf;

    return sum;
}
void match(int node,int l,int r,int lq,int rq,circle cur)
{
    if(tree[node].size()==0)return;

    if(tree[node][1]==0)return;

    if(l==lq&&r==rq)
    {
        int y1=pick(cur.y,-2*cur.r);
        int y2=pick(cur.y,2*cur.r);

        int start=lower_bound(nodes[node].begin(),nodes[node].end(),make_pair(y1,0))-nodes[node].begin();

        if(start==nodes[node].size())return;

        start=find_first(node,start);
        /*
        cout<<"circle "<<l<<" "<<r<<" "<<cur.id<<endl;

        for(auto k:nodes[node])
            cout<<k.first<<" "<<k.second<<"\t";cout<<endl;

        cout<<"y1= "<<y1<<" y2= "<<y2<<" start= "<<start<<endl;
        */
        if(start==-1)return;

        vector<int> to_pop={};

        while(start<nodes[node].size()&&nodes[node][start].first<=y2)
        {
            //cout<<"start= "<<start<<" "<<nodes[node][start].second<<" "<<inp[nodes[node][start].second].x<<" "<<inp[nodes[node][start].second].y<<" "<<inp[nodes[node][start].second].r<<endl;

            if(outp[nodes[node][start].second]==0&&touch(inp[nodes[node][start].second],cur))
            {
                //cout<<"touch "<<nodes[node][start].second<<" "<<cur.id<<" "<<endl;

                outp[nodes[node][start].second]=cur.id;
            }

            if(outp[nodes[node][start].second])to_pop.push_back(start);

            start=nxt[node][start];
        }

        for(auto k:to_pop)
        {
            if(0<=prv[node][k])
            {
                nxt[node][prv[node][k]]=nxt[node][k];
            }
            if(nxt[node][k]<nodes[node].size())
            {
                prv[node][nxt[node][k]]=prv[node][k];
            }

            sub_tree(node,1,0,nodes[node].size()-1,k);
        }

        return;
    }

    int av=(l+r)/2;

    if(lq<=av)match(node*2,l,av,lq,min(av,rq),cur);
    if(av<rq)match(node*2+1,av+1,r,max(av+1,lq),rq,cur);
}

int help_x[nmax];

void push(int which)
{
    outp[sorted_inp[which].id]=sorted_inp[which].id;

    int x1=pick(sorted_inp[which].x,-2*sorted_inp[which].r);
    int x2=pick(sorted_inp[which].x,+2*sorted_inp[which].r);

    int le=lower_bound(help_x+1,help_x+n+1,x1)-help_x;
    int ri=upper_bound(help_x+1,help_x+n+1,x2)-help_x;
    ri--;

    //cout<<"le= "<<le<<" ri= "<<ri<<endl;

    if(le>ri)return;

    match(1,1,n,le,ri,sorted_inp[which]);
}
signed main()
{
    scanf("%i",&n);
    for(int i=1;i<=n;i++)
    {
        scanf("%i%i%i",&inp[i].x,&inp[i].y,&inp[i].r);
        inp[i].id=i;
    }

    for(int i=1;i<=n;i++)sorted_inp[i]=inp[i],other_sorted[i]=inp[i];

    sort(other_sorted+1,other_sorted+n+1,cmp_y);

    sort(sorted_inp+1,sorted_inp+n+1,cmp);

    for(int i=1;i<=n;i++)help_x[i]=sorted_inp[i].x;

    sort(help_x+1,help_x+n+1);

    for(int i=1;i<=n;i++)
    {
        int pos=lower_bound(help_x+1,help_x+n+1,other_sorted[i].x)-help_x;

        update(1,1,n,pos,other_sorted[i]);
    }

    for(int i=1;i<MX;i++)
        if(nodes[i].size())
        {
            for(int j=0;j<nodes[i].size();j++)
            {
                prv[i].push_back(j-1);
                nxt[i].push_back(j+1);
            }
            build(i,1,0,nodes[i].size()-1);
        }

    for(int i=1;i<=n;i++)
        if(outp[sorted_inp[i].id]==0)
            push(i);

    for(int i=1;i<=n;i++)printf("%i ",outp[i]);printf("\n");
    return 0;
}

Compilation message

circle_selection.cpp: In function 'void build(int, int, int, int)':
circle_selection.cpp:56:29: warning: comparison of integer expressions of different signedness: 'std::vector<int>::size_type' {aka 'long unsigned int'} and 'int' [-Wsign-compare]
   56 |     while(tree[where].size()<=node)tree[where].push_back(0);
      |           ~~~~~~~~~~~~~~~~~~^~~~~~
circle_selection.cpp: In function 'void match(int, int, int, int, int, circle)':
circle_selection.cpp:147:17: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<std::pair<int, int> >::size_type' {aka 'long unsigned int'} [-Wsign-compare]
  147 |         if(start==nodes[node].size())return;
      |            ~~~~~^~~~~~~~~~~~~~~~~~~~
circle_selection.cpp:162:20: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<std::pair<int, int> >::size_type' {aka 'long unsigned int'} [-Wsign-compare]
  162 |         while(start<nodes[node].size()&&nodes[node][start].first<=y2)
      |               ~~~~~^~~~~~~~~~~~~~~~~~~
circle_selection.cpp:184:28: warning: comparison of integer expressions of different signedness: '__gnu_cxx::__alloc_traits<std::allocator<int>, int>::value_type' {aka 'int'} and 'std::vector<std::pair<int, int> >::size_type' {aka 'long unsigned int'} [-Wsign-compare]
  184 |             if(nxt[node][k]<nodes[node].size())
circle_selection.cpp: In function 'int main()':
circle_selection.cpp:249:26: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<std::pair<int, int> >::size_type' {aka 'long unsigned int'} [-Wsign-compare]
  249 |             for(int j=0;j<nodes[i].size();j++)
      |                         ~^~~~~~~~~~~~~~~~
circle_selection.cpp:261:5: warning: this 'for' clause does not guard... [-Wmisleading-indentation]
  261 |     for(int i=1;i<=n;i++)printf("%i ",outp[i]);printf("\n");
      |     ^~~
circle_selection.cpp:261:48: note: ...this statement, but the latter is misleadingly indented as if it were guarded by the 'for'
  261 |     for(int i=1;i<=n;i++)printf("%i ",outp[i]);printf("\n");
      |                                                ^~~~~~
circle_selection.cpp:222:10: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
  222 |     scanf("%i",&n);
      |     ~~~~~^~~~~~~~~
circle_selection.cpp:225:14: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
  225 |         scanf("%i%i%i",&inp[i].x,&inp[i].y,&inp[i].r);
      |         ~~~~~^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

Subtask #1
7.0 / 7.0

# Verdict Execution time Memory Grader output
1 Correct 69 ms 98808 KB Output is correct
2 Correct 69 ms 98808 KB Output is correct
3 Correct 76 ms 98936 KB Output is correct
4 Correct 70 ms 98808 KB Output is correct
5 Correct 75 ms 98808 KB Output is correct
6 Correct 70 ms 98936 KB Output is correct
7 Correct 71 ms 98936 KB Output is correct
8 Correct 71 ms 98936 KB Output is correct
9 Correct 69 ms 99064 KB Output is correct
10 Correct 76 ms 98936 KB Output is correct
11 Correct 69 ms 98936 KB Output is correct
12 Correct 68 ms 98936 KB Output is correct
13 Correct 71 ms 98868 KB Output is correct
14 Correct 68 ms 98936 KB Output is correct
15 Correct 71 ms 98936 KB Output is correct
16 Correct 72 ms 99320 KB Output is correct
17 Correct 73 ms 99336 KB Output is correct
18 Correct 72 ms 99320 KB Output is correct
19 Correct 92 ms 102520 KB Output is correct
20 Correct 84 ms 102520 KB Output is correct
21 Correct 86 ms 102396 KB Output is correct
22 Correct 92 ms 102392 KB Output is correct
23 Correct 92 ms 102392 KB Output is correct
24 Correct 90 ms 102392 KB Output is correct
25 Correct 87 ms 102392 KB Output is correct
26 Correct 90 ms 102372 KB Output is correct

Subtask #2
12.0 / 12.0

# Verdict Execution time Memory Grader output
1 Correct 1569 ms 376536 KB Output is correct
2 Correct 1630 ms 376912 KB Output is correct
3 Correct 1598 ms 377056 KB Output is correct
4 Correct 1655 ms 376028 KB Output is correct
5 Correct 1574 ms 360656 KB Output is correct
6 Correct 1783 ms 374860 KB Output is correct
7 Correct 1617 ms 372744 KB Output is correct
8 Correct 1664 ms 374156 KB Output is correct

Subtask #3
15.0 / 15.0

# Verdict Execution time Memory Grader output
1 Correct 73 ms 98808 KB Output is correct
2 Correct 776 ms 173792 KB Output is correct
3 Correct 2513 ms 375192 KB Output is correct
4 Correct 2533 ms 377144 KB Output is correct
5 Correct 2418 ms 369360 KB Output is correct
6 Correct 1136 ms 235356 KB Output is correct
7 Correct 542 ms 166500 KB Output is correct
8 Correct 139 ms 112304 KB Output is correct
9 Correct 2467 ms 373232 KB Output is correct
10 Correct 2421 ms 375248 KB Output is correct

Subtask #4
23.0 / 23.0

# Verdict Execution time Memory Grader output
1 Correct 2426 ms 377164 KB Output is correct
2 Correct 2172 ms 375128 KB Output is correct
3 Correct 1767 ms 375344 KB Output is correct
4 Correct 2228 ms 375784 KB Output is correct
5 Correct 2164 ms 375160 KB Output is correct
6 Correct 1766 ms 374612 KB Output is correct

Subtask #5
30.0 / 30.0

# Verdict Execution time Memory Grader output
1 Correct 69 ms 98808 KB Output is correct
2 Correct 69 ms 98808 KB Output is correct
3 Correct 76 ms 98936 KB Output is correct
4 Correct 70 ms 98808 KB Output is correct
5 Correct 75 ms 98808 KB Output is correct
6 Correct 70 ms 98936 KB Output is correct
7 Correct 71 ms 98936 KB Output is correct
8 Correct 71 ms 98936 KB Output is correct
9 Correct 69 ms 99064 KB Output is correct
10 Correct 76 ms 98936 KB Output is correct
11 Correct 69 ms 98936 KB Output is correct
12 Correct 68 ms 98936 KB Output is correct
13 Correct 71 ms 98868 KB Output is correct
14 Correct 68 ms 98936 KB Output is correct
15 Correct 71 ms 98936 KB Output is correct
16 Correct 72 ms 99320 KB Output is correct
17 Correct 73 ms 99336 KB Output is correct
18 Correct 72 ms 99320 KB Output is correct
19 Correct 92 ms 102520 KB Output is correct
20 Correct 84 ms 102520 KB Output is correct
21 Correct 86 ms 102396 KB Output is correct
22 Correct 92 ms 102392 KB Output is correct
23 Correct 92 ms 102392 KB Output is correct
24 Correct 90 ms 102392 KB Output is correct
25 Correct 87 ms 102392 KB Output is correct
26 Correct 90 ms 102372 KB Output is correct
27 Correct 103 ms 106344 KB Output is correct
28 Correct 101 ms 106232 KB Output is correct
29 Correct 103 ms 106232 KB Output is correct
30 Correct 116 ms 106104 KB Output is correct
31 Correct 110 ms 106232 KB Output is correct
32 Correct 115 ms 106236 KB Output is correct
33 Correct 506 ms 173924 KB Output is correct
34 Correct 521 ms 173928 KB Output is correct
35 Correct 520 ms 174220 KB Output is correct
36 Correct 689 ms 173536 KB Output is correct
37 Correct 704 ms 173628 KB Output is correct
38 Correct 722 ms 173608 KB Output is correct
39 Correct 570 ms 160260 KB Output is correct
40 Correct 593 ms 160348 KB Output is correct
41 Correct 563 ms 160016 KB Output is correct
42 Correct 469 ms 156808 KB Output is correct
43 Correct 561 ms 173540 KB Output is correct
44 Correct 627 ms 173540 KB Output is correct
45 Correct 568 ms 173580 KB Output is correct
46 Correct 573 ms 173540 KB Output is correct
47 Correct 604 ms 173540 KB Output is correct
48 Correct 641 ms 173576 KB Output is correct
49 Correct 576 ms 173540 KB Output is correct
50 Correct 574 ms 173612 KB Output is correct

Subtask #6
0 / 13.0

# Verdict Execution time Memory Grader output
1 Correct 69 ms 98808 KB Output is correct
2 Correct 69 ms 98808 KB Output is correct
3 Correct 76 ms 98936 KB Output is correct
4 Correct 70 ms 98808 KB Output is correct
5 Correct 75 ms 98808 KB Output is correct
6 Correct 70 ms 98936 KB Output is correct
7 Correct 71 ms 98936 KB Output is correct
8 Correct 71 ms 98936 KB Output is correct
9 Correct 69 ms 99064 KB Output is correct
10 Correct 76 ms 98936 KB Output is correct
11 Correct 69 ms 98936 KB Output is correct
12 Correct 68 ms 98936 KB Output is correct
13 Correct 71 ms 98868 KB Output is correct
14 Correct 68 ms 98936 KB Output is correct
15 Correct 71 ms 98936 KB Output is correct
16 Correct 72 ms 99320 KB Output is correct
17 Correct 73 ms 99336 KB Output is correct
18 Correct 72 ms 99320 KB Output is correct
19 Correct 92 ms 102520 KB Output is correct
20 Correct 84 ms 102520 KB Output is correct
21 Correct 86 ms 102396 KB Output is correct
22 Correct 92 ms 102392 KB Output is correct
23 Correct 92 ms 102392 KB Output is correct
24 Correct 90 ms 102392 KB Output is correct
25 Correct 87 ms 102392 KB Output is correct
26 Correct 90 ms 102372 KB Output is correct
27 Correct 1569 ms 376536 KB Output is correct
28 Correct 1630 ms 376912 KB Output is correct
29 Correct 1598 ms 377056 KB Output is correct
30 Correct 1655 ms 376028 KB Output is correct
31 Correct 1574 ms 360656 KB Output is correct
32 Correct 1783 ms 374860 KB Output is correct
33 Correct 1617 ms 372744 KB Output is correct
34 Correct 1664 ms 374156 KB Output is correct
35 Correct 73 ms 98808 KB Output is correct
36 Correct 776 ms 173792 KB Output is correct
37 Correct 2513 ms 375192 KB Output is correct
38 Correct 2533 ms 377144 KB Output is correct
39 Correct 2418 ms 369360 KB Output is correct
40 Correct 1136 ms 235356 KB Output is correct
41 Correct 542 ms 166500 KB Output is correct
42 Correct 139 ms 112304 KB Output is correct
43 Correct 2467 ms 373232 KB Output is correct
44 Correct 2421 ms 375248 KB Output is correct
45 Correct 2426 ms 377164 KB Output is correct
46 Correct 2172 ms 375128 KB Output is correct
47 Correct 1767 ms 375344 KB Output is correct
48 Correct 2228 ms 375784 KB Output is correct
49 Correct 2164 ms 375160 KB Output is correct
50 Correct 1766 ms 374612 KB Output is correct
51 Correct 103 ms 106344 KB Output is correct
52 Correct 101 ms 106232 KB Output is correct
53 Correct 103 ms 106232 KB Output is correct
54 Correct 116 ms 106104 KB Output is correct
55 Correct 110 ms 106232 KB Output is correct
56 Correct 115 ms 106236 KB Output is correct
57 Correct 506 ms 173924 KB Output is correct
58 Correct 521 ms 173928 KB Output is correct
59 Correct 520 ms 174220 KB Output is correct
60 Correct 689 ms 173536 KB Output is correct
61 Correct 704 ms 173628 KB Output is correct
62 Correct 722 ms 173608 KB Output is correct
63 Correct 570 ms 160260 KB Output is correct
64 Correct 593 ms 160348 KB Output is correct
65 Correct 563 ms 160016 KB Output is correct
66 Correct 469 ms 156808 KB Output is correct
67 Correct 561 ms 173540 KB Output is correct
68 Correct 627 ms 173540 KB Output is correct
69 Correct 568 ms 173580 KB Output is correct
70 Correct 573 ms 173540 KB Output is correct
71 Correct 604 ms 173540 KB Output is correct
72 Correct 641 ms 173576 KB Output is correct
73 Correct 576 ms 173540 KB Output is correct
74 Correct 574 ms 173612 KB Output is correct
75 Correct 1798 ms 377948 KB Output is correct
76 Correct 1811 ms 376984 KB Output is correct
77 Correct 1821 ms 376152 KB Output is correct
78 Correct 1748 ms 376500 KB Output is correct
79 Correct 1734 ms 375708 KB Output is correct
80 Correct 1732 ms 376260 KB Output is correct
81 Correct 2687 ms 375408 KB Output is correct
82 Correct 2558 ms 375356 KB Output is correct
83 Correct 2468 ms 375176 KB Output is correct
84 Correct 2834 ms 375220 KB Output is correct
85 Correct 2452 ms 375232 KB Output is correct
86 Correct 2400 ms 375492 KB Output is correct
87 Execution timed out 3123 ms 375096 KB Time limit exceeded
88 Halted 0 ms 0 KB -