BZOJ 4821 [Sdoi2017]相关分析 ——线段树
打开题面,看到许多$\sum$
woc,好神啊,SDOI好强啊
然后展开之后,woc,SDOI好弱啊,怎么T3出个线段树裸题啊。
最后写代码的时候,woc,SDOI怎么出个这么码农的题啊,怎么调啊。
想想自己加上考场$Debuff$是肯定写不出来的。
#include <map>
#include <cmath>
#include <queue>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
#define F(i,j,k) for (int i=j;i<=k;++i)
#define D(i,j,k) for (int i=j;i>=k;--i)
//#define double long double
#define ll long long
#define mp make_pair
#define maxn 400005 #define sum(x) (1LL*(x)*((x)+1)/2)
#define sum2(x) (1LL*(x)*((x)+1)*(2*(x)+1)/6) double sum_x[maxn],sum_y[maxn],sum_xy[maxn],sum_xx[maxn];
double tag_x[maxn],tag_y[maxn];
bool tag_c[maxn];
int x[maxn],y[maxn];
double sumx,sumy,sumxy,sumxx; int n,m; void update(int o,int l,int r)
{
sum_x[o]=sum_x[o<<1]+sum_x[o<<1|1];
sum_y[o]=sum_y[o<<1]+sum_y[o<<1|1];
sum_xy[o]=sum_xy[o<<1]+sum_xy[o<<1|1];
sum_xx[o]=sum_xx[o<<1]+sum_xx[o<<1|1];
} void add_tag_x(int o,int l,int r,double s)
{
tag_x[o]+=s;
sum_xx[o]+=1LL*s*s*(r-l+1)+2*s*sum_x[o];
sum_xy[o]+=1LL*sum_y[o]*s;
sum_x[o]+=1LL*(r-l+1)*s;
} void add_tag_y(int o,int l,int r,double t)
{
tag_y[o]+=t;
sum_xy[o]+=1LL*sum_x[o]*t;
sum_y[o]+=1LL*(r-l+1)*t;
} void add_cov(int o,int l,int r)
{
sum_y[o]=sum_x[o]=sum(r)-sum(l-1);
sum_xy[o]=sum_xx[o]=sum2(r)-sum2(l-1);
tag_x[o]=tag_y[o]=0; tag_c[o]=true;
} void pushdown(int o,int l,int r)
{
if (o==0||l==r) return;
int mid=l+r>>1;
if (tag_c[o])
{
add_cov(o<<1,l,mid);
add_cov(o<<1|1,mid+1,r);
tag_c[o]=false;
}
if (tag_x[o])
{
add_tag_x(o<<1,l,mid,tag_x[o]);
add_tag_x(o<<1|1,mid+1,r,tag_x[o]);
tag_x[o]=0;
}
if (tag_y[o])
{
add_tag_y(o<<1,l,mid,tag_y[o]);
add_tag_y(o<<1|1,mid+1,r,tag_y[o]);
tag_y[o]=0;
}
} void st_modify(int o,int l,int r,int L,int R,double s,double t)
{
if (L<=l&&r<=R){add_tag_x(o,l,r,s);add_tag_y(o,l,r,t);return;}
int mid=l+r>>1; pushdown(o,l,r);
if (L<=mid) st_modify(o<<1,l,mid,L,R,s,t);
if (R>mid) st_modify(o<<1|1,mid+1,r,L,R,s,t);
update(o,l,r);
} void modify(int l,int r,double s,double t)
{st_modify(1,1,n,l,r,s,t);} void st_cov(int o,int l,int r,int L,int R)
{
if (L<=l&&r<=R){add_cov(o,l,r);return;}
int mid=l+r>>1; pushdown(o,l,r);
if (L<=mid) st_cov(o<<1,l,mid,L,R);
if (R>mid) st_cov(o<<1|1,mid+1,r,L,R);
update(o,l,r);
} void cover(int l,int r)
{st_cov(1,1,n,l,r);} void st_query(int o,int l,int r,int L,int R)
{
if (L<=l&&r<=R){sumx+=sum_x[o],sumy+=sum_y[o],sumxx+=sum_xx[o];sumxy+=sum_xy[o];return;}
int mid=l+r>>1; pushdown(o,l,r);
if (R<=mid) st_query(o<<1,l,mid,L,R);
else if (L>mid) st_query(o<<1|1,mid+1,r,L,R);
else st_query(o<<1,l,mid,L,R),st_query(o<<1|1,mid+1,r,L,R);
} double query(int l,int r)
{
sumx=0; sumy=0;sumxx=0; sumxy=0;
double len=r-l+1; st_query(1,1,n,l,r);
printf("%.10lf\n",(len*sumxy-sumx*sumy)/(len*sumxx-sumx*sumx));
} void build(int o,int l,int r)
{
if (l==r)
{
sum_x[o]=x[l]; sum_y[o]=y[l];
sum_xy[o]=1LL*x[l]*y[l]; sum_xx[o]=1LL*x[l]*x[l];
return ;
}
int mid=l+r>>1;
build(o<<1,l,mid);
build(o<<1|1,mid+1,r);
update(o,l,r);
} int main()
{
scanf("%d%d",&n,&m);
F(i,1,n) scanf("%d",&x[i]);
F(i,1,n) scanf("%d",&y[i]);
build(1,1,n);
F(i,1,m)
{
int opt,l,r,s,t;
scanf("%d",&opt);
switch(opt)
{
case 1:
scanf("%d%d",&l,&r);
query(l,r);
break;
case 2:
scanf("%d%d%d%d",&l,&r,&s,&t);
modify(l,r,(double)s,(double)t);
break;
case 3:
scanf("%d%d%d%d",&l,&r,&s,&t);
cover(l,r);
modify(l,r,(double)s,(double)t);
}
}
}
BZOJ 4821 [Sdoi2017]相关分析 ——线段树的更多相关文章
- 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\). 操作三可以看成初始化一遍,然后同操作二. ...
- (WA)BZOJ 4821: [Sdoi2017]相关分析
二次联通门 : BZOJ 4821: [Sdoi2017]相关分析 2017.8.23 Updata 妈妈!!这道题卡我!!!就是不然我过!!!!! #include <cstdio> # ...
- ●BZOJ 4821 [Sdoi2017]相关分析
题链: http://www.lydsy.com/JudgeOnline/problem.php?id=4821 题解: 线段树是真的恶心,(也许是我的方法麻烦了一些吧)首先那个式子可以做如下化简: ...
- [Sdoi2017]相关分析 [线段树]
[Sdoi2017]相关分析 题意:沙茶线段树 md其实我考场上还剩一个多小时写了40分 其实当时写正解也可以吧1h也就写完了不过还要拍一下 正解代码比40分短2333 #include <io ...
- 【BZOJ4821】[Sdoi2017]相关分析 线段树
[BZOJ4821][Sdoi2017]相关分析 Description Frank对天文学非常感兴趣,他经常用望远镜看星星,同时记录下它们的信息,比如亮度.颜色等等,进而估算出星星的距离,半径等等. ...
- bzoj 4821 [Sdoi2017]相关分析
题面 https://www.lydsy.com/JudgeOnline/problem.php?id=4821 题解 做法显然 就是维护一颗线段树 里面装4个东西 区间x的和 区间y的和 区间$x^ ...
- 洛谷P3707 [SDOI2017]相关分析(线段树)
题目描述 Frank对天文学非常感兴趣,他经常用望远镜看星星,同时记录下它们的信息,比如亮度.颜色等等,进而估算出星星的距离,半径等等. Frank不仅喜欢观测,还喜欢分析观测到的数据.他经常分析两个 ...
- SDOI2017相关分析 线段树
题目 https://loj.ac/problem/2005 思路 \[ \sum_{L}^{R}{(x_i-x)^{2}} \] \[ \sum_{L}^{R}{(x_i^2-2*x_i*x+x^{ ...
随机推荐
- PHP生成类似类似优酷、腾讯视频等其他视频链的ID
不知道你注意了没有,类似优酷.腾讯视频等其他视频链接似乎类似这样的 http://v.youku.com/v_show/id_XNjA5MjE5OTM2.html 注意id_xxx那段,是不是看不懂了 ...
- Ubuntu18.04偏好设置
以下只是我个人的一些偏好设置,会持续更新 一.安装中文输入法和换用中文界面 Ubuntu18.04下安装搜狗输入法 Ubuntu18.04如何从英文界面更改为中文界面 二.更改系统自带设置 1.设置- ...
- Google 出品的 Java 编码规范,强烈推荐,权威又科学!
原文:google.github.io/styleguide/javaguide.html 译者:Hawstein 来源:hawstein.com/2014/01/20/google-java-sty ...
- javaweb基础(5)_servlet原理
一.Servlet简介 Servlet是sun公司提供的一门用于开发动态web资源的技术. Sun公司在其API中提供了一个servlet接口,用户若想用发一个动态web资源(即开发一个Java程序向 ...
- [LUOGU] P4251 [SCOI2015]小凸玩矩阵
行列看成点,格子看成边,二分一个边权,删去大于它的边,新图上的最大流>k则答案可以更优,小于k则调整左边界. #include<algorithm> #include<iost ...
- Docker 学习基本操作与守护式容器
Docker 学习基本操作与守护式容器 容器操作 运行容器 docker run --name指定名字 -istdin -ttty虚拟终端 在终端中用 exit 即可退出容器,并结束运行 查看容器 p ...
- 关于上传文件 非ajax提交 得到后台数据问题
<form name="configForm" id="configForm" method="post" action=" ...
- 安装ElasticSearch 6.1.1 head插件
https://blog.csdn.net/zoubf/article/details/79007908 主要参考了这个blog 才完成所有的配置,很好的参考资料
- Python学习笔记:time模块和datetime模块(时间和日期)
time模块 time模块通常用来操作时间戳信息(各种“秒”),常用的方法有: time.sleep(seconds):将当前程序阻塞指定秒数,然后继续运行程序. time.time():返回当前时间 ...
- SQL防止重复提交和Filter
/class User package com.neuedu.bean; import java.io.Serializable; public class User implements Seria ...