多维前缀和 学习笔记 模板及例题

最新推荐文章于 2024-08-11 17:15:11 发布
最新推荐文章于 2024-08-11 17:15:11 发布
阅读量182 收藏 1
点赞数 1
分类专栏： 模板 学习笔记 例题 文章标签： 算法
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_52028570/article/details/119651155
版权
模板 同时被 3 个专栏收录
21 篇文章 1 订阅
订阅专栏
例题
14 篇文章 0 订阅
订阅专栏
学习笔记
6 篇文章 1 订阅
订阅专栏

一维前缀和大家都会，二维就更好理解了。
一维初始化：

	for(i=1;i<=n;i++)
	{
		q[i]=(q[i-1]+a[i]);
	}

二维初始化其实花个图就可以了，这里直接上模板：

	q[1][1]=a[1][1];
	for(i=2;i<=m;i++)
	{
		q[1][i]=(q[1][i-1]+a[1][i]);
	}
	for(i=2;i<=n;i++)
	{
		q[i][1]=(q[i-1][1]+a[i][1]);
	}
	for(i=2;i<=n;i++)
	{
		for(j=2;j<=m;j++)
		{
			q[i][j]=(a[i][j]+q[i-1][j]+q[i][j-1]-q[i-1][j-1]);
		}
	}

这么写更简单：

for(i=1;i<=n;i++)
{
	for(j=1;j<=n;j++)
	{
		a[i][j]+=a[i][j-1];
	}
}
for(i=1;i<=n;i++)
{
	for(j=1;j<=n;j++)
	{
		a[i][j]+=a[i-1][j];
	}
}

三维推荐用后者：

for(i=1;i<=n;i++)
{
	for(j=1;j<=n;j++)
	{
		for(k=1;k<=n;k++)
		{
			a[i][j][k]+=a[i][j][k-1];
		}
	}
}
for(i=1;i<=n;i++)
{
	for(j=1;j<=n;j++)
	{
		for(k=1;k<=n;k++)
		{
			a[i][j][k]+=a[i][j-1][k];
		}
	}
}
for(i=1;i<=n;i++)
{
	for(j=1;j<=n;j++)
	{
		for(k=1;k<=n;k++)
		{
			a[i][j][k]+=a[i-1][j][k];
		}
	}
}

一维查询：


long long ywqzh(int x,int y)
{
	return (q[y]-q[x-1]);
}

二维则还要加上一个公共部分：
x,y是左上角坐标，xx,yy是右下角坐标。


long long ewqzh(int x,int y,int xx,int yy)
{
	return (q[xx][yy]-q[xx][y-1]-q[x-1][yy]+q[x-1][y-1]);
}

三维就更简单了 ：

inline long long swqzh(int aa,int bb,int cc,int i,int j,int k)
{
	return (q[aa][bb][cc]-q[i-1][bb][cc]-q[aa][j-1][cc]-q[aa][bb][k-1]+q[i-1][j-1][cc]+q[aa][j-1][k-1]+q[i-1][bb][k-1]-q[i-1][j-1][k-1]);
}

在这里插入图片描述

#include<bits/stdc++.h>
using namespace std;
int n,m,ma,mb,i,j,ii,jj,k,a[1111][1111],ax,ay,axx,ayy,hx[1111111],hy[1111111],hxx[1111111],hyy[1111111],kkk;
bool pc[12255872],pcc[12255872],pccc[12255872],pcccc[21788234];
long long q[1111][1111],h[1111111],qzh=0,has,hass,hasss,hassss;
long long ewqzh(int x,int y,int xx,int yy)
{
	return (q[xx][yy]-q[xx][y-1]-q[x-1][yy]+q[x-1][y-1]);
}
void up()
{
	int kkg=kkk;
	while(kkg>1&&h[kkg]<h[kkg/2])
	{
		h[0]=h[kkg];
		h[kkg]=h[kkg/2];
		h[kkg/2]=h[0];
		hx[0]=hx[kkg];
		hx[kkg]=hx[kkg/2];
		hx[kkg/2]=hx[0];
		hy[0]=hy[kkg];
		hy[kkg]=hy[kkg/2];
		hy[kkg/2]=hy[0];
		hxx[0]=hxx[kkg];
		hxx[kkg]=hxx[kkg/2];
		hxx[kkg/2]=hxx[0];
		hyy[0]=hyy[kkg];
		hyy[kkg]=hyy[kkg/2];
		hyy[kkg/2]=hyy[0];
		kkg/=2;
	}
}
void down()
{
	int kkgx=1,kkgy=2;
	while(((kkgy)<=kkk&&h[kkgy]<h[kkgx])||((kkgy+1)<=kkk&&h[kkgy+1]<h[kkgx]))
	{
		if(kkgy+1<=kkk&&h[kkgy]>h[kkgy+1])
		{
			kkgy++;
		}
		h[0]=h[kkgx];
		h[kkgx]=h[kkgy];
		h[kkgy]=h[0];
		hx[0]=hx[kkgx];
		hx[kkgx]=hx[kkgy];
		hx[kkgy]=hx[0];
		hy[0]=hy[kkgx];
		hy[kkgx]=hy[kkgy];
		hy[kkgy]=hy[0];
		hxx[0]=hxx[kkgx];
		hxx[kkgx]=hxx[kkgy];
		hxx[kkgy]=hxx[0];
		hyy[0]=hyy[kkgx];
		hyy[kkgx]=hyy[kkgy];
		hyy[kkgy]=hyy[0];
		kkgx=kkgy;
		kkgy*=2;
	}
}
void add(long long xxx,int x,int y,int xx,int yy)
{
	kkk++;
	h[kkk]=xxx;
	hx[kkk]=x;
	hy[kkk]=y;
	hxx[kkk]=xx;
	hyy[kkk]=yy;
	up();
}
void pop(int xxx)
{
	h[0]=h[xxx];
	h[xxx]=h[kkk];
	h[kkk]=h[0];
	hx[0]=hx[xxx];
	hx[xxx]=hx[kkk];
	hx[kkk]=hx[0];
	hy[0]=hy[xxx];
	hy[xxx]=hy[kkk];
	hy[kkk]=hy[0];
	hxx[0]=hxx[xxx];
	hxx[xxx]=hxx[kkk];
	hxx[kkk]=hxx[0];
	hyy[0]=hyy[xxx];
	hyy[xxx]=hyy[kkk];
	hyy[kkk]=hyy[0];
	kkk--;
	down();
}
int main()
{
	scanf("%d%d%d%d%d",&n,&m,&ma,&mb,&k);
	for(i=1;i<=n;i++)
	{
		for(j=1;j<=m;j++)
		{
			scanf("%d",&a[i][j]);
		}
	}
	q[1][1]=a[1][1];
	for(i=2;i<=m;i++)
	{
		q[1][i]=(q[1][i-1]+a[1][i]);
	}
	for(i=2;i<=n;i++)
	{
		q[i][1]=(q[i-1][1]+a[i][1]);
	}
	for(i=2;i<=n;i++)
	{
		for(j=2;j<=m;j++)
		{
			q[i][j]=(a[i][j]+q[i-1][j]+q[i][j-1]-q[i-1][j-1]);
		}
	}
	for(i=1;i<=(n-ma+1);i++)
	{
		for(j=1;j<=(m-mb+1);j++)
		{
			ax=i,ay=j,axx=i+ma-1,ayy=j+mb-1;
			qzh=ewqzh(ax,ay,axx,ayy);
			has=(qzh*11+ax*291143%12255871+ay*38833%12255871+axx*999983%12255871+ayy*122777%12255871)%12255871;
			pc[has]=1;
			hass=(qzh*13+ax*38833%12255871+ay*291143%12255871+axx*122777%12255871+ayy*999983%12255871)%12255871;
			pcc[has]=1;
			hasss=(qzh*13+ax*83%91815541+ay*89%91815541+axx*97%91815541+ayy*29%91815541)%91815541;
			pccc[hasss]=1;
			hasss=(qzh*2179+ax*29123%21788233+ay*113%21788233+axx*97%21788233+ayy*61%21788233)%21788233;
			pcccc[hassss]=1;
			add(qzh,ax,ay,axx,ayy);
		}
	}
	for(i=1;i<=k-1;i++)
	{
		ax=hx[1]-1,ay=hy[1],axx=hxx[1],ayy=hyy[1];
		if(ax>=1)
		{
			qzh=ewqzh(ax,ay,axx,ayy);
			has=(qzh*11+ax*291143%12255871+ay*38833%12255871+axx*999983%12255871+ayy*122777%12255871)%12255871;
			hass=(qzh*13+ax*38833%12255871+ay*291143%12255871+axx*122777%12255871+ayy*999983%12255871)%12255871;
			hasss=(qzh*13+ax*83%91815541+ay*89%91815541+axx*97%91815541+ayy*29%91815541)%91815541;
			hassss=(qzh*2179+ax*29123%21788233+ay*113%21788233+axx*97%21788233+ayy*61%21788233)%21788233;
			if(pc[has]!=1||pcc[hass]!=1||pccc[hasss]!=1||pcccc[hassss]!=1)
			{
				pc[has]=1;
				pcc[hass]=1;
				pccc[hasss]=1;
				pcccc[hassss]=1;
				add(qzh,ax,ay,axx,ayy);
			}
		}
		ax=hx[1],ay=hy[1]-1,axx=hxx[1],ayy=hyy[1];
		if(ay>=1)
		{
			qzh=ewqzh(ax,ay,axx,ayy);
			has=(qzh*11+ax*291143%12255871+ay*38833%12255871+axx*999983%12255871+ayy*122777%12255871)%12255871;
			hass=(qzh*13+ax*38833%12255871+ay*291143%12255871+axx*122777%12255871+ayy*999983%12255871)%12255871;
			hasss=(qzh*13+ax*83%91815541+ay*89%91815541+axx*97%91815541+ayy*29%91815541)%91815541;
			hassss=(qzh*2179+ax*29123%21788233+ay*113%21788233+axx*97%21788233+ayy*61%21788233)%21788233;
			if(pc[has]!=1||pcc[hass]!=1||pccc[hasss]!=1||pcccc[hassss]!=1)
			{
				pc[has]=1;
				pcc[hass]=1;
				pccc[hasss]=1;
				pcccc[hassss]=1;
				add(qzh,ax,ay,axx,ayy);
			}
		}
		ax=hx[1],ay=hy[1],axx=hxx[1]+1,ayy=hyy[1];
		if(axx<=n)
		{
			qzh=ewqzh(ax,ay,axx,ayy);
			has=(qzh*11+ax*291143%12255871+ay*38833%12255871+axx*999983%12255871+ayy*122777%12255871)%12255871;
			hass=(qzh*13+ax*38833%12255871+ay*291143%12255871+axx*122777%12255871+ayy*999983%12255871)%12255871;
			hasss=(qzh*13+ax*83%91815541+ay*89%91815541+axx*97%91815541+ayy*29%91815541)%91815541;
			hassss=(qzh*2179+ax*29123%21788233+ay*113%21788233+axx*97%21788233+ayy*61%21788233)%21788233;
			if(pc[has]!=1||pcc[hass]!=1||pccc[hasss]!=1||pcccc[hassss]!=1)
			{
				pc[has]=1;
				pcc[hass]=1;
				pccc[hasss]=1;
				pcccc[hassss]=1;
				add(qzh,ax,ay,axx,ayy);
			}
		}
		ax=hx[1],ay=hy[1],axx=hxx[1],ayy=hyy[1]+1;
		if(ayy<=m)
		{
			qzh=ewqzh(ax,ay,axx,ayy);
			has=(qzh*11+ax*291143%12255871+ay*38833%12255871+axx*999983%12255871+ayy*122777%12255871)%12255871;
			hass=(qzh*13+ax*38833%12255871+ay*291143%12255871+axx*122777%12255871+ayy*999983%12255871)%12255871;
			hasss=(qzh*13+ax*83%91815541+ay*89%91815541+axx*97%91815541+ayy*29%91815541)%91815541;
			hassss=(qzh*2179+ax*29123%21788233+ay*113%21788233+axx*97%21788233+ayy*61%21788233)%21788233;
			if(pc[has]!=1||pcc[hass]!=1||pccc[hasss]!=1||pcccc[hassss]!=1)
			{
				pc[has]=1;
				pcc[hass]=1;
				pccc[hasss]=1;
				pcccc[hassss]=1;
				add(qzh,ax,ay,axx,ayy);
			}
		}
		pop(1);
	}
	printf("%d\n",h[1]);
}

4hash爆70，以后有时间找一个更好的判重方法！
休息一下吧，去看耶稣好不好？

再来一题：
3186. 蜡笔 (Standard IO)
Time Limits: 5000 ms Memory Limits: 262144 KB Detailed Limits
Goto ProblemSet

Description
ABC生日收到N支蜡笔，每支蜡笔的颜色是三原色红绿蓝的组合，第i个蜡笔的颜色用Ri表示红色，Gi表示绿色，Bi表示蓝色。

蜡笔i和蜡笔j的颜色差异定义为max(|Ri-Rj|,|Gi-Gj|,|Bi-Bj|)，多支蜡笔的颜色差异定义为其中任意两个蜡笔颜色的最大差异值。

ABC想从N支蜡笔中选出K支出来，要求这K支蜡笔的颜色差异值最小。

Input
第一行输入N和K(2<=K<=N<=100,000)

接下来N行输入每支蜡笔的颜色值Ri,Gi,Bi(0<=Ri,Gi,Bi<=255)。

Output
输出最小颜色差异值。

Sample Input
输入1：

2 2

1 3 2

2 6 4

输入2：

3 2

3 3 4

1 6 4

1 1 2

输入3：

5 3

6 6 4

6 2 7

3 1 3

4 1 5

6 2 6

Sample Output
输出1：

3

输出2：

2

输出3：

2

Data Constraint
50%的数据满足0<=Ri,Gi,Bi<=20；

80%的数据满足0<=Ri,Gi,Bi<=50。

二分+三维前缀和：

#include<bits/stdc++.h>
using namespace std;
long long n,kkk,i,j,k,a[1111111],b[1111111],c[1111111],ans=99999999,star=0,q[301][301][301],kkg=260;
inline long long swqzh(int aa,int bb,int cc,int i,int j,int k)
{
	return (q[aa][bb][cc]-q[i-1][bb][cc]-q[aa][j-1][cc]-q[aa][bb][k-1]+q[i-1][j-1][cc]+q[aa][j-1][k-1]+q[i-1][bb][k-1]-q[i-1][j-1][k-1]);
}
bool check(int xd)
{
	int kl=0;
	for(int i=1;i<=(255-xd+1);i++)
	{
		int aa=i+xd;
		for(int j=1;j<=(255-xd+1);j++)
		{
			int bb=j+xd;
			for(int k=1;k<=(255-xd+1);k++)
			{
				int cc=k+xd;
				int kkgg=swqzh(aa,bb,cc,i,j,k);
//				cout<<i<<" "<<j<<" "<<k<<" "<<aa<<" "<<bb<<" "<<cc<<" "<<kkgg<<endl;
				if(kkgg>kl)
				{
					kl=kkgg;
				}
			}
		}
	}
	if(kl>=kkk) return true;
	return false;
}
int main()
{
	scanf("%lld%lld",&n,&kkk);
	for(i=1;i<=n;i++)
	{
		scanf("%lld%lld%lld",&a[i],&b[i],&c[i]);
		q[a[i]][b[i]][c[i]]++;
	}
	for(i=1;i<=kkg;i++)
	{
		for(j=1;j<=kkg;j++)
		{
			for(k=1;k<=kkg;k++)
			{
				q[i][j][k]+=q[i][j][k-1];
			}
		}
	}
	for(i=1;i<=kkg;i++)
	{
		for(j=1;j<=kkg;j++)
		{
			for(k=1;k<=kkg;k++)
			{
				q[i][j][k]+=q[i][j-1][k];
			}
		}
	}
	for(i=1;i<=kkg;i++)
	{
		for(j=1;j<=kkg;j++)
		{
			for(k=1;k<=kkg;k++)
			{
				q[i][j][k]+=q[i-1][j][k];
			}
		}
	}
	int l=1,r=256,mid;
	mid=(l+r)/2;
	while(l<r)
	{
		if(check(mid))
		{
			r=mid;
		}
		else
		{
			l=mid+1;
		}
		mid=(l+r)/2;
	}
//	cout<<check(3)<<endl;
	ans=l;
	printf("%lld\n",ans);
}
Starformula_神诘
关注
专栏目录
多维空间扁平响应式网页模板
01-22
在现代网页设计中，"多维空间扁平响应...总的来说，多维空间扁平响应式网页模板代表了现代网页设计的前沿，它结合了视觉吸引力、用户体验和跨平台兼容性的最佳实践，是设计师在创建引人入胜且实用的网页时的重要工具。
多维前缀和
09-19 515
问题描述 二维平面上，给定N（大约200000）个点，这些点的x和y的取值范围都是[1,3000]之间的整数，给定M（大约200000）个查询，每个查询输入一个点P（px,py）。对于每个查询，求N个点到点P的距离之和。 输入 第一行两个整数N和M 接下来N行表示N个点的x和y坐标 接下来M行表示M个查询的x和y坐标 3 3 5 5 5 10 10 5 1 1 5 5 10 10 输...
【算法笔记】前缀和，多维前缀和，一维差分
qq_46289793的博客
01-02 205
算法
前缀和（简讲）
cyyyyds857的博客
11-04 216
前缀和可以简单理解为「数列的前 项的和」，是一种重要的预处理方式，能大大降低查询的时间复杂度。C++ 标准库中实现了前缀和函数，定义于头文件中。
前缀和总结（一到多维）
最新发布
2302_77782043的博客
08-11 479
有的，读者可能发现我在写上面前缀和公式的时候都标记上了二进制数，实际上，我们可以利用二进制数来记忆前缀和的公式，这样不仅记忆轻松，在处理高维的前缀和时我们也不会漏。以上前缀和我们都可以很好的通过画图来理解和记忆，但是三维的就是立方体，我们每次需要通过画立方体来理解和记忆，实在是过于困难，有没有什么巧记的方法呢？比方说区间【l,r】,这也是很好理解的我们这段区间的部分和，就会等区间的最远端的和减去不属于区间的这部分和。我们在预处理时通常写下这样的公式，这个公式是很好理解的，（0，i）的前缀和就会等于。
高维前缀和 学习笔记
sophilex的博客
03-25 2311
显然一种暴力的方法就是枚举i的所有子集，考虑优化的话，我们可以尝试前缀和，因为显然如果有,我们可以将j的贡献都先算在z上面再传给i，而不必一个个来但是这样的话，就涉及到了高维的前缀和处理先来看看一维前缀和是怎么写的。
数据结构学习笔记.docx
11-26
数据结构是计算机科学中的核心概念，...以上就是数据结构学习笔记的主要内容，涵盖了时间复杂度、空间复杂度、线性表、栈、队列、字符串、数组和矩阵等核心概念。理解并掌握这些知识对于深入学习计算机科学至关重要。
MATLAB学习个人笔记总结.7z
04-27
通过阅读这份MATLAB学习笔记，读者将能够掌握MATLAB的基础操作和高级特性，为科研或工程实践打下坚实基础。在实际应用中，结合MATLAB的文档和示例，不断实践和探索，能进一步提升MATLAB的使用技能。
机器学习+深度学习笔记.pdf
07-25
7. Python基础语法及实操：Python的基础语法是人工智能学习的基础，Pandas、Numpy和Matplotlib是进行数据处理和可视化的必备工具。通过实际操作可以更深入地理解和掌握这些工具的使用方法，例如创建图对象、指定图像...
北大tensorflow学习笔记
11-09
在"北大tensorflow学习笔记"中，我们可以深入探索这个广泛应用于学术研究和工业实践的平台。 一、TensorFlow基础知识 1.1 **张量（Tensor）**：TensorFlow的核心概念，是多维数据数组，可以是标量、向量、矩阵甚至...
2021-08-17
weixin_53182424的博客
08-17 189
2021/8/17打卡算法前缀和&差分 前缀和 前缀和是一种重要的预处理，能大大降低查询的时间复杂度。可以简单理解为“数列的前 n 项的和”。 C++ 标准库中实现了前缀和函数 partial_sum，定义于头文件 中。 例题 有 N 个的正整数放到数组 A 里，现在要求一个新的数组B ，新数组的第 i 个数 B[i] 是原数组 A 第 0到第 i 个数的和。 输入： 5 1 2 3 4 5 输出： 1 3 6 10 15 解题思路 递推：B[0] = A[0]，对于 i>=1
【算法】前缀和(一维前缀和与二维前缀和)
kuangd_1992的博客
12-09 5721
前缀和是一种重要的预处理，能大大降低查询的时间复杂度。 【一维前缀和】 给定一个数组A[1,2,……n]，则它的前缀和数组为PrefixSum[1..n]。定义为：PrefixSum[i] = A[0]+A[1]+...+A[i-1]； 【例子】 A[5,6,7,8] --> PrefixSum[5,11,18,26] PrefixSum[0] =A[0] ; PrefixSum...
二维前缀和
Yee_mo的博客
02-14 724
C++ 二维前缀和、子矩阵的和
一维前缀和与二维前缀和
weixin_46310680的博客
11-05 241
P1719 最大加权矩形 #题目描述 为了更好的备战NOIP2013，电脑组的几个女孩子LYQ,ZSC,ZHQ认为，我们不光需要机房，我们还需要运动，于是就决定找校长申请一块电脑组的课余运动场地，听说她们都是电脑组的高手，校长没有马上答应他们，而是先给她们出了一道数学题，并且告诉她们：你们能获得的运动场地的面积就是你们能找到的这个最大的数字。 校长先给他们一个N*N矩阵。要求矩阵中最大加权矩形，即矩阵的每一个元素都有一权值，权值定义在整数集上。从中找一矩形，矩形大小无限制，是其中包含的所有元素的和最大 。矩
知识点 - 高维前缀和
qq_41848675的博客
09-03 907
知识点 - 高维前缀和 解决问题类型： 求解多维数组的前缀和，或者是求一个集合中所有有效子集的累加。 前置知识 前缀和 一维前缀和： for(int i=1;i<=n;i++) b[i]=b[i-1]+a[i]; 二维前缀和： for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) b[i][j]=b[i-1][j...
前缀和、二维前缀和与差分的小总结
热门推荐
Zookkk的博客
08-17 2万+
在了解二维前缀和之前，我们首先需要了解一下什么是前缀和。 如果我给你一串长度为n的数列a1,a2,a3......an,再给出m个询问，每次询问给出L，R两个数，要求给出区间[L,R]里的数的和，你会怎么做，若是没有了解过前缀和的人看到这道题的想法可能是对于m次询问，我每次都遍历一遍它给的区间，计算出答案，这样子的方法固然没错，但是其时间复杂度达到了O(n*m)，如果数据量稍微大一点就有可能超时...
高维前缀和总结
dft539533的博客
09-25 4115
前言 今天中午不知怎么的对这个东西产生了兴趣，感觉很神奇，结果花了一个中午多的时间来看QAQ 下面说下自己的理解。 高维前缀和一般解决这类问题： 对于所有的\(i,0\leq i\leq 2^n-1\)，求解\(\sum_{j\subset i}a_j\)。 显然，这类问题可以直接枚举子集求解，但复杂度为\(O(3^n)\)。如果我们施展高维前缀和的话，复杂度可以到\(O(n\cdo...
高维前缀和
继续学吧
11-22 2440
高维前缀和就是求关于一个集合子集(或超集)的状态的和 牛客有一道题写的很好 传送门 题面就已经说明了高维前缀和的原理 链接：https://ac.nowcoder.com/acm/contest/167/C 来源：牛客网 对于一个一维数组求部分和，可以使用如下代码 for (int i = 1; i &lt;= n; i++) { a[i] += a[i - 1]; } 对于一个二维数组...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值