POJ1118 Lining Up

快弄死我了 最后的原因是abs和fabs的区别。。。

说点收获：
1.cmp函数返回的是int，所以不要直接返回double相减的结果
2.define inf 1e9和eps 1e-9
3.在整数相除得到double时要1.0*
4.加上<cmath> ，用fabs取代abs，abs是用来整数取绝对值，应该用fabs
5.思路上，本题可从线的角度考虑转化为从点的角度考虑，这样变为N^2，然后不用hash，然后对斜率，用排序，又有logN。最后为O（N^2*logN）
6.注意边界情况
7.最后找错的时候采用了随机数数据和传说中正确的代码做比较，发现它们也还是有错，只是数据弱而已
8.下面的Discuss还是很有用的
9. GCD的思路也不错，应当会写GCD
此题做得真伤元气

#include <iostream>

#include <cmath>



using namespace std;



#define inf 1e9

#define eps 1e-9



int cmp(const void * a, const void * b)

{

	double * x = (double *) a;

	double * y = (double *) b;

	if (fabs( *x - *y) < eps) return 0 ;

	else if (*x < *y) return -1;

	else return 1;

}



int main()

{

	int count = 0;

	while (true)

	{

		int n;

		cin >> n;

		if (n == 0)

		{

			return 0;

		}

		count++;

		int max = 0;

		int* x = new int[n];

		int* y = new int[n];

		for (int i = 0; i < n; i++)

		{

			cin >> x[i];

			cin >> y[i];

		}

		if (n == 1)

		{

			max = 2;

		}

		else if (n == 2)

		{

			max = 2;

		}

		else

		{			

			for (int i = 0; i < n; i++)

			{

				double * tmp = new double[n];

				int index = 0;

				for (int j = i+1; j < n; j++)

				{

					if (x[i] == x[j]) tmp[index] = inf;

					else 

					{

						tmp[index] = 1.0 * (y[i] - y[j]) / (x[i] - x[j]);

					}

					index++;

				}

				qsort(tmp, index, sizeof (double), cmp);

				int k = 1;

				int prev = 0;

				int localMax = 2;

				while (k < index)

				{

					if (fabs(tmp[prev] - tmp[k]) < eps)

					{

						localMax++;

					}

					else

					{

						localMax = 2;

					}

					if (localMax > max) 

					{

							max = localMax;

					}

					prev = k;

					k++;

				}

			}		

		}

		cout << max << endl;

	}

	return 0;

}