Codeforces 138C Mushroom Gnomes - 2 线段树

Mushroom Gnomes - 2

感觉没啥东西， 用线段树算算每个被覆盖的概率， 坑点是有很多个在同一个点。

#include<bits/stdc++.h>
#define LL long long
#define fi first
#define se second
#define mk make_pair
#define PLL pair<LL, LL>
#define PLI pair<LL, int>
#define PII pair<int, int>
#define SZ(x) ((int)x.size())
#define ull unsigned long long

using namespace std;

const int N = 4e5 + ;
const int inf = 0x3f3f3f3f;
const LL INF = 0x3f3f3f3f3f3f3f3f;
const int mod = 1e9 + ;
const double eps = 1e-;
const double PI = acos(-);

int A[N], h[N], b[N], X[N], n, m, tot;
double pl[N], pr[N], z[N], val[N];

#define lson l, mid, rt << 1
#define rson mid + 1, r, rt << 1 | 1
double a[N << ], lazy[N << ];
void push(int rt) {
    if(fabs(lazy[rt] - ) > eps) {
        a[rt << ] *= lazy[rt];
        a[rt <<  | ] *= lazy[rt];
        lazy[rt << ] *= lazy[rt];
        lazy[rt <<  | ] *= lazy[rt];
        lazy[rt] = ;
    }
}
void build(int l, int r, int rt) {
    lazy[rt] = ;
    if(l == r) {
        a[rt] = ;
        return;
    }
    int mid = l + r >> ;
    build(lson); build(rson);
}
void update(int L, int R, double val, int l, int r, int rt) {
    if(L > R) return;
    if(l >= L && r <= R) {
        a[rt] *= val;
        lazy[rt] *= val;
        return;
    }
    int mid = l + r >> ;
    push(rt);
    if(L <= mid) update(L, R, val, lson);
    if(R > mid)  update(L, R, val, rson);
}
double query(int p, int l, int r, int rt) {
    if(l == r) return a[rt];
    int mid = l + r >> ;
    push(rt);
    if(p <= mid) return query(p, lson);
    else return query(p, rson);
}

int getId(int x) {
    return lower_bound(X, X + tot, x) - X + ;
}

int main() {
    scanf("%d%d", &n, &m);
    for(int i = ; i <= n; i++) {
        scanf("%d%d", &A[i], &h[i]);
        scanf("%lf%lf", &pl[i], &pr[i]);
        pl[i] /= ; pr[i] /= ;
        X[tot++] = A[i];
        X[tot++] = A[i] + h[i];
        X[tot++] = A[i] - h[i];
    }
    for(int i = ; i <= m; i++) {
        scanf("%d%lf", &b[i], &z[i]);
        X[tot++] = b[i];
    }
    sort(X, X + tot);
    tot = unique(X, X + tot) - X;
    build(, tot, );
    for(int i = ; i <= m; i++) {
        val[getId(b[i])] += z[i];
    }
    for(int i = ; i <= tot; i++) {
        update(i, i, val[i], , tot, );
    }
    for(int i = ; i <= n; i++) {
        update(getId(A[i]) + , getId(A[i] + h[i]),  - pr[i], , tot, );
        update(getId(A[i] - h[i]), getId(A[i]) - ,  - pl[i], , tot, );
    }
    double ans = ;
    for(int i = ; i <= tot; i++) {
        ans += query(i, , tot, );
    }
    printf("%.12f\n", ans);
    return ;
}

/*
*/