luogu3707 相关分析 (线段树)

把式子展开以后会发现，可以用线段树维护$x,y,x*y,x^2$分别的区间和

然后操作有区间加和区间修改

这个pushdown的时候，如果改和加的标记同时存在，那一定是先改再加，要不然加的标记已经被清掉了

所以在pushdown的时候，如果有改的标记，要把孩子的加的标记清掉

然后注意细节就行了（用*传数组 然后在函数里改了的话 它真的会改的 怎么就意识不到呢...）

 #include<bits/stdc++.h>
 #define pa pair<int,int>
 #define CLR(a,x) memset(a,x,sizeof(a))
 using namespace std;
 typedef long long ll;
 const int maxn=1e5+,inf=1e9;

 inline ll rd(){
     ll x=;char c=getchar();int neg=;
     while(c<''||c>''){if(c=='-') neg=-;c=getchar();}
     while(c>=''&&c<='') x=x*+c-'',c=getchar();
     return x*neg;
 }

 struct Node{
     double x,y,xy,x2;
     int l,r;
     Node(double a=,double b=,double c=,double d=,int e=,int f=){
         x=a,y=b,xy=c,x2=d,l=e,r=f;
     }
 }tr[maxn*];
 double laz[maxn*][];
 int ch[maxn*][],pct,X[maxn],Y[maxn];
 int N,M;

 Node operator + (Node a,Node b){
     Node p;
     p.l=a.l,p.r=b.r;
     p.x=a.x+b.x,p.y=a.y+b.y;
     p.xy=a.xy+b.xy,p.x2=a.x2+b.x2;
     return p;
 }

 inline void deal(Node &p,double *v){
     int r=p.r,l=p.l;
     if(v[]!=-inf){
         p.x=p.y=1ll*(p.r+p.l)*(p.r-p.l+)/;
         p.xy=p.x2=(1ll*r*(r+)*(*r+))/-(1ll*(l-)*l*(*l-))/;
         v[]+=v[],v[]+=v[];
     }

     p.x2+=(p.r-p.l+)*v[]*v[]+*v[]*p.x;
     p.xy+=v[]*p.y+v[]*p.x+(p.r-p.l+)*v[]*v[];
     p.x+=(p.r-p.l+)*v[],p.y+=(p.r-p.l+)*v[];
     if(v[]!=-inf) v[]-=v[],v[]-=v[];
 }

 inline void pushdown(int p){
     // return;
     if(!ch[p][]) return;
     if(laz[p][]==&&laz[p][]==&&laz[p][]==-inf) return;
     int a=ch[p][],b=ch[p][];
     if(laz[p][]!=-inf){
         laz[a][]=laz[a][]=laz[b][]=laz[b][]=;
         laz[a][]=laz[p][],laz[a][]=laz[p][];
         laz[b][]=laz[p][],laz[b][]=laz[p][];
     }
     laz[a][]+=laz[p][],laz[a][]+=laz[p][];
     laz[b][]+=laz[p][],laz[b][]+=laz[p][];

     deal(tr[a],laz[p]);
     deal(tr[b],laz[p]);
     laz[p][]=laz[p][]=,laz[p][]=laz[p][]=-inf;
 }

 void build(int &p,int l,int r){
     p=++pct;
     laz[p][]=laz[p][]=-inf;
     if(l==r){
         tr[p]=Node(X[l],Y[l],1ll*X[l]*Y[l],1ll*X[l]*X[l],l,r);
     }else{
         int m=l+r>>;
         build(ch[p][],l,m);
         build(ch[p][],m+,r);
         tr[p]=tr[ch[p][]]+tr[ch[p][]];
     }
 }

 void query(int p,int l,int r,int x,int y,Node &q){
     pushdown(p);
     if(x<=l&&r<=y){
         if(!q.l) q=tr[p];
         else q=q+tr[p];
     }else{
         int m=l+r>>;
         if(x<=m) query(ch[p][],l,m,x,y,q);
         if(y>=m+) query(ch[p][],m+,r,x,y,q);
     }
 }

 void add(int p,int l,int r,int x,int y,int s,int t){
     pushdown(p);
     if(x<=l&&r<=y){
         double v[];
         v[]=s,v[]=t;v[]=v[]=-inf;
         deal(tr[p],v);
         laz[p][]+=s,laz[p][]+=t;
         pushdown(p);
     }else{
         int m=l+r>>;
         if(x<=m) add(ch[p][],l,m,x,y,s,t);
         if(y>=m+) add(ch[p][],m+,r,x,y,s,t);
         tr[p]=tr[ch[p][]]+tr[ch[p][]];
     }
 }

 void change(int p,int l,int r,int x,int y,int s,int t){
     if(x<=l&&r<=y){
         double v[];
         v[]=v[]=,v[]=s,v[]=t;
         deal(tr[p],v);
         laz[p][]=laz[p][]=,laz[p][]=s,laz[p][]=t;
         pushdown(p);
     }else{
         pushdown(p);
         int m=l+r>>;
         if(x<=m) change(ch[p][],l,m,x,y,s,t);
         if(y>=m+) change(ch[p][],m+,r,x,y,s,t);
         tr[p]=tr[ch[p][]]+tr[ch[p][]];
     }
 }

 int main(){
     int i,j,k;
     N=rd(),M=rd();
     for(i=;i<=N;i++) X[i]=rd();
     for(i=;i<=N;i++) Y[i]=rd();
     build(i,,N);
     for(i=;i<=M;i++){
         int a=rd(),b=rd(),c=rd();
         if(a==){
             Node p;
             query(,,N,b,c,p);
             double xb=p.x/(c-b+),yb=p.y/(c-b+);
             double ans=;
             ans=p.xy-xb*p.y-yb*p.x+xb*yb*(c-b+);
             ans/=p.x2-*xb*p.x+xb*xb*(c-b+);
             printf("%.10lf\n",ans);
         }else{
             int d=rd(),e=rd();
             if(a==){
                 add(,,N,b,c,d,e);
             }else if(a==){
                 change(,,N,b,c,d,e);
             }
         }
     }
     return ;
 }