欢迎您访问程序员文章站本站旨在为大家提供分享程序员计算机编程知识!
您现在的位置是: 首页  >  IT编程

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

程序员文章站 2022-04-06 13:12:54
题目传送门题意: 给你n个点(n<=2000),问你最多有多少个点共圆且这个圆经过原点(0,0)。思路: 第一层枚举点a,第二层枚举点b,如果这两个点斜率相同,则没有交点,否则用map对交点坐标计数,每次维护最大值就行。(这题卡精度,我吐了 )时间复杂度是O(n2log(n))。**代码 **#include#pragma GCC optimize("Ofast")#define endl '\n'#define null NULL#defin...

题目传送门

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

思路: 第一层枚举点a,第二层枚举点b,如果这两个点斜率相同,则没有交点,否则用map对交点坐标计数,每次维护最大值就行。(这题卡精度,我吐了 )时间复杂度是O(n2log(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;
}

本文地址:https://blog.csdn.net/Joker_He/article/details/107324619