CEOI2020 道路（Roads） Solution

直接来构造。

考虑扫描线。从左到右扫，考虑当前扫到了一个左端点，我们把这个左端点连到其他点上。

我们可以找到这个点下方离他最近的线段，并且记下每条线段上方在扫描线左侧且最靠右，与这条线段中间没有其他线段的点，然后直接把左端点连到这样的点上就行了。容易证明这样的连发一定是对的。

找线段的过程用一个\(\text{set}\)来维护，对于每条线段记一个最后出现的在这条线段。

对于一些不好分类讨论的特殊情况，我们可以在所有线段的最下方加上一条直线，然后因为可能有斜率不存在的线段，我们可以把所有点随机旋转一下。就做完了

#include <cstdio>
#include <iostream>
#include <cmath>
#include <set>
#include <algorithm>
#define LL long long
#define db long double
using namespace std;
inline int read() {
	int res = 0, flag = 0; char ch = getchar();
	for(; !isdigit(ch); ch = getchar()) if(ch == '-') flag = 1;
	for(; isdigit(ch); ch = getchar()) res = (res << 1) + (res << 3) + (ch ^ 48);
	if(flag) res = ~res + 1; return res;
}
const db t = 1007.0/1009.0, cost = cos(t), sint = sin(t);
const db inf = 1e10;
const int N = 1e5 + 10;
int n;
struct Point {
	int id;
	db x, y;int rx, ry;
	inline void get() {x = rx * cost - ry * sint, y = ry * cost + rx * sint;}
	inline void print() {printf("%d %d ",rx,ry);}
	inline void print1() {printf("%Lf %Lf ",x,y);}
	inline bool operator < (const Point y) const{return (x < y.x);}
}p[N << 1], q[N << 1];
db nx;
struct Segment {
	int l, r; mutable int lst; db k, b;
	inline void get() {k = (p[l].y - p[r].y) / (p[l].x - p[r].x), b = p[l].y - p[l].x * k;}
	inline void print() {printf("seg %Lf %Lf\n",k,b);}
	inline bool operator < (const Segment y) const{return nx * k + b < nx * y.k + y.b;}
};
set <Segment> S; set <Segment>::iterator it;
int main() {
	n = read();
	for(register int i = 1; i <= n + n; i += 2) {
		p[i].rx = read(), p[i].ry = read(), p[i].get();
		p[i + 1].rx = read(), p[i + 1].ry = read(), p[i + 1].get();
		if(p[i + 1] < p[i]) swap(p[i], p[i + 1]); p[i].id = i, p[i + 1].id = i + 1;
		q[i] = p[i], q[i + 1] = p[i + 1];
	}
	sort(q + 1, q + n + n + 1);
	Segment tmp;
	tmp.k = 0, tmp.b = -inf, tmp.lst = 0, S.insert(tmp);
	for(register int i = 1, id; i <= n + n; ++i) {
		nx = q[i].x, id = q[i].id;
		if(id & 1) { // left
			tmp.l = id, tmp.r = id + 1, tmp.lst = id, tmp.get();
			it = S.insert(tmp).first; --it;
			if(it->lst) p[it->lst].print(), p[id].print(), printf("\n");
			it->lst = id;
		}
		else { // right
			tmp.l = id - 1, tmp.r = id, tmp.lst = id, tmp.get();
			it = S.find(tmp); --it;
			it->lst = id, S.erase(tmp);
		}
	}
}