2020牛客暑期多校第二场 B Boundary（计算几何）

最新推荐文章于 2020-08-11 10:10:15 发布

ACfunhsl

最新推荐文章于 2020-08-11 10:10:15 发布

阅读量186

点赞数

分类专栏：题解

本文链接：https://blog.csdn.net/Joker_He/article/details/107324619

版权

题解专栏收录该内容

279 篇文章 1 订阅

订阅专栏

题目传送门

题意： 给你n个点（n<=2000），问你最多有多少个点共圆且这个圆经过原点(0,0)。

思路： 第一层枚举点a，第二层枚举点b，如果这两个点斜率相同，则没有交点，否则用map对交点坐标计数，每次维护最大值就行。（~~这题卡精度，我吐了~~ ）时间复杂度是O(n²log(n))。

**代码 **

#include<bits/stdc++.h>
#pragma GCC optimize("Ofast")
#define endl '\n'
#define null NULL
#define ls p<<1
#define rs p<<1|1
#define fi first
#define se second
#define mp make_pair
#define pb push_back
#define ll long long
#define int long long
#define pii pair<int,int>
#define pdd pair<double,double>
#define ull unsigned long long
#define all(x) x.begin(),x.end()
#define IOS ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
#define ct cerr<<"Time elapsed:"<<1.0*clock()/CLOCKS_PER_SEC<<"s.\n";
char *fs,*ft,buf[1<<20];
#define gc() (fs==ft&&(ft=(fs=buf)+fread(buf,1,1<<20,stdin),fs==ft))?0:*fs++;
inline int read(){int x=0,f=1;char ch=gc();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=gc();}
while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=gc();}
return x*f;}
using namespace std;
const int N=2e3+5;
const int inf=0x3f3f3f3f;
const int mod=998244353;
const double eps=1e-7;
const double PI=acos(-1);
int x[N],y[N];
signed main()
{
    int n;
    cin>>n;
    for(int i=1;i<=n;i++)
        cin>>x[i]>>y[i];
    int res=0;
    for(int i=1;i<=n;i++)
    {
        map<pdd,int>ma;
        for(int j=i+1;j<=n;j++)
        {
            int f=x[i]*y[j]-x[j]*y[i];
            if(f==0)
                continue;
            double xx=1.0*(y[i]*y[j]*y[j]+y[i]*x[j]*x[j]-y[j]*y[i]*y[i]-x[i]*x[i]*y[j])/f;
            double yy=1.0*(x[i]*x[j]*x[j]+x[i]*y[j]*y[j]-x[j]*x[i]*x[i]-x[j]*y[i]*y[i])/f;
            ma[mp(xx,yy)]++;
            res=max(res,ma[mp(xx,yy)]);
        }
    }
    cout<<res+1<<endl;
}