标签:ted += 时间 ons img 现在 高度 tle max
Time Limit: 3000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 3105 Accepted Submission(s): 1103
1 /** 2 Finished time:9.15.2018 3 Author:hyacinth 4 Algorithm:计算几何判断点是否在多边形内 5 */ 6 #include<iostream> 7 #include<cstdio> 8 #include<cmath> 9 using namespace std; 10 const double esp=1e-8; 11 12 struct point 13 { 14 double x,y; 15 }; 16 point p[105]; 17 int n,m; 18 19 double k(point p1,point p2)//斜率 20 { 21 return (p2.y-p1.y)/(p2.x-p1.x); 22 } 23 bool onEdge(point q,point p1,point p2)//判断q是否在p1p2上 24 { 25 if(q.x >= min(p1.x,p2.x) && q.x <= max(p1.x,p2.x))//去除q在p1p2延长线上的情况 26 if(fabs(k(q,p1)-k(p2,p1)) <= esp)//通过斜率是否相等来判断q是否在p1p2上 27 return true; 28 return false; 29 } 30 bool inPolaygon(point q) 31 { 32 int cnt=0; 33 point p1,p2; 34 p1=p[0]; 35 for(int i=1;i <= n;++i) 36 { 37 p2=p[i%n]; 38 if(onEdge(q,p1,p2)) 39 return true; 40 if(p1.y != p2.y){//判断p1p2是否为水平边,水平边不作考虑 41 if(q.y > min(p1.y,p2.y) && q.y <= max(p1.y,p2.y) && q.x <= max(p1.x,p2.x)){//注意:此处q.y是严格>min(),此处是用来处理射线与顶点相交的情况的,当相交的顶点为凸顶点时,两条边都加上,当相交的顶点是凹顶点时,只加右边那条边 42 if(p1.x == p2.x)//当p1p2为竖线时,由条件q.x <= max(p1.x,p2.x)可得由p向右发出的射线与p1p2有交点 43 cnt++; 44 else{ 45 double x=p1.x+(q.y-p1.y)/k(p2,p1);//x : 经过q点且平行于x轴的直线与p1p2交点的横坐标 46 cnt += (q.x <= x ? 1:0);//如果q.x <= x,说明q向右发出的射线与p1p2有交点 47 } 48 } 49 } 50 p1=p2; 51 } 52 return cnt&1; 53 } 54 55 int main() 56 { 57 while(~scanf("%d",&n)) 58 { 59 for(int i=0;i < n;++i) 60 scanf("%lf%lf",&p[i].x,&p[i].y); 61 scanf("%d",&m); 62 for(int i=1;i <= m;++i) 63 { 64 point q; 65 scanf("%lf%lf",&q.x,&q.y); 66 if(inPolaygon(q)) 67 printf("Yes\n"); 68 else 69 printf("No\n"); 70 } 71 } 72 return 0; 73 }
标签:ted += 时间 ons img 现在 高度 tle max
原文地址:https://www.cnblogs.com/violet-acmer/p/9651799.html