(Relax ST1.13)POJ 2780 Linearity(给出若干个点,求最多有多少个点共线,不能使用n^3算法)

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cmath>

using namespace std;

const int maxn = 1010;
const int inf = INT_MIN;
struct Point{
int x;
int y;

bool operator<(const Point& b)const{//排序规则这里怎么定义都无所谓,但是一定要定义。。。
if(x == b.x){
return y < b.y;
}

return x < b.x;
}
}p[maxn];

double cal(const Point& a,const Point& b){
if(b.x == a.x){
return inf;
}

return (b.y - a.y)*1.0/(b.x - a.x);
}

int main(){
int n;
while(scanf("%d",&n)!=EOF){
int maxm = -10;

double xielv[maxn*2];

int i;
for(i = 0 ; i < n ; ++i){
scanf("%d%d",&p[i].x,&p[i].y);
}

sort(p,p+n);
int j;
for(i = 0 ; i < n-1 ; ++i){
int index = 0;
for(j = i+1 ; j < n ; ++j){
xielv[index++] = cal(p[i],p[j]);
}

sort(xielv,xielv+index);//对斜率进行排序....

/**
* 将斜率数组排序以后,这时候可能就会出现这么一种情况,
* 前后可能有好几个斜率是一样的...
* 以下就是数这样的斜率线每段有多少个...
*/
for(j = 0 ; j < index ; ++j){
int count = 0;
while(fabs(xielv[j] - xielv[j+1]) == 0){
count++;
j++;
}

maxm = (maxm>count)?maxm:count;
}

}

//		printf("%d\n",maxm+2);
cout<<maxm+2<<endl;
}

return 0;
}