（2020牛客暑期多校训练营）[二]Boundary

题目描述
给定平面上n个点，考虑所有原点(0,0)在边界上的圆，找出最多有几个点可以在同一个圆的边界上。输出数量。

例：考虑圆(x-1)^2+(y-1)2=2，我们会发现原点在圆的边界上，另外有三个给定的点(0,2),(2,0),(2,2)也在圆的边界上。
样例
输入
4
1 1
0 2
2 0
2 2
输出
3
思路
直接枚举，时间复杂度为O(n³)，显然超时。先枚举两个点，再算出圆心，统计有多少个圆心
在同一个点上。n个圆心代表有n+1个点在圆上
代码如下
#include<bits/stdc++.h>
using namespace std;
const int MAXN=2005;
int n,ans;
double a[MAXN],b[MAXN],x,y;
typedef pair<double,double> p;
map<p,int> mx;
int main()
{
scanf("%d",&n);
for(int i=1;i<=n;i++) scanf("%lf%lf",a+i,b+i);
for(int i=1;i<=n;i++)
{
mx.clear();
for(int j=i+1;j<=n;j++)
if(a[i]*b[j]!=a[j]*b[i])
{
x=((b[j]-b[i])*b[i]*b[j]-a[i]*a[i]*b[j]+a[j]*a[j]*b[i])/(a[j]*b[i]-a[i]*b[j]);
y=((a[j]-a[i])*a[i]*a[j]-b[i]*b[i]*a[j]+b[j]*b[j]*a[i])/(b[j]*a[i]-b[i]*a[j]);
mx[p(x,y)]++; ans=max(mx[p(x,y)],ans);
}
}
printf("%d\n",ans+1);
}