#include<iostream>
#include<cstring>
#include<iomanip>
#include<string>
#include<cstdio>
#include<queue>
#include<cmath>
#include<algorithm>
using namespace std;
long double t,ll,rr,kk,anss=0;
struct node{
	long double s;
	long long l,r;
	long double lan=0;
}tree[1000000];
long double ans=0;
long double da;
void down(long long n){
	tree[n*2].lan+=tree[n].lan;
	tree[n*2+1].lan+=tree[n].lan;
	tree[n*2].s+=(tree[n*2].r-tree[n*2].l+1)*tree[n].lan;
	tree[n*2+1].s+=(tree[n*2+1].r-tree[n*2+1].l+1)*tree[n].lan;
	tree[n].lan=0;
}
void build(long long l,long long r,long long n){
	tree[n].l=l;
	tree[n].r=r;
	if(l==r){
		cin >>tree[n].s;
		return;
	}
	build(l,(l+r)/2,n*2);
	build((l+r)/2+1,r,n*2+1);
	tree[n].s=tree[n*2].s+tree[n*2+1].s;
}
void qjcx(long long l, long long r,long long n){
	if(ll<=tree[n].l&&tree[n].r<=rr){
		ans+=tree[n].s;
		return;
	}
	down(n);
	if(ll<=(l+r)/2)qjcx(l,(l+r)/2,n*2);
	if(rr>(l+r)/2)qjcx((l+r)/2+1,r,n*2+1);
	tree[n].s=tree[n*2].s+tree[n*2+1].s;
}
void qjxg(long long l, long long r,long long n){
	if(ll<=tree[n].l&&tree[n].r<=rr){
		tree[n].s+=((r-l)+1)*kk;
	//	cout<<"------"<<((r-l)+1)*kk<<endl;
		tree[n].lan+=kk;
		return;
	}
	if(tree[n].lan)down(n);
	if(ll<=(l+r)/2)qjxg(l,(l+r)/2,n*2);
	if(rr>(l+r)/2)qjxg((l+r)/2+1,r,n*2+1);
	tree[n].s=tree[n*2].s+tree[n*2+1].s;
}
void qjcxx(long long l, long long r,long long n){
	//if(l>rr||r<ll)return; 
	if(l==r){
		anss+=(tree[n].s-da)*(tree[n].s-da);
//		cout <<"["<<l<<"] ";
		return;
	}
	down(n);
	qjcxx(l,(l+r)/2,n*2);
	qjcxx((l+r)/2+1,r,n*2+1);
	tree[n].s=tree[n*2].s+tree[n*2+1].s;
}
long long n,m; 
int  main(){
	cin >>n>>m;
	build(1,n,1);
//	for(long long i=1;i<=n*4;i++){
//		cout <<i<<":["<<tree[i].l<<" "<<tree[i].r<<"]"<<tree[i].s<<endl;
//	}
	for(long long i=1;i<=m;i++){
		ans=0;
		anss=0;
		cin >>t;
		if(t==1){
			cin >>ll>>rr>>kk;
			qjxg(1,n,1);
		}
		if(t==2){
			cin >>ll>>rr;
			qjcx(1,n,1);
			cout <<fixed<<setprecision(4)<<ans/(rr-ll+1)<<endl; 
		}
		if(t==3){
			cin >>ll>>rr;
			qjcx(1,n,1);
			da=ans/(rr-ll+1);
//			cout <<"da"<<da<<endl;
			qjcxx(1,n,1);
			cout <<fixed<<setprecision(4)<<anss/(rr-ll+1)<<endl;
		}
//		long long t=1,t1=1;
//		for(long long j=1;j<=9;j++){
//			cout <<j<<":["<<tree[j].l<<" "<<tree[j].r<<"]"<<"("<<tree[j].s<<") ";
//			if(t==j){
//				cout <<endl;
//				t+=2*t1;
//				t1++;
//			}
//		}
//		cout <<endl; 
	}
	return 0;
}
// 4 4 1 1 0