hdu1154(点是否在多边形内部、两线段交点)
来源:互联网 发布:淘宝买东西怎么收货 编辑:程序博客网 时间:2024/05/17 22:54
point intersection(line &u,line &v)//两线段之间的交点
{
point ret=u.a;
double t=((u.a.x-v.a.x)*(v.a.y-v.b.y) - (u.a.y-v.a.y)*(v.a.x-v.b.x))/((u.a.x-u.b.x)*(v.a.y-v.b.y)-(u.a.y-u.b.y)*(v.a.x-v.b.x));
ret.x+=(u.b.x-u.a.x)*t;
ret.y+=(u.b.y-u.a.y)*t;
return ret;
}
题意是一条直线与多边形相交,求相交的总长,最后没怎么搞懂
#include <iostream>
#include <string.h>
#include <algorithm>
#include <cmath>
#include <stdio.h>
using namespacestd;
#define maxn 1200
#define eps 1e-8
struct point
{
double x;
double y;
}po[maxn],rec[maxn];
struct line
{
point a;
point b;
}temp,temp1;
double xmulit(point &a,point &b,point &c)//差乘
{
return (a.x-b.x)*(a.y-c.y)-(a.y-b.y)*(a.x-c.x);
}
bool across(point &a,point &b,point &c,point &d)//直线ab和线段cd是否相交
{
double p=xmulit(a,b,c),p1=xmulit(a,b,d);
if( fabs(p1) <= eps ||fabs(p) <= eps )//c或者d在ab上
returntrue;
if( p*p1 < -eps )//c和d在ab两侧,不用判断a和b是不是在cd两侧,因为ab是直线line
returntrue;
returnfalse;
}
bool one_line(point &a,point &b,point &c,point &d)//直线ab和线段cd是否重合,相当于p==0&&p1==0
{
double p=xmulit(a,b,c),p1=xmulit(a,b,d);
if( fabs(p1) < eps &&fabs(p) < eps )
returntrue;
returnfalse;
}
bool is_equal(point &a,point &b)//判断点a和点b是否相等
{
return (fabs(a.x-b.x) <= eps) && (fabs(a.y-b.y) <=eps);
}
point intersection(line &u,line &v)//两线段之间的交点
{
point ret=u.a;
double t=((u.a.x-v.a.x)*(v.a.y-v.b.y) - (u.a.y-v.a.y)*(v.a.x-v.b.x))/((u.a.x-u.b.x)*(v.a.y-v.b.y)-(u.a.y-u.b.y)*(v.a.x-v.b.x));
ret.x+=(u.b.x-u.a.x)*t;
ret.y+=(u.b.y-u.a.y)*t;
return ret;
}
int n,m;
double dis(point &a,point &b)
{
returnsqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
int cmp(point a,point b)
{
if(fabs(a.x-b.x)<eps)
return a.y<b.y;
return a.x<b.x;
}
bool on_segment(point pi,point pj,point pk)//判断点pk时候在线段pi, pj上
{
if(xmulit(pi, pj, pk)==0)
{
if(pk.x>=min(pi.x,pj.x)&&pk.x<=max(pi.x,pj.x)&&pk.y>=min(pi.y,pj.y)&&pk.y<=max(pi.y,pj.y))
returntrue;
}
returnfalse;
}
bool segments_intersect(point p1,point p2,point p3,point p4)//判断线段是否相交
{
double d1=xmulit(p3,p4,p1);
double d2=xmulit(p3,p4,p2);
double d3=xmulit(p1,p2,p3);
double d4=xmulit(p1,p2,p4);
if(d1*d2<0&&d3*d4<0)
returntrue;
else if(d1==0&&on_segment(p3,p4,p1))
returntrue;
else if(d2==0&&on_segment(p3,p4,p2))
returntrue;
else if(d3==0&&on_segment(p1,p2,p3))
returntrue;
else if(d4==0&&on_segment(p1,p2,p4))
returntrue;
returnfalse;
}
int inpoto(point a)//判断点是否在多边形的内部
{
int i;
point b,c,d;
b.y=a.y;
b.x=1e15;//定义射线
int count=0;
for(i=0;i<n;i++)
{
c = po[i];
d = po[i +1];
if(on_segment(c,d,a))//该点在多边形的一条边上
return1;
if(abs(c.y-d.y)<eps)
continue;
if(on_segment(a,b,c))//和顶点相交的情况,如果y值较大则取
{
if(c.y>d.y)
count++;
}
elseif(on_segment(a,b,d))//和顶点相交的情况,如果y值较大则取
{
if(d.y>c.y)
count++;
}
elseif(segments_intersect(a,b,c,d))//和边相交
count++;
}
return count%2;//当L和多边形的交点数目C是奇数的时候,P在多边形内,是偶数的话P在多边形外。
}
point mid(point &a,point &b)
{
point c;
c.x=(a.x+b.x)/2;
c.y=(a.y+b.y)/2;
return c;
}
double find_ans()
{
int i,pos=0;
double ans=0;
po[n]=po[0];
for(i=0;i<n;i++)
{
if(one_line(temp.a,temp.b,po[i],po[i+1]))//线段在直线上
{
rec[pos++]=po[i];
rec[pos++]=po[i+1];
//printf("pos=%d\n",pos);
continue;
}
if(across(temp.a,temp.b,po[i],po[i+1]))//线段与直线相交
{
temp1.a=po[i],temp1.b=po[i+1];
rec[pos++]=intersection(temp,temp1);//两条线的交点
}
}
for(int i=0;i<pos;i++)
printf("%lf %lf\n",rec[i].x,rec[i].y);
printf("*****\n");
sort(rec,rec+pos,cmp);
for(int i=0;i<pos;i++)
printf("%lf %lf\n",rec[i].x,rec[i].y);
printf("~~~~~\n");
for(i=0;i<pos-1;i++)
{
if(inpoto(mid(rec[i],rec[i+1])))
ans+=dis(rec[i],rec[i+1]);
}
return ans;
}
int main()
{
int i;
while(scanf("%d %d",&n,&m))
{
if(m==0 &&m==0)
return0;
for(i=0;i<n;i++)
scanf("%lf%lf",&po[i].x,&po[i].y);
for(i=0;i<m;i++)
{
scanf("%lf%lf%lf%lf",&temp.a.x,&temp.a.y,&temp.b.x,&temp.b.y);
printf("%.3lf\n",find_ans());
}
}
return0;
}
- hdu1154(点是否在多边形内部、两线段交点)
- HDU1154求线段在多边形内长度
- ZOJ1081 Points Within(判断点是否在多边形内部)
- 判断一个点是否在多边形内部
- 判断一个点是否在多边形内部
- 判定一个点是否在多边形内部
- C判断点是否在多边形内部
- 判断一个点是否在多边形内部
- 判断点是否在多边形内部
- 判断点是否在多边形内部
- 判断点是否在多边形内部
- 判断点是否在多边形内部
- hdu1756判断点是否在多边形内部
- 判断点是否在多边形内部
- 判断点是否在多边形内部
- 17.8.19 校内赛 解题报告【求线段交点+凸包+求多边形面积】【判定点是否在多边形内】【二分答案+半平面交】
- 判断一个点是否在一个多边形内部
- 判断一个点是否在多边形内部的问题
- 三、ArcGIS API for JavaScript之地图与图层(2)
- maven pom.xml文件详解
- 【Java编程】调用方法判断一个数是否是素数
- POJ 3073 Spam 笔记
- HDU 1963 && POJ 2063
- hdu1154(点是否在多边形内部、两线段交点)
- Effective C++(第3版)pdf
- CountDownLatch使用、源码浅析
- 使用QSerialport读取串口数据并处理
- UESTC
- 文章标题
- 队列
- Java优先级队列简单实现
- HDU