hdu4082 求相似三角形，模拟题

A了4081顺便看了一下4082，感觉题目不好翻译，前几段都是没用的，意思就是给你n个点(0<n<18)，数据范围很小…

求这n个点能组成最多有多少个相似三角形。

暴力枚举，但是要注意很多细节，首先：

①去重点，在图论里总是会有去掉重边，可是计算几何里很少遇到去掉重点的问题，所以容易被遗忘，so……记住！！！！

②判断三个点能否组成三角形，即是否共线。

后台数据应该水了，一开始的代码，对于

4

1 1

2 2

3 3

4 4

这组数据输出1也AC了~

#include <stdio.h>
#include <iostream>
#include <math.h>
#include <algorithm>
#define eps 1e-8
using namespace std;
struct point{
    int x,y;
}p[25];
struct triangle{
    double a,b,c;
}t[410];
int cross(point p0,point p1,point p2){
    return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);
}
double dist(point p1,point p2){
    return sqrt((double)(p2.x-p1.x)*(p2.x-p1.x)+(p2.y-p1.y)*(p2.y-p1.y));
}
bool similar(triangle t1,triangle t2){
    if( (fabs(t1.a*t2.b-t1.b*t2.a)<eps)&&
        (fabs(t1.b*t2.c-t1.c*t2.b)<eps)&&
        (fabs(t1.a*t2.c-t1.c*t2.a)<eps) )
    return 1;
    return 0;
}
int main(){
    int n;
    while(scanf("%d",&n)!=EOF&&n){
        for(int i=0;i<n;i++){
            scanf("%d%d",&p[i].x,&p[i].y);
            for(int j=0;j<i;j++){
                if(p[i].x==p[j].x&&p[i].y==p[j].y){
                    i--;n--;
                    break;
                }
            }
        }
        if(n<3){
            printf("0\n");
            continue;
        }
        int top=0;
        double d[3];
        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[k])!=0){
                        d[0]=dist(p[i],p[j]);
                        d[1]=dist(p[i],p[k]);
                        d[2]=dist(p[k],p[j]);
                        sort(d,d+3);
                        t[top].a=d[0];
                        t[top].b=d[1];
                        t[top].c=d[2];
                        top++;
                    }
                }
            }
        }
        if(top==0){
            printf("0\n");
            continue;
        }
        int _max=0;
        for(int i=0;i<top;i++){
            int res=0;
            for(int j=i+1;j<top;j++){
                if(similar(t[i],t[j])) res++;
            }
            _max=max(res,_max);
        }
        printf("%d\n",_max+1);
    }
    return 0;
}