luogu P5142 区间方差 十分优美的线段树

又来了个维护方差的线段树。。。。

大致推导过程（字丑多包涵QAQ）

注意取模时要加一些100000007防止出现负数。。

#include<cstdio>
#include<iostream>
#define ll long long
#define R register long long
#define ls (tr<<1)
#define rs (tr<<1|1)
using namespace std;
const long long M=;
int n,m;
inline ll g() {
    R ret=,fix=; register char ch; while(!isdigit(ch=getchar())) fix=ch=='-'?-:fix;
    do ret=ret*+(ch^); while(isdigit(ch=getchar())); return ret*fix;
}
ll sum[],po[];
void build(int tr,int l,int r) {
    if(l==r) {sum[tr]=g(),po[tr]=sum[tr]*sum[tr]; return;}
    R md=(l+r)>>;
    build(ls,l,md),build(rs,md+,r);
    sum[tr]=sum[ls]+sum[rs],po[tr]=po[ls]+po[rs];
}
void add(int tr,int l,int r,int pos,ll d) {
    if(l==r) {sum[tr]=d%M; po[tr]=d*d%M; return ;}
    R md=(l+r)>>;
    if(pos>md) add(rs,md+,r,pos,d);
    else if(pos<md+) add(ls,l,md,pos,d);
    sum[tr]=(sum[ls]+sum[rs])%M,po[tr]=(po[ls]+po[rs])%M;
}
ll query(int tr,int l,int r,int LL,int RR) {
    if(l==LL&&r==RR) return sum[tr];
    R md=(l+r)>>;
    if(LL>md) return query(rs,md+,r,LL,RR);
    else if(RR<md+) return query(ls,l,md,LL,RR);
    else return (query(ls,l,md,LL,md)+query(rs,md+,r,md+,RR))%M;
}
ll _query(int tr,int l,int r,int LL,int RR) {
    if(l==LL&&r==RR) return po[tr];
    R md=(l+r)>>;
    if(LL>md) return _query(rs,md+,r,LL,RR);
    else if(RR<md+) return _query(ls,l,md,LL,RR);
    else return (_query(ls,l,md,LL,md)+_query(rs,md+,r,md+,RR))%M;
}
inline ll qpow(int a) {
    R ret=,p=M-; a%=M;
    for(;p;p>>=,(a*=a)%=M) if(p&) (ret*=a)%=M;
    return ret;
}
inline ll Inv(int x) {
    if(x==) return ; if(x<) return ;
    return (M-M/x)*Inv(M%x)%M;
}
signed main() {
    n=g(),m=g();
    build(,,n);
    for(R i=;i<=m;++i) {
        R k=g(),l=g(),r=g();
        if(k&) add(,,n,l,r);
        else {
            R inv=Inv(r-l+),ave=query(,,n,l,r)*inv%M;
            printf("%lld\n",((_query(,,n,l,r)%M-(ave<<)*query(,,n,l,r)%M+M)*inv%M+ave*ave%M+M)%M);
        }
    }
}

也可树状数组，就不多说了（自己没写）（怕是树状数组更易写？）

---

2019.04.12