BZOJ4817 SDOI2017 相关分析
4821: [Sdoi2017]相关分析
Time Limit: 10 Sec Memory Limit: 128 MBSec Special Judge
Description
Input
Output
Sample Input
1 2 3
1 2 3
1 1 3
2 2 3 -3 2
1 1 2
3 1 2 2 1
1 1 3
Sample Output
-1.5000000000
-0.6153846154

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <vector>
#include <cstring>
#include <queue>
#include <complex>
#include <stack>
#define LL long long int
#define dob long double
#define ls (x<<1)
#define rs (x<<1|1)
using namespace std; const int N = ;
struct Tree{
dob x,y,xx,xy;
Tree operator +(const Tree &t){
return (Tree){x+t.x,y+t.y,xx+t.xx,xy+t.xy};
}
}Tr[N*];
int n,m,lazy_vis[N];
dob X[N/],Y[N/],lazy_add1[N],lazy_add2[N],lazy_set1[N],lazy_set2[N]; inline int gi(){
int x=,res=;char ch=getchar();
while(ch>''||ch<''){if(ch=='-')res*=-;ch=getchar();}
while(ch<=''&&ch>='')x=x*+ch-,ch=getchar();
return x*res;
} inline void build(int x,int l,int r){
if(l==r){
Tr[x]=(Tree){X[l],Y[l],1.0*X[l]*X[l],1.0*X[l]*Y[l]};
return;
}
int mid=(l+r)>>;
build(ls,l,mid);build(rs,mid+,r);
Tr[x]=Tr[ls]+Tr[rs];
} inline dob calc(dob l,dob r){
return 0.5*(l+r)*(r-l+);
} inline dob calcpow(dob l,dob r){
l-=;
dob p1=1.0*(r)*(r+)*(*r+)/6.0;
dob p2=1.0*(l)*(l+)*(*l+)/6.0;
return p1-p2;
} inline void down(int x,int l,int r){
int mid=(l+r)>>,sl=mid-l+,sr=r-mid;
if(lazy_vis[x]){
lazy_add1[ls]=lazy_add1[rs]=lazy_add2[ls]=lazy_add2[rs]=;
lazy_vis[ls]=lazy_vis[rs]=;
dob S=lazy_set1[x],T=lazy_set2[x];
lazy_set1[ls]=lazy_set1[rs]=lazy_set1[x];
lazy_set2[ls]=lazy_set2[rs]=lazy_set2[x];
Tr[ls].xx=1.0*sl*S*S+2.0*S*calc(l,mid)+calcpow(l,mid);
Tr[rs].xx=1.0*sr*S*S+2.0*S*calc(mid+,r)+calcpow(mid+,r);
Tr[ls].xy=1.0*sl*S*T+1.0*(S+T)*calc(l,mid)+calcpow(l,mid);
Tr[rs].xy=1.0*sr*S*T+1.0*(S+T)*calc(mid+,r)+calcpow(mid+,r);
Tr[ls].x=1.0*sl*S+calc(l,mid);Tr[rs].x=1.0*sr*S+calc(mid+,r);
Tr[ls].y=1.0*sl*T+calc(l,mid);Tr[rs].y=1.0*sr*T+calc(mid+,r);
lazy_vis[x]=;
}
if(lazy_add1[x] || lazy_add2[x]){
dob S=lazy_add1[x],T=lazy_add2[x];
lazy_add1[ls]+=S;lazy_add1[rs]+=S;
lazy_add2[ls]+=T;lazy_add2[rs]+=T;
Tr[ls].xx+=2.0*Tr[ls].x*S+1.0*sl*S*S;
Tr[rs].xx+=2.0*Tr[rs].x*S+1.0*sr*S*S;
Tr[ls].xy+=1.0*Tr[ls].x*T+1.0*Tr[ls].y*S+1.0*sl*S*T;
Tr[rs].xy+=1.0*Tr[rs].x*T+1.0*Tr[rs].y*S+1.0*sr*S*T;
Tr[ls].x+=1.0*sl*S;Tr[rs].x+=1.0*sr*S;
Tr[ls].y+=1.0*sl*T;Tr[rs].y+=1.0*sr*T;
lazy_add1[x]=lazy_add2[x]=;
}
} inline Tree query_1(int x,int l,int r,int xl,int xr){
if(xl<=l && r<=xr)return Tr[x];
down(x,l,r);int mid=(l+r)>>;
if(xr<=mid)return query_1(ls,l,mid,xl,xr);
else if(xl>mid)return query_1(rs,mid+,r,xl,xr);
return query_1(ls,l,mid,xl,mid)+query_1(rs,mid+,r,mid+,xr);
} inline void update_2(int x,int l,int r,int xl,int xr,dob S,dob T){
if(xl<=l && r<=xr){
lazy_add1[x]+=S;lazy_add2[x]+=T;
Tr[x].xx+=2.0*Tr[x].x*S+1.0*(r-l+)*S*S;
Tr[x].xy+=1.0*Tr[x].x*T+1.0*Tr[x].y*S+1.0*(r-l+)*S*T;
Tr[x].x+=1.0*(r-l+)*S;Tr[x].y+=1.0*(r-l+)*T;
return;
}
down(x,l,r);int mid=(l+r)>>;
if(xr<=mid)update_2(ls,l,mid,xl,xr,S,T);
else if(xl>mid)update_2(rs,mid+,r,xl,xr,S,T);
else update_2(ls,l,mid,xl,mid,S,T),update_2(rs,mid+,r,mid+,xr,S,T);
Tr[x]=Tr[ls]+Tr[rs];
} inline void update_3(int x,int l,int r,int xl,int xr,dob S,dob T){
if(xl<=l && r<=xr){
lazy_add1[x]=lazy_add2[x]=;
lazy_vis[x]=;lazy_set1[x]=S;lazy_set2[x]=T;
Tr[x].xx=1.0*(r-l+)*S*S+2.0*S*calc(1.0*l,1.0*r)+calcpow(1.0*l,1.0*r);
Tr[x].xy=1.0*(r-l+)*S*T+1.0*(S+T)*calc(l,r)+calcpow(l,r);
Tr[x].x=1.0*(r-l+)*S+calc(l,r);Tr[x].y=1.0*(r-l+)*T+calc(l,r);
return;
}
down(x,l,r);int mid=(l+r)>>;
if(xr<=mid)update_3(ls,l,mid,xl,xr,S,T);
else if(xl>mid)update_3(rs,mid+,r,xl,xr,S,T);
else update_3(ls,l,mid,xl,mid,S,T),update_3(rs,mid+,r,mid+,xr,S,T);
Tr[x]=Tr[ls]+Tr[rs];
} int main()
{
/*freopen(".in","r",stdin);
freopen(".out","w",stdout);*/
n=gi();m=gi();
for(int i=;i<=n;++i)X[i]=gi();
for(int i=;i<=n;++i)Y[i]=gi();
build(,,n);
for(int i=;i<=m;++i){
int type=gi();
if(type==){
int l=gi(),r=gi();
Tree ans=query_1(,,n,l,r);
dob fz=ans.xy-ans.x*ans.y/(r-l+);
dob fm=ans.xx-ans.x*ans.x/(r-l+);
printf("%.10Lf\n",fz/fm);
}
if(type==){
int l=gi(),r=gi(),S=gi(),T=gi();
update_2(,,n,l,r,1.0*S,1.0*T);
}
if(type==){
int l=gi(),r=gi(),S=gi(),T=gi();
update_3(,,n,l,r,1.0*S,1.0*T);
}
} /*fclose(stdin);
fclose(stdout);*/
return ;
}
相关分析
BZOJ4817 SDOI2017 相关分析的更多相关文章
- [Sdoi2017]相关分析 [线段树]
[Sdoi2017]相关分析 题意:沙茶线段树 md其实我考场上还剩一个多小时写了40分 其实当时写正解也可以吧1h也就写完了不过还要拍一下 正解代码比40分短2333 #include <io ...
- [题目] Luogu P3707 [SDOI2017]相关分析
参考资料:[Luogu 3707] SDOI2017 相关分析 P3707 [SDOI2017]相关分析 TFRAC FRAC DFRAC \(\tfrac{\sum}{1}\) \(\frac{\s ...
- 【BZOJ4821】[SDOI2017]相关分析(线段树)
[BZOJ4821][SDOI2017]相关分析(线段树) 题面 BZOJ 洛谷 题解 看看询问要求的东西是什么.把所有的括号拆开,不难发现要求的就是\(\sum x,\sum y,\sum xy,\ ...
- [Luogu 3707] SDOI2017 相关分析
[Luogu 3707] SDOI2017 相关分析 前言 Capella 和 Frank 一样爱好天文学. 她常在冬季的夜晚,若有所思地望着东北方上空的五边形中,最为耀眼的一个顶点. 那一抹金黄曾带 ...
- 4821: [Sdoi2017]相关分析
4821: [Sdoi2017]相关分析 链接 分析: 大力拆式子,化简,然后线段树.注意精度问题与爆longlong问题. 代码: #include<cstdio> #include&l ...
- P3707 [SDOI2017]相关分析
P3707 [SDOI2017]相关分析 线段树裸题?但是真的很麻烦QAQ 题目给的式子是什么不用管,大力拆开,就是\(\frac{\sum x_iy_i-\overline xy_i-\overli ...
- AC日记——[SDOI2017]相关分析 洛谷 P3707
[SDOI2017]相关分析 思路: 裸线段树: (玄学ac): 代码: #include <bits/stdc++.h> using namespace std; #define max ...
- 【BZOJ4821】[Sdoi2017]相关分析 线段树
[BZOJ4821][Sdoi2017]相关分析 Description Frank对天文学非常感兴趣,他经常用望远镜看星星,同时记录下它们的信息,比如亮度.颜色等等,进而估算出星星的距离,半径等等. ...
- (WA)BZOJ 4821: [Sdoi2017]相关分析
二次联通门 : BZOJ 4821: [Sdoi2017]相关分析 2017.8.23 Updata 妈妈!!这道题卡我!!!就是不然我过!!!!! #include <cstdio> # ...
随机推荐
- spring boot自定义log4j2日志文件
背景:因为从 spring boot 1.4开始的版本就要用log4j2 了,支持的格式有json和xml两种格式,此次实践主要使用的是xml的格式定义日志说明. spring boot 1.5.8. ...
- memcached内存管理
前言 memcached默认情况下采用了名为Slab Allocator的机制来管理内存.在该机制出现以前,内存的分配是通过对所有记录简单地进行malloc和free来进行的.但是,这种方式会导致内存 ...
- 08-图8 How Long Does It Take
原题: Given the relations of all the activities of a project, you are supposed to find the earliest co ...
- km算法入门
本文知识均由笔者自学,文章有错误之处请不吝指出. 笔者刷数模题的时候有一道题考到了"二分图最大权分配",需要用到KM算法,但是书上对KM算法的介绍又臭又长,更何况有些同学" ...
- MySql基础入门-mysql的结构层次
了解MySql必须牢牢记住其体系结构图,Mysql是由SQL接口,解析器,优化器,缓存,存储引擎组成的. 1.Connectors指的是不同语言中与SQL的交互 2.Management Servei ...
- linux上mysql安装与卸载
以下步骤运行环境是centos6.5 1.查找以前是否装有mysql命令:rpm -qa|grep -i mysql2.删除mysql删除命令:rpm -e --nodeps 包名3.删除老版本m ...
- windows配置git
每次要使用git指令的时候都要去打开git bash 操作,太麻烦,要想直接在dos窗口下使用git指令需要再进行如下环境变量配置. 1.系统环境变量path添加:D://programFiles/g ...
- WPF水珠效果按钮组
效果图 相关知识 这部分基本就是废话,网上都能找到,我只不过是整理了以下.建议先不看,用到的时候可以回来看看 贝塞尔曲线 先来看两组图,有助于理解什么是贝塞尔曲线(图片取自维基百科,参考链接1) 二次 ...
- 使用hiredis实现pipeline方式访问
1.介绍 hiredis: 一个c/c++的访问redis的api库 地址:https://github.com/redis/hiredis pipeline方式: redis中的pipeline方式 ...
- Ext表单验证
//大家在很多的extjs代码中都看到了这两个,他们都起提示作用的Ext.QuickTips.init();//支持tips提示Ext.form.Field.prototype.msgTarget=' ...