如题，代码：

#include <bits/stdc++.h>
using namespace std;

#define int long long
#define maxn 100005
struct SegmentTree {
    int l, r;
    double add;
    double sum, squ;
} t[maxn * 4];
int n, m;
int a[maxn];

void push_up(int p) {
    t[p].sum = t[p << 1].sum + t[(p << 1) | 1].sum;
    t[p].squ = t[p << 1].squ + t[(p << 1) | 1].squ;
}

void build(int p, int l, int r) {
    t[p].l = l, t[p].r = r;
    if (l == r) { t[p].sum = a[l], t[p].squ = a[l] * a[l]; return; }
    int mid = (t[p].l + t[p].r) >> 1;
    build(p << 1, l, mid), build((p << 1) | 1, mid + 1, r);
    push_up(p);
}

void spread(int p) {
    if (t[p].add) {
        t[p << 1].squ += 2 * t[p].add * t[p << 1].sum + t[p].add * t[p].add * (t[p << 1].r - t[p << 1].l + 1);
        t[(p << 1) | 1].squ += 2 * t[p].add * t[(p << 1) | 1].sum + t[p].add * t[p].add * (t[(p << 1) | 1].r - t[(p << 1) | 1].l + 1);
        t[p << 1].sum += (t[p << 1].r - t[p << 1].l + 1) * t[p].add;
        t[(p << 1) | 1].sum += (t[(p << 1) | 1].r - t[(p << 1) | 1].l + 1) * t[p].add;
        t[p << 1].add += t[p].add, t[(p << 1) | 1].add += t[p].add;
        t[p].add = 0;
    }
}

void update(int p, int l, int r, int v) {
    if (l <= t[p].l && r >= t[p].r) {
        t[p].squ += 2 * v * t[p].sum + v * v * (t[p].r - t[p].l + 1);
        t[p].sum += (t[p].r - t[p].l + 1) * v;
        t[p].add += v;
    }
    spread(p);
    int mid = (t[p].l + t[p].r) >> 1;
    if (l <= mid) update(p << 1, l, r, v);
    if (r > mid) update((p << 1) | 1, l, r, v);
    push_up(p);
}

double ask(int p, int l, int r) {
    if (l <= t[p].r && r >= t[p].r) return t[p].sum;
    spread(p);
    int mid = (t[p].l + t[p].r) >> 1, val = 0;
    if (l <= mid) val += ask(p << 1, l, r);
    if (r > mid) val += ask((p << 1) | 1, l, r);
    return val;
}

double query(int p, int l, int r) {
    if (l <= t[p].l && r >= t[p].r) return t[p].squ;
    spread(p);
    int mid = (t[p].l + t[p].r) >> 1, val = 0;
    if (l <= mid) val += ask(p << 1, l, r);
    if (r > mid) val += ask((p << 1) | 1, l, r);
    return val;
}

signed main() {
    cin >> n >> m;
    for (int i = 1; i <= n; i++) cin >> a[i];
    build(1, 1, n);
    while (m--) {
        int op, x, y, k;
        cin >> op >> x >> y;
        if (op == 1) {
            cin >> k;
            update(1, x, y, k);
        }
        else if (op == 2) printf("%4lf\n", (ask(1, x, y) / (y - x + 1)));
        else printf("%4lf\n", (query(1, x, y) / (y - x + 1) - ask(1, x, y) * ask(1, x, y)));
    }
    return 0;
}