HDU 4082 Hou Yi's secret（相似三角形最大数量）

题目

题意：

后羿射了一些箭形成一些坑，这些坑可以组成三角形，求相似三角形最大的数量。

什么是相似三角形最大数量？

相似三角形可能形成很多个，我们把不同样子的三角形定义为不同类的三角形，那么在同一类的三角形中，数量最大的就是答案。

注意：数据中的点可能用重复，要去重。

#include <bits/stdc++.h>

using namespace std;
double eps=1e-6;
bool is[205][205];
struct Triangle{
    double a,b,c;
    bool v;
}t[1000];
struct Point
{
    double x;
    double y;
}a[20];
double dis(Point a,Point b)
{
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
bool line(Point a,Point b,Point c)
{
    if(fabs((b.x-a.x)*(c.y-a.y)-(c.x-a.x)*(b.y-a.y))<eps)return true;
    else return false;
}
void add(Point a,Point b,Point c,int J)///按照顺序添加三角形的三边
{
    //cout<<J<<" "<<a.x<<" "<<a.y<<" "<<b.x<<" "<<b.y<<" "<<c.x<<" "<<c.y<<"            ";
    double d1=dis(a,b);
    double d2=dis(a,c);
    double d3=dis(b,c);
    double sum=d1+d2+d3;
    t[J].a=min(d1,min(d2,d3));
    t[J].c=max(d1,max(d2,d3));
    t[J].b=sum-t[J].a-t[J].c;
   //cout<<t[J].a<<" "<<t[J].b<<" "<<t[J].c<<endl;
    t[J].v=false;
}
bool judge(Triangle x,Triangle y)///判断两个三角形是否相似，用乘法
{
    if((fabs(y.a*x.b-x.a*y.b)<eps)&&(fabs(y.b*x.c-x.b*y.c)<eps)&&(fabs(y.a*x.c-x.a*y.c)<eps))
        return true;
    else return false;
}
int main()
{
    ios::sync_with_stdio(false);///加个速
        int n;
        int  x,y;
        while(cin>>n&&n)
        {   memset(is,false,sizeof(is));
            int Y=0;
            for(int i=0;i<n;i++)
            {
                cin>>x>>y;
                x+=100;y+=100;
                if(!is[x][y]){is[x][y]=true;a[Y].x=x,a[Y].y=y;Y++;}///记忆化，防止数据是重复的点
            }
            int J=0;
            for(int i=0;i<Y;i++)
                for(int j=i+1;j<Y;j++)
                    for(int k=j+1;k<Y;k++)
            {
                if(!line(a[i],a[j],a[k]))///三点不共线构成三角形
                {
                    add(a[i],a[j],a[k],J);
                    J++;
                }
            }
            int ans=0;
            for(int i=0;i<J;i++)
            {   if(t[i].v)continue;
                int temp=1;
                t[i].v=true;
                for(int j=i+1;j<J;j++)
                {
                    if(t[j].v)continue;

                    if(judge(t[i],t[j])){/*cout<<"i="<<i<<" j="<<j<<endl;*/temp++;t[j].v=true;}///相似的一类三角形全标记

                }
                ans=max(ans,temp);//cout<<endl;
            }
            cout<<ans<<endl;

        }
    return 0;
}