codevs 5967 [SDOI2017]相关分析
/*
WA://50分
last:(r-l+1)<-- (r-mid)
(r-l+1)<-- (mid-l+1)
now:int mid=l+r>>1;
tr[lc].sxy+=tr[lc].sx*addy[k]+tr[lc].sy*addx[k]+(real)(mid-l+1)*addx[k]*addy[k];
tr[lc].sx2+=(real)(mid-l+1)*addx[k]*addx[k]+2.0*tr[lc].sx*addx[k];
tr[lc].sx+=(real)(mid-l+1)*addx[k];
tr[lc].sy+=(real)(mid-l+1)*addy[k]; tr[rc].sxy+=tr[rc].sx*addy[k]+tr[rc].sy*addx[k]+(real)(r-mid)*addx[k]*addy[k];
tr[rc].sx2+=(real)(r-mid)*addx[k]*addx[k]+2.0*tr[rc].sx*addx[k];
tr[rc].sx+=(real)(r-mid)*addx[k];
tr[rc].sy+=(real)(r-mid)*addy[k];
addx[k]=addy[k]=0;
*/
#include<cmath>
#include<cstdio>
#include<cstdlib>
#include<iostream>
#define pf(x)((x)*(x))
#define lc k<<1
#define rc k<<1|1
#define setfile(x) freopen(#x".in","r",stdin);freopen(#x".out","w",stdout);
using namespace std;
typedef double real;
const real inf=2e9;
const int N=2e5+;
const int M=1e3+;
int n,m,bsize;
real ax[N],ay[N];
struct sgt{real sx,sy,sxy,sx2;}tr[N<<];
real addx[N<<],addy[N<<];
real setx[N<<],sety[N<<];
bool tag[N<<];
inline int read(){
int x=,f=;char ch=getchar();
while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}
while(ch>=''&&ch<=''){x=x*+ch-'';ch=getchar();}
return x*f;
}
sgt operator +(const sgt &a,const sgt &b){
sgt c;
c.sx=a.sx+b.sx;
c.sy=a.sy+b.sy;
c.sxy=a.sxy+b.sxy;
c.sx2=a.sx2+b.sx2;
return c;
}
inline real sumsq1(real l,real r){return (l+r)*(r-l+)/;}
inline real sumsq2(real l,real r){return (r*(r+)*(r*+)-(l-)*l*(l*-))/;}
void updata(int k){
tr[k].sx=tr[lc].sx+tr[rc].sx;
tr[k].sy=tr[lc].sy+tr[rc].sy;
tr[k].sxy=tr[lc].sxy+tr[rc].sxy;
tr[k].sx2=tr[lc].sx2+tr[rc].sx2;
}
void build(int k,int l,int r){
// tag[k]=0;
if(l==r){
tr[k].sx=ax[l];
tr[k].sy=ay[l];
tr[k].sxy=ax[l]*ay[l];
tr[k].sx2=ax[l]*ax[l];
return ;
}
int mid=l+r>>;
build(lc,l,mid);
build(rc,mid+,r);
updata(k);
}
void pushdown(int k,int l,int r){
if(l==r) return ;
int mid=l+r>>;
if(addx[k]!=||addy[k]!=){
if(tag[lc]){
setx[lc]+=addx[k];
sety[lc]+=addy[k];
}
else{
addx[lc]+=addx[k];
addy[lc]+=addy[k];
}
if(tag[rc]){
setx[rc]+=addx[k];
sety[rc]+=addy[k];
}
else{
addx[rc]+=addx[k];
addy[rc]+=addy[k];
}
int mid=l+r>>;
tr[lc].sxy+=tr[lc].sx*addy[k]+tr[lc].sy*addx[k]+(real)(mid-l+)*addx[k]*addy[k];
tr[lc].sx2+=(real)(mid-l+)*addx[k]*addx[k]+2.0*tr[lc].sx*addx[k];
tr[lc].sx+=(real)(mid-l+)*addx[k];
tr[lc].sy+=(real)(mid-l+)*addy[k]; tr[rc].sxy+=tr[rc].sx*addy[k]+tr[rc].sy*addx[k]+(real)(r-mid)*addx[k]*addy[k];
tr[rc].sx2+=(real)(r-mid)*addx[k]*addx[k]+2.0*tr[rc].sx*addx[k];
tr[rc].sx+=(real)(r-mid)*addx[k];
tr[rc].sy+=(real)(r-mid)*addy[k];
addx[k]=addy[k]=;
}
if(tag[k]){
tag[lc]=tag[rc]=;
setx[lc]=setx[k];
setx[rc]=setx[k];
sety[lc]=sety[k];
sety[rc]=sety[k]; real li1,li2,ri1,ri2;
li1=sumsq1(l,mid);
li2=sumsq2(l,mid);
ri1=sumsq1(mid+,r);
ri2=sumsq2(mid+,r);
tr[lc].sxy=setx[k]*sety[k]*(real)(mid-l+);
tr[lc].sx2=setx[k]*setx[k]*(real)(mid-l+); tr[lc].sx=setx[k]*(real)(mid-l+)+li1;
tr[lc].sy=sety[k]*(real)(mid-l+)+li1;
tr[lc].sxy+=li1*(setx[k]+sety[k])+li2;
tr[lc].sx2+=*li1*setx[k]+li2; tr[rc].sxy=setx[k]*sety[k]*(real)(r-mid);
tr[rc].sx2=setx[k]*setx[k]*(real)(r-mid);
tr[rc].sx=setx[k]*(real)(r-mid)+ri1;
tr[rc].sy=sety[k]*(real)(r-mid)+ri1;
tr[rc].sxy+=ri1*(setx[k]+sety[k])+ri2;
tr[rc].sx2+=*ri1*setx[k]+ri2;
tag[k]=;
} }
void change(int k,int l,int r,int x,int y,real s,real t){
if(l==x&&r==y){
tr[k].sxy+=tr[k].sx*t+tr[k].sy*s+(real)(r-l+)*s*t;
tr[k].sx2+=(real)(r-l+)*s*s+2.0*tr[k].sx*s;
tr[k].sx+=(real)(r-l+)*s;
tr[k].sy+=(real)(r-l+)*t;
if(tag[k]){
setx[k]+=s;
sety[k]+=t;
}
else{
addx[k]+=s;
addy[k]+=t;
}
return ;
}
pushdown(k,l,r);
int mid=l+r>>;
if(y<=mid) change(lc,l,mid,x,y,s,t);
else if(x>mid) change(rc,mid+,r,x,y,s,t);
else change(lc,l,mid,x,mid,s,t),change(rc,mid+,r,mid+,y,s,t);
updata(k);
}
void cover(int k,int l,int r,int x,int y,real s,real t){ if(l==x&&r==y){
real li1,li2;
li1=sumsq1(l,r);
li2=sumsq2(l,r);
tr[k].sxy=s*t*(real)(r-l+);
tr[k].sx2=s*s*(real)(r-l+);
tr[k].sx=s*(real)(r-l+)+li1;
tr[k].sy=t*(real)(r-l+)+li1;
tr[k].sxy+=li1*(s+t)+li2;
tr[k].sx2+=*li1*s+li2; tag[k]=;
setx[k]=s;
sety[k]=t;
addx[k]=addy[k]=;
return ;
}
pushdown(k,l,r);
int mid=l+r>>;
if(y<=mid) cover(lc,l,mid,x,y,s,t);
else if(x>mid) cover(rc,mid+,r,x,y,s,t);
else cover(lc,l,mid,x,mid,s,t),cover(rc,mid+,r,mid+,y,s,t);
updata(k);
}
sgt query(int k,int l,int r,int x,int y){
if(l==x&&r==y) return tr[k];
pushdown(k,l,r);
int mid=l+r>>;
if(y<=mid) return query(lc,l,mid,x,y);
else if(x>mid) return query(rc,mid+,r,x,y);
else return query(lc,l,mid,x,mid)+query(rc,mid+,r,mid+,y);
}
real find(int x,int y){
real X=,Y=,XY=,X2=,Xp,Yp,fz=,fm=,A;
sgt ans=query(,,n,x,y);int l=x,r=y;
X=ans.sx;
Y=ans.sy;
XY=ans.sxy;
X2=ans.sx2;
Xp=X/(real)(r-l+);
Yp=Y/(real)(r-l+);
fz=XY-X*Yp-Xp*Y+Xp*Yp*(r-l+);
fm=X2-2.0*X*Xp+Xp*Xp*(r-l+);
// if(fz==0||fm==0) return 0;
return fz/fm;
}
void work(){
build(,,n);
real s,t;
for(int i=,opt,x,y;i<=m;i++){
opt=read();x=read();y=read();
if(opt==) printf("%.7lf\n",find(x,y));
if(opt==) scanf("%lf%lf",&s,&t),change(,,n,x,y,s,t);
if(opt==) scanf("%lf%lf",&s,&t),cover(,,n,x,y,s,t);
}
}
int main(){
setfile(relative);
n=read();m=read();
for(int i=;i<=n;i++) scanf("%lf",&ax[i]);
for(int i=;i<=n;i++) scanf("%lf",&ay[i]);
work();
fclose(stdin);
fclose(stdout);
return ;
}
codevs 5967 [SDOI2017]相关分析的更多相关文章
- BZOJ4817 SDOI2017 相关分析
4821: [Sdoi2017]相关分析 Time Limit: 10 Sec Memory Limit: 128 MBSec Special Judge Description Frank对天文 ...
- [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对天文学非常感兴趣,他经常用望远镜看星星,同时记录下它们的信息,比如亮度.颜色等等,进而估算出星星的距离,半径等等. ...
随机推荐
- MS-SQL 删除数据库所有的表
godeclare @tbname varchar(250)declare #tb cursor for select name from sysobjects where objectpropert ...
- JavaScript(一):JavaScript简介
一.什么是JavaScript JavaScript是一种具有面向对象能力的.解释性的程序设计语言.更具体一点,它是基于对象和事件驱动并具有相对安全性的客户端脚本语言.因为他不需要在一个语言环境下运行 ...
- MATLAB 安装使用libsvm详细步骤
根据本文后面部分博友提出的在配置过程中出现的问题,其中需要特别强调的一点:整个过程,都是在 libsvm-3.12\matlab目录下操作的.如果这一点你忽视了,你不可能解决配置中报的Bug,即使重新 ...
- delete指针以后应赋值为NULL——QT deletelater指针以后也同样要马上赋值为NULL
delete p后,只是释放了指针指向的内存空间.p并不会自动被置为NULL,而且指针还在,同时还指向了之前的地址 delete NULL编译器不会报错(因为delete空指针是合法的) 例: 对一个 ...
- 网络协议之bt---bt协议详解 DHT篇(下)
-------------------------author:pkf -------------------------------qq:1327706646 ------------------- ...
- CentOS查看登录用户以及踢出用户
查看登录用户,使用w命令 [root@lnmp ~]# w 18:51:18 up 35 min, 2 users, load average: 0.00, 0.00, 0.00 USER ...
- 今天搞log4net插入错误日志去mysql数据库的时候出现了点问题,已解决。记录下解决方案
先上图 配置log4net的时候要填这项,可是这个value我不知道啊.....上图里的value是我用下面的方法获取的 MySqlConnection con = new MySqlConnecti ...
- jQuery.ajax实现根据不同的Content-Type做出不同的响应
使用H5+ASP.NET General Handler开发项目,使用ajax进行前后端的通讯.有一个场景需求是根据服务器返回的不同数据类型,前端进行不同的响应,这里记录下如何使用$.ajax实现该需 ...
- [转]wcout输出中文却不显示出来
准备使用UNICODE来写个控制台测试程序发现,cout无法输出UNICODE的中文字符.查找C++标准看到,其提供了wcin.wcout.wcerr.wclog用于处理wchar_t字符的输入输出. ...
- css两列等高布局
布局方案 等高布局有几种不同的方法,但目前为止我认为浏览器兼容最好最简便的应该是padding补偿法.首先把列的padding-bottom设为一个足够大的值,再把列的margin-bottom设一个 ...