C++求平面上不重合的n个点最多构成多少条两两互不平行（包括重合）的直线

UPC NO.78场 问题 E: 阅兵队形 plane
题目描述
70 周年阅兵的时候，飞机在空中排练着队形，Yyx 很好奇，他想知道这么训练有素的队形到底是如何造就的呢？他记录下了飞行路径上的各个端点。
他发现：把整个天空看做一个平面直角坐标系，飞行路径是所有过任意两个端点的直线。
如果这些飞机可能会撞在一起，或者说只要这些直线有交点，就可能发生事故。
在所有直线中应该最少删除多少条直线使得剩下的直线两两都不相互平行（重合也是平行）。
求出最多可以构成多少条两两互不平行的直线。

输入
第一行，整数 n。
接下来 n 行，每行两个整数 x,y，表示这个端点的横坐标与纵坐标。
输出
一行，一个整数 ans，表示答案。

样例输入
4
-1 1
-2 0
0 0
1 1
样例输出
4

提示
对于所有数据，2≤N≤200，-1000≤Xi≤1000，-1000≤Yi≤1000。
在所有直线中应该最少删除多少条直线使得剩下的直线两两都不相互平行。

下面2种算法思路的直线都是通过n^2枚举出来的

#include<bits/stdc++.h>
using namespace std;
int n,ans,x1[520],y1[520];
map<double,bool>f;

int main()
{
    cin>>n;
    for(int i=1; i<=n; i++)  cin>>x1[i]>>y1[i];
    for(int i=1; i<=n; i++)
        for(int j=i+1; j<=n; j++)
        {
            int lx = x1[i]-x1[j];
            int ly = y1[i]-y1[j];
            double k;
            if(lx==0)  k=5.25;
            else k = (double)ly/lx;
            if(!f[k]) {f[k]=1;    ans++;}
        }
    cout<<ans<<endl;
    return 0;
}

/*
另一写法：求出每两个点间的斜率，排序后判断
#include<bits/stdc++.h>
#define N 201
#define eps 1e-6
#define PI acos(-1.0)
using namespace std;
struct Point {double x, y;} p[N];
double a[N*N];
int n;

double FindSlewRate(Point p1, Point p2)
{
    Point p;
    p.x = p2.x-p1.x;
    p.y = p2.y-p1.y;
    if(abs(p.x)<eps)  return PI/2;     //斜率不存在
    double tmp = atan(p.y/p.x);
    if(tmp<0)  return PI+tmp;
    return tmp;
}
 
int cmp(const void *c, const void *d){
    return *(double *)c > *(double *)d ? 1:-1;
}

int main()
{
    while(~scanf("%d", &n))
    {
        for(int i=0; i<n; i++)
            scanf("%lf%lf", &p[i].x, &p[i].y);
        
        int rt=0;
        for(int i=0; i<n; i++)
            for(int j=i+1; j<n; j++)
                a[rt++] = FindSlewRate(p[i], p[j]);
        
        qsort(a,rt,sizeof(a[0]),cmp);
        
        int ans=1;
        for(int i=1; i<rt; i++)
            if(a[i] != a[i-1])    ans++;
        
        printf("%d\n", ans);
    }
    return 0;
}
*/

#include<bits/stdc++.h>
using namespace std;
int x[520],y[520];
bool flag[5200][5200];
int n,ans;
int main()
{
    cin>>n;
    for(int i=1;i<=n;i++)    cin>>x[i]>>y[i];
    
    for(int i=1;i<=n;i++)
        for(int j=i+1;j<=n;j++)
        {
            int lx=x[i]-x[j];
            int ly=y[i]-y[j];
            if(lx<0) 
            {
                lx = -lx;
                ly = -ly;
            }
            else if(lx==0)    ly=abs(ly);
            int gcd = __gcd(lx,abs(ly));
            lx /= gcd;
            ly /= gcd;
            if(!flag[lx][ly+2005])
            {
                ans++;
                flag[lx][ly+2005] = 1;
            }
        }
        
    cout<<ans<<endl;
    return 0;
}