[题解]

/*
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]相关分析的更多相关文章

  1. BZOJ4817 SDOI2017 相关分析

    4821: [Sdoi2017]相关分析 Time Limit: 10 Sec  Memory Limit: 128 MBSec  Special Judge Description Frank对天文 ...

  2. [Sdoi2017]相关分析 [线段树]

    [Sdoi2017]相关分析 题意:沙茶线段树 md其实我考场上还剩一个多小时写了40分 其实当时写正解也可以吧1h也就写完了不过还要拍一下 正解代码比40分短2333 #include <io ...

  3. [题目] Luogu P3707 [SDOI2017]相关分析

    参考资料:[Luogu 3707] SDOI2017 相关分析 P3707 [SDOI2017]相关分析 TFRAC FRAC DFRAC \(\tfrac{\sum}{1}\) \(\frac{\s ...

  4. 【BZOJ4821】[SDOI2017]相关分析(线段树)

    [BZOJ4821][SDOI2017]相关分析(线段树) 题面 BZOJ 洛谷 题解 看看询问要求的东西是什么.把所有的括号拆开,不难发现要求的就是\(\sum x,\sum y,\sum xy,\ ...

  5. [Luogu 3707] SDOI2017 相关分析

    [Luogu 3707] SDOI2017 相关分析 前言 Capella 和 Frank 一样爱好天文学. 她常在冬季的夜晚,若有所思地望着东北方上空的五边形中,最为耀眼的一个顶点. 那一抹金黄曾带 ...

  6. 4821: [Sdoi2017]相关分析

    4821: [Sdoi2017]相关分析 链接 分析: 大力拆式子,化简,然后线段树.注意精度问题与爆longlong问题. 代码: #include<cstdio> #include&l ...

  7. P3707 [SDOI2017]相关分析

    P3707 [SDOI2017]相关分析 线段树裸题?但是真的很麻烦QAQ 题目给的式子是什么不用管,大力拆开,就是\(\frac{\sum x_iy_i-\overline xy_i-\overli ...

  8. AC日记——[SDOI2017]相关分析 洛谷 P3707

    [SDOI2017]相关分析 思路: 裸线段树: (玄学ac): 代码: #include <bits/stdc++.h> using namespace std; #define max ...

  9. 【BZOJ4821】[Sdoi2017]相关分析 线段树

    [BZOJ4821][Sdoi2017]相关分析 Description Frank对天文学非常感兴趣,他经常用望远镜看星星,同时记录下它们的信息,比如亮度.颜色等等,进而估算出星星的距离,半径等等. ...

随机推荐

  1. SQL Server查询某个字段存在哪些表中

    一.查询SQL Server中所有的表 SQL语句:SELECT * FROM sys.tables name列表示所有的表名. 二.查询SQL Server中所有的列 SQL语句:SELECT * ...

  2. libRTMP文档

    https://rtmpdump.mplayerhq.hu/ 原文地址:http://rtmpdump.mplayerhq.hu/librtmp.3.html git clone git://git. ...

  3. if、for、while、do 等语句自占一行

    if.for.while.do 等语句自占一行,执行语句不得紧跟其后.不论 执行语句有多少都要加{}.这样可以防止书写失误. #include <iostream> /* run this ...

  4. python小工具

    http://blog.csdn.net/pipisorry/article/details/46754515 python复制.删除文件代码.python代码出错重新启动 python遍历和删除指定 ...

  5. <img/>标签onerror事件在IE下的bug和解决方法

    IE下打开网页时,会弹出“Stack overflow at line: 0”的弹框.经分析,这个bug是由于img标签的onerror事件引起的.程序中用到的代码片段如下:正常情况下显示src所指路 ...

  6. android基站定位程序获取地理位置

    目录 一.设置界面 二.为按钮绑定事件 三.获取基站信息 四.获取经纬度 五.获取物理位置 六.显示结果 七.运行程序 八.总结 九.程序代码 正文 在Android操作系统下,基站定位其实很简单,先 ...

  7. android SDK manager 快速更新【转】

  8. ChemDraw教程:如何查看和删除俗名

    化学范畴里,允许用俗名表示ChemDraw原子标记或原子标记的一部分,可以定义俗名的快捷键也可以自由查看或删除俗名,熟练掌握可以提高ChemDraw软件使用效率,下面将具体介绍此部分内容. 一.查看俗 ...

  9. java jdk-awt.font在centos上中文乱码的问题, 安装中文字体

    有需求生成一个二维码,并且有一段文本说明,但是使用awt.font来生成中文时,一直存在乱码的问题.网上的解决办法有几种,但是在centos上亲测有用的就是如下的方法. Java代码如下:new ja ...

  10. idea-java项目配置

    导入项目后,工程结构配置: 如果不加入tomcat 运行库,项目会报servlet jar 找不到的异常 tomcat服务器配置