#include <bits/stdc++.h>
using namespace std;
template <typename T1, typename T2>
bool Max(T1& a, T2 b) {
  return a < b ? a = b, 1 : 0;
}
template <typename T1, typename T2>
bool Min(T1& a, T2 b) {
  return a > b ? a = b, 1 : 0;
}
template <typename T, typename AD>
struct ST {
  vector<T> st;
  vector<AD> lazy;
  int n, h;
  static T combine(T a, T b) { return a + b; }
  static T defaultT() { return 0; }
  void add(int i, int sz, AD change) {
    st[i] = change.apply(st[i], sz);
    if (sz > 1) {
      lazy[i].combine(change);
    }
  }
  int hsb_loc(int n) { return 31 - __builtin_clz(n); }
  void init(int _n) {
    n = _n;
    h = hsb_loc(n) + 1;
    st.assign(n << 1, defaultT());
    lazy.assign(n, AD());
  }
  void init(T* a, int _n) {
    init(_n);
    for (int i = 0; i < n; i++) {
      st[n + i] = a[i];
    }
    for (int i = n - 1; i > 0; i--) {
      st[i] = combine(st[i << 1], st[i << 1 | 1]);
    }
  }
  void push(int p) {
    for (int s = h; s > 0; --s) {
      int i = p >> s;
      if (lazy[i]) {
        add(i << 1, 1 << (s - 1), lazy[i]);
        add(i << 1 | 1, 1 << (s - 1), lazy[i]);
        lazy[i] = AD();
      }
    }
  }
  T query(int l, int r) {
    l += n;
    r += n;
    push(l);
    push(r);
    T ret = defaultT();
    for (; l <= r; l >>= 1, r >>= 1) {
      if (l & 1) ret = combine(ret, st[l++]);
      if (!(r & 1)) ret = combine(ret, st[r--]);
    }
    return ret;
  }
  void update(int l, int r, AD update) {
    l += n;
    r += n;
    push(l);
    push(r);
    int l0 = l, r0 = r;
    for (int k = 1; l <= r; l >>= 1, r >>= 1, k <<= 1) {
      if (l & 1) add(l++, k, update);
      if (!(r & 1)) add(r--, k, update);
    }
    while (!(l0 & 1)) l0 >>= 1;
    for (int i = l0; i > 1; i >>= 1) {
      st[i >> 1] = combine(st[i], st[i ^ 1]);
    }
    while (r0 & 1) r0 >>= 1;
    for (int i = r0; i > 1; i >>= 1) {
      st[i >> 1] = combine(st[i], st[i ^ 1]);
    }
  }
};
struct upd {
  double m, a;
  upd() : m(1), a(0) {}
  upd(double _m, double _a) : m(_m), a(_a) {}
  double apply(double num, int sz) { return num * m + a * sz; }
  double combine(upd rhs) {
    m *= rhs.m;
    a = a * rhs.m + rhs.a;
  }
  operator bool() { return m != 1 || a != 0; }
};
const int N = 2000001, inf = 1e9;
double a[N];
void solve() {
  int n, q;
  cin >> n >> q;
  for (int i = 0; i < (n); ++i) cin >> a[i];
  ST<double, upd> st;
  st.init(a, n);
  while (q--) {
    int t, l1, r1, l2, r2;
    cin >> t >> l1 >> r1;
    l1--, r1--;
    if (t == 1) {
      cin >> l2 >> r2;
      l2--, r2--;
      double n1 = r1 - l1 + 1, n2 = r2 - l2 + 1;
      double exp1 = st.query(l1, r1) / n1;
      double exp2 = st.query(l2, r2) / n2;
      st.update(l1, r1, upd((n1 - 1) / n1, exp2 / n1));
      st.update(l2, r2, upd((n2 - 1) / n2, exp1 / n2));
    } else {
      double ans = st.query(l1, r1);
      cout << ans << '\n';
    }
  }
}
void init() {}
void gen() {}
int main(int argc, char** argv) {
  ios::sync_with_stdio(false);
  cin.tie(0);
  cout.setf(ios::fixed);
  cout.precision(20);
  if (argc > 1 && strcmp(argv[1], "gen") == 0) {
    gen();
    return 0;
  }
  init();
  int TC = 1;
  for (int TI = 1; TI <= (TC); ++TI) {
    solve();
  }
  return 0;
}
