hrbustoj 1142:围困（计算几何基础题，判断点是否在三角形内）

围困

Time Limit: 1000 MS     Memory Limit: 65536 K

Total Submit: 360(138 users) Total Accepted: 157(129 users)   Rating:  Special Judge: No

Description

Leyni是一名猎人，一天他在草原中散步，遇到了三只老狼，三只老狼围成了一个三角形，如果Leyni在这个三角形中，那么他必死无疑，否则他还有一线生机。

现在已知三只老狼的坐标和Leyni的坐标，请问，Leyni还有活下去的机会吗？

 

Input

本题有多组测试数据，对于每组测试数据：

第一行，含有4对坐标，x1, y1, x2, y2, x3, y3, xLeyni, yLeyni。前3对坐标表示三只老狼的位置，最后一对表示Leyni的坐标。坐标的范围在[-10000, 10000]，且都是整数。

注：输入数据保证Leyni不会在三角形的边上。

 

Output

对于每组测试数据：

第一行，如果Leyni还有一线生机，那么请输出Live，否则请输出Die

 

Sample Input

-1 -1 1 -1 0 1 0 0

-1 -1 1 -1 0 1 10 10

Sample Output

Die

Live

Author

齐达拉图

---

 

　　计算几何，基础题，判断点是否在三角形内。

　　这道题可以用来练习叉积的应用，思路是如果Leyni位置的点与三个顶点构成的向量与三个边代表的向量（必须是一个顺序，都是逆时针，或都是顺时针）的叉积都大于0或者都小于0，则Leyni在三角形内。需要注意的是不能只考虑一种情况，例如只考虑>0的情况，而不考虑<0的情况，这样会WA。

　　复杂代码：

 1 #include <stdio.h>

 2 #define eps 1e-10

 3 typedef struct{    //定义点

 4     double x,y;  5 }Point;  6 typedef Point Vector;    //定义向量

 7 double Cross(Vector a,Vector b)    //叉积

 8 {  9     return a.x*b.y - a.y*b.x; 10 } 11 bool isn(Point a,Point b,Point c,Point p)    //判断点p是否在三角形abc内

12 { 13     Vector ab,bc,ca;    //创建向量

14  Vector ap,bp,cp; 15     ab.x = b.x - a.x;ab.y = b.y - a.y;    //向量赋值

16     bc.x = c.x - b.x;bc.y = c.y - b.y; 17     ca.x = a.x - c.x;ca.y = a.y - c.y; 18     ap.x = p.x - a.x;ap.y = p.y - a.y; 19     bp.x = p.x - b.x;bp.y = p.y - b.y; 20     cp.x = p.x - c.x;cp.y = p.y - c.y; 21 

22     if( Cross(ab,ap)>eps && Cross(bc,bp)>eps && Cross(ca,cp)>eps || 

23         Cross(ab,ap)<eps && Cross(bc,bp)<eps && Cross(ca,cp)<eps)    //核心代码。注意不能单考虑一种情况（>0 或 <0）

24         return true; 25     else 

26         return false; 27 } 28 int main() 29 { 30     int x1,y1,x2,y2,x3,y3,xLeyni,yLeyni; 31     while(scanf("%d%d%d%d%d%d%d%d",&x1,&y1,&x2,&y2,&x3,&y3,&xLeyni,&yLeyni)!=EOF){    //输入四对坐标

32         if( (xLeyni==x1 && yLeyni==y1) ||

33             (xLeyni==x2 && yLeyni==y2) ||

34             (xLeyni==x3 && yLeyni==y3) ){ 35             printf("Die\n");    //如果Leyni正好在三角形顶点上，则直接死掉

36             continue; 37  } 38         Point a,b,c,p;    //将输入的坐标赋给点

39         a.x = double(x1),a.y = double(y1); 40         b.x = double(x2),b.y = double(y2); 41         c.x = double(x3),c.y = double(y3); 42         p.x = double(xLeyni),p.y = double(yLeyni); 43 

44         if(isn(a,b,c,p))    //判断是否在三角形内

45             printf("Die\n"); 46         else

47             printf("Live\n"); 48  } 49     return 0; 50 }

　　优化代码：

 1 #include <stdio.h>

 2 typedef struct {    //定义点

 3     double x,y;  4 }Point;  5 double Cross(Point a,Point b,Point c)    //三点求叉积

 6 {  7     return (b.x-a.x)*(c.y-a.y) - (b.y-a.y)*(c.x-a.x);  8 }  9 int main() 10 { 11     Point a,b,c,p;    //定义四点

12     while(scanf("%lf%lf%lf%lf%lf%lf%lf%lf",&a.x,&a.y,&b.x,&b.y,&c.x,&c.y,&p.x,&p.y)!=EOF){    //输入四对坐标

13         if( (Cross(a,b,p)>0 && Cross(b,c,p)>0 && Cross(c,a,p)>0) || 

14             (Cross(a,b,p)<0 && Cross(b,c,p)<0 && Cross(c,a,p)<0) )    //判断点p是否在三角形abc内。注意不能只考虑一种情况(<0 或 >0)

15             printf("Die\n"); 16         else 

17             printf("Live\n"); 18  } 19     return 0; 20 }

 

Freecode : www.cnblogs.com/yym2013