hdu 4082 计算几何 相似三角形

题意 给出n个点，判断组成的三角形中最多的一组相似三角形的个数。
比较菜，没写过这种题，不过不难。由于没写过，不知道会不会有精度问题，于是决定不开根号。
首先这题的坑 有可能给重复的点，重复的点要先去掉，然后枚举，判断三点共线，将能组成三角形的存到一个三角形结构体中。
然后对于每个三角形保存三个角的cos值，但由于担心精度问题 因为 cos = ( a^ 2 +b ^ 2 - c ^ 2 ) / ( 2 * a * b )
将 2 * a * b 乘到等式左边 再两边平方 得到 4 * a ^ 2 * b ^ 2 * cos ^ 2 = ( a ^ 2 + b ^ 2 - c ^ 2 ) ^ 2 
将等式中cos的系数设为x ， 等式右边设为y  ，然后取g = gcd （ x , y ) , x /= g , y /= g, 将x,y 存在结构体Cos中
保存在三角形中的三个Cos先排序再存，方便比较。
最后枚举三角形寻找最多的相似三角形即可。
( 代码写的比较丑，训练时候也懒得改了。。。。 而且发现好像不会出现精度问题。。。 ）

#include <bits/stdc++.h>
using namespace std;

struct Cos{
    long long x,y;
    bool operator < (Cos A)
    {
        double x1,x2,y1,y2;
        x1=x;
        x2=A.x;
        y1=y;
        y2=A.y;
        return y1/x1<y2/x2;
    }
    bool judge(Cos A)
    {
        if(x==A.x&&y==A.y)
            return true;
        else return false;
    }
};

struct Point{
    long long x,y;
    Point operator - (Point A)
    {
        Point p1;
        p1.x=x-A.x;
        p1.y=y-A.y;
        return p1;
    }
    bool operator < (Point A)
    {
        if(x<A.x)
            return true;
        else if(x==A.x&&y<A.y)
            return true;
        else return false;
    }
    bool judge(Point A)
    {
        if(x==A.x&&y==A.y)
            return true;
        else return false;
    }
};

struct Triangle{
    Point p1,p2,p3;
    Cos cos1,cos2,cos3;
    bool judge (Triangle A)
    {
        if(cos1.judge(A.cos1)&&cos2.judge(A.cos2)&&cos3.judge(A.cos3))
            return true;
        else return false;
    }
};

Point p_oint[25],point[25];
Triangle triangle[1005];

bool judge_seg(Point p1,Point p2,Point p3)
{
    Point t1,t2;
    t1=p1-p2;
    t2=p1-p3;
    if(t1.x*t2.y==t1.y*t2.x)
        return true;
    else return false;
}

long long seg_len(Point p1,Point p2)
{
    return ((p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y));
}

long long gcd(long long a,long long b)
{
    return b?gcd(b,a%b):a;
}

int main()
{
    int a,n,m,x,y,i,j,k;
    while(cin>>n)
    {
        if(n==0)
            break;
        for(i=0;i<n;i++)
            cin>>p_oint[i].x>>p_oint[i].y;
        int now=0;
        sort(p_oint,p_oint+n);
        for(i=0;i<n;i++)
        {
            if(i<n-1)
            {
                if(!p_oint[i].judge(p_oint[i+1]))
                {
                    point[now]=p_oint[i];
                    now++;
                }
            }
            else
            {
                point[now]=p_oint[i];
                now++;
            }
        }
        m=0;
        for(i=0;i<now;i++)
        {
            for(j=i+1;j<now;j++)
            {
                for(k=j+1;k<now;k++)
                {
                    if(!judge_seg(point[i],point[j],point[k]))
                    {
                        triangle[m].p1=point[i];
                        triangle[m].p2=point[j];
                        triangle[m++].p3=point[k];
                    }
                }
            }
        }
        long long ans=0;
        for(i=0;i<m;i++)
        {
            Cos test[3];
            test[0].x=4*(seg_len(triangle[i].p1,triangle[i].p2))*(seg_len(triangle[i].p1,triangle[i].p3));
            test[0].y=(seg_len(triangle[i].p1,triangle[i].p2)+seg_len(triangle[i].p1,triangle[i].p3)-seg_len(triangle[i].p2,triangle[i].p3))*(seg_len(triangle[i].p1,triangle[i].p2)+seg_len(triangle[i].p1,triangle[i].p3)-seg_len(triangle[i].p2,triangle[i].p3));
            test[1].x=4*(seg_len(triangle[i].p2,triangle[i].p1))*(seg_len(triangle[i].p2,triangle[i].p3));
            test[1].y=(seg_len(triangle[i].p2,triangle[i].p1)+seg_len(triangle[i].p2,triangle[i].p3)-seg_len(triangle[i].p1,triangle[i].p3))*(seg_len(triangle[i].p2,triangle[i].p1)+seg_len(triangle[i].p2,triangle[i].p3)-seg_len(triangle[i].p1,triangle[i].p3));
            test[2].x=4*(seg_len(triangle[i].p3,triangle[i].p1))*(seg_len(triangle[i].p3,triangle[i].p2));
            test[2].y=(seg_len(triangle[i].p3,triangle[i].p1)+seg_len(triangle[i].p3,triangle[i].p2)-seg_len(triangle[i].p1,triangle[i].p2))*(seg_len(triangle[i].p3,triangle[i].p1)+seg_len(triangle[i].p3,triangle[i].p2)-seg_len(triangle[i].p1,triangle[i].p2));
            long long g1,g2,g3;
            g1=gcd(test[0].x,test[0].y);
            g2=gcd(test[1].x,test[1].y);
            g3=gcd(test[2].x,test[2].y);
            test[0].x/=g1;
            test[0].y/=g1;
            test[1].x/=g2;
            test[1].y/=g2;
            test[2].x/=g3;
            test[2].y/=g3;
            sort(test,test+3);
            triangle[i].cos1=test[0];
            triangle[i].cos2=test[1];
            triangle[i].cos3=test[2];
        }
        for(i=0;i<m;i++)
        {
            long long sum=1;
            for(j=i+1;j<m;j++)
            {
                if(triangle[i].judge(triangle[j]))
                {
                    sum++;
                }
            }
            ans=max(sum,ans);
        }
        cout<<ans<<endl;
    }
    return 0;
}