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 ;
}
luogu3707 相关分析 (线段树)的更多相关文章
- [Sdoi2017]相关分析 [线段树]
[Sdoi2017]相关分析 题意:沙茶线段树 md其实我考场上还剩一个多小时写了40分 其实当时写正解也可以吧1h也就写完了不过还要拍一下 正解代码比40分短2333 #include <io ...
- 【BZOJ4821】【SDOI2017】相关分析 [线段树]
相关分析 Time Limit: 10 Sec Memory Limit: 128 MB[Submit][Status][Discuss] Description Frank对天文学非常感兴趣,他经 ...
- 【BZOJ4821】[Sdoi2017]相关分析 线段树
[BZOJ4821][Sdoi2017]相关分析 Description Frank对天文学非常感兴趣,他经常用望远镜看星星,同时记录下它们的信息,比如亮度.颜色等等,进而估算出星星的距离,半径等等. ...
- BZOJ 4821 [Sdoi2017]相关分析 ——线段树
打开题面,看到许多$\sum$ woc,好神啊,SDOI好强啊 然后展开之后,woc,SDOI好弱啊,怎么T3出个线段树裸题啊. 最后写代码的时候,woc,SDOI怎么出个这么码农的题啊,怎么调啊. ...
- 洛谷P3707 [SDOI2017]相关分析(线段树)
题目描述 Frank对天文学非常感兴趣,他经常用望远镜看星星,同时记录下它们的信息,比如亮度.颜色等等,进而估算出星星的距离,半径等等. Frank不仅喜欢观测,还喜欢分析观测到的数据.他经常分析两个 ...
- BZOJ 4821: [Sdoi2017]相关分析 线段树 + 卡精
考试的时候切掉了,然而卡精 + 有一个地方忘开 $long long$,完美挂掉 $50$pts. 把式子化简一下,然后直接拿线段树来维护即可. Code: // luogu-judger-enabl ...
- BZOJ.4821.[SDOI2017]相关分析(线段树)
BZOJ LOJ 洛谷 恶心的拆式子..然后就是要维护\(\sum x_i,\ \sum y_i,\ \sum x_iy_i,\ \sum x_i^2\). 操作三可以看成初始化一遍,然后同操作二. ...
- SDOI2017相关分析 线段树
题目 https://loj.ac/problem/2005 思路 \[ \sum_{L}^{R}{(x_i-x)^{2}} \] \[ \sum_{L}^{R}{(x_i^2-2*x_i*x+x^{ ...
- LOJ #2005. 「SDOI2017」相关分析 线段树维护回归直线方程
题目描述 \(Frank\) 对天文学非常感兴趣,他经常用望远镜看星星,同时记录下它们的信息,比如亮度.颜色等等,进而估算出星星的距离,半径等等. \(Frank\) 不仅喜欢观测,还喜欢分析观测到的 ...
随机推荐
- 【持续更新中···】Linux下的小技巧
1.Linux回到上级文件的命令: cd ..回到上一级目录(注意:cd 和..中间有空格) cd ~回到home目录 cd -回到某一目录
- .net core实践系列之短信服务-Api的SDK的实现与测试
前言 上一篇<.net core实践系列之短信服务-Sikiro.SMS.Api服务的实现>讲解了API的设计与实现,本篇主要讲解编写接口的SDK编写还有API的测试. 或许有些人会认为, ...
- REST-framework快速构建API--源码解析
一.APIView 通过APIView实现API的过程如下: urls.py url(r'^books/$', views.BookView.as_view(),name="books&qu ...
- Quartz.net 定时任务之储存与持久化和集群(源码)
一.界面 1.这篇博客不上教程.直接看结果(包括把quartz任务转换成Windows服务) (1).主界面 (2).添加任务(默认执行) (3).编辑(默认开启) (4).关闭和开启 2.代码说明 ...
- ActiveMQ 填坑记
前言 MQ是现在大型系统架构中必不可少的一个重要中间件,之前有偏文章<MQ(消息队列)常见的应用场景解析>介绍过MQ的应用场景,现在流行的几个MQ是rabbitmq,rocketma,ka ...
- ExtJS初探:在项目中使用ExtJS
注意:本文写作时间是 2013 年,所讲的 ExtJS 如今早已过时,请勿学习! -------------------------------- 今天ExtJS官网发布了ExtJS最新正式版4.2. ...
- 关于git的一些体会:
周忠贤github链接:https://github.com/zhouzhongxian git学习心得:通过这次的学习,体会到了许多东西只要你用心去做,就没有什么做不成,,这次体会到了网上学习的重要 ...
- 第七组团队项目——专业课程资源共享平台——需求分析&原型设计
一.项目目标.定位需求: (1)目标:在教师.学生之间建立一个综合的.全面的.快捷的.高效的免费课程和学习资源共享.交流与推荐的开放性平台,实现多维和动态的推荐与分类检索服务. (2)定位:学生与教师 ...
- ubuntu16.04下载安装navicate
1.下载试用版本地址: https://www.navicat.com.cn/download/navicat-premium 2.解压缩 tar -zxvf /home/rain/download ...
- FreeMaker使用HashMap
private Map<String, Object> variables; <input type="hidden" id="tongzhisbm&q ...