#include <bits/stdc++.h>
using namespace std;
constexpr int N=1e5+7;
using lf=double;
int n,m,x[N],y[N];
inline lf ss(lf l,lf r){
	return (l+r)*(r-l+1)/2;
}
inline lf sqs(lf l,lf r){
	--l;
	return r*(r+1)*(2*r+1)/6.0-l*(l+1)*(2*l+1)/6.0;
}
struct Res{
	lf x,y,x2,xy;
};
Res operator+(const Res &a,const Res &b){
	return Res{a.x+b.x,a.y+b.y,a.x2+b.x2,a.xy+b.xy};
}
struct Seg{
	struct Tag{
		bool e;
		lf x,y;
	};
	struct{
		int l,r;
		lf xs,ys,x2s,xys;
		Tag lz3,lz2;
		inline int len()const{return r-l+1;}
		void init(int x,int y){
			xs=x,ys=y,
			x2s=1ll*x*x,
			xys=1ll*x*y;
			lz2=lz3={false,0.0,0.0};
		}
	}tr[N<<2];
	void doasn(int k,Tag tg){
		if(!tg.e) return ;
		tr[k].lz2.e=false;
		tr[k].lz2.x=tr[k].lz2.y=0;
		tr[k].lz3=tg;
		tr[k].xs=ss(tr[k].l+tg.x,tr[k].r+tg.x);
		tr[k].ys=ss(tr[k].l+tg.y,tr[k].r+tg.y);
		tr[k].x2s=sqs(tr[k].l+tg.x,tr[k].r+tg.x);
		tr[k].xys=tr[k].len()*tg.x*tg.y+ss(tr[k].l,tr[k].r)*(tg.x+tg.y)+sqs(tr[k].l,tr[k].r);
	}
	void doadd(int k,Tag d){
		if(!d.e) return ;
		tr[k].lz2.e = true;
		tr[k].lz2.x += d.x,
		tr[k].lz2.y += d.y;
		tr[k].x2s += 2*d.x*tr[k].xs+tr[k].len()*d.x*d.x;
		tr[k].xys += d.x*tr[k].ys+d.y*tr[k].xs+d.x*d.y*tr[k].len();
		tr[k].xs += tr[k].len()*d.x;
		tr[k].ys += tr[k].len()*d.y;
	}
	void pup(int k){
		tr[k].xs = tr[k<<1].xs+tr[k<<1|1].xs;
		tr[k].ys = tr[k<<1].ys+tr[k<<1|1].ys;
		tr[k].x2s = tr[k<<1].x2s+tr[k<<1|1].x2s;
		tr[k].xys = tr[k<<1].xys+tr[k<<1|1].xys;
	}
	void pdn(int k){
		if(tr[k].lz3.e){
			doasn(k<<1,tr[k].lz3);
			doasn(k<<1|1,tr[k].lz3);
			tr[k].lz3={false,0.0,0.0};
		}
		if(tr[k].lz2.e){
			doadd(k<<1,tr[k].lz2);
			doadd(k<<1|1,tr[k].lz2);
			tr[k].lz2={false,0.0,0.0};
		}
	}
	void build(int k,int l,int r){
		tr[k].l=l,tr[k].r=r;
		tr[k].lz2=tr[k].lz3={false,0.0,0.0};
		if(l==r){
			tr[k].init(x[l],y[l]);
			return ;
		}
		int mi((l+r)>>1);
		build(k<<1,l,mi),build(k<<1|1,mi+1,r),
		pup(k);
	}
	void modify(int k,int l,int r,int ty,int s,int t){
		if(tr[k].l>=l&&tr[k].r<=r){
			if(ty==2) doadd(k,Tag{true,1.0*s,1.0*t});
			else doasn(k,Tag{true,1.0*s,1.0*t});
			return ;
		}
		pdn(k); int mi((tr[k].l+tr[k].r)>>1);
		if(l<=mi) modify(k<<1,l,r,ty,s,t);
		if(r>mi) modify(k<<1|1,l,r,ty,s,t);
		pup(k);
	}
	Res query(int k,int l,int r){
		if(tr[k].l>=l&&tr[k].r<=r){
			return Res{tr[k].xs,tr[k].ys,tr[k].x2s,tr[k].xys};
		}
		pdn(k); int mi((tr[k].l+tr[k].r)>>1);
		if(r<=mi) return query(k<<1,l,r);
		else if(l>mi) return query(k<<1|1,l,r);
		else return query(k<<1,l,r)+query(k<<1|1,l,r);
	}
}seg;
int main(){
	cin.tie(0)->sync_with_stdio(0);
	cout<<fixed<<setprecision(10);
	cin>>n>>m;
	for(int i{1};i<=n;++i){
		cin>>x[i];
	}
	for(int i{1};i<=n;++i){
		cin>>y[i];
	}
	seg.build(1,1,n);
	while(m--){
		int t,l,r,x{0},y{0};
		cin>>t>>l>>r;
		if(t==1){
			auto res=seg.query(1,l,r);
			int len{r-l+1};
			lf x0=res.x/len,y0=res.y/len;
			cout<< (res.xy-res.y*x0-res.x*y0+x0*y0*len) / (res.x2-2*res.x*x0+x0*x0*len) <<'\n';
		}
		else{
			cin>>x>>y;
			seg.modify(1,l,r,t,x,y);
		}
	}
	return 0;
}

/*
 (x-x0)(y-y0)
= xy-x0y-y0x+x0y0
 (x-x0)^2
= x^2 -2x*x0+x0^2

sum xy, sum x, sum y, sum x^2
*/