P3707 [SDOI2017]相关分析

最新推荐文章于 2024-10-12 16:35:22 发布
最新推荐文章于 2024-10-12 16:35:22 发布
阅读量108 收藏
点赞数
分类专栏： 线段树
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/lycccccccc/article/details/120336443
版权

线段树 回归方程 区间操作 维护系数 算法实现
关键词由CSDN通过智能技术生成

题目链接： P3707 [SDOI2017]相关分析

题目大意： 你有一个数列，你需要维护某一段数的回归方程的系数 a a a a = ∑ L R ( x i − x ˉ ) ( y i − y ˉ ) ∑ L R ( x i − x ˉ ) 2 a=\frac{\sum_{L}^R{(x_i-\bar{x})(y_i-\bar{y})}}{\sum_L^R(x_i-\bar{x})^2} a=∑LR​(xi​−xˉ)2∑LR​(xi​−xˉ)(yi​−yˉ​)​有三个操作分别是求区间 [ L , R ] [L,R] [L,R]上的 a a a，将区间 [ L , R ] [L,R] [L,R]上的所有 x + s x+s x+s,所有 y + t y+t y+t，以及将区间 [ L , R ] [L,R] [L,R]上的 x x x值改为 i + s i+s i+s, y y y的值改为 i + t i+t i+t其中 ( L ≤ i ≤ R ) (L\leq{i}\leq{R}) (L≤i≤R)

题目分析： 这个题看起来不难，但是一看就麻烦的要死；我们先把 a a a的那个式子拆开，可以写为 a = ∑ L R x i y i − ∑ L R x i y ˉ ∑ L R ( x i 2 − 2 x i x ˉ + x ˉ 2 ) a=\frac{\sum_L^Rx_iy_i-\sum_L^Rx_i\bar{y}}{\sum_L^R(x_i^2-2x_i\bar{x}+\bar{x}^2)} a=∑LR​(xi2​−2xi​xˉ+xˉ2)∑LR​xi​yi​−∑LR​xi​yˉ​​然后就可以考虑通过线段树维护这个式子了，我们需要维护的量有 ∑ x y , ∑ x 2 , ∑ y , ∑ x \sum xy,\sum x^2,\sum y,\sum x ∑xy,∑x2,∑y,∑x
先看比较简单的操作2，对于 x x x和 y y y的修改都是简单的区间操作，考虑对于 ∑ L R x y \sum_L^R{xy} ∑LR​xy，修改之后变为了 ∑ L R ( x + s ) ( y + t ) \sum_L^R{(x+s)(y+t)} ∑LR​(x+s)(y+t),因此我们要在原来的值的基础上加上 t ∑ L R x + s ∑ L R y + ( R − L + 1 ) ∗ s t t\sum_L^R{x}+s\sum_L^R{y}+(R-L+1)*st t∑LR​x+s∑LR​y+(R−L+1)∗st，对于维护 x 2 x^2 x2的操作与之类似。需要注意的一点是2操作与3操作的次序问题，我们假设总是先执行3操作，后执行2操作，因此需要在执行3操作时将之前的2操作清零。
然后我们看一下比较复杂的三操作，对于 ∑ L R x y \sum_L^R xy ∑LR​xy 将其修改为 ∑ L R ( i + s ) ( i + t ) \sum_L^R{(i+s)(i+t)} ∑LR​(i+s)(i+t)我们只需要求出 ∑ L R ( i 2 + ( s + t ) i + s t ) \sum_L^R{(i^2+(s+t)i+st)} ∑LR​(i2+(s+t)i+st)，以这个值来替换原来的值，同样 x 2 x^2 x2可以用类似的方法求出来。下面放题目代码

题目代码：

#include<stdio.h>
#include<algorithm>
#include<string.h>
#define mxn 100005
using namespace std;
int n,m;
struct Tree{
	double xsum,ysum,xysum,x2sum;
	double x1lazy,y1lazy,x2lazy,y2lazy;
	int lazy2,lazy1;
}tree[mxn<<2];
double x[mxn],y[mxn];
void pushup(int loc){
	tree[loc].xsum=tree[loc<<1].xsum+tree[loc<<1|1].xsum;
	tree[loc].ysum=tree[loc<<1].ysum+tree[loc<<1|1].ysum;
	tree[loc].x2sum=tree[loc<<1].x2sum+tree[loc<<1|1].x2sum;
	tree[loc].xysum=tree[loc<<1].xysum+tree[loc<<1|1].xysum;
}
void Build(int loc,int L,int R){
	if(L==R){
		tree[loc].xsum=x[L];
		tree[loc].ysum=y[L];
		tree[loc].x2sum=x[L]*x[L];
		tree[loc].xysum=x[L]*y[L];
		return;
	}
	int mid=(L+R)>>1;
	Build(loc<<1,L,mid);
	Build(loc<<1|1,mid+1,R);
	pushup(loc);
}
double calc(double x){
	return x*(x+1)*(2*x+1)/6;
}
void pushdown(int loc,int L,int R){
	int mid=(L+R)>>1;
	if(tree[loc].lazy2){
		double s=tree[loc].x2lazy,t=tree[loc].y2lazy;
		tree[loc<<1].lazy2=tree[loc<<1|1].lazy2=1;
		tree[loc<<1].x2lazy=tree[loc<<1|1].x2lazy=s;
		tree[loc<<1].y2lazy=tree[loc<<1|1].y2lazy=t;
		tree[loc<<1].x1lazy=tree[loc<<1|1].x1lazy=0;
		tree[loc<<1].y1lazy=tree[loc<<1|1].y1lazy=0;
		tree[loc<<1].lazy1=tree[loc<<1|1].lazy1=0;
		tree[loc].lazy2=0;tree[loc].x2lazy=tree[loc].y2lazy=0;
		tree[loc<<1].xsum=(mid-L+1)*1.0*(s+(mid+L)*0.5);
		tree[loc<<1].ysum=(mid-L+1)*1.0*(t+(mid+L)*0.5);
		tree[loc<<1].x2sum=calc(mid)-calc(L-1)+2*s*((mid-L+1)*1.0)*((mid+L)*0.5)+(mid-L+1)*s*s;
		tree[loc<<1].xysum=calc(mid)-calc(L-1)+(s+t)*((mid-L+1)*1.0)*((mid+L)*0.5)+(mid-L+1)*s*t;
		tree[loc<<1|1].xsum=(R-mid)*1.0*(s+(mid+R+1)*0.5);
		tree[loc<<1|1].ysum=(R-mid)*1.0*(t+(mid+R+1)*0.5);
		tree[loc<<1|1].x2sum=calc(R)-calc(mid)+2*s*((R-mid)*1.0)*((mid+R+1)*0.5)+(R-mid)*s*s;
		tree[loc<<1|1].xysum=calc(R)-calc(mid)+(s+t)*((R-mid)*1.0)*((mid+R+1)*0.5)+(R-mid)*s*t;
	}
	if(tree[loc].lazy1){
		double s=tree[loc].x1lazy,t=tree[loc].y1lazy;
		tree[loc<<1].x2sum+=2*s*tree[loc<<1].xsum+(mid-L+1)*s*s;
		tree[loc<<1].xysum+=t*tree[loc<<1].xsum+s*tree[loc<<1].ysum+(mid-L+1)*s*t;
		tree[loc<<1].xsum+=(mid-L+1)*s;
		tree[loc<<1].ysum+=(mid-L+1)*t;
		tree[loc<<1].x1lazy+=s;tree[loc<<1].y1lazy+=t;tree[loc<<1].lazy1=1;
		tree[loc<<1|1].x2sum+=2*s*tree[loc<<1|1].xsum+(R-mid)*s*s;
		tree[loc<<1|1].xysum+=t*tree[loc<<1|1].xsum+s*tree[loc<<1|1].ysum+(R-mid)*s*t;
		tree[loc<<1|1].xsum+=(R-mid)*s;
		tree[loc<<1|1].ysum+=(R-mid)*t;
		tree[loc<<1|1].x1lazy+=s;tree[loc<<1|1].y1lazy+=t;tree[loc<<1|1].lazy1=1;
		tree[loc].lazy1=0;tree[loc].x1lazy=0;tree[loc].y1lazy=0;
	}
}
void Update1(int loc,int L,int R,int l,int r,double s,double t){
	if(l<=L&&r>=R){
		tree[loc].x2sum+=2*s*tree[loc].xsum+(R-L+1)*s*s;
		tree[loc].xysum+=t*tree[loc].xsum+s*tree[loc].ysum+(R-L+1)*s*t;
		tree[loc].xsum+=(R-L+1)*s;
		tree[loc].ysum+=(R-L+1)*t;
		tree[loc].x1lazy+=s;
		tree[loc].y1lazy+=t;
		tree[loc].lazy1=1;
		return;
	}
	int mid=(L+R)>>1;
	pushdown(loc,L,R);
	if(r<=mid)Update1(loc<<1,L,mid,l,r,s,t);
	else if(l>mid)Update1(loc<<1|1,mid+1,R,l,r,s,t);
	else Update1(loc<<1,L,mid,l,mid,s,t),Update1(loc<<1|1,mid+1,R,mid+1,r,s,t);
	pushup(loc);
}
void Update2(int loc,int L,int R,int l,int r,double s,double t){
	if(l<=L&&r>=R){
		tree[loc].xsum=(R-L+1)*1.0*(s+(R+L)*0.5);
		tree[loc].ysum=(R-L+1)*1.0*(t+(R+L)*0.5);
		tree[loc].x2sum=calc(R)-calc(L-1)+2*s*((R-L+1)*1.0)*((R+L)*0.5)+(R-L+1)*s*s;
		tree[loc].xysum=calc(R)-calc(L-1)+(s+t)*((R-L+1)*1.0)*((R+L)*0.5)+(R-L+1)*s*t;
		tree[loc].lazy2=1;tree[loc].x2lazy=s;tree[loc].y2lazy=t;
		tree[loc].x1lazy=tree[loc].y1lazy=0;tree[loc].lazy1=0;
		return;
	}
	int mid=(L+R)>>1;
	pushdown(loc,L,R);
	if(r<=mid)Update2(loc<<1,L,mid,l,r,s,t);
	else if(l>mid)Update2(loc<<1|1,mid+1,R,l,r,s,t);
	else Update2(loc<<1,L,mid,l,mid,s,t),Update2(loc<<1|1,mid+1,R,mid+1,r,s,t);
	pushup(loc);
}
double Query_x(int loc,int L,int R,int l,int r){
	if(l<=L&&r>=R)return tree[loc].xsum;
	int mid=(L+R)>>1;
	pushdown(loc,L,R);
	if(r<=mid)return Query_x(loc<<1,L,mid,l,r);
	else if(l>mid)return Query_x(loc<<1|1,mid+1,R,l,r);
	else return Query_x(loc<<1,L,mid,l,mid)+Query_x(loc<<1|1,mid+1,R,mid+1,r);
}
double Query_y(int loc,int L,int R,int l,int r){
	if(l<=L&&r>=R)return tree[loc].ysum;
	int mid=(L+R)>>1;
	pushdown(loc,L,R);
	if(r<=mid)return Query_y(loc<<1,L,mid,l,r);
	else if(l>mid)return Query_y(loc<<1|1,mid+1,R,l,r);
	else return Query_y(loc<<1,L,mid,l,mid)+Query_y(loc<<1|1,mid+1,R,mid+1,r);
}
double Query_x2(int loc,int L,int R,int l,int r){
	if(l<=L&&r>=R)return tree[loc].x2sum;
	int mid=(L+R)>>1;
	pushdown(loc,L,R);
	if(r<=mid)return Query_x2(loc<<1,L,mid,l,r);
	else if(l>mid)return Query_x2(loc<<1|1,mid+1,R,l,r);
	else return Query_x2(loc<<1,L,mid,l,mid)+Query_x2(loc<<1|1,mid+1,R,mid+1,r);
}
double Query_xy(int loc,int L,int R,int l,int r){
	if(l<=L&&r>=R)return tree[loc].xysum;
	int mid=(L+R)>>1;
	pushdown(loc,L,R);
	if(r<=mid)return Query_xy(loc<<1,L,mid,l,r);
	else if(l>mid)return Query_xy(loc<<1|1,mid+1,R,l,r);
	else return Query_xy(loc<<1,L,mid,l,mid)+Query_xy(loc<<1|1,mid+1,R,mid+1,r);
}
void work(int L,int R){
	double xb=Query_x(1,1,n,L,R)/(1.0*(R-L+1)),yb=Query_y(1,1,n,L,R)/(1.0*(R-L+1));
	double xs=Query_x(1,1,n,L,R),ys=Query_y(1,1,n,L,R),x2s=Query_x2(1,1,n,L,R),xys=Query_xy(1,1,n,L,R);
	//printf("%.2f %.2f %.2f %.2f\n",xb,yb,x2s,xys);
	double ans=xys-yb*xs-xb*ys+xb*yb*(R-L+1);
	ans/=(x2s-2*xb*xs+(R-L+1)*xb*xb);
	printf("%.12f\n",ans);
}
int main()
{
	scanf("%d %d",&n,&m);
	for(int i=1;i<=n;i++)
		scanf("%lf",&x[i]);
	for(int i=1;i<=n;i++)
		scanf("%lf",&y[i]);
	Build(1,1,n);
	while(m--){
		int opt,L,R;
		double S,T;
		scanf("%d",&opt);
		//for(int i=1;i<=n;i++)
		//	printf("%.2f %.2f\n",Query_x(1,1,n,i,i),Query_y(1,1,n,i,i));
		//putchar(10);
		if(m==0)
			int a=1;
		if(opt==1)scanf("%d %d",&L,&R),work(L,R);
		if(opt==2)scanf("%d %d %lf %lf",&L,&R,&S,&T),Update1(1,1,n,L,R,S,T);
		if(opt==3)scanf("%d %d %lf %lf",&L,&R,&S,&T),Update2(1,1,n,L,R,S,T);
	}
	return 0;
}
lycccccccc
关注
专栏目录
bzoj 4818: [Sdoi2017]序列计数 动态规划+矩阵乘法
beginend
04-11 1743
题意Alice想要得到一个长度为nn的序列，序列中的数都是不超过mm的正整数，而且这nn个数的和是pp的倍数。 Alice还希望，这nn个数中，至少有一个数是质数。 Alice想知道，有多少个序列满足她的要求。 对100%100\%的数据，1≤n≤109,1≤m≤2×107,1≤p≤1001\leq n \leq 10^9,1\leq m \leq 2×10^7,1\leq p\leq 100
水题-sdoi2017-相关分析
lym01803的博客
05-30 360
线段树.维护四个值:∑x,∑y,∑xy,∑x2\sum x,\sum y,\sum xy,\sum x^2 3操作拆分成一个赋值操作外加一个2操作就行了。预处理一下平方和之类的就好了。 注意几个问题： 1.不要用long long（会爆long long）,要用double. 2.最后是把2,3操作的lazy标记分开写。处理3标记时，直接清空2标记.贴个代码：#include<cstdio
[Luogu 3707] SDOI2017 相关分析
weixin_30673715的博客
02-28 64
[Luogu 3707] SDOI2017 相关分析 前言 Capella 和 Frank 一样爱好天文学。 她常在冬季的夜晚，若有所思地望着东北方上空的五边形中，最为耀眼的一个顶点。 那一抹金黄曾带给她多少美好的遐想！ 直到有一天，她做了这一题，并饱受线段树的折磨。 概述 这是一道区间操作的题目，解法有线段树与分块等。 题主选择线段树进行讲解（因为没用分块写这题）。 准备与计算答案 ...
luogu P3707 [SDOI2017]相关分析
weixin_30920513的博客
11-06 73
传送门 对于题目要求的东西,考虑拆开懒得拆了 ,可以发现有\(\sum x\sum y\sum x^2\sum xy\)四个变量影响最终结果,考虑维护这些值 下面记\(l,r\)为区间两端点 首先是区间加操作,可以这样维护\[\sum(x+s)=\sum x+(r-l+1)s\]\[\sum(y+s)=\sum y+(r-l+1)s\]\[\sum(x+s)^2=\sum x^2+2s\sum x...
bzoj4821 && luogu3707 SDOI2017相关分析（线段树，数学）
weixin_30484247的博客
12-22 78
题目大意 给定n个元素的数列，每一个元素有x和y两种元素，现在有三种操作： \(1\ L\ R\) 设\(xx\)为\([l,r]\)的元素的\(x_i\)的平均值，\(yy\)同理 求 \(\frac{\sum_{i=L}^R(x_i-xx)(y_i-yy)} {\sum_{i=L}^R(xi-xx)^2}\) \(2\ L\ R\ S\ T\\) 将\([L,R]\)中的每个元素的\(x...
SDOI2017 相关分析
weixin_34198881的博客
03-11 108
SDOI2017 相关分析 题意： 题目传送门 题解： 看了\(loj\) 300多份代码，似乎我是唯二写分块的？刚开始感觉线段树比较难写，就开始码分块了。现在知道为什么没人写分块了…… 很容易想到将式子进行拆分，然后维护各种东西： 操作1：\(a = \frac{\sum_{i = L}^{R}(x_i * y_i + x_i * \overline{y} + \overline{x} * y_...
[题目] Luogu P3707 [SDOI2017]相关分析
weixin_30505225的博客
03-29 76
参考资料：[Luogu 3707] SDOI2017 相关分析 P3707 [SDOI2017]相关分析 TFRAC FRAC DFRAC \(\tfrac{\sum}{1}\) \(\frac{\sum}{1}\) \(\dfrac{\sum}{1}\) \[\bar{x}=\frac{1}{R-L+1}\sum x_i​\] \[\bar{y}=\frac{1}{R-L...
洛谷 P3707 [SDOI2017]相关分析 线段树
liangzihao1的博客
06-23 326
题目描述 Frank对天文学非常感兴趣，他经常用望远镜看星星，同时记录下它们的信息，比如亮度、颜色等等，进而估算出星星的距离，半径等等。 Frank不仅喜欢观测，还喜欢分析观测到的数据。他经常分析两个参数之间（比如亮度和半径）是否存在某种关系。 现在Frank要分析参数 X 与 Y 之间的关系。他有 n 组观测数据，第 i 组观测数据记录了 xixix_i和 yiyiy_i。 他需要一下几...
[SDOI2017]相关分析
araw94333的博客
08-28 152
把求\(a\)的式子拆开，然后把所有的需要求的东西用线段树维护即可。 // luogu-judger-enable-o2 /* mail: mleautomaton@foxmail.com author: MLEAutoMaton This Code is made by MLEAutoMaton */ #include<stdio.h> #include<...
洛谷P3707 [SDOI2017]相关分析(线段树)
weixin_30526593的博客
03-03 63
题目描述 Frank对天文学非常感兴趣，他经常用望远镜看星星，同时记录下它们的信息，比如亮度、颜色等等，进而估算出星星的距离，半径等等。 Frank不仅喜欢观测，还喜欢分析观测到的数据。他经常分析两个参数之间（比如亮度和半径）是否存在某种关系。 现在Frank要分析参数XX与YY之间的关系。他有nn组观测数据，第ii组观测数据记录了x_ixi​和y_iyi​。他需要一下几种操作 ...
bzoj 4820: [Sdoi2017]硬币游戏 概率dp+高斯消元+KMP
beginend
01-05 516
题意 周末同学们非常无聊，有人提议，咱们扔硬币玩吧，谁扔的硬币正面次数多谁胜利。 大家纷纷觉得这个游戏非常符合同学们的特色，但只是扔硬币实在是太单调了。 同学们觉得要加强趣味性，所以要找一个同学扔很多很多次硬币，其他同学记录下正反面情况。 用 H 表示正面朝上， 用 T 表示反面朝上，扔很多次硬币后，会得到一个硬币序列。比如HTT 表示第一次正面朝上，后两次反面朝上。 但扔到什么时候停止
P2484 [SDOI2011]打地鼠 题解
if_OF的博客
11-04 303
锤子每次只能将每个洞里打掉一只地鼠，所以对于每次击打，都会打掉r*c只地鼠。而每次击打的范围内的每个洞里必须都要有地鼠，所以地鼠的总数必须被r*c整除。
P5566 [SDOI2008]红黑树 贪心|动态规划
youth518的博客
09-20 350
题意： 求一颗nnn个节点的红黑树红色节点数目最多和最少 红黑树是满足如下性质的的染色二叉搜索树： 每个结点被染成红色或黑色； 每个前端结点为黑色结点； 任一红结点的子结点均为黑结点； 在从任一结点到其子孙前端结点的所有路径上具有相同的黑结点数。 二叉搜索树结点中的空指针看作是指向一个空结点，则称这类空结点为二叉搜索树的前端结点。并规定所有前端结点的高度为 -1。 范围&性质：1≤n≤50001\le n\le 50001≤n≤5000 分析： 贪心 首先对于一颗nnn个点的红黑树，他必定存
[SDOI2006]保安站岗 洛谷p2458
MrTinTin的博客
10-30 333
题目描述 五一来临，某地下超市为了便于疏通和指挥密集的人员和车辆，以免造成超市内的混乱和拥挤，准备临时从外单位调用部分保安来维持交通秩序。 已知整个地下超市的所有通道呈一棵树的形状；某些通道之间可以互相望见。总经理要求所有通道的每个端点（树的顶点）都要有人全天候看守，在不同的通道端点安排保安所需的费用不同。 一个保安一旦站在某个通道的其中一个端点，那么他除了能看守住他所站的那个端点，也能看到...
每日学习一个数据结构-图
wendao76的专栏
10-08 1179
在数学中，图是描述于一组对象的结构，其中某些对象对在某种意义上是“相关的”。这些对象对应于称为顶点的数学抽象（也称为节点或点），并且每个相关的顶点对都称为边（也称为链接或线）。通常，图形以图解形式描绘为顶点的一组点或环，并通过边的线或曲线连接。图（Graph）是由顶点的有穷非空集合和顶点之间的连线（边）的集合组成，通常表示为G=(V, E)，其中G表示一个图，V(G)代表图G中的顶点集合，E(G)代表图G中的边集合。
数据结构-5.1.树的定义和基本术语
ADCvbV的博客
10-08 210
1.祖先结点：从所指节点上方开始，沿一个方向一直到根结点之间(包括根结点)所遇到的结点都是祖先结点：如：结点"你"的祖先结点有"父亲"和"爷爷";2.子孙结点：从一个结点出发，他的所有分支以及他的分支的分支，一直到没有分支时的所有结点:如：结点"父亲"的子孙结点有"你","K","L"和"F";3.双亲结点(父结点/前驱结点)：如结点"你"的双亲结点为"父亲";4.孩子结点(后继结点)：如结点"你"的孩子结点为"K"和"L";
利用编程思维做题之求节点 x 在二叉树中的双亲节点算法
跟着信仰走
10-12 541
通过递归遍历二叉树，可以有效地查找目标节点的双亲节点。如果树结构较为复杂，递归的方法也能较为清晰地处理问题。
浙大数据结构：07-图5 Saving James Bond - Hard Version
qq_74924951的博客
10-08 1389
bfs难题
numexpr-2.8.3-cp38-cp38-win_amd64.whl
最新发布
10-14
numexpr-2.8.3-cp38-cp38-win_amd64.whl
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值