#include <cmath>
#include <cstdio>
using namespace std;
const double eps=1e-12;
const int N=100003,P=15;
int n; double a[N];
namespace SGT{
	struct Node{
		int l,r; double mul,sum[P+1];
	}node[N*3];
	Node operator+(Node l,Node r){
		Node ans={l.l,r.r,1};
		for(int i=1;i<=P;i++)
			ans.sum[i]=l.sum[i]+r.sum[i];
		return ans;
	}
	void update(int u,double x){
		int i; double j; node[u].mul*=x;
		for(i=1,j=x;i<=P;i++,j*=x)
			node[u].sum[i]*=j;
	}
	void push_up(int u){
		node[u]=node[u*2]+node[u*2+1];
	}
	void push_down(int u){
		if(fabs(node[u].mul-1)<eps) return ;
		update(u*2,node[u].mul),update(u*2+1,node[u].mul),node[u].mul=1;
	}
	void build(int u,int l,int r){
		if(l==r){
			int i; double j; node[u]={l,r,1};
			for(i=1,j=a[l];i<=P;i++,j*=a[l])
				node[u].sum[i]=j;
		}else{
			int mid=(l+r)/2;
			build(u*2  ,l,mid  );
			build(u*2+1,mid+1,r);
			push_up(u);
		}
	}
	void modify(int u,int l,int r,double x){
		if(node[u].l>=l&&node[u].r<=r) update(u,x);
		else{
			push_down(u);
			int mid=(node[u].l+node[u].r)/2;
			if(l<=mid) modify(u*2  ,l,r,x);
			if(r> mid) modify(u*2+1,l,r,x);
			push_up(u);
		}
	}
	Node query(int u,int l,int r){
		if(node[u].l>=l&&node[u].r<=r) return node[u];
		int mid=(node[u].l+node[u].r)/2; push_down(u);
		if(r<=mid) return query(u*2  ,l,r);
		if(l> mid) return query(u*2+1,l,r);
		return query(u*2,l,r)+query(u*2+1,l,r);
	}
}
int main(){
//	freopen("event.in","r",stdin);
//	freopen("event.out","w",stdout);
	double x,ans; SGT::Node s1;
	int i,j,q,opt,l,r;
	while(scanf("%d%d",&n,&q)!=EOF){
		for(i=1;i<=n;i++)
			scanf("%lf",&a[i]);
		SGT::build(1,1,n);
		for(i=1;i<=q;i++){
			scanf("%d%d%d",&opt,&l,&r);
			switch(opt){
				case 0:
					s1=SGT::query(1,l,r),ans=0;
					for(j=1;j<=P;j++) ans-=s1.sum[j]/j;
					printf("%.12f\n",ans); break;
				case 1:
					scanf("%lf",&x);
					SGT::modify(1,l,r,x);
			}
		}
	}
//	fclose(stdin);
//	fclose(stdout);
	return 0;
}