#include <bits/stdc++.h>
#define int long long
#define double long double
using namespace std;
const int NMAX=2e5+5;
int n, r;
map<int, int> v;
/*
Delta=4R^2( (a_i-c_i)^2+(b_i-d_i)^2 ) -4(a_i_d_i-b_ic_i)^2
k_{1,2}=(-2(a_ic_i+b_id_i-a_i^2-b_i^2) \pm sqrt(Delta) ) / (2( (a_i-c_i)^2+(b_i-d_i)^2 ) )
*/
pair<int, int> secs(int a, int b, int c, int d)
{
    int dlt=r*r*((a-c)*(a-c)+(b-d)*(b-d))-(a*d-b*c)*(a*d-b*c);
    if(dlt<-1e-9)
        return {-1, -1};
    double k1=(a*a+b*b-a*c-b*d-sqrtl(dlt))/((a-c)*(a-c)+(b-d)*(b-d)),
           k2=(a*a+b*b-a*c-b*d+sqrtl(dlt))/((a-c)*(a-c)+(b-d)*(b-d));
    if(max(k1, k2)<-1e-9)
        return {-1, -1};
    return {ceil(max(min(k1, k2), (double)0)), floor(max(max(k1, k2), (double)0))};
}
signed main()
{
    cin>>n>>r;
    for(int i=1;i<=n;i++)
    {
        int a, b, c, d;
        cin>>a>>b>>c>>d;
        pair<int, int> sus=secs(a, b, c, d);
        if(sus.second!=-1)
        {
            v[sus.first]++;
            v[sus.second+1]--;
        }
    }
    int ans=0, cnt=0;
    for(auto x:v)
    {
        cnt+=x.second;
        ans=max(ans, cnt);
    }
    cout<<ans;
    return 0;
}
| # | Verdict  | Execution time | Memory | Grader output | 
|---|
| Fetching results... | 
| # | Verdict  | Execution time | Memory | Grader output | 
|---|
| Fetching results... | 
| # | Verdict  | Execution time | Memory | Grader output | 
|---|
| Fetching results... | 
| # | Verdict  | Execution time | Memory | Grader output | 
|---|
| Fetching results... | 
| # | Verdict  | Execution time | Memory | Grader output | 
|---|
| Fetching results... |