不过样例但40pts求调

#include<bits/stdc++.h>
using namespace std;
double n,m;
struct tree{
	long long l,r;
	double sum,add,psum;
}t[400000];
double a[100001];
void push_up(long long p){
	t[p].sum=t[p<<1].sum+t[p<<1|1].sum;
	t[p].psum=t[p<<1].psum+t[p<<1|1].psum;
}
void push_down(long long p){
	if(t[p].add){
		long long len=t[p].r-t[p].l+1;
		t[p<<1].psum+=2*t[p].add*t[p<<1].psum+(len-len/2)*t[p].add*t[p].add;
		t[p<<1|1].psum+=2*t[p].add*t[p<<1|1].psum+(len/2)*t[p].add*t[p].add;
		t[p<<1].sum+=t[p].add*(len-len/2);
		t[p<<1|1].sum+=t[p].add*(len/2);
		t[p<<1].add+=t[p].add;
		t[p<<1|1].add+=t[p].add;
		t[p].add=0;
	}
}
void build(long long p,long long l,long long r){
	t[p].l=l,t[p].r=r;
	if(l==r){
		t[p].sum=a[l],t[p].psum=a[l]*a[l];
		return;
	}
	long long mid=(l+r)>>1;
	build(p<<1,l,mid);
	build(p<<1|1,mid+1,r);
	push_up(p);
}
void updata(long long p,long long l,long long r,double x){
	if(l<=t[p].l&&t[p].r<=r){
		t[p].sum+=x*(t[p].r-t[p].l+1);
		t[p].add+=x;
		return;
	}
	push_down(p);
	long long mid=(t[p].l+t[p].r)>>1;
	if(l<=mid) updata(p<<1,l,r,x);
	if(r>mid) updata(p<<1|1,l,r,x);
	push_up(p);
}
double query(long long p,long long l,long long r){
	if(l<=t[p].l&&t[p].r<=r){
		return t[p].sum;
	}
	push_down(p);
	double sum=0,mid=(t[p].l+t[p].r)/2;
	if(l<=mid) sum+=query(p<<1,l,r);
	if(r>mid) sum+=query(p<<1|1,l,r);
	return sum;
}
double pquery(long long p,long long l,long long r){
	if(l<=t[p].l&&t[p].r<=r){
		return t[p].psum;
	}
	push_down(p);
	double sum=0,mid=(t[p].l+t[p].r)>>1;
	if(l<=mid) sum+=pquery(p<<1,l,r);
	if(r>mid) sum+=pquery(p<<1|1,l,r);
	return sum;
}
int main(){
	cin>>n>>m;
	for(int i=1;i<=n;i++) cin>>a[i];
	build(1,1,n);
	while(m--){
		int op;
		cin>>op;
		if(op==1){
			long long x,y;
			double k;
			cin>>x>>y>>k;
			updata(1,x,y,k);
		}
		else if(op==2){
			long long x,y;
			cin>>x>>y;
			printf("%.4lf\n",query(1,x,y)/(y-x+1));
		}
		else{
			long long x,y;
			cin>>x>>y;
			printf("%.4lf\n",pquery(1,x,y)/(y-x+1)-((query(1,x,y)/(y-x+1))*(query(1,x,y)/(y-x+1))));
		}
	}
	return 0;
}