这题觉得似乎不难的样子……
但硬是没有做出来，真是不知说什么好喵~
注意到没有两条共线的线段具有公共点，没有重合的线段
说明每个十字形最多涉及一个水平线段和一个竖直线段
这说明我们可以二分了:每条线段两端砍掉delta长度后，有没有线段有公共点？
很水的扫描线吧~ 按 X 轴扫描，先把横的线段扫进去，再用竖的线段去查询
但是写法是很重要的，我说我不用线段树你信喵？我说我连离散化都不用你信喵？
约大爷教你做人！@Ruchiose
1 #include <cstdio> 2 #include <algorithm> 3 #include <set> 4 const int size=100000; 5 6 namespace IOspace 7 { 8 inline int getint() 9 { 10 register int num=0; 11 register char ch=0, last; 12 do last=ch, ch=getchar(); while (ch<'0' || ch>'9'); 13 do num=num*10+ch-'0', ch=getchar(); while (ch>='0' && ch<='9'); 14 if (last=='-') num=-num; 15 return num; 16 } 17 inline void putint(int num, char ch=' ') 18 { 19 char stack[10]; 20 register int top=0; 21 if (num==0) stack[top=1]='0'; 22 for ( ;num;num/=10) stack[++top]=num%10+'0'; 23 for ( ;top;top--) putchar(stack[top]); 24 if (ch) putchar(ch); 25 } 26 } 27 28 struct line 29 { 30 int x, y, l; 31 inline line() {} 32 inline line(int _x, int _y, int _l):x(_x), y(_y), l(_l) {} 33 }; 34 struct node 35 { 36 int type; 37 int x, y1, y2; 38 inline node() {} 39 inline node(int _type, int _x, int _y1, int _y2):type(_type), x(_x), y1(_y1), y2(_y2) {} 40 inline bool operator < (node a) const {return x!=a.x?x<a.x:type*y2>a.type*a.y2;} 41 }; 42 43 int n, m, k; 44 line h[size], s[size]; 45 node q[size]; 46 std::multiset<int> t; 47 inline int max(int x, int y) {return x>y?x:y;} 48 inline void swap(int & a, int & b) {int t=a; a=b; b=t;} 49 inline bool query(int, int); 50 inline bool check(int); 51 52 int main() 53 { 54 n=IOspace::getint(); m=IOspace::getint(); 55 for (int i=0;i<n;i++) s[i].x=IOspace::getint(), s[i].y=IOspace::getint(), s[i].l=IOspace::getint(); 56 for (int i=0;i<m;i++) h[i].x=IOspace::getint(), h[i].y=IOspace::getint(), h[i].l=IOspace::getint(); 57 58 int left=0, right=0, mid; 59 for (int i=0;i<n;i++) right=max(right, s[i].l); 60 while (left+1<right) 61 { 62 mid=(left+right)>>1; 63 if (check(mid)) left=mid; 64 else right=mid; 65 } 66 67 IOspace::putint(left); 68 69 return 0; 70 } 71 72 inline bool query(int l, int r) 73 { 74 return t.lower_bound(l)!=t.upper_bound(r); 75 } 76 inline bool check(int d) 77 { 78 k=0; 79 for (int i=0;i<m;i++) if (h[i].l>=(d<<1)) 80 { 81 q[k++]=node(1, h[i].x+d, h[i].y, 1); 82 q[k++]=node(1, h[i].x+h[i].l-d, h[i].y, -1); 83 } 84 for (int i=0;i<n;i++) if (s[i].l>=(d<<1)) 85 q[k++]=node(0, s[i].x, s[i].y+d, s[i].y+s[i].l-d); 86 std::sort(q, q+k); 87 t.clear(); 88 for (int i=0;i<k;i++) 89 { 90 if (q[i].type) 91 { 92 if (q[i].y2==1) t.insert(q[i].y1); 93 else t.erase(t.find(q[i].y1)); 94 } 95 else 96 if (query(q[i].y1, q[i].y2)) 97 return true; 98 } 99 return false; 100 }
只存操作(包括插入、删除、询问)真是无比的高大上，跪跪跪
还有那鬼畜的查询写法
(set).lower_bound(l)!=(set).upper_bound(r)
可以查询 [l..r] 中是否有数存在，至于为什么这样是对的而不是别的样子了，大家在纸上画画就知道了吧
其实更好理解的是
当 (set).lower_bound(l)==(set).upper_bound(r) 时， [l..r] 中一定没有数
lower_bound是为了将 l 算入 [l..r] 中，upper_bound(r) 也是为了将 r 算入 [l..r] 中
真是跪跪跪跪跪跪跪跪