线段交点

2024-01-24 13:58
文章标签 交点 线段

本文主要是介绍线段交点,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

线段交点

题目描述:
给定N个线段。求线段交点数。
保证没有两条线段共线
输入:
一行一个整数N,表示线段的个数
第2~N+1行,每行四个实数,x1,y1,x2,y2,表示线段的两个端点(x1,y1)和(x2,y2)
输出:
一行一个整数,表示线段交点数。
样例输入1:
3
0.00 0.00 1.00 1.00
0.00 1.00 1.00 0.00
0.00 0.00 1.00 0.00
样例输出1:
3
样例输入2:
3
0.00 0.00 1.00 1.00
0.00 0.00 1.00 0.00
0.00 0.00 0.00 1.00
样例输出2:
1
(没错这就是上一个题改了一下,其实是我把之前题看错了做的解法,然后,emmm,改成这个题意,就OK了)
先把坐标都换成向量,然后计算交点。

double getCrossPointx(STU s1, STU s2)
{Point a;a.y=s2.y2-s2.y1;a.x=s2.x2-s2.x1;Point b;b.y=s1.y1-s2.y1;b.x=s1.x1-s2.x1;Point c;c.y=s1.y2-s2.y1;c.x=s1.x2-s2.x1;double d1 = fabs(waiji(a, b));double d2 = fabs(waiji(a, c));double t = d1 / (d1 + d2);return 1.0*s1.x1 + (s1.x2 - s1.x1) * t;
}
double getCrossPointy(STU s1, STU s2)
{Point a;a.y=s2.y2-s2.y1;a.x=s2.x2-s2.x1;Point b;b.y=s1.y1-s2.y1;b.x=s1.x1-s2.x1;Point c;c.y=s1.y2-s2.y1;c.x=s1.x2-s2.x1;double d1 = fabs(waiji(a, b));double d2 = fabs(waiji(a, c));double t = d1 / (d1 + d2);return 1.0*s1.y1 + (s1.y2 - s1.y1) * t;
}

跟前面那个结合一下,就变成了这样:

#include<bits/stdc++.h>
using namespace std;
struct Point//向量
{double x;double y;
};
typedef struct
{double x1;double y1;double x2;double y2;Point p1;Point p2;
}STU;
double maxx(double a,double b)
{if(a>b) return a;else return b;
}
double minn(double a,double b)
{if(a>b) return b;else return a;
}
bool judge1(STU a,STU b)
{if(maxx(a.x1,a.x2)<minn(b.x1,b.x2)||maxx(a.y1,a.y2)<minn(b.y1,b.y2)||minn(a.x1,a.x2)>maxx(b.x1,b.x2)||minn(a.y1,a.y2)>maxx(b.y1,b.y2)){return false;}else return true;
}
bool judge2(STU a,STU b)
{if((1.0*(a.x1-b.x1)*(b.y2-b.y1)-1.0*(a.y1-b.y1)*(b.x2-b.x1))*(1.0*(a.x2-b.x1)*(b.y2-b.y1)-1.0*(a.y2-b.y1)*(b.x2-b.x1))>0.000001){return false;}else if((1.0*(b.x1-a.x1)*(a.y2-a.y1)-1.0*(b.y1-a.y1)*(a.x2-a.x1))*(1.0*(b.x2-a.x1)*(a.y2-a.y1)-1.0*(b.y2-a.y1)*(a.x2-a.x1))>0.000001){return false;}else return true;
}
bool xiangjiao(STU a,STU b)
{if(!judge1(a,b)) return false;else if(!judge2(a,b)) return false;else return true;
}
double waiji(Point a, Point b)
{return 1.0*a.x*b.y - a.y*b.x;
}
double getCrossPointx(STU s1, STU s2)//求交点横坐标
{Point a;a.y=s2.y2-s2.y1;a.x=s2.x2-s2.x1;Point b;b.y=s1.y1-s2.y1;b.x=s1.x1-s2.x1;Point c;c.y=s1.y2-s2.y1;c.x=s1.x2-s2.x1;double d1 = fabs(waiji(a, b));double d2 = fabs(waiji(a, c));double t = d1 / (d1 + d2);return 1.0*s1.x1 + (s1.x2 - s1.x1) * t;
}
double getCrossPointy(STU s1, STU s2)//求交点纵坐标
{Point a;a.y=s2.y2-s2.y1;a.x=s2.x2-s2.x1;Point b;b.y=s1.y1-s2.y1;b.x=s1.x1-s2.x1;Point c;c.y=s1.y2-s2.y1;c.x=s1.x2-s2.x1;double d1 = fabs(waiji(a, b));double d2 = fabs(waiji(a, c));double t = d1 / (d1 + d2);return 1.0*s1.y1 + (s1.y2 - s1.y1) * t;
}
int main()
{int n,sum=0;cin>>n;STU stu[n+1];for(int i=1;i<=n;i++){scanf("%lf%lf%lf%lf",&stu[i].x1,&stu[i].y1,&stu[i].x2,&stu[i].y2);stu[i].p1.x=stu[i].x1;stu[i].p1.y=stu[i].y1;stu[i].p2.x=stu[i].x2;stu[i].p2.y=stu[i].y2;}int u=0,v=0;Point po[100000];for(int i=1;i<n;i++){for(int j=i+1;j<=n;j++){if(xiangjiao(stu[i],stu[j])) sum++;Point p;p.x=getCrossPointx(stu[i],stu[j]);p.y=getCrossPointy(stu[i],stu[j]);for(int k=0;k<u;k++)//判断是否重复{if(p.x==po[k].x&&p.y==po[k].y){v=1;break;}}if(v==0){po[u].x=p.x;po[u].y=p.y;u++;}else sum--;//重复,则不应该+1,但前面已经加了,故在此处减去}}printf("%d\n",sum);return 0;
}

这篇关于线段交点的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/639911

相关文章

poj3468(线段树成段更新模板题)

题意:包括两个操作:1、将[a.b]上的数字加上v;2、查询区间[a,b]上的和 下面的介绍是下解题思路: 首先介绍  lazy-tag思想:用一个变量记录每一个线段树节点的变化值,当这部分线段的一致性被破坏我们就将这个变化值传递给子区间,大大增加了线段树的效率。 比如现在需要对[a,b]区间值进行加c操作,那么就从根节点[1,n]开始调用update函数进行操作,如果刚好执行到一个子节点,

hdu1394(线段树点更新的应用)

题意:求一个序列经过一定的操作得到的序列的最小逆序数 这题会用到逆序数的一个性质,在0到n-1这些数字组成的乱序排列,将第一个数字A移到最后一位,得到的逆序数为res-a+(n-a-1) 知道上面的知识点后,可以用暴力来解 代码如下: #include<iostream>#include<algorithm>#include<cstring>#include<stack>#in

hdu1689(线段树成段更新)

两种操作:1、set区间[a,b]上数字为v;2、查询[ 1 , n ]上的sum 代码如下: #include<iostream>#include<algorithm>#include<cstring>#include<stack>#include<queue>#include<set>#include<map>#include<stdio.h>#include<stdl

hdu 1754 I Hate It(线段树,单点更新,区间最值)

题意是求一个线段中的最大数。 线段树的模板题,试用了一下交大的模板。效率有点略低。 代码: #include <stdio.h>#include <string.h>#define TREE_SIZE (1 << (20))//const int TREE_SIZE = 200000 + 10;int max(int a, int b){return a > b ? a :

hdu 1166 敌兵布阵(树状数组 or 线段树)

题意是求一个线段的和,在线段上可以进行加减的修改。 树状数组的模板题。 代码: #include <stdio.h>#include <string.h>const int maxn = 50000 + 1;int c[maxn];int n;int lowbit(int x){return x & -x;}void add(int x, int num){while

poj 1127 线段相交的判定

题意: 有n根木棍,每根的端点坐标分别是 px, py, qx, qy。 判断每对木棍是否相连,当他们之间有公共点时,就认为他们相连。 并且通过相连的木棍相连的木棍也是相连的。 解析: 线段相交的判定。 首先,模板中的线段相交是不判端点的,所以要加一个端点在直线上的判定; 然后,端点在直线上的判定这个函数是不判定两个端点是同一个端点的情况的,所以要加是否端点相等的判断。 最后

HDU4737线段树

题目大意:给定一系列数,F(i,j)表示对从ai到aj连续求或运算,(i<=j)求F(i,j)<=m的总数。 const int Max_N = 100008 ;int sum[Max_N<<2] , x[Max_N] ;int n , m ;void push_up(int t){sum[t] = sum[t<<1] | sum[t<<1|1] ;}void upd

zoj 1721 判断2条线段(完全)相交

给出起点,终点,与一些障碍线段。 求起点到终点的最短路。 枚举2点的距离,然后最短路。 2点可达条件:没有线段与这2点所构成的线段(完全)相交。 const double eps = 1e-8 ;double add(double x , double y){if(fabs(x+y) < eps*(fabs(x) + fabs(y))) return 0 ;return x + y ;

圆与线段的交点

poj 3819  给出一条线段的两个端点,再给出n个圆,求出这条线段被所有圆覆盖的部分占了整条线段的百分比。 圆与线段的交点 : 向量AB 的参数方程  P = A + t * (B - A)      0<=t<=1 ; 将点带入圆的方程即可。  注意: 有交点 0 <= t <= 1 ; 此题求覆盖的部分。 则 若求得 t  满足 ; double ask(d

Codeforces 482B 线段树

求是否存在这样的n个数; m次操作,每次操作就是三个数 l ,r,val          a[l] & a[l+1] &......&a[r] = val 就是区间l---r上的与的值为val 。 也就是意味着区间[L , R] 每个数要执行 | val 操作  最后判断  a[l] & a[l+1] &......&a[r] 是否= val import ja