CF895 E. Eyes Closed(线段树 期望)

题意

sol

今天考试的t3，，我本来留了一个小时去写。但是t2一刚就刚了两个小时

最后也没来的及写。。

然后考完 开始写，，25min就a了。。

感觉自己太高估自己的思维，太低估自己的码力了。。。

这题比较简单吧

期望的和等于和的期望

然后线段树维护每个节点的值就可以了

交换的时候分自己不变和变成哪个数讨论一下

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<map>
#include<vector>
#include<set>
#include<queue>
#include<cmath>
#include<iostream>
#define pair pair<int, int>
#define mp(x, y) make_pair(x, y)
#define fi first
#define se second
#define ls k << 1
#define rs k << 1 | 1
#define int long long 
#define ll long long 
//#define getchar() (p1 == p2 && (p2 = (p1 = buf) + fread(buf, 1, 1<<22, stdin), p1 == p2) ? eof : *p1++)
//char buf[(1 << 22)], *p1 = buf, *p2 = buf;
using namespace std;
const int maxn = 1e6 + 10, inf = 1e9 + 10, mod = 1e9 + 7;
const double eps = 1e-9;
inline int read() {
    char c = getchar(); int x = 0, f = 1;
    while(c < '0' || c > '9') {if(c == '-') f = -1; c = getchar();}
    while(c >= '0' && c <= '9') x = x * 10 + c - '0', c = getchar();
    return x * f;
}
int n, q;
struct node {
    int l, r, siz;
    double sum, mul, add;
}t[maxn];
int a[maxn];
void pr(int k, int son) {
    t[son].sum *= t[k].mul;
    t[son].sum += t[son].siz * (t[k].add);
    t[son].mul *= t[k].mul;
    t[son].add *= t[k].mul;
    t[son].add += t[k].add;
}
void pushdown(int k) {
    if(t[k].mul == 1 && t[k].add == 0) return ;
    pr(k, ls); pr(k, rs); 
    t[k].mul = 1; t[k].add = 0; 
}
void update(int k) {
    t[k].sum = t[ls].sum + t[rs].sum;
}
void build(int k, int ll, int rr) {
    t[k].siz = rr - ll + 1; t[k].l = ll; t[k].r = rr; t[k].mul = 1;
    if(ll == rr) {
        t[k].sum = a[ll];
        return ;
    }
    int mid = (ll + rr) >> 1;
    build(ls, ll, mid); build(rs, mid + 1, rr);
    update(k);
}
void intervalmul(int k, int ll, int rr, double val) {
    if(ll <= t[k].l && t[k].r <= rr) {
        t[k].sum *= val; t[k].mul *= val; t[k].add *= val;
        return ;
    }
    pushdown(k);
    int mid = t[k].l + t[k].r >> 1;
    if(ll <= mid) intervalmul(ls, ll, rr, val); 
    if(rr >  mid) intervalmul(rs, ll, rr, val);
    update(k);
}
void intervaladd(int k, int ll, int rr, double val) {
    if(ll <= t[k].l && t[k].r <= rr) {
        t[k].sum += t[k].siz * val; t[k].add += val;
        return ;
    }
    pushdown(k);
    int mid = t[k].l + t[k].r >> 1;
    if(ll <= mid) intervaladd(ls, ll, rr, val); 
    if(rr > mid)  intervaladd(rs, ll, rr, val);
    update(k);
}
double query(int k, int ll, int rr) {
    double ans = 0;
    if(ll <= t[k].l && t[k].r <= rr) {
        ans += t[k].sum; return ans;
    }
    pushdown(k);
    int mid = t[k].l + t[k].r >> 1;
    if(ll <= mid) ans += query(ls, ll, rr); 
    if(rr > mid)  ans += query(rs, ll, rr);
    return ans;
}
main() {
    n = read(); q = read();
    for(int i = 1; i <= n; i++) a[i] = read();
    build(1, 1, n);
    while(q--) {
        int opt = read();
        if(opt == 1) {
            int l1 = read(), r1 = read(), l2 = read(), r2 = read();
            double d1 = query(1, l1, r1);
            double d2 = query(1, l2, r2), len1 = r1 - l1 + 1, len2 = r2 - l2 + 1;
            intervalmul(1, l1, r1, (len1 - 1) / len1 ); intervalmul(1, l2, r2, (len2 - 1) / len2);
            intervaladd(1, l1, r1, d2 / len1 / len2); intervaladd(1, l2, r2, d1 / len1 / len2);
        } else {
            int l = read(), r = read();
            printf("%.10lf\n", query(1, l, r));
        }
    }
    return 0;
}