SDOI2017相关分析 线段树
题目
思路
\]
\]
\]
\]
\]
\]
\]
\]
\]
\]
\]
\]
其实不用这么麻烦的、、
好了,剩下的去维护吧
好吧,我太菜了
要开long doule
代码
#include <bits/stdc++.h>
#define ll long double
#define ls rt<<1
#define rs rt<<1|1
using namespace std;
const int N=2e5+7;
int read() {
int x=0,f=1;char s=getchar();
for (;s>'9'||s<'0';s=getchar()) if(s=='-') f=-1;
for (;s>='0'&&s<='9';s=getchar()) x=x*10+s-'0';
return x*f;
}
ll x[N],y[N];
struct node {
int l,r,siz;
ll tot[2],pingfang,chengji;
ll lazy,S,T;
}e[N<<2];
void pushup(int rt) {
e[rt].tot[0]=e[ls].tot[0]+e[rs].tot[0];
e[rt].tot[1]=e[ls].tot[1]+e[rs].tot[1];
e[rt].pingfang=e[ls].pingfang+e[rs].pingfang;
e[rt].chengji=e[ls].chengji+e[rs].chengji;
}
ll calc(int x) {return (ll)x*(x+1)/2;};
ll calc2(int x) {return (ll)x*(x+1)/2*(2*x+1)/3;};
void pushdown(int rt) {
if(e[rt].lazy) {
e[ls].tot[0]=e[ls].tot[1]=calc(e[ls].r)-calc(e[ls].l-1);
e[ls].pingfang=e[ls].chengji=calc2(e[ls].r)-calc2(e[ls].l-1);
e[ls].S=e[ls].T=0;
e[ls].lazy=1;
e[rs].tot[0]=e[rs].tot[1]=calc(e[rs].r)-calc(e[rs].l-1);
e[rs].pingfang=e[rs].chengji=calc2(e[rs].r)-calc2(e[rs].l-1);
e[rs].S=e[rs].T=0;
e[rs].lazy=1;
e[rt].lazy=0;
}
if(e[rt].S||e[rt].T) {
e[ls].chengji+=e[ls].tot[0]*e[rt].T+e[ls].tot[1]*e[rt].S+e[rt].S*e[rt].T*e[ls].siz;
e[ls].pingfang+=e[ls].tot[0]*2*e[rt].S+e[rt].S*e[rt].S*e[ls].siz;
e[ls].tot[0]+=e[ls].siz*e[rt].S;
e[ls].tot[1]+=e[ls].siz*e[rt].T;
e[ls].S+=e[rt].S;
e[ls].T+=e[rt].T;
e[rs].chengji+=e[rs].tot[0]*e[rt].T+e[rs].tot[1]*e[rt].S+e[rt].S*e[rt].T*e[rs].siz;
e[rs].pingfang+=e[rs].tot[0]*2*e[rt].S+e[rt].S*e[rt].S*e[rs].siz;
e[rs].tot[0]+=e[rs].siz*e[rt].S;
e[rs].tot[1]+=e[rs].siz*e[rt].T;
e[rs].S+=e[rt].S;
e[rs].T+=e[rt].T;
e[rt].S=e[rt].T=0;
}
}
void build(int l,int r,int rt) {
e[rt].l=l,e[rt].r=r,e[rt].siz=r-l+1;
if(l==r) {
e[rt].tot[0]=x[l];
e[rt].tot[1]=y[l];
e[rt].pingfang=x[l]*x[l];
e[rt].chengji=x[l]*y[l];
return;
}
int mid=(l+r)>>1;
build(l,mid,ls);
build(mid+1,r,rs);
pushup(rt);
}
void modify_1(int L,int R,ll S,ll T,int rt) {
if(L<=e[rt].l&&e[rt].r<=R) {
e[rt].chengji+=e[rt].tot[0]*T+e[rt].tot[1]*S+S*T*e[rt].siz;
e[rt].pingfang+=e[rt].tot[0]*2*S+S*S*e[rt].siz;
e[rt].tot[0]+=e[rt].siz*S;
e[rt].tot[1]+=e[rt].siz*T;
e[rt].S+=S;
e[rt].T+=T;
return;
}
pushdown(rt);
int mid=(e[rt].l+e[rt].r)>>1;
if(L<=mid) modify_1(L,R,S,T,ls);
if(R>mid) modify_1(L,R,S,T,rs);
pushup(rt);
}
void modify_2(int L,int R,int rt) {
if(L<=e[rt].l&&e[rt].r<=R) {
e[rt].tot[0]=e[rt].tot[1]=calc(e[rt].r)-calc(e[rt].l-1);
e[rt].pingfang=e[rt].chengji=calc2(e[rt].r)-calc2(e[rt].l-1);
e[rt].S=e[rt].T=0;
e[rt].lazy=1;
return;
}
pushdown(rt);
int mid=(e[rt].l+e[rt].r)>>1;
if(L<=mid) modify_2(L,R,ls);
if(R>mid) modify_2(L,R,rs);
pushup(rt);
}
ll query(int L,int R,int opt,int rt) {
if(L<=e[rt].l&&e[rt].r<=R) {
if(opt==0) return e[rt].tot[0];
if(opt==1) return e[rt].tot[1];
if(opt==2) return e[rt].pingfang;
if(opt==3) return e[rt].chengji;
}
pushdown(rt);
int mid=(e[rt].l+e[rt].r)>>1;
ll ans=0;
if(L<=mid) ans+=query(L,R,opt,ls);
if(R>mid) ans+=query(L,R,opt,rs);
pushup(rt);
return ans;
}
int main() {
int n=read(),m=read();
for(int i=1;i<=n;++i) x[i]=read();
for(int i=1;i<=n;++i) y[i]=read();
build(1,n,1);
while(m--) {
int opt=read(),L=read(),R=read();
if(opt==1) {
ll tot_x=query(L,R,0,1);
ll tot_y=query(L,R,1,1);
ll tot_x_x=query(L,R,2,1);
ll tot_x_y=query(L,R,3,1);
long double a=tot_x_y-tot_x*tot_y/(R-L+1);
long double b=tot_x_x-tot_x*tot_x/(R-L+1);
printf("%.10Lf\n",(long double)a/b);
} else if(opt==2) {
ll S=read(),T=read();
modify_1(L,R,(ll)S,(ll)T,1);
} else {
ll S=read(),T=read();
modify_2(L,R,1);
modify_1(L,R,(ll)S,(ll)T,1);
}
}
return 0;
}
SDOI2017相关分析 线段树的更多相关文章
- [Sdoi2017]相关分析 [线段树]
[Sdoi2017]相关分析 题意:沙茶线段树 md其实我考场上还剩一个多小时写了40分 其实当时写正解也可以吧1h也就写完了不过还要拍一下 正解代码比40分短2333 #include <io ...
- 【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\). 操作三可以看成初始化一遍,然后同操作二. ...
- 【BZOJ4821】【SDOI2017】相关分析 [线段树]
相关分析 Time Limit: 10 Sec Memory Limit: 128 MB[Submit][Status][Discuss] Description Frank对天文学非常感兴趣,他经 ...
- luogu3707 相关分析 (线段树)
把式子展开以后会发现,可以用线段树维护$x,y,x*y,x^2$分别的区间和 然后操作有区间加和区间修改 这个pushdown的时候,如果改和加的标记同时存在,那一定是先改再加,要不然加的标记已经被清 ...
- LOJ #2005. 「SDOI2017」相关分析 线段树维护回归直线方程
题目描述 \(Frank\) 对天文学非常感兴趣,他经常用望远镜看星星,同时记录下它们的信息,比如亮度.颜色等等,进而估算出星星的距离,半径等等. \(Frank\) 不仅喜欢观测,还喜欢分析观测到的 ...
随机推荐
- 软工网络15团队作业4——Alpha阶段敏捷冲刺7.0
1.每天举行站立式会议,提供当天站立式会议照片一张. 2.项目每个成员的昨天进展.存在问题.今天安排. 成员 昨天已完成 今天计划完成 郭炜埕 实现前端各界面的跳转连接 学习后端相关知识 郑晓丽 完善 ...
- FSDB Dumper
FSDB:Fast Signal Database 相比较于VCD文件,FSDB文件的大小比VCD波形小5-50倍. 各家的仿真器都支持在simulation的过程中,直接生成FSDB文件 将VCD文 ...
- 五 js对象简介
对象简介 js中没有"类"的概念,只有对象. A:对象声明方式有三种 ------------1.调用Object函数创建对象: var person = new Object; ...
- python -- 判断给定的参数是否是地理位置的经度和纬度
例子: coordinates = "2.3,-8ui" 无效地理位置 coordinates = "2.3,-8i" 有效地理位置 方法一:使用try...e ...
- Gardener Bo (树剖 + 问题分解)
一开始没看懂计算答案的第四部和update2,很迷.然后一直推敲之后才发下我计算的时候漏掉一个关键点.没有把加值的影响放到父节点上. #include<bits/stdc++.h> #de ...
- 即时通信系统中实现全局系统通知,并与Web后台集成【附C#开源即时通讯系统(支持广域网)——QQ高仿版IM最新源码】
像QQ这样的即时通信软件,时不时就会从桌面的右下角弹出一个小窗口,或是显示一个广告.或是一个新闻.或是一个公告等.在这里,我们将其统称为“全局系统通知”.很多使用C#开源即时通讯系统——GGTalk的 ...
- FastDFS:Java客户都实现文件的上传、下载、修改、删除
客户端版本:fastdfs_client_v1.24.jar 配置文件 connect_timeout = 200 network_timeout = 3000 charset = UTF-8 htt ...
- Java8函数式编程探秘
引子 将行为作为数据传递 怎样在一行代码里同时计算一个列表的和.最大值.最小值.平均值.元素个数.奇偶分组.指数.排序呢? 答案是思维反转!将行为作为数据传递. 文艺青年的代码如下所示: public ...
- python 文件路径名,文件名,后缀名的操作
需要使用路径名来获取文件名,目录名,绝对路径等等. 使用os.path 模块中的函数来操作路径名.下面是一个交互式例子来演示一些关键的特性: >>> import os >&g ...
- 转:【专题十二】实现一个简单的FTP服务器
引言: 休息一个国庆节后好久没有更新文章了,主要是刚开始休息完心态还没有调整过来的, 现在差不多进入状态了, 所以继续和大家分享下网络编程的知识,在本专题中将和大家分享如何自己实现一个简单的FTP服务 ...