线段树——Luogu1471 方差

题面：Luogu1471
一眼数据结构题
然后我们要维护区间方差？带区间修改？
听起来是不是很麻烦？
首先题目给的方差公式：

这个是平均数
我们可以化开来看：

然后因为中间一项的a[1]+a[2]+…+a[n]=*n
是不是？
所以继续化：

然后我们只要用线段树维护一下区间和还有平方和就可以了
但是区间修改怎么下传平方和呢？
我们把要加上的数记为v，那么：

然后就好了，lazy标记下传的时候搞一下就好了
据说分块比这个好写。。。我不会
代码有点长。。。

#include <cstdio>
#include <algorithm>
#include <cmath>
#include <cstring>
#include <iostream>
#include <ctime>
#include <map>
#include <queue>
#include <cstdlib>
#include <string>
#include <climits>
#include <set>
#include <vector>
using namespace std;
int n,m,lt[400001],rt[400001];
double t[400001],pt[400001],a[100001],add[400001];
inline double sqr(double x){return x*x;}
inline void clean(int nod){
    if(!add[nod])return;
    if(lt[nod]!=rt[nod]){
        pt[nod*2]+=2*add[nod]*t[nod*2]+sqr(add[nod])*(rt[nod*2]-lt[nod*2]+1);
        pt[nod*2+1]+=2*add[nod]*t[nod*2+1]+sqr(add[nod])*(rt[nod*2+1]-lt[nod*2+1]+1);
        t[nod*2]+=(rt[nod*2]-lt[nod*2]+1)*add[nod];
        t[nod*2+1]+=(rt[nod*2+1]-lt[nod*2+1]+1)*add[nod];
        add[nod*2]+=add[nod];add[nod*2+1]+=add[nod];
    }add[nod]=0;
}
inline void build(int l,int r,int nod){
    lt[nod]=l;rt[nod]=r;add[nod]=0;
    if(l==r){t[nod]=a[l];pt[nod]=sqr(a[l]);return;}
    int mid=l+r>>1;
    build(l,mid,nod*2);build(mid+1,r,nod*2+1);
    t[nod]=t[nod*2]+t[nod*2+1];pt[nod]=pt[nod*2]+pt[nod*2+1];
}
inline void xg(int i,int j,double v,int nod){
    clean(nod);
    if(lt[nod]>=i&&rt[nod]<=j){
        add[nod]+=v;pt[nod]+=2*v*t[nod]+sqr(v)*(rt[nod]-lt[nod]+1);
        t[nod]+=(rt[nod]-lt[nod]+1)*v;
        return;
    }
    int mid=lt[nod]+rt[nod]>>1;
    if(i<=mid)xg(i,j,v,nod*2);
    if(j>mid)xg(i,j,v,nod*2+1);
    t[nod]=t[nod*2]+t[nod*2+1];pt[nod]=pt[nod*2]+pt[nod*2+1];   
}
inline double ssum(int i,int j,int nod){
    clean(nod);
    if(lt[nod]>=i&&rt[nod]<=j)return t[nod];
    int mid=lt[nod]+rt[nod]>>1;double ans=0;
    if(i<=mid)ans+=ssum(i,j,nod*2);
    if(j>mid)ans+=ssum(i,j,nod*2+1);
    return ans;
} 
inline double spsum(int i,int j,int nod){
    clean(nod);
    if(lt[nod]>=i&&rt[nod]<=j)return pt[nod];
    int mid=lt[nod]+rt[nod]>>1;double ans=0;
    if(i<=mid)ans+=spsum(i,j,nod*2);
    if(j>mid)ans+=spsum(i,j,nod*2+1);
    return ans;
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)scanf("%lf",&a[i]);
    build(1,n,1);
    for(int i=1;i<=m;i++){
        int p,x,y;scanf("%d%d%d",&p,&x,&y);
        if(p==1){
            double z;scanf("%lf",&z);
            xg(x,y,z,1);
        }else{
            double a1=ssum(x,y,1)/(double)(y-x+1);
            if(p==2)printf("%.4lf\n",a1);
            else{
                double a2=spsum(x,y,1)/(double)(y-x+1);
                printf("%.4lf\n",a2-sqr(a1));
            }
        }
    }
    return 0;
}