对空防御的训练 改编自BZOJ3165 （线段树永久化标记 李超线段树）

对空防御的训练

256MB / 1s ; defense.cpp / c / pas / in / out

【题目描述】
秋月十分擅长对空作战。尽管如此，必要的训练也不能懈怠。
在一次训练中，会有m个事件发生。
1. 出现一架飞机在(x1,y1)到(x2,y2)的线段上出现并来回飞行。
2. 秋月接到指令，要在x坐标为a处击坠一架飞行路线经过该位置的飞机。
为了确保击坠的可能性最大，秋月将选择在x位置时高度最低的飞机。
由于只是训练任务，秋月并不会真正地击坠飞机。你只需要在每次接到指令时告诉她高度最低的飞机的y坐标就可以了。

【输入格式】
第一行一个整数m，表示事件数。
接下来m行，每行有若干正整数，第一个数D表示事件类型。
若D=1，表示出现了一架飞机，接下来四个正整数x1,y1,x2,y2。
若D=2，接下来一个正整数a，表示秋月询问在x=a处能击坠的飞机的最低高度是多少。

【输出格式】
对于每个询问操作，输出一行一个实数，表示该飞机的y坐标。保证存在这样的飞机。

【样例数据】
defense1.in
4
1 1 1 3 3
2 3
1 1 3 3 1
2 3
defense1.out
3
1
defense2.in
7
1 3 6 5 7
1 6 1 9 7
1 1 3 7 7
2 3
2 7
1 5 2 10 4
2 7
defense2.out
4.333333
3
2.8

【数据范围】
对于30%的数据，m≤5000。
对于另30%的数据，1≤x,y,a≤100000。
对于100%的数据，1≤m≤200000，1≤x,y,a≤10^9，x1＜x2。
【评分方式】
本题有自定义校验器。当你的答案与标准答案的绝对误差不超过1e-4时，便被认为正确。

思路：
线段树维护线段，对于可以判断整体最小线段的区间直接更新，否则进行递归处理。处理到长度为1 的节点即可停止。永久化标记，查询时需要考虑路径上所有节点的标记，取最小值。
其实和BZOJ3165差不多，只是求最小的y，不强制在线而已。

#include <cstdio>#include <cstring>#include <iostream>#include <algorithm>#include <cmath>#define LL long long#define inf 1000000000#define N 200010using namespace std;int m, root, tot=1;const double eps=1e-8;struct Edge{    int lson, rson;    double ly, ry;}tree[N*30];inline int read(){    int x=0,f=1;char ch=getchar();    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}    return x*f;}inline double get(int x, int X1, double Y1, int X2, double Y2){//线段上横坐标为x时的y     return (Y2-Y1) * (x-X1) / (X2-X1) + Y1;}int modify(int pos, int l, int r, int x1, double y1, int x2, double y2){    if ( !pos ){        pos = ++tot;        tree[pos].ly = tree[pos].ry = inf;    }    int mid = (l + r) >> 1;    if (x1 <= l && x2 >= r){//线段覆盖整个区间         double ly = get(l, x1, y1, x2, y2), ry = get(r, x1, y1, x2, y2);        double sl = ly - tree[pos].ly, sr = ry - tree[pos].ry;        if (sl<=eps && sr<=eps){//维护更低的y             tree[pos].ly = ly;            tree[pos].ry = ry;            return pos;        }        if (sl>=-eps && sr>=-eps)            return pos;        //-----------------------------------------------------交错线段        if (fabs(sl) < fabs(sr)){//由相似三角形可得交点在左区间             if (sr <= eps){                tree[pos].lson = modify(tree[pos].lson, l, mid, l, tree[pos].ly, r, tree[pos].ry);//递归处理左区间                 tree[pos].ly = ly;                tree[pos].ry = ry;//当前区间相当于右区间（更新） 那么就不需要新增右区间了             }             else                tree[pos].lson = modify(tree[pos].lson, l, mid, x1, y1, x2, y2);        }         else {//同上             if (sl <= eps){                tree[pos].rson = modify(tree[pos].rson, mid+1, r, l, tree[pos].ly, r, tree[pos].ry);                tree[pos].ly = ly;                tree[pos].ry = ry;            }             else                tree[pos].rson = modify(tree[pos].rson, mid+1, r, x1, y1, x2, y2);        }    }    else {        if (x1 <= mid)            tree[pos].lson = modify(tree[pos].lson, l, mid, x1, y1, x2, y2);        if (x2 > mid)            tree[pos].rson = modify(tree[pos].rson, mid+1, r, x1, y1, x2, y2);    }    return pos;}double query(int pos, int l, int r, int x){    if ( !pos ) return inf;    if ( l == r ) return tree[pos].ly;    int mid = (l + r) >> 1;    double ans = get(x, l, tree[pos].ly, r, tree[pos].ry);    if (x <= mid)//递归左右区间         return min(ans, query(tree[pos].lson, l, mid, x));    else        return min(ans, query(tree[pos].rson, mid+1, r, x));}int main(){    freopen ("defense.in", "r", stdin);    freopen ("defense.out", "w", stdout);    m = read();    int x1, x2;    double y1, y2;    root = 1;    tree[root].ly = tree[root].ry = inf;    while ( m-- ){        x1 = read();        if (x1 == 1){            x1 = read(); y1 = read(); x2 = read(); y2 = read();            modify(root, 1, inf, x1, y1, x2, y2);        }        else {            x1 = read();            printf("%.10f\n", query(root, 1, inf, x1));        }    }    return 0;}