关于几何算法的一个简单问题 [问题点数：20分，结帖人jwd_1_cool]
-
本版专家分：0结帖率 100% -
回复次数7楼主谁能给我贴个
判断点是否在多边形内或上的c/c++程序啊？
另：如何将输入的多边形点，（可能是顺时针，也可能是逆时针顺序）转化到一个数组里按逆时针排列。
参考 zju1081 ,这道题过的人不少，若能相告，3x问题点数：20分
-
-
CSDN今日推荐
-
本版专家分：0
-
- 蓝花 2005年12月 专题开发/技术/项目大版内专家分月排行榜第三
-
-
#1 得分：15给你
//#include <iostream>
#include <fstream>
#include <cmath>
using namespace std;
#define INF 1E200 //无穷大
#define EP 1E-10 //无穷小
ifstream cin("1081.in");
ofstream cout("1081.out");
struct POINT
{
double x;
double y;
}poly[100];
struct LINESEG
{
POINT s;
POINT e;
};
double max(double a,double b)
{
return a>b?a:b;
}
double min(double a,double b)
{
return a<b?a:b;
}
double multiply(POINT sp,POINT ep,POINT op)
{
return((sp.x-op.x)*(ep.y-op.y)-(ep.x-op.x)*(sp.y-op.y));
}
bool online(LINESEG l,POINT p)
{
return( (multiply(l.e,p,l.s)==0) &&( ( (p.x-l.s.x)*(p.x-l.e.x)<=0 )&&( (p.y-l.s.y)*(p.y-l.e.y)<=0 ) ) );
}
/* 如果线段u和v相交(包括相交在端点处)时，返回true */
bool intersect(LINESEG u,LINESEG v)
{
return( (max(u.s.x,u.e.x)>=min(v.s.x,v.e.x))&& //排斥实验
(max(v.s.x,v.e.x)>=min(u.s.x,u.e.x))&&
(max(u.s.y,u.e.y)>=min(v.s.y,v.e.y))&&
(max(v.s.y,v.e.y)>=min(u.s.y,u.e.y))&&
(multiply(v.s,u.e,u.s)*multiply(u.e,v.e,u.s)>=0)&& //跨立实验
(multiply(u.s,v.e,v.s)*multiply(v.e,u.e,v.s)>=0));
}
/* (线段u和v相交)&&(交点不是双方的端点) 时返回true */
bool intersect_A(LINESEG u,LINESEG v)
{
return((intersect(u,v))&&
(!online(u,v.s))&&
(!online(u,v.e))&&
(!online(v,u.e))&&
(!online(v,u.s)));
}
/* 射线法判断点q与多边形polygon的位置关系，要求polygon为简单多边形，顶点按逆时针排列
如果点在多边形内： 返回0
如果点在多边形边上：返回1
如果点在多边形外： 返回2 */
int insidepolygon(int vcount,POINT Polygon[],POINT q)
{
int c=0,i,n;
LINESEG l1,l2;
bool bintersect_a,bonline1,bonline2,bonline3;
double r1,r2;
l1.s=q;
l1.e=q;
l1.e.x=double(INF);
n=vcount;
for (i=0;i<vcount;i++)
{
l2.s=Polygon[i];
l2.e=Polygon[(i+1)%n];
if(online(l2,q))return 1; //如果点在边上，返回1
if ( (bintersect_a=intersect_A(l1,l2))|| //相交且不在端点
(
(bonline1=online(l1,Polygon[(i+1)%n]))&& //第二个端点在射线上
(
(!(bonline2=online(l1,Polygon[(i+2)%n])))&& //前一个端点和后一个端点在射线两侧
((r1=multiply(Polygon[i],Polygon[(i+1)%n],l1.s)*multiply(Polygon[(i+1)%n],Polygon[(i+2)%n],l1.s))>0)
||
(bonline3=online(l1,Polygon[(i+2)%n]))&& //下一条边是水平线，前一个端点和后一个端点在射线两侧
((r2=multiply(Polygon[i],Polygon[(i+2)%n],l1.s)*multiply(Polygon[(i+2)%n],Polygon[(i+3)%n],l1.s))>0)
)
)
) c++;
}
if(c%2 == 1)
return 0;
else
return 2;
}
int N,M;
int main()
{
int cn=1;
while(cin>> N && N)
{
if(cn > 1)cout<<endl;
cout<<"Problem "<<cn<<':'<<endl;cn++;
cin>>M;
for(int i=0;i<N;i++)
{
cin>>poly[i].x>>poly[i].y;
}
POINT q;
for(int i=0;i<M;i++)
{
cin>>q.x>>q.y;
if(insidepolygon(N,poly,q) < 2)
{
cout<<"Within\n";
}
else
{
cout<<"Outside\n";
}
}
}
return 0;
}
-
本版专家分：0 -
#2 得分：5超超的速度好快啊。刚要re就回复了....
ZOJ 1081以前讨论过了
http://search.csdn.net/expert/topic/6/603/2002/9/11/1013280.htm
-
本版专家分：0
-
#3 得分：0哈哈，mmmcd(超超) 动作真的好快啊
判断点是否在多边行内，简单得说就是划射线，从这个点向左或又划一条平行的与x轴的射线，用射线和多边行的边去判断交点，如果交点个数是奇数的话，说明点在多边行内了
-
本版专家分：0
-
#4 得分：0晕翻，kbsoft(景乐) 你给的那个贴子是什么时候的啊，小虎竟然说自己是第一次做计算几何的题目？？
这道题目 典型的很，楼上这位说的我也知道。
我WA了，以后判断了一下，可能又两个原因：
1。 使用int 后面叉乘精度出错，改成double 了
2。 也是我疑惑的，题目输入数据没说是按什么顺序啊，可是象 超超 写的函数里 Polygon[],必须是逆时针的， 这怎么说？？？？
-
本版专家分：0
-
#5 得分：0我怎么也看不出我和超超的程序有什么区别，但就是wa
我也试过copy超超的小函数，但还是不行，现在将的未改的最初程序，贴上来，实在是痛苦啊
#include<iostream.h>
#include<math.h>
int n,m;
#define infinity 1e200
#define EP 1e-10
struct TPoint{
double x,y;
};
struct TLineSeg{
TPoint a,b;
};
double max(double a,double b)
{
return a>b?a:b;
}
double min(double a,double b)
{
return a<b?a:b;
}
double multiply(TPoint p1,TPoint p2,TPoint p0)
{
return((p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y));
}
//确定两条线段是否相交
bool intersect(TLineSeg u,TLineSeg v)
{
return( (max(u.a.x,u.b.x)>=min(v.a.x,v.b.x))&&
(max(v.a.x,v.b.x)>=min(u.a.x,u.b.x))&&
(max(u.a.y,u.b.y)>=min(v.a.y,v.b.y))&&
(max(v.a.y,v.b.y)>=min(u.a.y,u.b.y))&&
(multiply(v.a,u.b,u.a)*multiply(u.b,v.b,u.a)>=0)&&
(multiply(u.a,v.b,v.a)*multiply(v.b,u.b,v.a)>=0));
}
//判断点p是否在线段l上 不包括端点
bool online(TLineSeg l,TPoint p)
{
return( (multiply(l.b,p,l.a)==0)&&( ((p.x-l.a.x)*(p.x-l.b.x)<0 )||( (p.y-l.a.y)*(p.y-l.b.y)<0 )) );
}
//包括端点
bool online2(TLineSeg l,TPoint p)
{
return( (multiply(l.b,p,l.a)==0)&&( ((p.x-l.a.x)*(p.x-l.b.x)<=0 )||( (p.y-l.a.y)*(p.y-l.b.y)<=0 )) );
}
//判断两个点是否相等
bool Euqal_Point(TPoint p1,TPoint p2)
{
return((fabs(p1.x-p2.x)<EP)&&(fabs(p1.y-p2.y)<EP));
}
//一种线段相交判断函数，当且仅当u,v相交并且交点不是u,v的端点时函数为true;
bool intersect_A(TLineSeg u,TLineSeg v)
{
return((intersect(u,v))&&
(!Euqal_Point(u.a,v.a))&&
(!Euqal_Point(u.a,v.b))&&
(!Euqal_Point(u.b,v.a))&&
(!Euqal_Point(u.b,v.b)));
}
/*===============================================
判断点q是否在多边形Polygon内，
其中多边形是任意的凸或凹多边形，
Polygon中存放多边形的逆时针顶点序列
================================================*/
int InsidePolygon(int vcount,TPoint Polygon[],TPoint q)
{
int c=0,i,n;
TLineSeg l1,l2;
l1.a=q;
l1.b=q;
l1.b.x=infinity;
n=vcount;
for (i=0;i<vcount;i++)
{
l2.a=Polygon[i];
l2.b=Polygon[(i+1)%n];
if(online2(l2,q))return 1;
if ( (intersect_A(l1,l2))||
(
(online(l1,Polygon[(i+1)%n]))&&
(
(!online(l1,Polygon[(i+2)%n]))&&
(multiply(Polygon[i],Polygon[(i+1)%n],l1.a)*multiply(Polygon[(i+1)%n],Polygon[(i+2)%n],l1.a)>0)
||
(online(l1,Polygon[(i+2)%n]))&&
(multiply(Polygon[i],Polygon[(i+2)%n],l1.a)*multiply(Polygon[(i+2)%n],Polygon[(i+3)%n],l1.a)>0)
)
)
) c++;
}
return(c%2!=0);
}
TPoint p[100];
int main()
{
int cs(1);
bool first=true;
while(cin>>n){
if(n==0)return 0;
if(!first) cout<<endl;
cout<<"Problem "<<cs++<<':'<<endl;
cin>>m;
TPoint t;
for(int i=0;i<n;i++){
cin>>t.x>>t.y;
p[i]=t;
}
for(int i=0;i<m;i++){
cin>>t.x>>t.y;
if(InsidePolygon(n,p,t))
cout<<"Within"<<endl;
else cout<<"Outside"<<endl;
}
first=false;
}
return 1;
}
-
本版专家分：0
-
- 铜牌 2002年7月 总版技术专家分月排行榜第三
-
-
红花
2002年7月 PHP大版内专家分月排行榜第一
2002年5月 PHP大版内专家分月排行榜第一
2002年7月 C/C++大版内专家分月排行榜第一
-
-
蓝花
2003年1月 PHP大版内专家分月排行榜第三
2002年7月 专题开发/技术/项目大版内专家分月排行榜第三
-
-
#6 得分：0呵呵, 好久没来转转了, 1081没写总结, 不然可以帖出来给大家看看:)
前段时间又写了一个角度和法的, 感觉比射线法简单得多了. ^^
-
本版专家分：0
-
#7 得分：0继续
- 其他相关推荐
- 简单的图论问题(Dijkstra)
- O - 简单的图论问题? Time Limit:5000MS Memory Limit:65535KB 64bit IO Format: Submit Status Description 给一个 n 行 m 列的迷宫,每个格子要么是障碍物要么是空地。每个空地里都有一个权值。你的 任务是从找一条(r1,c1)到(r2,c2
- 计算几何较全模板（有包含三维的）
- 转自:http://xiaobm.ycool.com/post.2715020.html #include #include #include #include #define MAX_N 100 using namespace std; ///////////////////////////////////////////////////////////////////
- 一些简单的java编程题(3) ————小球落地问题
- 版权声明：本文为博主原创文章，未经博主允许不得转载。 /* (程序头部注释开始) 版权声明：保留个人权利。 如程序有不当之处，敬请指正。 * 文件名称： 《一些简单的java编程题(3) ————小球落地问题》 * 作 者： 那就去巴黎 * 完成日期： 2018 年 01 月 30 日 * 版 本 号： J
- 关于求余数问题的一个简单方法.doc
- 关于求余数问题的一个简单方法
- 计算几何算法C++实现
- 用C++语言实现的各种计算几何算法, 如判断点是否在多边形内, 多边形凸凹性判断, 寻找点集凸包的卷包裹法,求简单多边形的重心,求从多边形外一点出发到该多边形的切线,判断线段是否相交等等!
- 计算几何题目集锦
- 其实也谈不上推荐，只是自己做过的题目而已，甚至有的题目尚未AC，让在挣扎中。之所以推荐计算几何题，是因为，本人感觉ACM各种算法中计算几何算是比较实际的算法，在很多领域有着重要的用途 计算几何题的特点与做题要领： 1.大部分不会很难，少部分题目思路很巧妙 2.做计算几何题目，模板很重要，模板必须高度可靠。 3.要注意代码的组织，因为计算几何的题目很容易上两百行代码，里面大部分是模板。如果代码
- 计算机几何常用算法
- 1、矢量减法 设二维矢量 P = （x1,y1） ，Q = (x2,y2) 则矢量减法定义为： P - Q = ( x1 - x2 , y1 - y2 ) 显然有性质 P - Q = - ( Q - P ) 如不加说明，下面所有的点都看作矢量，两点的减法就是矢量相减； 2、矢量叉积 设矢量P = （x1,y1） ，Q = (x2,y2) 则矢量叉积定义为： P × Q =
- 开源几何计算数学库
- a) Geos Geos的全称就是Geometry Engine Open Source，参考网站就是：http://geos.refractions.net，从名称就可以看出，Geos就是判断几何体的过程。 Geos是用C++语言写的处理计算几何的开源库。 b) CGAL CGAL的全称就是Computational Geometry Algorithms Library
- 一些三维计算几何算法
- //三维几何函数库 #include #define eps 1e-8 #define zero(x) (((x)>0?(x):-(x))<eps) struct point3{double x,y,z;}; struct line3{point3 a,b;}; struct plane3{point3 a,b,c;}; //计算cross product U x V point3 xmult
- 多态的一个简单例子
- package com; public class Animal { public int a; public void move() { System.out.println("移动"); } public static void main(String[] args) { // 多态调用时，父类变量指向子类对象，所能调用的方法和属性仅限于父类变量中定义