#include <bits/stdc++.h>
using namespace ::std;
const long long maxn = 1e5 + 500;
const long long mod = 1e9 + 7;
const long long inf = 1e9 + 500;
long double lazyz[maxn * 8];
long double lazys[maxn * 8];
long double sum[maxn * 8];
void shift(long long node, long long l, long long r) {
  long long mid = (l + r) / 2;
  lazyz[2 * node] *= lazyz[node];
  lazys[2 * node] *= lazyz[node];
  sum[2 * node] *= lazyz[node];
  lazyz[2 * node + 1] *= lazyz[node];
  lazys[2 * node + 1] *= lazyz[node];
  sum[2 * node + 1] *= lazyz[node];
  lazyz[node] = 1;
  lazys[2 * node] += lazys[node];
  sum[2 * node] += (mid - l) * lazys[node];
  lazys[2 * node + 1] += lazys[node];
  sum[2 * node + 1] += (r - mid) * lazys[node];
  lazys[node] = 0;
}
void updatez(long long L, long long R, long long node, long long l, long long r,
             long double z) {
  if (l <= L && R <= r) {
    shift(node, L, R);
    lazyz[node] *= z;
    sum[node] *= z;
    return;
  }
  if (R <= l || r <= L) {
    return;
  }
  long long mid = (L + R) / 2;
  shift(node, L, R);
  updatez(L, mid, 2 * node, l, r, z);
  updatez(mid, R, 2 * node + 1, l, r, z);
  sum[node] = sum[2 * node] + sum[2 * node + 1];
}
void updates(long long L, long long R, long long node, long long l, long long r,
             long double z) {
  if (l <= L && R <= r) {
    shift(node, L, R);
    lazys[node] += z;
    sum[node] += z * (R - L);
    return;
  }
  if (R <= l || r <= L) {
    return;
  }
  long long mid = (L + R) / 2;
  shift(node, L, R);
  updates(L, mid, 2 * node, l, r, z);
  updates(mid, R, 2 * node + 1, l, r, z);
  sum[node] = sum[2 * node] + sum[2 * node + 1];
}
long double finds(long long L, long long R, long long node, long long l,
                  long long r) {
  if (l <= L && R <= r) {
    return sum[node];
  }
  if (R <= l || r <= L) {
    return 0;
  }
  long long mid = (L + R) / 2;
  shift(node, L, R);
  return finds(L, mid, 2 * node, l, r) + finds(mid, R, 2 * node + 1, l, r);
}
int main() {
  ios_base::sync_with_stdio(0);
  cin.tie(0);
  cout.tie(0);
  fill(lazyz, lazyz + maxn * 4, 1);
  long long n, m;
  cin >> n >> m;
  for (long long i = 0; i < n; i++) {
    long long v;
    cin >> v;
    updates(0, n, 1, i, i + 1, v);
  }
  for (long long i = 0; i < m; i++) {
    long long typ;
    cin >> typ;
    if (typ == 1) {
      long long l1, l2, r1, r2;
      cin >> l1 >> r1 >> l2 >> r2;
      l1--;
      l2--;
      r1--;
      r2--;
      long long f = (r1 - l1 + 1) * (r2 - l2 + 1);
      long double s2 = finds(0, n, 1, l2, r2 + 1) / f;
      long double s1 = finds(0, n, 1, l1, r1 + 1) / f;
      long double p = (r1 - l1) / ((long double)(r1 - l1 + 1));
      updatez(0, n, 1, l1, r1 + 1, p);
      updates(0, n, 1, l1, r1 + 1, s2);
      p = (r2 - l2) / ((long double)(r2 - l2 + 1));
      updatez(0, n, 1, l2, r2 + 1, p);
      updates(0, n, 1, l2, r2 + 1, s1);
    } else {
      long long l, r;
      cin >> l >> r;
      l--;
      r--;
      cout << fixed << setprecision(13) << finds(0, n, 1, l, r + 1) << endl;
    }
  }
}
