fzu 2273 几何学

Problem 2273 Triangles

Accept: 207    Submit: 701
Time Limit: 1000 mSec    Memory Limit : 262144 KB

 Problem Description

This is a simple problem. Given two triangles A and B, you should determine they are intersect, contain or disjoint. (Public edge or point are treated as intersect.)

 Input

First line contains an integer T (1 ≤ T ≤ 10), represents there are T test cases.

For each test case: X1 Y1 X2 Y2 X3 Y3 X4 Y4 X5 Y5 X6 Y6. All the coordinate are integer. (X1,Y1) , (X2,Y2), (X3,Y3) forms triangles A ; (X4,Y4) , (X5,Y5), (X6,Y6) forms triangles B.

-10000<=All the coordinate <=10000

 Output

For each test case, output “intersect”, “contain” or “disjoint”.

 Sample Input

2

0 0 0 1 1 0 10 10 9 9 9 10

0 0 1 1 1 0 0 0 1 1 0 1

 Sample Output

disjoint

intersect

题意： 给你 6 个点 前三个点为一个三角形，后三个点为一个三角形。 问你两个三角形是相交 相离 还是包含关系

思路：我们怎么判断三角形相含呢 肯定是一个三角形的三个点 都在另一个三角形的内部。我们怎么判断一

个点是否在另一个三角形的内部呢？ 有两种方法，第一种就是 判断这个点和三角形的任意一条边组成

的三角形的面积和 如果大于原三角形的面积 那么肯定就是在外部 不然就是在内部或者在边上，

第二种方法 就是做叉乘 假设一个三角形的三个点为a，b，c 点d 是否在内部 我们将向量

ab 与向量ad 得到ans1 向量ab 与向量ac 做叉乘得到ans2 如果正负相

同 表示 c 和d 在同一边 这样对于 bc 判断 ad 是否在同一边 对于ca 判断 bd 是否在一边。

那我们怎么判断三角形是相交呢？ 可以通过判断这个三角形的任意一条边 是否和 另一个三角形的

边相交。判断两个线段相交： 线段是否相交

其他的情况肯定就是 相离。

代码：

#include<stdio.h>
#include<string.h>
#include<math.h>
#include<iostream>
#include<algorithm>

using namespace std;

const double eps=1e-8;
/*
typedef struct point
{
	double x,y;
	point(double x=0,double y=0 ):x(x),y(y){}
}vvecotr;
*/
typedef struct point
{
    double x,y;
    point(double x=0,double y=0):x(x),y(y){}
}vvector;

vvector operator + (const vvector &a,const vvector &b)
{
	return vvector(a.x+b.x,a.y+b.y);
}

vvector operator - (const vvector &a,const vvector &b)
{
	return vvector(a.x-b.x,a.y-b.y);
}

vvector operator * (const vvector &a,double p)
{
	return vvector(a.x*p,a.y*p);
}

vvector operator / (const vvector &a,double p)
{
	return vvector(a.x/p,a.y/p);
}

int dcmp(double x)
{
	if(fabs(x)<eps) return 0;
	return x<0?-1:1;
}

double dot(const vvector &a,const vvector &b)
{
	return a.x*b.x+a.y*b.y;
}

double corss(const vvector &a,const vvector &b)
{
	return a.x*b.y-b.x*a.y;
}

bool judxianduanjiao(const point &a1,const point &a2,const point &b1,const point &b2)
{
	double c1=corss(a2-a1,b1-a1);  double c2=corss(a2-a1,b2-a1);
	double c3=corss(b2-b1,a1-b1);  double c4=corss(b2-b1,a2-b1);
	return dcmp(c1)*dcmp(c2)<=0&&dcmp(c3)*dcmp(c4)<=0;
}

bool sameside(point a,point b,point c,point d)
{
	vvector ab=b-a,ac=c-a,ad=d-a;
	double jud1=corss(ab,ac);
	double jud2=corss(ab,ad);
	if(jud1*jud2>=0) return 1;
	return 0;
}

bool intri(point a,point b,point c,point d)
{
	return sameside(a,b,c,d)&&sameside(b,c,a,d)&&sameside(c,a,b,d);
}

point p[7];

int main()
{
	int cas;
	cin>>cas;
	while(cas--)
	{
		for(int i=0;i<6;i++)
		{
			cin>>p[i].x>>p[i].y;
		}
		
		int flag=0;
		for(int i = 0; i<3; i++){  
            if(judxianduanjiao(p[i],p[(i+1)%3],p[3],p[4]) || judxianduanjiao(p[i],p[(i+1)%3],p[3],p[5]) || judxianduanjiao(p[i],p[(i+1)%3],p[5],p[4]))  
            {  
                flag = 1;  
                break;  
            }  
        }
        
        
        //printf("flag : %d \n",flag); 
        
        if(flag)
        {
        	printf("intersect\n");
        	continue;
		}
		
		if(intri(p[0],p[1],p[2],p[3])&&intri(p[0],p[1],p[2],p[4])&&intri(p[0],p[1],p[2],p[5]))
		{
			printf("contain\n");
			continue;
		}
		
		if(intri(p[3],p[4],p[5],p[0])&&intri(p[3],p[4],p[5],p[0])&&intri(p[3],p[4],p[5],p[0]))
		{
			printf("contain\n");
			continue;
		}
		
		printf("disjoint\n");
	}
	return 0;
}