HDU 4082 Hou Yi's secret 伪相似三角形 && 2011 Asia Beijing Regional Contest

题意是：给你很多点，每三个不共线的点可以连成三角形，相似的三角形可以归为一个集合，问最大的集合三角形的数量。

Trick：有三点共线以及重合的点，所以需要点去重和判定三点共线

解法：水暴力，去重后枚举每三个点，如果不共线就把三条边排序，除以最大公约数，放入结构体，再加一重for循环 判定这个三角形是否存在。O(n^4) = 15ms  数据太水了！！

#include<stdio.h>
#include<algorithm>
#include<string.h>
#include<math.h>
#include<stdlib.h>
#include<iostream>
using namespace std;
const int maxn=500000;
const double eps=1e-8;
const double pi = acos(-1.0);

inline int dcmp(const double& x) {return (x>eps)-(x<-eps);}

struct Point
{
    int x, y;
    Point(int x=0, int y=0): x(x),y(y) {}
    Point operator-(const Point& p) const {return Point(x-p.x,y-p.y);}
    bool operator<(const Point& p) const {if(x!=p.x) return x<p.x; return y<p.y;}
    bool operator==(const Point& p) const {return x==p.x&&y==p.y;}
    bool read(){scanf("%d%d",&x,&y);return true;}
};

inline int gcd(int a, int b)
{
    if(b==0) return a;
    return gcd(b,a%b);
}

inline int dist(const Point& a,const Point& b) {return (a.x-b.x)*(a.x-b.x) + (a.y-b.y)*(a.y-b.y);}
inline int Cross (const Point& a, const Point& b) {return a.x*b.y-a.y*b.x;}
Point p[maxn];

struct node
{
    int a, b, c, d;
    node(int a=0, int b=0, int c=0, int d=0):a(a),b(b),c(c),d(d) {}
    bool operator==(const node& n) {return a==n.a&&b==n.b&&c==n.c&&d==n.d ;}
};

pair<node,int> v[maxn];
int tnum = 0;
int main()
{
    //freopen("in.txt","r",stdin);
    int n;
    while(~scanf("%d", &n) && n)
    {
        for(int i = 0 ; i < n ; ++ i)
        {
            p[i].read();
        }
        sort(p,p+n);
        n = unique(p,p+n)-p;
        tnum = 0;
        for(int i = 0 ; i < n ; ++ i)
        {
            for(int j = i + 1; j < n ; ++ j)
            {
                for(int k = j + 1; k < n ; ++ k)
                {
                    if(Cross(p[i]-p[j], p[i]-p[k])==0) continue;
                    int d1 = dist(p[i],p[j]), d2 = dist(p[j],p[k]), d3 = dist(p[i],p[k]);
                    if(d1>d2) swap(d1, d2); if(d1>d3) swap(d1, d3); if(d2>d3) swap(d2, d3);
                    int d4 = d1; int temp;
                    temp = gcd(d1, d2); d1/=temp;d2/=temp;
                    temp = gcd(d3, d4); d3/=temp;d4/=temp;
                    node tmp = node(d1,d2,d3,d4);
                    bool flag = false;
                    for(int i = 0 ; i < tnum ; ++i )
                    {
                        if(v[i].first==tmp)
                        {
                            v[i].second ++;
                            flag = true;
                            break;
                        }
                    }
                    if(!flag)
                    {
                        v[tnum].first = tmp;
                        v[tnum++].second = 1;
                    }
                }
            }
        }
        int ans = 0;
        for(int i = 0 ; i < tnum ; ++ i)
        {
            ans = max(ans, v[i].second);
        }
        printf("%d\n",ans);
    }
    return 0;
}