链接:https://ac.nowcoder.com/acm/contest/3570/I

来源:牛客网

题目描述

最近ZSC和他的女朋友NULL吵架了。因为ZSC是一个直男,他不知道该怎么办,于是他寻求恋爱之子磊子哥的帮助。“比起磊子,我更需要女朋友”/doge。

矛盾就如一条条线,纠在一起,越来越乱,ZSC和NULL的矛盾可以看作二维直角坐标系,在平面中有N条线段,线段的两个端点的坐标分别是 pi_x,pi_y和qi_x,qi_y。当两个线段有至少一个公共点时,就认为它们是相连的。事情往往都有关联,所以通过相连的线段连在一起的两条线段也是相连的。

恋爱之子磊子哥每次祈祷都可以消除一条线段以及所有与它相连的线段,自从磊子哥成为恋爱之子后特别忙,现在他向你求助“那你能帮帮我吗,我想知道我最少要祈祷几次才能消除所有线段”。学过算法的你,能帮助磊子哥求出最少的次数吗?

输入描述:

输入第一行是一个整数N(1≤N≤4∗103),表示线段的个数,接下来N行,每行四个整数pi_x,pi_y和qi_x,qi_y表示该线段的两个端点坐标。

输出描述: 输出最少需要的祈祷数。

  1. 示例1
  2. 输入
  3. 复制
  4. 4
  5. 0 0 6 0
  6. 6 -4 6 4
  7. 0 4 4 4
  8. 2 2 2 6
  9. 输出
  10. 2

说明

  1. 示例2
  2. 输入
  3. 复制
  4. 3
  5. 0 1 1 0
  6. 0 2 2 0
  7. 0 3 3 0
  8. 输出
  9. 复制
  10. 3
  11. 说明

  1. 示例3
  2. 输入
  3. 复制
  4. 29
  5. 1 8 2 9
  6. 2 9 11 9
  7. 11 9 9 8
  8. 9 8 1 8
  9. 1 6 2 7
  10. 2 7 2 8
  11. 3 6 5 8
  12. 5 6 6 7
  13. 6 7 6 8
  14. 7 6 9 8
  15. 7 7 8 8
  16. 1 3 1 5
  17. 1 5 3 5
  18. 3 5 3 3
  19. 3 3 1 3
  20. 2 1 2 3
  21. 0 2 2 2
  22. 2 2 4 2
  23. 1 -1 2 1
  24. 3 -1 2 1
  25. 5 3 5 5
  26. 5 5 7 5
  27. 7 5 7 3
  28. 7 3 5 3
  29. 6 1 6 3
  30. 4 2 6 2
  31. 6 2 8 2
  32. 5 -1 6 1
  33. 7 -1 6 1
  34. 输出
  35. 复制
  36. 2

说明

思路如下

题意:题目给出 n 条直线的端点,通过每次次的祈祷,可以把其中相交的线段都去掉,问我们需要祈祷多少次才可以,把所有线段都消除,那么我们就要考虑两个问题:怎么判断直线相交、怎么一次把所有相交的直线都消去。对于前一个问题我也没看懂到大佬的题解,但对于第二个问题,我们可以明确的看出可以并查集来解决,通过判断线段相交,把所有的相交的线段视为一个集合,这样有几个集合我们就需要祈祷 几次

并查集总结传送门

题解如下(贴出大佬的代码一起共赏)

  1. #include<set>
  2. #include<map>
  3. #include<cmath>
  4. #include<ctime>
  5. #include<queue>
  6. #include<stack>
  7. #include<cstdio>
  8. #include<string>
  9. #include<vector>
  10. #include<cstdlib>
  11. #include<cstring>
  12. #include<iomanip>
  13. #include<iostream>
  14. #include<algorithm>
  15. #define fi first
  16. #define se second
  17. #define db double
  18. #define gcd __gcd
  19. #define pb push_back
  20. #define lowbit(x) (x & (-x))
  21. #define PII pair<int, int>
  22. #define all(x) x.begin(), x.end()
  23. #define debug(x) cout << #x << " = " << x << endl
  24. #define rep(i, a, b) for(__typeof(b) i = a; i <= (b); i++)
  25. #define Rep(i, a, b) for(__typeof(a) i = a; i >= (b); i--)
  26. #define Mod(a,b) a<b?a:a%b+b
  27. template<class T> inline T qmin(T a, T b) { return a < b ? a : b; }
  28. template<class T> inline T qmax(T a, T b) { return a > b ? a : b; }
  29. typedef long long ll;
  30. typedef unsigned long long ull;
  31. const db PI = acos(-1);
  32. const int inf = 0x3f3f3f3f;
  33. const int mod = (int)1e9 + 1;
  34. const int maxn = (int)1e5 + 5;//remember to modify it, No RE&nbs***bsp;MLE
  35. const ll INF = 0x3f3f3f3f3f3f3f3f;
  36. using namespace std;
  37. #define cross(p1, p2, p3) ((p2.x - p1.x) * (p3.y - p1.y) - (p3.x - p1.x) * (p2.y - p1.y))
  38. #define cross0p(p1, p2, p3) sign(cross(p1,p2,p3))
  39. const db eps = 1e-9;
  40. inline int sign(db a) { return a < -eps ? -1: a > eps; }
  41. inline int cmp(db a, db b) { return sign(a-b); }
  42. struct P{
  43. db x, y;
  44. P() {}
  45. P(db _x, db _y): x(_x), y(_y) {}
  46. P operator+(P p) { return {x + p.x, y + p.y}; }
  47. P operator-(P p) { return {x - p.x, y - p.y}; }
  48. P operator*(db d) { return {x * d, y * d}; }
  49. P operator/(db d) { return {x / d, y / d}; }
  50. bool operator<(P p)const {
  51. int c = cmp(x, p.x);
  52. if(c) return c == -1;
  53. return cmp(y, p.y) == -1;
  54. }
  55. bool operator==(P o) const {
  56. return cmp(x, o.x) == 0 && cmp(y, o.y) == 0;
  57. }
  58. db distTo(P p) { return (*this - p).abs(); }
  59. db abs() { return sqrt(abs2()); }
  60. db abs2() { return x * x + y * y; }
  61. db dot(P p) { return x * p.x + y * p.y; }
  62. db det(P p) { return x * p.y - y * p.x; }
  63. };
  64. bool intersect(db l1, db r1, db l2, db r2)
  65. {
  66. if(l1 > r1)swap(l1, r1); if(l2 > r2) swap(l2, r2);
  67. return !( cmp(r1, l2) == -1 || cmp(r2, l1) == -1);
  68. }
  69. bool isSS(P p1, P p2, P q1, P q2){
  70. return intersect(p1.x, p2.x, q1.x, q2.x) && intersect(p1.y, p2.y, q1.y, q2.y) && cross0p(p1, p2, q1) * cross0p(p1, p2, q2) <= 0 && cross0p(q1, q2, p1) * cross0p(q1, q2, p2) <= 0;
  71. }
  72. //并查集开始
  73. int f[maxn];
  74. int find(int x) //递归查找x元素的父节点
  75. {
  76. if(f[x] == x)
  77. return x;
  78. return f[x] = find(f[x]); //注意这里在查找父节点的时候 已经把节点的路径进行了压缩(进行了优化)
  79. }
  80. void Union(int x,int y)
  81. {
  82. int fx = find(x),fy = find(y);
  83. if(fx != fy) f[fx] = fy; //x,y是是有关系的,但是他们的父节点不同所以要 把其中的一个赋接待你作为子节点 拼接到临一个父节点上
  84. } //这样x、 y 就有了相同的父节点,就属于同一个集合了
  85. int main()
  86. {
  87. int n;
  88. P p[10005], q[10005];
  89. scanf("%d", &n);
  90. for(int i = 1 ; i <= n ; i++){
  91. f[i] = i; //并查集元素的初始化
  92. scanf("%lf %lf %lf %lf", &p[i].x, &p[i].y, &q[i].x, &q[i].y);
  93. }
  94. for(int i = 1 ; i <= n ; i++){
  95. for(int j = i + 1 ; j <= n ; j++){
  96. if(isSS(p[i], q[i], p[j], q[j]))Union(i, j);//如果为真i,j 是有关系的,所以他们应该在同一个集合中
  97. }
  98. }
  99. int ans = 0;
  100. for(int i = 1 ; i <= n ; i++)if(f[i] == i)ans++;
  101. printf("%d\n", ans);
  102. return 0;
  103. }

I、恋爱之子的更多相关文章

  1. XXOOJL

    她的鞋子放在外面没拿进来很显眼,我们俩正抱在床上,刚做完什么也没穿,她也没擦.听到门外的脚步声.赶紧穿.她在我房间里没出声,但我父亲肯定看出来了.于是问我:母亲去哪了,然后他去找她. 太尴尬了,那晚我 ...

  2. 【0门槛】PR稿的自我修养

    本文来自网易云社区 作者:巩爽 十一过完,离2018年结束就只剩下85天啦!是不是2016年许下的2017年的梦想,在2018年还没有实现? 做过的项目仿佛都小有成就,可惜只是内部自嗨,想做域外宣传却 ...

  3. 二狗子 、初恋及HTTPS

    最近二狗子宅在老家,最悠闲的就是泡壶茶看着院子的风景发呆一下午.今天,二狗子看到了对面自己暗恋的小翠花,看着美好的小翠花二狗子不禁想起了自己美好的初恋. 二狗子的初恋在初中,那个时候学校禁止带手机.上 ...

  4. 深入理解MySql子查询IN的执行和优化

    IN为什么慢? 在应用程序中使用子查询后,SQL语句的查询性能变得非常糟糕.例如: SELECT driver_id FROM driver where driver_id in (SELECT dr ...

  5. Hawk 6. 高级话题:子流程系统

    子流程的定义 当流程设计的越来越复杂,越来越长时,就难以进行管理了.因此,采用模块化的设计才会更加合理.本节我们介绍子流程的原理和使用. 所谓子流程,就是能先构造出一个流程,然后被其他流程调用.被调用 ...

  6. margin折叠-从子元素margin-top影响父元素引出的问题

    正在做一个手机端电商项目,顶部导航栈的布局是一个div包含一个子div,如果给在正常文档流中的子div一个垂直margin-top,神奇的现象出现了,两父子元素的边距没变,但父div跟着一起往下走了! ...

  7. DevExpress第三方控件使用实例之ASPxPopupControl弹出子窗体

    弹出页面控件:ASPxPopupControl, <dxpc:ASPxPopupControl ID="popubCtr" runat="server" ...

  8. Mysql - 性能优化之子查询

    记得在做项目的时候, 听到过一句话, 尽量不要使用子查询, 那么这一篇就来看一下, 这句话是否是正确的. 那在这之前, 需要介绍一些概念性东西和mysql对语句的大致处理. 当Mysql Server ...

  9. 面向组合子设计Coder

    面向组合子 面向组合子(Combanitor-Oriented),是最近帮我打开新世界大门的一种pattern.缘起haskell,又见monad与ParseC,终于ajoo前辈的几篇文章. 自去年9 ...

随机推荐

  1. 怎么用一行HTML代码搭建一个谷歌镜像?

    <iframe src="google.com"></iframe>

  2. Echarts 自定义legend图片,修改点击之后的颜色图解

    第一个问题:echarts 可以自定义图例的图标,百度上很多回答都是引用的相对路径,但是不知道为啥,我的vue项目就是引用不显示,在network里面找不到相应图片 后来我想了个法子,就是先获取到这个 ...

  3. href="#"和href=“javascript:void(0)”的区别

    void是javascript中的关键字,该操作符指定要计算一个表达式但是不返回值. <a href="javascript:void(0);">点我没有反应的!< ...

  4. js数组方法全

    js数组方法大全 一:前言          转载   作者:九夏       出处:https://www.cnblogs.com/jiuxia/ 我们在学到js中数组的时候,我们会接触到js中数组 ...

  5. windows7免费永久激活方法分享

    前言 我相信,这里肯定有看过我上一篇博客的同学. 我说了,为解决windows7激活问题,我会找一个比较好的方法. 首先先让大家看一看激活前windows7的计算机属性: 显示是未激活的.下面就是方法 ...

  6. 回想笔记 瞎比比 域名注册 解析绑定ip 下载证书 设置证书 重定向http到https请求

    2019.7.27 回想笔记 拥有腾讯云服务器一台 阿里云注册5元域名,进行备案 完成之后 使用解析 绑定服务器ip地址 ,使用域名可以访问到web服务器而不是通过直接暴露ip地址进行访问 证书购买 ...

  7. vue 接口 对象数组的传入传出

    数组传入 对象的话用{...obj} <selectFile :tableData="[...tableData]" @uploadSuccess="data =& ...

  8. Java第一节课考试

    1 package kaoshi; import java.util.Scanner; public class ScoreInformation { Scanner input=new Scanne ...

  9. Linux 中useradd命令的使用

    Linux 系统中通常都是root用户具有超级权限,超级用户root一般是不需要创建的,然而很多时候root用户不是任何人都可以使用的,毕竟最高权限的用户,任意使用的话,会对系统造成很多不必要的破坏. ...

  10. 图解汉诺塔问题( Java 递归实现)

    汉诺塔简介 最近在看数据结构和算法,遇到了一个非常有意思的问题--汉诺塔问题. 先看下百度百科是怎么定义汉诺塔的规则的: 汉诺塔(又称河内塔)问题是源于印度一个古老传说的益智玩具.大梵天创造世界的时候 ...