XDU 1008 - 数星星

思路：若比较三点是否在同一条直线上直接不好比，将其化为斜率，若斜率相等则在同一直线上。

            从前往后枚举每一个点作为起点，然后每次计算出它之后的点到它的斜率，对斜率进行排序，找出斜率相等的最多点数即为所求。

#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<algorithm>
using namespace std;
#define eps 1e-10

int n;
struct point
{
    double x,y;
};
struct point star[21000+10];

bool cmp(double x,double y)
{
    return x<y;
}

int main()
{
    int i,j;
    while(scanf("%d",&n)!=EOF)
    {
        int ans=0,max_num;
        if(n==0)
            printf("0\n");
        else
        {
            for(i=1;i<=n;i++)
                scanf("%lf%lf",&star[i].x,&star[i].y);
            for(i=1;i<n;i++)//以每一个点作为斜率中的一个点
            {
                int num=0;
                max_num=0;
                double gradient[21005];//定义当前点的斜率
                for(j=i+1;j<=n;j++)//求出所有点的斜率
                {
                    if(star[i].x==star[j].x)//x坐标相同，斜率无穷大
                        max_num++;
                    else                         //一般情况
                        gradient[++num]=(double)(star[j].y-star[i].y)/(star[j].x-star[i].x);
                }
                if(max_num>ans)
                    ans=max_num;
                max_num=1;

                sort(gradient+1,gradient+1+num,cmp);
                //从现在开始寻找斜率最多相同的为止
                double pre=gradient[1];
                for(j=2;j<=num;j++)
                {
                    if(fabs(pre-gradient[j])<eps)//当前斜率和上一条斜率相等
                        max_num++;
                    else                        //当前斜率与上一条斜率不等
                    {
                        if(max_num>ans)
                            ans=max_num;
                        max_num=1;
                        pre=gradient[j];
                    }
                }
                if(max_num>ans)//必须要再加上一个
                    ans=max_num;
            }
            printf("%d\n",ans+1);
        }
    }
    return 0;
}