Hou Yi's secret(UVALive - 5714)枚举 精度处理

一、题目大意

给出n个点，看最多能凑成多少个相似三角形。

注意：坑点挺多。

二、题目思路

记录下来点，三层循环枚举所有可能组成的三角形，二层循环找出哪一组的相似三角形最多。

坑点：

(1)给的点会有重复的。

解决方法：设置一个vis二维数组，记录下点，注意记录的是vis[100+x][100+y]。

(2)如何判断三个点是否会组成一个三角形。

解决方法：此处不能用“任意两边之和大于第三边，任意两边之差小于第三边”。

用到的事 如果三个点在一条直线上，那么就不能构成三角形，并且注意由于之后需要判断相似，边的大小应该排序。

(3)如何判断三角形相似。

解决方法：并不是单纯的三个边对应成比例，而是要有精度的处理。

三、附代码

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<stack>
#include<queue>
#include<cstring>
#include<string>
#include<set>
#include<cmath>
#include<map>
#include<sstream>
using namespace std;
#define inf 0x3f3f3f3f
typedef long long LL;
const int maxn = 18 + 8;

struct Node
{
    int x,y;
}node[maxn];

struct Tri
{
    double x,y,z;
}tri[maxn * maxn];

int n,xx,yy;
double x1,x2,x3;
double Cal(int x1,int x2,int y1,int y2)
{
    return sqrt( (x2-x1) * (x2-x1) + (y2-y1) * (y2-y1) );
}
bool Judge(Node a, Node b, Node c)
{
    x1 = Cal(a.x,b.x,a.y,b.y);
    x2 = Cal(a.x,c.x,a.y,c.y);
    x3 = Cal(b.x,c.x,b.y,c.y);
    if(x1 > x2) swap(x1,x2);
    if(x1 > x3) swap(x1,x3);
    if(x2 > x3) swap(x2,x3);
    if( (a.y - b.y) * (a.x - c.x)  == (a.x - b.x)  * (a.y - c.y) )  return false;
    return true;
}

bool JudgeTri(Tri a,Tri b)
{
    if(fabs(a.x * b.y - b.x * a.y) <= 1e-6 && fabs(a.x * b.z -a.z * b.x) <= 1e-6 && fabs(a.z * b.y - a.y * b.z) <= 1e-6) return true;
    return false;
}
int main()
{
    while(cin >> n)
    {
        if(!n) break;
        int cnt = 0;
        int num = 0;
        int ans = 0;
        int vis[205][205];
        memset(vis,0,sizeof vis);
        for(int i = 0; i < n; ++i)
        {
            cin >> xx >> yy;
            if(!vis[xx+100][yy+100])
            {
                node[cnt].x = xx;
                node[cnt++].y = yy;
                vis[xx+100][yy+100] = true;
            }
        }

        for(int i = 0; i < cnt; ++i){
            for(int j = i+1; j < cnt; ++j) {
                for(int k = j+1; k < cnt; ++k){
                    if(Judge(node[i],node[j],node[k])){
                        tri[num].x = x1;
                        tri[num].y = x2;
                        tri[num++].z = x3;
                    }
                }
            }
        }

        //cout << num << endl;

        int res;
        for(int i = 0; i < num; ++i)
        {
            res = 0;
            for(int j = 0; j < num; ++j)
            {
                if(JudgeTri(tri[i],tri[j]))
                {
                    res++;
                }

            }
            ans = max(ans,res);
        }
        cout << ans << endl;
    }
}