POJ 1569 Myacm Triangles【判定点在三角形内】

                                     POJ 1569 Myacm Triangles

                                                         http://poj.org/problem?id=1569

 

题意

给你n个点的坐标,要你找出面积最大的由3个点构成的三角形,且该三角形的内部或边上不能有点. 点数<=15。

输入

多组输入，输入整数N（4<=N<=15 ，N=0表示结束），N表示有N个点，每个点的编号为‘A’、‘B’...'Z'。接下来N行，每行由一个字符和两个整数组成，分别表示点的编号和点的二维坐标。

输出

每组输出题目要求的一个三角形（按升序输出三角形顶点的编号）。

样例输入

6
A 1 0
B 4 0
C 0 3
D 1 3
E 4 4
F 0 6
4
A 0 0
B 1 0
C 99 0
D 99 99
0

样例输出

BEF
BCD

分析

枚举每个三角形，然后判断其他点是否在三角形的内部或边上。判断一个点是否在一个三角形的边上或内部的方法如下：

方法一: 如果一个点P在三角形ABC内部(或边上),那么三角形PAB+三角形PBC+三角形PAC的面积 == ABC的面积。否则前者面积和应该更大。

方法二: 如果一个点P在三角形ABC内部(或边上), Cross(A-B,P-B)与Cross(B-C,P-C),与Cross(C-A,P-A) 一定要么同时>=0,要么同时<=0. (注意Cross内ABC三点的书写顺序不能改)。

本程序使用方法一，具体看程序。

C++程序

#include<iostream>
#include<cmath>

using namespace std;

const int N=20;
const double EPS=1e-8;

struct Point{
	double x,y;
	Point(){}
	Point(double x,double y):x(x),y(y){}
	Point operator -(const Point &a)const 
	{
		return Point(x-a.x,y-a.y);
	}
	double operator ^(const Point &a)const
	{
		return x*a.y-y*a.x;
	}
}; 

//存储结果 
struct Ans{
	Ans()
	{
		area=-1.0;
		i=j=k=-1;
	}
	double area;
	int i,j,k;
};

Point p[N];

//判断其他点是否在ijk构成的三角形中 
bool check(int i,int j,int k,double area,int n)
{
	for(int t=0;t<n;t++)
	{
		if(t==i||t==j||t==k) continue;
		double s1=fabs((p[i]-p[t])^(p[j]-p[t]));
		double s2=fabs((p[i]-p[t])^(p[k]-p[t]));
		double s3=fabs((p[j]-p[t])^(p[k]-p[t]));
		double sum=(s1+s2+s3)/2.0;
		if(fabs(sum-area)<EPS) 
		  return true;
	}
	return false;
}

int main()
{
	int n;
	while(~scanf("%d",&n)&&n)
	{
		char c;
		for(int i=0;i<n;i++)
		  scanf(" %c%lf%lf",&c,&p[i].x,&p[i].y);
		Ans ans;//答案 
		for(int i=0;i<n;i++)
		  for(int j=i+1;j<n;j++)
		    for(int k=j+1;k<n;k++)
		    {
		    	//计算i,j,k三点构成的三角形的面积
				double area=fabs((p[i]-p[j])^(p[k]-p[j]))/2.0;
				if(area<ans.area) continue;
				if(!check(i,j,k,area,n))
				{
					ans.area=area;
					ans.i=i;
					ans.j=j;
					ans.k=k;
				}
			}
		printf("%c%c%c\n",ans.i+'A',ans.j+'A',ans.k+'A');
	}
	return 0;
}