#include "bits/stdc++.h"
#define st first
#define nd second
using lint = int64_t;
constexpr int MOD = int(1e9) + 7;
constexpr int INF = 0x3f3f3f3f;
constexpr int NINF = 0xcfcfcfcf;
constexpr lint LINF = 0x3f3f3f3f3f3f3f3f;
const long double PI = acosl(-1.0);
// Returns -1 if a < b, 0 if a = b and 1 if a > b.
int cmp_double(double a, double b = 0, double eps = 1e-9) {
return a + eps > b ? b + eps > a ? 0 : 1 : -1;
}
using namespace std;
using ld = long double;
bool inCircle(lint x, lint y, lint dx, lint dy, lint r, lint t){
return ( (x + dx*t)*(x + dx*t) + (y + dy*t)*(y + dy*t) ) <= r*r;
}
pair<bool, lint> firstInCircle(lint x, lint y, lint dx, lint dy, lint r, ld dbt){
lint s = lint(dbt) - 3;
lint e = s + 6;
for(lint t = s; t<=e; t++)
if(inCircle(x,y,dx,dy,r, t)) return {true,t};
return {false, -1};
}
pair<bool, lint> lastInCircle(lint x, lint y, lint dx, lint dy, lint r, ld dbt){
lint s = lint(dbt) + 3;
lint e = s - 6;
for(lint t = s; t>=e; t--)
if(inCircle(x,y,dx,dy,r, t)) return {true,t};
return {false, -1};
}
vector<lint> intersection(lint x, lint y, lint dx, lint dy, lint r){
vector<lint> answer;
lint a = (dx*dx + dy*dy);
lint b = (2ll*x*dx + 2ll*y*dy);
lint c = -(r*r - x*x - y*y);
lint delta = b*b - 4ll*a*c;
if(delta < 0) return answer;
if(delta == 0){
ld t = ld(-b) / ld(2*a);
auto left = firstInCircle(x,y,dx,dy,r,t);
if(!left.first) return {};
auto right = lastInCircle(x,y,dx,dy,r,t);
if(!right.first) return {};
answer = {left.second, right.second};
}else{
auto left = firstInCircle(x,y,dx,dy,r, ( ld(-b) - sqrtl((ld)delta))/ ld(2*a) );
if(!left.first) return {};
auto right = lastInCircle(x,y,dx,dy,r, (ld(-b) + sqrtl((ld)delta))/ ld(2*a) );
if(!right.first) return {};
answer = {left.second, right.second};
}
if(answer[0] < 0) answer[0] = 0;
if(answer[0] > answer[1]) answer = {};
else if(answer.back() < 0) answer = {};
return answer;
}
vector<lint> comprTime;
int main() {
cin.tie(nullptr)->sync_with_stdio(false);
comprTime = {-1};
int n;
lint r; cin>>n>>r;
vector<vector<lint>> stars;
for(int i=0; i<n; i++){
lint x1, y1, x2, y2;
cin>>x1>>y1>>x2>>y2;
auto v = intersection(x1, y1, x2-x1, y2-y1, r);
if(!v.empty()){
stars.push_back(v);
comprTime.push_back(v[0]);
comprTime.push_back(v[1]);
}
}
sort(comprTime.begin(), comprTime.end());
comprTime.resize(distance(
comprTime.begin(),
unique(comprTime.begin(), comprTime.end())
));
comprTime.push_back(comprTime.back()+1);
vector<lint> pref(comprTime.size());
for(auto star : stars){
int l = lower_bound(comprTime.begin(),comprTime.end(), star[0]) - comprTime.begin();
int r = upper_bound(comprTime.begin(),comprTime.end(), star[1]) - comprTime.begin();
pref[l]++;
pref[r]--;
}
for(int i=1; i<int(comprTime.size()); i++)
pref[i] += pref[i-1];
cout<<*max_element(pref.begin(), pref.end())<<"\n";
return 0;
}
/*
[ ]Leu o problema certo???
[ ]Ver se precisa de long long
[ ]Viu o limite dos fors (é n? é m?)
[ ]Tamanho do vetor, será que é 2e5 em vez de 1e5??
[ ]Testar sample
[ ]Testar casos de borda
[ ]1LL no 1LL << i
[ ]Testar mod (é 1e9+7, mesmo?, será que o mod não ficou negativo?)
*/
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
0 ms |
212 KB |
Output is correct |
2 |
Correct |
0 ms |
212 KB |
Output is correct |
3 |
Correct |
1 ms |
340 KB |
Output is correct |
4 |
Correct |
1 ms |
340 KB |
Output is correct |
5 |
Correct |
1 ms |
340 KB |
Output is correct |
6 |
Correct |
1 ms |
340 KB |
Output is correct |
7 |
Correct |
1 ms |
340 KB |
Output is correct |
8 |
Correct |
0 ms |
212 KB |
Output is correct |
9 |
Correct |
0 ms |
212 KB |
Output is correct |
10 |
Correct |
0 ms |
212 KB |
Output is correct |
11 |
Correct |
0 ms |
212 KB |
Output is correct |
12 |
Correct |
0 ms |
212 KB |
Output is correct |
13 |
Correct |
0 ms |
212 KB |
Output is correct |
14 |
Correct |
1 ms |
340 KB |
Output is correct |
15 |
Correct |
1 ms |
340 KB |
Output is correct |
16 |
Correct |
1 ms |
340 KB |
Output is correct |
17 |
Correct |
1 ms |
340 KB |
Output is correct |
18 |
Correct |
1 ms |
340 KB |
Output is correct |
19 |
Correct |
1 ms |
340 KB |
Output is correct |
20 |
Correct |
1 ms |
340 KB |
Output is correct |
21 |
Correct |
1 ms |
340 KB |
Output is correct |
22 |
Correct |
1 ms |
340 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
0 ms |
212 KB |
Output is correct |
2 |
Correct |
0 ms |
212 KB |
Output is correct |
3 |
Correct |
1 ms |
340 KB |
Output is correct |
4 |
Correct |
1 ms |
340 KB |
Output is correct |
5 |
Correct |
1 ms |
340 KB |
Output is correct |
6 |
Correct |
1 ms |
340 KB |
Output is correct |
7 |
Correct |
1 ms |
340 KB |
Output is correct |
8 |
Correct |
0 ms |
212 KB |
Output is correct |
9 |
Correct |
0 ms |
212 KB |
Output is correct |
10 |
Correct |
0 ms |
212 KB |
Output is correct |
11 |
Correct |
0 ms |
212 KB |
Output is correct |
12 |
Correct |
0 ms |
212 KB |
Output is correct |
13 |
Correct |
0 ms |
212 KB |
Output is correct |
14 |
Correct |
1 ms |
340 KB |
Output is correct |
15 |
Correct |
1 ms |
340 KB |
Output is correct |
16 |
Correct |
1 ms |
340 KB |
Output is correct |
17 |
Correct |
1 ms |
340 KB |
Output is correct |
18 |
Correct |
1 ms |
340 KB |
Output is correct |
19 |
Correct |
1 ms |
340 KB |
Output is correct |
20 |
Correct |
1 ms |
340 KB |
Output is correct |
21 |
Correct |
1 ms |
340 KB |
Output is correct |
22 |
Correct |
1 ms |
340 KB |
Output is correct |
23 |
Correct |
47 ms |
6968 KB |
Output is correct |
24 |
Correct |
46 ms |
7104 KB |
Output is correct |
25 |
Correct |
48 ms |
7084 KB |
Output is correct |
26 |
Correct |
48 ms |
7008 KB |
Output is correct |
27 |
Correct |
39 ms |
3796 KB |
Output is correct |
28 |
Correct |
43 ms |
6460 KB |
Output is correct |
29 |
Correct |
44 ms |
6440 KB |
Output is correct |
30 |
Correct |
44 ms |
6440 KB |
Output is correct |
31 |
Correct |
45 ms |
6448 KB |
Output is correct |
32 |
Correct |
60 ms |
7064 KB |
Output is correct |
33 |
Correct |
59 ms |
7028 KB |
Output is correct |
34 |
Correct |
59 ms |
7000 KB |
Output is correct |
35 |
Correct |
59 ms |
7096 KB |
Output is correct |
36 |
Correct |
59 ms |
7096 KB |
Output is correct |
37 |
Correct |
40 ms |
4792 KB |
Output is correct |
38 |
Correct |
41 ms |
4792 KB |
Output is correct |
39 |
Correct |
37 ms |
3792 KB |
Output is correct |
40 |
Correct |
39 ms |
4104 KB |
Output is correct |
41 |
Correct |
36 ms |
3756 KB |
Output is correct |
42 |
Correct |
38 ms |
3756 KB |
Output is correct |
43 |
Correct |
144 ms |
14544 KB |
Output is correct |
44 |
Correct |
148 ms |
14540 KB |
Output is correct |
45 |
Correct |
143 ms |
14664 KB |
Output is correct |
46 |
Correct |
148 ms |
14444 KB |
Output is correct |
47 |
Correct |
141 ms |
14444 KB |
Output is correct |
48 |
Correct |
122 ms |
13860 KB |
Output is correct |
49 |
Correct |
118 ms |
13880 KB |
Output is correct |
50 |
Correct |
124 ms |
13792 KB |
Output is correct |
51 |
Correct |
120 ms |
13832 KB |
Output is correct |
52 |
Correct |
120 ms |
13844 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
0 ms |
212 KB |
Output is correct |
2 |
Correct |
0 ms |
212 KB |
Output is correct |
3 |
Correct |
1 ms |
340 KB |
Output is correct |
4 |
Correct |
1 ms |
340 KB |
Output is correct |
5 |
Correct |
1 ms |
340 KB |
Output is correct |
6 |
Correct |
1 ms |
340 KB |
Output is correct |
7 |
Correct |
1 ms |
340 KB |
Output is correct |
8 |
Correct |
0 ms |
212 KB |
Output is correct |
9 |
Correct |
0 ms |
212 KB |
Output is correct |
10 |
Correct |
0 ms |
212 KB |
Output is correct |
11 |
Correct |
0 ms |
212 KB |
Output is correct |
12 |
Correct |
0 ms |
212 KB |
Output is correct |
13 |
Correct |
0 ms |
212 KB |
Output is correct |
14 |
Correct |
1 ms |
340 KB |
Output is correct |
15 |
Correct |
1 ms |
340 KB |
Output is correct |
16 |
Correct |
1 ms |
340 KB |
Output is correct |
17 |
Correct |
1 ms |
340 KB |
Output is correct |
18 |
Correct |
1 ms |
340 KB |
Output is correct |
19 |
Correct |
1 ms |
340 KB |
Output is correct |
20 |
Correct |
1 ms |
340 KB |
Output is correct |
21 |
Correct |
1 ms |
340 KB |
Output is correct |
22 |
Correct |
1 ms |
340 KB |
Output is correct |
23 |
Correct |
0 ms |
212 KB |
Output is correct |
24 |
Correct |
1 ms |
340 KB |
Output is correct |
25 |
Correct |
1 ms |
340 KB |
Output is correct |
26 |
Correct |
1 ms |
340 KB |
Output is correct |
27 |
Correct |
1 ms |
340 KB |
Output is correct |
28 |
Correct |
1 ms |
340 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
0 ms |
212 KB |
Output is correct |
2 |
Correct |
0 ms |
212 KB |
Output is correct |
3 |
Correct |
1 ms |
340 KB |
Output is correct |
4 |
Correct |
1 ms |
340 KB |
Output is correct |
5 |
Correct |
1 ms |
340 KB |
Output is correct |
6 |
Correct |
1 ms |
340 KB |
Output is correct |
7 |
Correct |
1 ms |
340 KB |
Output is correct |
8 |
Correct |
0 ms |
212 KB |
Output is correct |
9 |
Correct |
0 ms |
212 KB |
Output is correct |
10 |
Correct |
0 ms |
212 KB |
Output is correct |
11 |
Correct |
0 ms |
212 KB |
Output is correct |
12 |
Correct |
0 ms |
212 KB |
Output is correct |
13 |
Correct |
0 ms |
212 KB |
Output is correct |
14 |
Correct |
1 ms |
340 KB |
Output is correct |
15 |
Correct |
1 ms |
340 KB |
Output is correct |
16 |
Correct |
1 ms |
340 KB |
Output is correct |
17 |
Correct |
1 ms |
340 KB |
Output is correct |
18 |
Correct |
1 ms |
340 KB |
Output is correct |
19 |
Correct |
1 ms |
340 KB |
Output is correct |
20 |
Correct |
1 ms |
340 KB |
Output is correct |
21 |
Correct |
1 ms |
340 KB |
Output is correct |
22 |
Correct |
1 ms |
340 KB |
Output is correct |
23 |
Correct |
47 ms |
6968 KB |
Output is correct |
24 |
Correct |
46 ms |
7104 KB |
Output is correct |
25 |
Correct |
48 ms |
7084 KB |
Output is correct |
26 |
Correct |
48 ms |
7008 KB |
Output is correct |
27 |
Correct |
39 ms |
3796 KB |
Output is correct |
28 |
Correct |
43 ms |
6460 KB |
Output is correct |
29 |
Correct |
44 ms |
6440 KB |
Output is correct |
30 |
Correct |
44 ms |
6440 KB |
Output is correct |
31 |
Correct |
45 ms |
6448 KB |
Output is correct |
32 |
Correct |
60 ms |
7064 KB |
Output is correct |
33 |
Correct |
59 ms |
7028 KB |
Output is correct |
34 |
Correct |
59 ms |
7000 KB |
Output is correct |
35 |
Correct |
59 ms |
7096 KB |
Output is correct |
36 |
Correct |
59 ms |
7096 KB |
Output is correct |
37 |
Correct |
40 ms |
4792 KB |
Output is correct |
38 |
Correct |
41 ms |
4792 KB |
Output is correct |
39 |
Correct |
37 ms |
3792 KB |
Output is correct |
40 |
Correct |
39 ms |
4104 KB |
Output is correct |
41 |
Correct |
36 ms |
3756 KB |
Output is correct |
42 |
Correct |
38 ms |
3756 KB |
Output is correct |
43 |
Correct |
144 ms |
14544 KB |
Output is correct |
44 |
Correct |
148 ms |
14540 KB |
Output is correct |
45 |
Correct |
143 ms |
14664 KB |
Output is correct |
46 |
Correct |
148 ms |
14444 KB |
Output is correct |
47 |
Correct |
141 ms |
14444 KB |
Output is correct |
48 |
Correct |
122 ms |
13860 KB |
Output is correct |
49 |
Correct |
118 ms |
13880 KB |
Output is correct |
50 |
Correct |
124 ms |
13792 KB |
Output is correct |
51 |
Correct |
120 ms |
13832 KB |
Output is correct |
52 |
Correct |
120 ms |
13844 KB |
Output is correct |
53 |
Correct |
0 ms |
212 KB |
Output is correct |
54 |
Correct |
1 ms |
340 KB |
Output is correct |
55 |
Correct |
1 ms |
340 KB |
Output is correct |
56 |
Correct |
1 ms |
340 KB |
Output is correct |
57 |
Correct |
1 ms |
340 KB |
Output is correct |
58 |
Correct |
1 ms |
340 KB |
Output is correct |
59 |
Correct |
47 ms |
7024 KB |
Output is correct |
60 |
Correct |
44 ms |
3508 KB |
Output is correct |
61 |
Correct |
51 ms |
3820 KB |
Output is correct |
62 |
Correct |
43 ms |
4156 KB |
Output is correct |
63 |
Correct |
55 ms |
3640 KB |
Output is correct |
64 |
Correct |
52 ms |
3464 KB |
Output is correct |
65 |
Correct |
52 ms |
3508 KB |
Output is correct |
66 |
Correct |
53 ms |
3784 KB |
Output is correct |
67 |
Correct |
50 ms |
3432 KB |
Output is correct |
68 |
Correct |
64 ms |
7080 KB |
Output is correct |
69 |
Correct |
63 ms |
7084 KB |
Output is correct |
70 |
Correct |
64 ms |
7040 KB |
Output is correct |
71 |
Correct |
63 ms |
7028 KB |
Output is correct |
72 |
Correct |
65 ms |
7128 KB |
Output is correct |
73 |
Correct |
44 ms |
4544 KB |
Output is correct |
74 |
Correct |
45 ms |
6412 KB |
Output is correct |
75 |
Correct |
45 ms |
6388 KB |
Output is correct |
76 |
Correct |
40 ms |
4012 KB |
Output is correct |
77 |
Correct |
139 ms |
13072 KB |
Output is correct |
78 |
Correct |
138 ms |
12952 KB |
Output is correct |
79 |
Correct |
141 ms |
13104 KB |
Output is correct |
80 |
Correct |
154 ms |
12932 KB |
Output is correct |
81 |
Correct |
144 ms |
12948 KB |
Output is correct |
82 |
Correct |
125 ms |
13828 KB |
Output is correct |
83 |
Correct |
129 ms |
13840 KB |
Output is correct |
84 |
Correct |
131 ms |
13816 KB |
Output is correct |
85 |
Correct |
124 ms |
13872 KB |
Output is correct |
86 |
Correct |
129 ms |
13956 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
0 ms |
212 KB |
Output is correct |
2 |
Correct |
0 ms |
212 KB |
Output is correct |
3 |
Correct |
1 ms |
340 KB |
Output is correct |
4 |
Correct |
1 ms |
340 KB |
Output is correct |
5 |
Correct |
1 ms |
340 KB |
Output is correct |
6 |
Correct |
1 ms |
340 KB |
Output is correct |
7 |
Correct |
1 ms |
340 KB |
Output is correct |
8 |
Correct |
0 ms |
212 KB |
Output is correct |
9 |
Correct |
0 ms |
212 KB |
Output is correct |
10 |
Correct |
0 ms |
212 KB |
Output is correct |
11 |
Correct |
0 ms |
212 KB |
Output is correct |
12 |
Correct |
0 ms |
212 KB |
Output is correct |
13 |
Correct |
0 ms |
212 KB |
Output is correct |
14 |
Correct |
1 ms |
340 KB |
Output is correct |
15 |
Correct |
1 ms |
340 KB |
Output is correct |
16 |
Correct |
1 ms |
340 KB |
Output is correct |
17 |
Correct |
1 ms |
340 KB |
Output is correct |
18 |
Correct |
1 ms |
340 KB |
Output is correct |
19 |
Correct |
1 ms |
340 KB |
Output is correct |
20 |
Correct |
1 ms |
340 KB |
Output is correct |
21 |
Correct |
1 ms |
340 KB |
Output is correct |
22 |
Correct |
1 ms |
340 KB |
Output is correct |
23 |
Correct |
47 ms |
6968 KB |
Output is correct |
24 |
Correct |
46 ms |
7104 KB |
Output is correct |
25 |
Correct |
48 ms |
7084 KB |
Output is correct |
26 |
Correct |
48 ms |
7008 KB |
Output is correct |
27 |
Correct |
39 ms |
3796 KB |
Output is correct |
28 |
Correct |
43 ms |
6460 KB |
Output is correct |
29 |
Correct |
44 ms |
6440 KB |
Output is correct |
30 |
Correct |
44 ms |
6440 KB |
Output is correct |
31 |
Correct |
45 ms |
6448 KB |
Output is correct |
32 |
Correct |
60 ms |
7064 KB |
Output is correct |
33 |
Correct |
59 ms |
7028 KB |
Output is correct |
34 |
Correct |
59 ms |
7000 KB |
Output is correct |
35 |
Correct |
59 ms |
7096 KB |
Output is correct |
36 |
Correct |
59 ms |
7096 KB |
Output is correct |
37 |
Correct |
40 ms |
4792 KB |
Output is correct |
38 |
Correct |
41 ms |
4792 KB |
Output is correct |
39 |
Correct |
37 ms |
3792 KB |
Output is correct |
40 |
Correct |
39 ms |
4104 KB |
Output is correct |
41 |
Correct |
36 ms |
3756 KB |
Output is correct |
42 |
Correct |
38 ms |
3756 KB |
Output is correct |
43 |
Correct |
144 ms |
14544 KB |
Output is correct |
44 |
Correct |
148 ms |
14540 KB |
Output is correct |
45 |
Correct |
143 ms |
14664 KB |
Output is correct |
46 |
Correct |
148 ms |
14444 KB |
Output is correct |
47 |
Correct |
141 ms |
14444 KB |
Output is correct |
48 |
Correct |
122 ms |
13860 KB |
Output is correct |
49 |
Correct |
118 ms |
13880 KB |
Output is correct |
50 |
Correct |
124 ms |
13792 KB |
Output is correct |
51 |
Correct |
120 ms |
13832 KB |
Output is correct |
52 |
Correct |
120 ms |
13844 KB |
Output is correct |
53 |
Correct |
47 ms |
7024 KB |
Output is correct |
54 |
Correct |
44 ms |
3508 KB |
Output is correct |
55 |
Correct |
51 ms |
3820 KB |
Output is correct |
56 |
Correct |
43 ms |
4156 KB |
Output is correct |
57 |
Correct |
55 ms |
3640 KB |
Output is correct |
58 |
Correct |
52 ms |
3464 KB |
Output is correct |
59 |
Correct |
52 ms |
3508 KB |
Output is correct |
60 |
Correct |
53 ms |
3784 KB |
Output is correct |
61 |
Correct |
50 ms |
3432 KB |
Output is correct |
62 |
Correct |
64 ms |
7080 KB |
Output is correct |
63 |
Correct |
63 ms |
7084 KB |
Output is correct |
64 |
Correct |
64 ms |
7040 KB |
Output is correct |
65 |
Correct |
63 ms |
7028 KB |
Output is correct |
66 |
Correct |
65 ms |
7128 KB |
Output is correct |
67 |
Correct |
44 ms |
4544 KB |
Output is correct |
68 |
Correct |
45 ms |
6412 KB |
Output is correct |
69 |
Correct |
45 ms |
6388 KB |
Output is correct |
70 |
Correct |
40 ms |
4012 KB |
Output is correct |
71 |
Correct |
139 ms |
13072 KB |
Output is correct |
72 |
Correct |
138 ms |
12952 KB |
Output is correct |
73 |
Correct |
141 ms |
13104 KB |
Output is correct |
74 |
Correct |
154 ms |
12932 KB |
Output is correct |
75 |
Correct |
144 ms |
12948 KB |
Output is correct |
76 |
Correct |
125 ms |
13828 KB |
Output is correct |
77 |
Correct |
129 ms |
13840 KB |
Output is correct |
78 |
Correct |
131 ms |
13816 KB |
Output is correct |
79 |
Correct |
124 ms |
13872 KB |
Output is correct |
80 |
Correct |
129 ms |
13956 KB |
Output is correct |
81 |
Correct |
0 ms |
212 KB |
Output is correct |
82 |
Correct |
1 ms |
340 KB |
Output is correct |
83 |
Correct |
1 ms |
340 KB |
Output is correct |
84 |
Correct |
1 ms |
340 KB |
Output is correct |
85 |
Correct |
1 ms |
340 KB |
Output is correct |
86 |
Correct |
1 ms |
340 KB |
Output is correct |
87 |
Incorrect |
61 ms |
3384 KB |
Output isn't correct |
88 |
Halted |
0 ms |
0 KB |
- |