前言

略

題目相關(guān)

鏈接

題目大意

單點(diǎn)修改，區(qū)間方差

數(shù)據(jù)范圍

$n,q\le 100000$

題解

區(qū)間方差的定義
$$\frac{1}{r?l+1}\sum _{i=l}^r(a_i?a{)}^2$$
其中$$a=\frac{1}{r?l+1}\sum _{i=l}^r{a}_i$$
直接展開1式
$$\begin{array}{rl}Ans& =\frac{1}{r?l+1}\sum _{i=l}^r(a_i?a{)}^2\\ & =\frac{1}{r?l+1}\sum _{i=l}^r(a_i^2?2a_{i}a+a^2)\end{array}$$
第三項(xiàng)求$a$很方便，區(qū)間求和即可
我們發(fā)現(xiàn)，第二項(xiàng)可以直接套用區(qū)間求和的結(jié)果，乘起來即可
第一項(xiàng)是區(qū)間平方和，直接維護(hù)
兩個(gè)BIT維護(hù)區(qū)間和和平方和即可

代碼

#include<bits/stdc++.h> typedef long long ll; #define rg register template <typename T> inline void read(T&x){char cu=getchar();x=0;bool fla=0;while(!isdigit(cu)){if(cu=='-')fla=1;cu=getchar();}while(isdigit(cu))x=x*10+cu-'0',cu=getchar();if(fla)x=-x;} template <typename T> inline void printe(const T x){if(x>=10)printe(x/10);putchar(x%10+'0');} template <typename T> inline void print(const T x){if(x<0)putchar('-'),printe(-x);else printe(x);} inline int lowbit(const int x){return x&(-x);} const int mod=1000000007; void md(int &x){if(x>=mod)x-=mod;} int Md(const int x){return x>=mod?x-mod:x;} inline ll pow(ll x,ll y) {ll res=1;for(;y;y>>=1,x=x*x%mod)if(y&1)res=res*x%mod;return res; } int n,m,v[100001]; struct BIT {int a[100001];inline void add(int u,const int v){while(u<=n)md(a[u]+=v),u+=lowbit(u);}inline int query(int u){ll res=0;while(u)res+=a[u],u-=lowbit(u);return res%mod;} }A,B; int main() {read(n),read(m);for(rg int i=1;i<=n;i++)read(v[i]),A.add(i,v[i]),B.add(i,(ll)v[i]*v[i]%mod);for(rg int i=1;i<=m;i++){int opt,a,b;read(opt),read(a),read(b);if(opt==1)A.add(a,mod-v[a]),B.add(a,mod-(ll)v[a]*v[a]%mod),v[a]=b,A.add(a,v[a]),B.add(a,(ll)v[a]*v[a]%mod);else{a--;const int lenth=b-a,inv=pow((ll)lenth,mod-2);const int q=Md(mod+A.query(b)-A.query(a)),p=(ll)q*inv%mod;print((ll)(mod+mod+(ll)p*p%mod*lenth%mod-(ll)q*p%mod*2%mod+B.query(b)-B.query(a))%mod*inv%mod),putchar('\n');}}return 0; }

總結(jié)

模板題，增強(qiáng)信心（大霧）

總結(jié)

以上是生活随笔為你收集整理的[luogu5142]区间方差的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。